24 if (i<n-1) A(i,i+1)=-1;
28 LapackSolver.
solve(U,F);
29 double residual=
normi(A*U-F);
31 std::cout <<
"residual:" << residual << std::endl;
A::value_type normi(const A &a)
Maximum norm of array or expression.
void solve(TArray< T > &Sol, const TArray< T > &Rhs) const
Solve LU factorized system.
Main header of the library.
One dimensional array class.
Lapack LU factorization class.