4 : m(M), n(N), kl(KL), ku(KU), array(
Array), darray(Darray)
7 std::cerr <<
"# [MARK] dgbmatrix::dgbmatrix()"
20 std::cerr <<
"# [NOTE] dgbmatrix::dgbmatrix() "
21 <<
"A new 0x0:0,0 matrix at " <<
Array <<
" has been made."
29 : m(M), n(N), kl(KL), ku(KU), array(
Array), darray(Darray)
32 std::cerr <<
"# [MARK] dgbmatrix::dgbmatrix(const dgbmatrix&)"
41 Array =
new double[(KL+KU+1)*N];
42 Darray =
new double*[N];
43 for(
int i=0;
i<N;
i++){ Darray[
i] =&
Array[
i*(KL+KU+1)]; }
46 dcopy_((KL+KU+1)*N, mat.Array, 1,
Array, 1);
49 std::cerr <<
"# [NOTE] dgbmatrix::dgbmatrix(const dgbmatrix&) "
50 <<
"A new matrix at " <<
Array <<
" has been made." << std::endl;
57 : m(M), n(N), kl(KL), ku(KU), array(
Array), darray(Darray)
60 std::cerr <<
"# [MARK] dgbmatrix::dgbmatrix(const _dgbmatrix&)"
72 std::cerr <<
"# [NOTE] dgbmatrix::dgbmatrix(const dgbmatrix&) "
73 <<
"A new matrix at " <<
Array <<
" has been made." << std::endl;
80 const long& _kl,
const long& _ku)
81 : m(M), n(N), kl(KL), ku(KU), array(
Array), darray(Darray)
84 std::cerr <<
"# [MARK] dgbmatrix::dgbmatrix(const long&, const long&, const long&, const long&)"
89 if( _m<0 || _n<0 || _kl<0 || _ku<0 || _m<_kl || _n<_ku ){
90 std::cerr <<
"[ERROR] dgbmatrix::dgbmatrix"
91 <<
"(const long, const long, const long const long)"
93 <<
"It is impossible to make a matrix you ordered. "
96 << _m <<
"," << _n <<
","<< _ku <<
"," << _kl <<
")."
107 Array =
new double[(KL+KU+1)*N];
108 Darray =
new double*[N];
109 for(
int i=0;
i<N;
i++){ Darray[
i] =&
Array[
i*(KL+KU+1)]; }
112 std::cerr <<
"# [NOTE] dgbmatrix::dgbmatrix"
113 <<
"(const long&, const long&, const long&,const long&) "
114 <<
"A new matrix at " <<
Array <<
" has been made." << std::endl;
121 : m(M), n(N), kl(KL), ku(KU), array(
Array), darray(Darray)
124 std::cerr <<
"# [MARK] dgbmatrix::dgbmatrix(const char*)"
135 std::cerr <<
"# [NOTE] dgbmatrix::dgbmatrix(const char*) "
136 <<
"A new matrix at " <<
Array <<
" has been made." << std::endl;
149 std::cerr <<
"# [MARK] dgbmatrix::~dgbmatrix()"
154 std::cerr <<
"# [NOTE] dgbmatrix::~dgbmatrix() "
155 <<
"A matrix at " <<
Array
156 <<
" is going to be deleted." << std::endl;
friend _dgematrix i(const dgbmatrix &)
Definition: dgbmatrix-calc.hpp:22
long KL
lower band width
Definition: _dgbmatrix.hpp:9
long M
matrix row size
Definition: _dgbmatrix.hpp:7
long N
matrix column size
Definition: _dgbmatrix.hpp:8
double * Array
1D Array to store matrix data
Definition: _dgbmatrix.hpp:11
~dgbmatrix()
Definition: dgbmatrix-constructor.hpp:146
double ** Darray
array of pointers of column head addresses
Definition: _dgbmatrix.hpp:12
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
dgbmatrix()
Definition: dgbmatrix-constructor.hpp:3
double * Array
1D Array to store vector data
Definition: _drovector.hpp:8
long KU
upper band width
Definition: _dgbmatrix.hpp:10