summaryrefslogtreecommitdiff
path: root/packages/hmatrix/src/Numeric/GSL/gsl-vector.c
diff options
context:
space:
mode:
Diffstat (limited to 'packages/hmatrix/src/Numeric/GSL/gsl-vector.c')
-rw-r--r--packages/hmatrix/src/Numeric/GSL/gsl-vector.c82
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 }
135int 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
161int 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
188int mapCAux(int code, KGCVEC(x), GCVEC(r)) { 137int 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
299int 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
315int 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
331int mapValCAux(int code, gsl_complex* pval, KGCVEC(x), GCVEC(r)) { 249int mapValCAux(int code, gsl_complex* pval, KGCVEC(x), GCVEC(r)) {
332 int k; 250 int k;