8 std::cerr <<
"# [MARK] zgbmatrix::zgbsv(zgematrix&)"
14 std::cerr <<
"[ERROR] zgbmatrix::zgbsv(zgematrix&) " << 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 zgbsv_(N, KL, KU, NRHS, newmat.Array, LDAB, IPIV, mat.Array, LDB, INFO);
35 std::cerr <<
"[WARNING] zgbmatrix::zgbsv(zgematrix&) "
36 <<
"Serious trouble happend. INFO = "<< INFO <<
"." << std::endl;
48 std::cerr <<
"# [MARK] zgbmatrix::zgbsv(zcovector&)"
54 std::cerr <<
"[ERROR] zgbmatrix::zgbsv(zcovector&) " << 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 zgbsv_(N, KL, KU, NRHS, newmat.Array, LDAB, IPIV, vec.Array, LDB, INFO);
75 std::cerr <<
"[WARNING] zgbmatrix::zgbsv(zcovector&) "
76 <<
"Serious trouble happend. INFO = "<< INFO <<
"." << std::endl;
std::complex< double > & operator()(const long &, const long &)
Definition: zgbmatrix-io.hpp:3
friend _zgematrix i(const zgbmatrix &)
Definition: zgbmatrix-calc.hpp:22
long zgbsv(zgematrix &)
Definition: zgbmatrix-lapack.hpp:5
Complex Double-precision General Dence Matrix Class.
Definition: zgematrix.hpp:3
Complex Double-precision General Band Matrix Class.
Definition: zgbmatrix.hpp:3
Complex Double-precision Column Vector Class.
Definition: zcovector.hpp:3
long const & n
matrix column size (readable)
Definition: zgbmatrix.hpp:17
friend void swap(zgbmatrix &, zgbmatrix &)
Definition: zgbmatrix-misc.hpp:167