diff options
author | Alberto Ruiz <aruiz@um.es> | 2008-11-05 11:02:45 +0000 |
---|---|---|
committer | Alberto Ruiz <aruiz@um.es> | 2008-11-05 11:02:45 +0000 |
commit | 21e13ae0a13befb5cb8feb7c52bcd4b4e4cda953 (patch) | |
tree | dc41066b6fd545d76a6fa3114f1d5cf3473a642c /lib/Data/Packed/Internal/auxi.c | |
parent | 02805ad64715373347b34bac2f75cbb866563ba2 (diff) |
diag using ST
Diffstat (limited to 'lib/Data/Packed/Internal/auxi.c')
-rw-r--r-- | lib/Data/Packed/Internal/auxi.c | 26 |
1 files changed, 0 insertions, 26 deletions
diff --git a/lib/Data/Packed/Internal/auxi.c b/lib/Data/Packed/Internal/auxi.c index bbb8cd4..c449b9a 100644 --- a/lib/Data/Packed/Internal/auxi.c +++ b/lib/Data/Packed/Internal/auxi.c | |||
@@ -113,32 +113,6 @@ int constantC(gsl_complex* pval, CVEC(r)) { | |||
113 | } | 113 | } |
114 | 114 | ||
115 | 115 | ||
116 | int diagR(KRVEC(d),RMAT(r)) { | ||
117 | REQUIRES(dn==rr && rr==rc,BAD_SIZE); | ||
118 | DEBUGMSG("diagR"); | ||
119 | int i,j; | ||
120 | for (i=0;i<rr;i++) { | ||
121 | for(j=0;j<rc;j++) { | ||
122 | rp[i*rc+j] = i==j?dp[i]:0.; | ||
123 | } | ||
124 | } | ||
125 | OK | ||
126 | } | ||
127 | |||
128 | int diagC(KCVEC(d),CMAT(r)) { | ||
129 | REQUIRES(dn==rr && rr==rc,BAD_SIZE); | ||
130 | DEBUGMSG("diagC"); | ||
131 | int i,j; | ||
132 | gsl_complex zero; | ||
133 | GSL_SET_COMPLEX(&zero,0.,0.); | ||
134 | for (i=0;i<rr;i++) { | ||
135 | for(j=0;j<rc;j++) { | ||
136 | rp[i*rc+j] = i==j?dp[i]:zero; | ||
137 | } | ||
138 | } | ||
139 | OK | ||
140 | } | ||
141 | |||
142 | int conjugate(KCVEC(x),CVEC(t)) { | 116 | int conjugate(KCVEC(x),CVEC(t)) { |
143 | REQUIRES(xn==tn,BAD_SIZE); | 117 | REQUIRES(xn==tn,BAD_SIZE); |
144 | DEBUGMSG("conjugate"); | 118 | DEBUGMSG("conjugate"); |