diff options
Diffstat (limited to 'packages/hmatrix/src/Numeric/GSL/gsl-vector.c')
-rw-r--r-- | packages/hmatrix/src/Numeric/GSL/gsl-vector.c | 82 |
1 files changed, 0 insertions, 82 deletions
diff --git a/packages/hmatrix/src/Numeric/GSL/gsl-vector.c b/packages/hmatrix/src/Numeric/GSL/gsl-vector.c index 1eaf66c..99891d4 100644 --- a/packages/hmatrix/src/Numeric/GSL/gsl-vector.c +++ b/packages/hmatrix/src/Numeric/GSL/gsl-vector.c | |||
@@ -132,57 +132,6 @@ inline gsl_complex complex_signum(gsl_complex z) { | |||
132 | 132 | ||
133 | #define OP(C,F) case C: { for(k=0;k<xn;k++) rp[k] = F(xp[k]); OK } | 133 | #define OP(C,F) case C: { for(k=0;k<xn;k++) rp[k] = F(xp[k]); OK } |
134 | #define OPV(C,E) case C: { for(k=0;k<xn;k++) rp[k] = E; OK } | 134 | #define OPV(C,E) case C: { for(k=0;k<xn;k++) rp[k] = E; OK } |
135 | int mapR(int code, KRVEC(x), RVEC(r)) { | ||
136 | int k; | ||
137 | REQUIRES(xn == rn,BAD_SIZE); | ||
138 | DEBUGMSG("mapR"); | ||
139 | switch (code) { | ||
140 | OP(0,sin) | ||
141 | OP(1,cos) | ||
142 | OP(2,tan) | ||
143 | OP(3,fabs) | ||
144 | OP(4,asin) | ||
145 | OP(5,acos) | ||
146 | OP(6,atan) /* atan2 mediante vectorZip */ | ||
147 | OP(7,sinh) | ||
148 | OP(8,cosh) | ||
149 | OP(9,tanh) | ||
150 | OP(10,gsl_asinh) | ||
151 | OP(11,gsl_acosh) | ||
152 | OP(12,gsl_atanh) | ||
153 | OP(13,exp) | ||
154 | OP(14,log) | ||
155 | OP(15,sign) | ||
156 | OP(16,sqrt) | ||
157 | default: ERROR(BAD_CODE); | ||
158 | } | ||
159 | } | ||
160 | |||
161 | int mapF(int code, KFVEC(x), FVEC(r)) { | ||
162 | int k; | ||
163 | REQUIRES(xn == rn,BAD_SIZE); | ||
164 | DEBUGMSG("mapF"); | ||
165 | switch (code) { | ||
166 | OP(0,sin) | ||
167 | OP(1,cos) | ||
168 | OP(2,tan) | ||
169 | OP(3,fabs) | ||
170 | OP(4,asin) | ||
171 | OP(5,acos) | ||
172 | OP(6,atan) /* atan2 mediante vectorZip */ | ||
173 | OP(7,sinh) | ||
174 | OP(8,cosh) | ||
175 | OP(9,tanh) | ||
176 | OP(10,gsl_asinh) | ||
177 | OP(11,gsl_acosh) | ||
178 | OP(12,gsl_atanh) | ||
179 | OP(13,exp) | ||
180 | OP(14,log) | ||
181 | OP(15,sign) | ||
182 | OP(16,sqrt) | ||
183 | default: ERROR(BAD_CODE); | ||
184 | } | ||
185 | } | ||
186 | 135 | ||
187 | 136 | ||
188 | int mapCAux(int code, KGCVEC(x), GCVEC(r)) { | 137 | int mapCAux(int code, KGCVEC(x), GCVEC(r)) { |
@@ -296,37 +245,6 @@ int mapQ(int code, KQVEC(x), QVEC(r)) { | |||
296 | } | 245 | } |
297 | 246 | ||
298 | 247 | ||
299 | int mapValR(int code, double* pval, KRVEC(x), RVEC(r)) { | ||
300 | int k; | ||
301 | double val = *pval; | ||
302 | REQUIRES(xn == rn,BAD_SIZE); | ||
303 | DEBUGMSG("mapValR"); | ||
304 | switch (code) { | ||
305 | OPV(0,val*xp[k]) | ||
306 | OPV(1,val/xp[k]) | ||
307 | OPV(2,val+xp[k]) | ||
308 | OPV(3,val-xp[k]) | ||
309 | OPV(4,pow(val,xp[k])) | ||
310 | OPV(5,pow(xp[k],val)) | ||
311 | default: ERROR(BAD_CODE); | ||
312 | } | ||
313 | } | ||
314 | |||
315 | int mapValF(int code, float* pval, KFVEC(x), FVEC(r)) { | ||
316 | int k; | ||
317 | float val = *pval; | ||
318 | REQUIRES(xn == rn,BAD_SIZE); | ||
319 | DEBUGMSG("mapValF"); | ||
320 | switch (code) { | ||
321 | OPV(0,val*xp[k]) | ||
322 | OPV(1,val/xp[k]) | ||
323 | OPV(2,val+xp[k]) | ||
324 | OPV(3,val-xp[k]) | ||
325 | OPV(4,pow(val,xp[k])) | ||
326 | OPV(5,pow(xp[k],val)) | ||
327 | default: ERROR(BAD_CODE); | ||
328 | } | ||
329 | } | ||
330 | 248 | ||
331 | int mapValCAux(int code, gsl_complex* pval, KGCVEC(x), GCVEC(r)) { | 249 | int mapValCAux(int code, gsl_complex* pval, KGCVEC(x), GCVEC(r)) { |
332 | int k; | 250 | int k; |