26 if (i<n-1) A(i,i+1)=-1;
32 UmfpackSolver.
solve(U,F);
33 double residual=
normi(A*U-F);
35 std::cout <<
"residual:" << residual << std::endl;
Sparse matrix class using CRS storage scheme.
A::value_type normi(const A &a)
Maximum norm of array or expression.
void flush()
Re-create the internal data structure in order to accomodated all newly created elements.
void solve(TArray< T > &Sol, const TArray< T > &Rhs)
Solve LU factorized system.
Main header of the library.
One dimensional array class.
Bridge class for using umfpack as solver for vmatrix.
void update()
Perform actual computation of LU factorization.