Menu

[r136]: / trunk / include / _dssmatrix- / _dssmatrix-calc.hpp  Maximize  Restore  History

Download this file

47 lines (42 with data), 1.5 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
//=============================================================================
/*! return transposed _dssmatrix */
inline _dssmatrix t(const _dssmatrix& mat)
{VERBOSE_REPORT;
return mat;
}
///////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////
//=============================================================================
/*! search the index of element having the largest absolute value
in 0-based numbering system */
inline void idamax(long& i, long& j, const _dssmatrix& mat)
{VERBOSE_REPORT;
std::vector<dcomponent>::const_iterator itx(mat.data.begin());
double vmax(0);
for(std::vector<dcomponent>::const_iterator it=mat.data.begin(); it!=mat.data.end(); it++){
if(vmax < it->v){
vmax=fabs(it->v);
itx=it;
}
}
i=itx->i;
j=itx->j;
mat.destroy();
}
//=============================================================================
/*! return its largest absolute value */
inline double damax(const _dssmatrix& mat)
{VERBOSE_REPORT;
std::vector<dcomponent>::const_iterator itx(mat.data.begin());
double vmax(0);
for(std::vector<dcomponent>::const_iterator it=mat.data.begin(); it!=mat.data.end(); it++){
if(vmax < it->v){
vmax=fabs(it->v);
itx=it;
}
}
double val(itx->v);
mat.destroy();
return val;
}
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.