6 std::cerr <<
"# [MARK] t(const dgbmatrix&)"
10 dgbmatrix newmat(mat.N, mat.M, mat.KU, mat.KL);
11 for(
long i=0; i<newmat.M; i++){
12 for(
long j=
max(0,i-newmat.KL); j<min(newmat.N,i+newmat.KU+1); j++){
13 newmat(i,j) =mat(j,i);
25 std::cerr <<
"# [MARK] i(const dgbmatrix&)"
31 std::cerr <<
"[ERROR] i(dgbmatrix&) " << std::endl
32 <<
"This matrix is not square and has no inverse matrix."
34 <<
"Your input was (" << mat.M <<
"x" << mat.N <<
")."
43 mat_cp.
dgbsv(mat_inv);
double max(double v1, double v2)
Return maximum.
Definition: variousFunctions.cpp:355
Real Double-precision General Dence Matrix Class.
Definition: dgematrix.hpp:3
(DO NOT USE) Smart-temporary Real Double-precision General Dence Matrix Class
Definition: _dgematrix.hpp:3
void identity()
Definition: dgematrix-misc.hpp:38
friend _dcovector t(const drovector &)
Definition: drovector-calc.hpp:3
Real Double-precision General Band Matrix Class.
Definition: dgbmatrix.hpp:3
(DO NOT USE) Smart-temporary Real Double-precision General Band Matrix Class
Definition: _dgbmatrix.hpp:3
long dgbsv(dgematrix &)
Definition: dgbmatrix-lapack.hpp:5