30 if (i<n-1) A(i,i+1)=-1;
35 double residual=
normi(A*U-F);
37 std::cout <<
"residual:" << residual << std::endl;
A::value_type normi(const A &a)
Maximum norm of array or expression.
static std::shared_ptr< TSolverLapackLU< T > > create(const std::shared_ptr< TMatrix< T >> pMatrix)
Static wrapper around constructor.
static std::shared_ptr< TMatrix< T > > create(index n1, index n2)
Construct empty square matrix.
Main header of the library.
static std::shared_ptr< TArray1< T > > create(index n1)
Construct smart pointer empty 1D Array.