29 if (i<n-1) A(i,i+1)=-1;
36 double residual=
normi(A*U-F);
38 std::cout <<
"residual:" << residual << std::endl;
static std::shared_ptr< TSolverUMFPACK< T > > create(const std::shared_ptr< TSparseMatrix< T >> pA)
Create LU factorization class.
A::value_type normi(const A &a)
Maximum norm of array or expression.
Main header of the library.
static std::shared_ptr< TArray1< T > > create(index n1)
Construct smart pointer empty 1D Array.
static std::shared_ptr< TSparseMatrix< T > > create(index n1, index n2)
Static wrapper around corresponding constructor.