6 std::cerr <<
"# [MARK] dgematrix::clear()"
11 std::cerr <<
"# [NOTE] dgematrix::clear() "
12 <<
" An array at " <<
Array
13 <<
" is going to be cleared." << std::endl;
29 std::cerr <<
"# [MARK] dgematrix::zero()"
33 for(
long i=0;
i<M*N;
i++){
Array[
i] =0.0; }
41 std::cerr <<
"# [MARK] dgematrix::identity()"
47 std::cerr <<
"[ERROR] dgematrix::identity()" << std::endl
48 <<
"Only square matrix can be a identity matrix." << std::endl
49 <<
"The matrix size was " << M <<
"x" << N <<
"." << std::endl;
54 for(
long i=0;
i<M*N;
i++){
Array[
i] =0.0; }
63 std::cerr <<
"# [MARK] dgematrix::chsign()"
75 std::cerr <<
"# [MARK] dgematrix::copy(const dgematrix&)"
80 std::cerr <<
"# [NOTE] dgematrix::copy(const dgematrix&) "
81 <<
"A dgematrix at " <<
Array <<
" is going to be deleted.";
87 Array =
new double[mat.M*mat.N];
89 Darray =
new double*[N];
90 for(
int i=0;
i<N;
i++){ Darray[
i] =&Array[
i*M]; }
92 dcopy_(mat.M*mat.N, mat.Array, 1, Array, 1);
95 std::cerr <<
" Then, a COPY of a dgematrix has been cleated at "
96 << Array <<
"." << std::endl;
106 std::cerr <<
"# [MARK] dgematrix::shallow_copy(const _dgematrix&)"
111 std::cerr <<
"# [NOTE] dgematrix:shallow_copy(const _dgematrix&) "
112 <<
"A dgematrix at " << Array <<
" is going to be deleted, "
113 <<
"and point at " << mat.
Array <<
" instead." << std::endl;
129 std::cerr <<
"# [MARK] dgematrix::resize(const long&, const long&)"
135 std::cerr <<
"[ERROR] dgematrix::resize(const long&, const long&)"
137 <<
"Matrix sizes must be positive integers." << std::endl
138 <<
"Your input was (" << _m <<
"," << _n <<
")." << std::endl;
146 Array =
new double[M*N];
148 Darray =
new double*[N];
149 for(
int i=0;
i<N;
i++){ Darray[
i] =&Array[
i*M]; }
157 std::cerr <<
"# [MARK] swap(dgematrix&, dgematrix&)"
161 long A_m(A.M), A_n(A.N);
162 double* A_array(A.Array);
163 double** A_darray = A.Darray;
179 std::cerr <<
"# [MARK] _(dgematrix&)"
187 newmat.
Array =mat.Array;
188 newmat.
Darray =mat.Darray;
void clear()
Definition: dgematrix-misc.hpp:3
void zero()
Definition: dgematrix-misc.hpp:26
void resize(const long &, const long &)
Definition: dgematrix-misc.hpp:126
double * Array
1D Array to store vector data
Definition: _drovector.hpp:8
double ** Darray
array of pointers of column head addresses
Definition: _dgematrix.hpp:10
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
long N
matrix column size
Definition: _dgematrix.hpp:8
void copy(const dgematrix &)
Definition: dgematrix-misc.hpp:72
void shallow_copy(const _dgematrix &)
Definition: dgematrix-misc.hpp:103
double * Array
1D array to store matrix data
Definition: _dgematrix.hpp:9
void chsign()
Definition: dgematrix-misc.hpp:60
friend _dgematrix i(const dgematrix &)
Definition: dgematrix-calc.hpp:21
long M
matrix row size
Definition: _dgematrix.hpp:7
double & operator()(const long &, const long &)
Definition: dgematrix-io.hpp:3