Menu

[r137]: / trunk / include / dgematrix- / dgematrix-constructor.hpp  Maximize  Restore  History

Download this file

90 lines (77 with data), 2.7 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
//=============================================================================
/*! dgematrix constructor without arguments */
inline dgematrix::dgematrix()
{VERBOSE_REPORT;
//////// initialize ////////
m =0;
n =0;
array =NULL;
darray =NULL;
}
///////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////
//=============================================================================
/*! dgematrix copy constructor */
inline dgematrix::dgematrix(const dgematrix& mat)
{VERBOSE_REPORT;
//////// initialize ////////
m =mat.m;
n =mat.n;
array =new double[m*n];
darray =new double*[n];
for(int i=0; i<n; i++){ darray[i] =&array[i*m]; }
//////// copy ////////
dcopy_(m*n, mat.array, 1, array, 1);
}
//=============================================================================
/*! dgematrix constructor to cast _dgematrix */
inline dgematrix::dgematrix(const _dgematrix& mat)
{VERBOSE_REPORT;
m =mat.m;
n =mat.n;
array =mat.array;
darray =mat.darray;
mat.nullify();
}
///////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////
//=============================================================================
/*! dgematrix constructor with size specification */
inline dgematrix::dgematrix(const long& _m, const long& _n)
{VERBOSE_REPORT;
#ifdef CPPL_DEBUG
if( _m<0 || _n<0 ){
ERROR_REPORT;
std::cerr << "Matrix sizes must be positive integers. " << std::endl
<< "Your input was (" << _m << "," << _n << ")." << std::endl;
exit(1);
}
#endif//CPPL_DEBUG
//////// initialize ////////
m =_m;
n =_n;
array =new double[m*n];
darray =new double*[n];
for(int i=0; i<n; i++){ darray[i] =&array[i*m]; }
}
//=============================================================================
/*! dgematrix constructor with filename */
inline dgematrix::dgematrix(const char* filename)
{VERBOSE_REPORT;
array =NULL;
darray =NULL;
//// read ////
read(filename);
}
///////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////
//=============================================================================
/*! dgematrix destructor */
inline dgematrix::~dgematrix()
{VERBOSE_REPORT;
delete [] darray;
delete [] array;
}
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.