#include <iostream>
#include <array>
#include <complex>
#include "../logging.h"
Go to the source code of this file.
◆ solve()
template<typename M , typename X >
X chemmisol::gauss::solve |
( |
const M & |
m, |
|
|
const X & |
y |
|
) |
| |
Solves m*x=y using the Gaussian elimination algorithm.
- Parameters
-
m | Matrix of size n*n. |
y | Matrix of size n. |
- Returns
- Vector x such that m*x=y.
- Template Parameters
-
M | Matrix type. |
X | Vector type. |
◆ solve< double, double >()
Specialisation of the Gaussian elimination for scalar value.
The Gaussiant elimination is useless in this case, but this specialisation is necessary for the implementation of the Newton solver for scalar values.
- Parameters
-
- Returns
- Value x such that f*x=y.