Menu

[r31]: / trunk / include / _dgsmatrix- / _dgsmatrix-calc.hpp  Maximize  Restore  History

Download this file

55 lines (45 with data), 1.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
//=============================================================================
/*! return transposed _dgsmatrix */
inline _dgsmatrix t(const _dgsmatrix& mat)
{
#ifdef CPPL_VERBOSE
std::cerr << "# [MARK] t(const _dgsmatrix&)"
<< std::endl;
#endif//CPPL_VERBOSE
long* mat_Indx(mat.Indx);
mat.Indx =mat.Jndx;
mat.Jndx =mat_Indx;
std::vector< std::pair<long,long> >* mat_Row(mat.Row);
mat.Row =mat.Col;
mat.Col =mat_Row;
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 _dgsmatrix& mat)
{
#ifdef CPPL_VERBOSE
std::cerr << "# [MARK] idamax(long&, long&, const _dgsmatrix&)"
<< std::endl;
#endif//CPPL_VERBOSE
long index( idamax_(mat.VOL, mat.Array, 1) -1 );
i =mat.Indx[index];
j =mat.Jndx[index];
mat.destroy();
}
//=============================================================================
/*! return its largest absolute value */
inline double damax(const _dgsmatrix& mat)
{
#ifdef CPPL_VERBOSE
std::cerr << "# [MARK] damax(const _dgsmatrix&)"
<< std::endl;
#endif//CPPL_VERBOSE
double val( mat.Array[idamax_(mat.VOL, mat.Array, 1) -1] );
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.