4 : m(M), n(N), array(
Array), darray(Darray)
7 std::cerr <<
"# [MARK] dgematrix::dgematrix()"
18 std::cerr <<
"# [NOTE] dgematrix::dgematrix() "
19 <<
"A new 0x0 matrix at " <<
Array
20 <<
" has been made." << std::endl;
31 : m(M), n(N), array(
Array), darray(Darray)
34 std::cerr <<
"# [MARK] dgematrix::dgematrix(const dgematrix&)"
41 Array =
new double[M*N];
42 Darray =
new double*[N];
43 for(
int i=0;
i<N;
i++){ Darray[
i] =&
Array[
i*M]; }
46 dcopy_(M*N, mat.Array, 1,
Array, 1);
49 std::cerr <<
"# [NOTE] dgematrix::dgematrix(const dgematrix&) "
50 <<
"A new matrix at " <<
Array <<
" has been made." << std::endl;
57 : m(M), n(N), array(
Array), darray(Darray)
60 std::cerr <<
"# [MARK] dgematrix::dgematrix(const _dgematrix&)"
70 std::cerr <<
"# [NOTE] dgematrix::dgematrix(const _dgematrix&) "
71 <<
"A new matrix pointing at " <<
Array <<
" has been made."
83 : m(M), n(N), array(
Array), darray(Darray)
86 std::cerr <<
"# [MARK] dgematrix::dgematrix(const dgbmatrix&)"
93 Array =
new double[M*N];
94 Darray =
new double*[N];
95 for(
int i=0;
i<N;
i++){ Darray[
i] =&
Array[
i*M]; }
99 for(
long i=0;
i<mat.M;
i++){
100 for(
long j=max(0,
i-mat.KL); j<min(N,
i+mat.KU+1); j++){
105 std::cerr <<
"# [NOTE] dgematrix::dgematrix(const dgbmatrix&) "
106 <<
"A new matrix at " <<
Array <<
" has been made." << std::endl;
113 : m(M), n(N), array(
Array), darray(Darray)
116 std::cerr <<
"# [MARK] dgematrix::dgematrix(const _dgbmatrix&)"
123 Array =
new double[M*N];
124 Darray =
new double*[N];
125 for(
int i=0;
i<N;
i++){ Darray[
i] =&
Array[
i*M]; }
129 for(
long i=0;
i<mat.
M;
i++){
130 for(
long j=max(0,
i-mat.
KL); j<min(N,
i+mat.
KU+1); j++){
138 std::cerr <<
"# [NOTE] dgematrix::dgematrix(const dgbmatrix&) "
139 <<
"A new matrix at " <<
Array <<
" has been made." << std::endl;
150 : m(M), n(N), array(
Array), darray(Darray)
153 std::cerr <<
"# [MARK] dgematrix::dgematrix(const dsymatrix&)"
160 Array =
new double[N*N];
161 Darray =
new double*[N];
162 for(
int i=0;
i<N;
i++){ Darray[
i] =&
Array[
i*M]; }
166 for(
long i=0;
i<mat.N;
i++) {
167 for(
long j=0; j<mat.N; j++ ) {
172 std::cerr <<
"# [NOTE] dgematrix::dgematrix(const dsymatrix&) "
173 <<
"A new matrix at " <<
Array <<
" has been made." << std::endl;
180 : m(M), n(N), array(
Array), darray(Darray)
183 std::cerr <<
"# [MARK] dgematrix::dgematrix(const _dsymatrix&)"
196 for(
long i=0;
i<mat.
N;
i++){
for(
long j=0; j<
i; j++){
201 std::cerr <<
"# [NOTE] dgematrix::dgematrix(const _dsymatrix&) "
202 <<
"A matrix at " <<
Array <<
" has been transformed into "
203 <<
"dgematrix." << std::endl;
214 : m(M), n(N), array(
Array), darray(Darray)
217 std::cerr <<
"# [MARK] dgematrix::dgematrix(const dssmatrix&)"
224 Array =
new double[M*N];
225 Darray =
new double*[N];
226 for(
int i=0;
i<N;
i++){ Darray[
i] =&
Array[
i*M]; }
230 for(
long c=0; c<mat.VOL; c++) {
231 operator()(mat.Indx[c],mat.Jndx[c]) =mat.Array[c];
235 std::cerr <<
"# [NOTE] dgematrix::dgematrix(const dssmatrix&) "
236 <<
"A new matrix at " <<
Array <<
" has been made." << std::endl;
243 : m(M), n(N), array(
Array), darray(Darray)
246 std::cerr <<
"# [MARK] dgematrix::dgematrix(const _dssmatrix&)"
253 Array =
new double[M*N];
254 Darray =
new double*[N];
255 for(
int i=0;
i<N;
i++){ Darray[
i] =&
Array[
i*M]; }
259 for(
long c=0; c<mat.
VOL; c++){
266 std::cerr <<
"# [NOTE] dgematrix::dgematrix(const _dssmatrix&) "
267 <<
"A new matrix at " <<
Array <<
" has been made." << std::endl;
278 : m(M), n(N), array(
Array), darray(Darray)
281 std::cerr <<
"# [MARK] dgematrix::dgematrix(const long&, const long&)"
287 std::cerr <<
"[ERROR] dgematrix::dgematrix(const long, const long)"
289 <<
"Matrix sizes must be positive integers. " << std::endl
290 <<
"Your input was (" << _m <<
"," << _n <<
")." << std::endl;
298 Array =
new double[M*N];
299 Darray =
new double*[N];
300 for(
int i=0;
i<N;
i++){ Darray[
i] =&
Array[
i*M]; }
303 std::cerr <<
"# [NOTE] dgematrix(long,long) "
304 <<
"A new matrix at " <<
Array <<
" has been made." << std::endl;
311 : m(M), n(N), array(
Array), darray(Darray)
314 std::cerr <<
"# [MARK] dgematrix::dgematrix(const char*)"
325 std::cerr <<
"# [NOTE] dgematrix::dgematrix(const char*) "
326 <<
"A new matrix at " <<
Array <<
" has been made." << std::endl;
339 std::cerr <<
"# [MARK] dgematrix::~dgematrix()"
344 std::cerr <<
"# [NOTE] dgematrix::~dgematrix() "
345 <<
"A matrix at " <<
Array <<
" is going to be deleted." << std::endl;
double * Array
1D array to store non-zero matrix data
Definition: _dssmatrix.hpp:11
~dgematrix()
Definition: dgematrix-constructor.hpp:336
void zero()
Definition: dgematrix-misc.hpp:26
long KL
lower band width
Definition: _dgbmatrix.hpp:9
(DO NOT USE) Smart-temporary Real Double-precision Sparse Matrix Class
Definition: _dssmatrix.hpp:3
long M
matrix row size
Definition: _dgbmatrix.hpp:7
void destroy() const
Definition: _dgbmatrix-misc.hpp:3
double * Array
1D Array to store vector data
Definition: _drovector.hpp:8
long M
matrix row size
Definition: _dssmatrix.hpp:7
double ** Darray
array of pointers of column head addresses
Definition: _dgematrix.hpp:10
Real Double-precision General Dence Matrix Class.
Definition: dgematrix.hpp:3
long N
matrix column size
Definition: _dgbmatrix.hpp:8
double * Array
1D Array to store matrix data
Definition: _dsymatrix.hpp:8
Real Double-precision Symmetric Matrix Class [L-type (UPLO=L) Strage].
Definition: dsymatrix.hpp:3
double ** Darray
array of pointers of column head addresses
Definition: _dsymatrix.hpp:9
long * Jndx
1D array to store the j-index of non-zero matrix components
Definition: _dssmatrix.hpp:13
(DO NOT USE) Smart-temporary Real Double-precision General Dence Matrix Class
Definition: _dgematrix.hpp:3
long N
matrix column size
Definition: _dgematrix.hpp:8
Real Double-precision General Band Matrix Class.
Definition: dgbmatrix.hpp:3
long VOL
the number of non-zero components
Definition: _dssmatrix.hpp:10
double * Array
1D array to store matrix data
Definition: _dgematrix.hpp:9
long N
matrix column or row size
Definition: _dsymatrix.hpp:7
(DO NOT USE) Smart-temporary Real Double-precision General Band Matrix Class
Definition: _dgbmatrix.hpp:3
long N
matrix column size
Definition: _dssmatrix.hpp:8
(DO NOT USE) Smart-temporary Real Double-precision Symmetric Matrix Class
Definition: _dsymatrix.hpp:3
Real Double-precision Sparse Matrix Class.
Definition: dssmatrix.hpp:3
friend _dgematrix i(const dgematrix &)
Definition: dgematrix-calc.hpp:21
long * Indx
1D array to store the i-index of non-zero matrix components
Definition: _dssmatrix.hpp:12
void destroy() const
Definition: _dssmatrix-misc.hpp:3
long KU
upper band width
Definition: _dgbmatrix.hpp:10
dgematrix()
Definition: dgematrix-constructor.hpp:3
long M
matrix row size
Definition: _dgematrix.hpp:7
double & operator()(const long &, const long &)
Definition: dgematrix-io.hpp:3