Menu

[r37]: / trunk / include / dgbmatrix- / dgbmatrix-constructor.hpp  Maximize  Restore  History

Download this file

142 lines (120 with data), 4.6 kB

  1
  2
  3
  4
  5
  6
  7
  8
  9
 10
 11
 12
 13
 14
 15
 16
 17
 18
 19
 20
 21
 22
 23
 24
 25
 26
 27
 28
 29
 30
 31
 32
 33
 34
 35
 36
 37
 38
 39
 40
 41
 42
 43
 44
 45
 46
 47
 48
 49
 50
 51
 52
 53
 54
 55
 56
 57
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
//============================================================================
/*! dgbmatrix constructor */
inline dgbmatrix::dgbmatrix()
: m(M), n(N), kl(KL), ku(KU), array(Array), darray(Darray)
{
#ifdef CPPL_VERBOSE
std::cerr << "# [MARK] dgbmatrix::dgbmatrix()" << std::endl;
#endif//CPPL_VERBOSE
//////// initialize ////////
M =0;
N =0;
KL=0;
KU=0;
Array =NULL;
Darray =NULL;
#ifdef CPPL_VERBOSE
std::cerr << "# [NOTE] dgbmatrix::dgbmatrix()" << " A new 0x0:0,0 matrix at " << Array << " has been made." << std::endl;
#endif//CPPL_VERBOSE
}
//============================================================================
/*! dgbmatrix copy constructor */
inline dgbmatrix::dgbmatrix(const dgbmatrix& mat)
: m(M), n(N), kl(KL), ku(KU), array(Array), darray(Darray)
{
#ifdef CPPL_VERBOSE
std::cerr << "# [MARK] dgbmatrix::dgbmatrix(const dgbmatrix&)" << std::endl;
#endif//CPPL_VERBOSE
//////// initialize ////////
M =mat.M;
N =mat.N;
KL=mat.KL;
KU=mat.KU;
Array =new double[(KL+KU+1)*N];
Darray =new double*[N];
for(int i=0; i<N; i++){ Darray[i] =&Array[i*(KL+KU+1)]; }
//////// copy ////////
dcopy_((KL+KU+1)*N, mat.Array, 1, Array, 1);
#ifdef CPPL_VERBOSE
std::cerr << "# [NOTE] dgbmatrix::dgbmatrix(const dgbmatrix&) " << "A new matrix at " << Array << " has been made." << std::endl;
#endif//CPPL_VERBOSE
}
//============================================================================
/*! dgbmatrix constructor to cast _dgbmatrix */
inline dgbmatrix::dgbmatrix(const _dgbmatrix& mat)
: m(M), n(N), kl(KL), ku(KU), array(Array), darray(Darray)
{
#ifdef CPPL_VERBOSE
std::cerr << "# [MARK] dgbmatrix::dgbmatrix(const _dgbmatrix&)" << std::endl;
#endif//CPPL_VERBOSE
M =mat.M;
N =mat.N;
KL =mat.KL;
KU =mat.KU;
Array =mat.Array;
Darray =mat.Darray;
mat.nullify();
#ifdef CPPL_VERBOSE
std::cerr << "# [NOTE] dgbmatrix::dgbmatrix(const dgbmatrix&)" << " A new matrix at " << Array << " has been made." << std::endl;
#endif//CPPL_VERBOSE
}
//============================================================================
/*! dgbmatrix constructor with size specification */
inline dgbmatrix::dgbmatrix(const long& _m, const long& _n,
const long& _kl, const long& _ku)
: m(M), n(N), kl(KL), ku(KU), array(Array), darray(Darray)
{
#ifdef CPPL_VERBOSE
std::cerr << "# [MARK] dgbmatrix::dgbmatrix(const long&, const long&, const long&, const long&)" << std::endl;
#endif//CPPL_VERBOSE
#ifdef CPPL_DEBUG
if( _m<0 || _n<0 || _kl<0 || _ku<0 || _m<_kl || _n<_ku ){
std::cerr << "[ERROR] dgbmatrix::dgbmatrix(const long, const long, const long const long)" << std::endl
<< "It is impossible to make a matrix you ordered. " << std::endl
<< "Your input was (" << _m << "," << _n << ","<< _ku << "," << _kl << ")." << std::endl;
exit(1);
}
#endif//CPPL_DEBUG
//////// initialize ////////
M =_m;
N =_n;
KL =_kl;
KU =_ku;
Array =new double[(KL+KU+1)*N];
Darray =new double*[N];
for(int i=0; i<N; i++){ Darray[i] =&Array[i*(KL+KU+1)]; }
#ifdef CPPL_VERBOSE
std::cerr << "# [NOTE] dgbmatrix::dgbmatrix(const long&, const long&, const long&,const long&)" << " A new matrix at " << Array << " has been made." << std::endl;
#endif//CPPL_VERBOSE
}
//============================================================================
/*! dgbmatrix constructor with filename */
inline dgbmatrix::dgbmatrix(const char* filename)
: m(M), n(N), kl(KL), ku(KU), array(Array), darray(Darray)
{
#ifdef CPPL_VERBOSE
std::cerr << "# [MARK] dgbmatrix::dgbmatrix(const char*)" << std::endl;
#endif//CPPL_VERBOSE
Array =NULL;
Darray =NULL;
//// read ////
read(filename);
#ifdef CPPL_VERBOSE
std::cerr << "# [NOTE] dgbmatrix::dgbmatrix(const char*)" << " A new matrix at " << Array << " has been made." << std::endl;
#endif//CPPL_VERBOSE
}
///////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////
//============================================================================
/*! dgbmatrix destructor */
inline dgbmatrix::~dgbmatrix()
{
#ifdef CPPL_VERBOSE
std::cerr << "# [MARK] dgbmatrix::~dgbmatrix()" << std::endl;
std::cerr << "# [NOTE] dgbmatrix::~dgbmatrix()" << " A matrix at " << Array << " is going to be deleted." << std::endl;
#endif//CPPL_VERBOSE
//////// delete Array ////////
delete [] Array;
delete [] Darray;
}
Want the latest updates on software, tech news, and AI?
Get latest updates about software, tech news, and AI from SourceForge directly in your inbox once a month.