Menu

[r108]: / trunk / include / _dssmatrix- / _dssmatrix-cast.hpp  Maximize  Restore  History

Download this file

98 lines (78 with data), 3.1 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
//============================================================================
/*! cast to _zhsmatrix */
inline _zhsmatrix _dssmatrix::to_zhsmatrix() const
{
#ifdef CPPL_VERBOSE
std::cerr << "# [MARK] _dssmatrix::to_zhsmatrix() const" << std::endl;
#endif//CPPL_VERBOSE
zhsmatrix newmat(N,Data.size());
for(std::vector<dcomponent>::const_iterator it=Data.begin(); it!=Data.end(); it++){
newmat.fput(it->i, it->j, std::complex<double>(it->v,0.0));
}
#ifdef CPPL_DEBUG
std::cerr << "# [NOTE] _dssmatrix::operator _zhsmatrix() const" << " A new casted matrix at " << newmat.array << " has been made." << std::endl;
#endif//CPPL_DEBUG
destroy();
return _(newmat);
}
//////////////////////////////////////////////////////////////////////////////
//////////////////////////////////////////////////////////////////////////////
//////////////////////////////////////////////////////////////////////////////
//============================================================================
/*! convert to _dgematrix */
inline _dgematrix _dssmatrix::to_dgematrix() const
{
#ifdef CPPL_VERBOSE
std::cerr << "# [MARK] _dssmatrix::to_dgematrix() const" << std::endl;
#endif//CPPL_VERBOSE
dgematrix newmat(M,N);
newmat.zero();
for(std::vector<dcomponent>::const_iterator it=Data.begin(); it!=Data.end(); it++){
newmat(it->i, it->j) =it->v;
newmat(it->j, it->i) =it->v;
}
#ifdef CPPL_DEBUG
std::cerr << "# [NOTE] _dssmatrix::to_dgematrix() const" << " A new casted matrix at " << newmat.array << " has been made." << std::endl;
#endif//CPPL_DEBUG
destroy();
return _(newmat);
}
//============================================================================
/*! convert to _dsymatrix */
inline _dsymatrix _dssmatrix::to_dsymatrix() const
{
#ifdef CPPL_VERBOSE
std::cerr << "# [MARK] _dssmatrix::to_dsymatrix() const" << std::endl;
#endif//CPPL_VERBOSE
dsymatrix newmat(N);
newmat.zero();
for(std::vector<dcomponent>::const_iterator it=Data.begin(); it!=Data.end(); it++){
newmat(it->i, it->j) =it->v;
}
#ifdef CPPL_DEBUG
std::cerr << "# [NOTE] _dssmatrix::to_dsymatrix() const" << " A new casted matrix at " << newmat.array << " has been made." << std::endl;
#endif//CPPL_DEBUG
destroy();
return _(newmat);
}
//============================================================================
/*! convert to _dgsmatrix */
inline _dgsmatrix _dssmatrix::to_dgsmatrix() const
{
#ifdef CPPL_VERBOSE
std::cerr << "# [MARK] _dssmatrix::to_dgsmatrix() const" << std::endl;
#endif//CPPL_VERBOSE
dgsmatrix newmat(M,N,Data.size()*2);
newmat.zero();
for(std::vector<dcomponent>::const_iterator it=Data.begin(); it!=Data.end(); it++){
newmat.fput(it->i, it->j, it->v);
if(it->i!=it->j){
newmat.fput(it->j, it->i, it->v);
}
}
#ifdef CPPL_VERBOSE
std::cerr << "# [NOTE] _dssmatrix::to_dgsmatrix() const" << " A new casted matrix at " << newmat.array << " has been made." << std::endl;
#endif//CPPL_VERBOSE
destroy();
return _(newmat);
}
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.