8 std::cerr <<
"# [MARK] dgbmatrix::dgbsv(dgematrix&)"
14 std::cerr <<
"[ERROR] dgbmatrix::dgbsv(dgematrix&) " << std::endl
15 <<
"These matrix and vector cannot be solved." << std::endl
16 <<
"Your input was (" << M <<
"x" << N <<
") and ("
17 << mat.M <<
"x" << mat.N <<
")." << std::endl;
23 for(
long i=0;
i<M;
i++){
for(
long j=
max(0,
i-KL); j<min(N,
i+KU+1); j++){
27 long NRHS(mat.N), LDAB(2*KL+KU+1),
28 *IPIV(
new long[N]), LDB(mat.M), INFO(1);
29 dgbsv_(N, KL, KU, NRHS, newmat.Array, LDAB, IPIV, mat.Array, LDB, INFO);
35 std::cerr <<
"[WARNING] dgbmatrix::dgbsv(dgematrix&) "
36 <<
"Serious trouble happend. INFO = "<< INFO <<
"." << std::endl;
48 std::cerr <<
"# [MARK] dgbmatrix::dgbsv(dcovector&)"
54 std::cerr <<
"[ERROR] dgbmatrix::dgbsv(dcovector&) " << std::endl
55 <<
"These matrix and vector cannot be solved." << std::endl
56 <<
"Your input was (" << M <<
"x" << N <<
") and ("
57 << vec.L <<
")." << std::endl;
63 for(
long i=0;
i<M;
i++){
for(
long j=
max(0,
i-KL); j<min(N,
i+KU+1); j++){
67 long NRHS(1), LDAB(2*KL+KU+1),
68 *IPIV(
new long[N]), LDB(vec.L), INFO(1);
69 dgbsv_(N, KL, KU, NRHS, newmat.Array, LDAB, IPIV, vec.Array, LDB, INFO);
75 std::cerr <<
"[WARNING] dgbmatrix::dgbsv(dgematrix&) "
76 <<
"Serious trouble happend. INFO = "<< INFO <<
"." << std::endl;
friend _dgematrix i(const dgbmatrix &)
Definition: dgbmatrix-calc.hpp:22
double max(double v1, double v2)
Return maximum.
Definition: variousFunctions.cpp:355
double & operator()(const long &, const long &)
Definition: dgbmatrix-io.hpp:3
Real Double-precision General Dence Matrix Class.
Definition: dgematrix.hpp:3
Real Double-precision General Band Matrix Class.
Definition: dgbmatrix.hpp:3
Real Double-precision Column Vector Class.
Definition: dcovector.hpp:3
long dgbsv(dgematrix &)
Definition: dgbmatrix-lapack.hpp:5
friend void swap(dgbmatrix &, dgbmatrix &)
Definition: dgbmatrix-misc.hpp:167