summaryrefslogtreecommitdiff
path: root/packages/hmatrix/src/Numeric/GSL/gsl-vector.c
diff options
context:
space:
mode:
authorAlberto Ruiz <aruiz@um.es>2014-05-14 20:25:30 +0200
committerAlberto Ruiz <aruiz@um.es>2014-05-14 20:25:30 +0200
commitbdd718066eeb1dffc790d7f48398f24b9201c6ea (patch)
tree4a342f01ce372659c0e7757cecd821de157d411b /packages/hmatrix/src/Numeric/GSL/gsl-vector.c
parentce03cbdf416db2af2830b4b3b7d1722bb26f6546 (diff)
moved prod
Diffstat (limited to 'packages/hmatrix/src/Numeric/GSL/gsl-vector.c')
-rw-r--r--packages/hmatrix/src/Numeric/GSL/gsl-vector.c54
1 files changed, 0 insertions, 54 deletions
diff --git a/packages/hmatrix/src/Numeric/GSL/gsl-vector.c b/packages/hmatrix/src/Numeric/GSL/gsl-vector.c
index 3c2db4d..7e3c548 100644
--- a/packages/hmatrix/src/Numeric/GSL/gsl-vector.c
+++ b/packages/hmatrix/src/Numeric/GSL/gsl-vector.c
@@ -88,60 +88,6 @@
88#define BAD_FILE 2003 88#define BAD_FILE 2003
89 89
90 90
91int prodF(KFVEC(x),FVEC(r)) {
92 DEBUGMSG("prodF");
93 REQUIRES(rn==1,BAD_SIZE);
94 int i;
95 float res = 1;
96 for (i = 0; i < xn; i++) res *= xp[i];
97 rp[0] = res;
98 OK
99}
100
101int prodR(KRVEC(x),RVEC(r)) {
102 DEBUGMSG("prodR");
103 REQUIRES(rn==1,BAD_SIZE);
104 int i;
105 double res = 1;
106 for (i = 0; i < xn; i++) res *= xp[i];
107 rp[0] = res;
108 OK
109}
110
111int prodQ(KQVEC(x),QVEC(r)) {
112 DEBUGMSG("prodQ");
113 REQUIRES(rn==1,BAD_SIZE);
114 int i;
115 gsl_complex_float res;
116 float temp;
117 res.dat[0] = 1;
118 res.dat[1] = 0;
119 for (i = 0; i < xn; i++) {
120 temp = res.dat[0] * xp[i].dat[0] - res.dat[1] * xp[i].dat[1];
121 res.dat[1] = res.dat[0] * xp[i].dat[1] + res.dat[1] * xp[i].dat[0];
122 res.dat[0] = temp;
123 }
124 rp[0] = res;
125 OK
126}
127
128int prodC(KCVEC(x),CVEC(r)) {
129 DEBUGMSG("prodC");
130 REQUIRES(rn==1,BAD_SIZE);
131 int i;
132 gsl_complex res;
133 double temp;
134 res.dat[0] = 1;
135 res.dat[1] = 0;
136 for (i = 0; i < xn; i++) {
137 temp = res.dat[0] * xp[i].dat[0] - res.dat[1] * xp[i].dat[1];
138 res.dat[1] = res.dat[0] * xp[i].dat[1] + res.dat[1] * xp[i].dat[0];
139 res.dat[0] = temp;
140 }
141 rp[0] = res;
142 OK
143}
144
145 91
146int toScalarR(int code, KRVEC(x), RVEC(r)) { 92int toScalarR(int code, KRVEC(x), RVEC(r)) {
147 REQUIRES(rn==1,BAD_SIZE); 93 REQUIRES(rn==1,BAD_SIZE);