summaryrefslogtreecommitdiff
path: root/src/elementaryFunctions/pow/testDoublePow.c
diff options
context:
space:
mode:
authortorset2009-02-19 14:12:35 +0000
committertorset2009-02-19 14:12:35 +0000
commit6bec7e5598ad5cf78e56d82daf5a61f7d87839a8 (patch)
treed0951a78dd5be3730b9514c5b2d95386a2f133b3 /src/elementaryFunctions/pow/testDoublePow.c
parent31ebf3501c23dcd9d9d4b2452e91b6da5623060d (diff)
downloadscilab2c-6bec7e5598ad5cf78e56d82daf5a61f7d87839a8.tar.gz
scilab2c-6bec7e5598ad5cf78e56d82daf5a61f7d87839a8.tar.bz2
scilab2c-6bec7e5598ad5cf78e56d82daf5a61f7d87839a8.zip
Modify arrays pow : now computes A.^B instead of A.^b, A,B=matrices,b=scalar
Diffstat (limited to 'src/elementaryFunctions/pow/testDoublePow.c')
-rw-r--r--src/elementaryFunctions/pow/testDoublePow.c18
1 files changed, 12 insertions, 6 deletions
diff --git a/src/elementaryFunctions/pow/testDoublePow.c b/src/elementaryFunctions/pow/testDoublePow.c
index 4149b300..57707f2a 100644
--- a/src/elementaryFunctions/pow/testDoublePow.c
+++ b/src/elementaryFunctions/pow/testDoublePow.c
@@ -463,14 +463,20 @@ void zpowsTest(void) {
void dpowaTest(void) {
double in1[]=SOURCE;
- double in2=EXPAND;
+ double in2[200];
+ double ZERO[200]={0};
+ double ONE[200];
double res[]=RESULT;
double out1[200],out2[200],out3[200];
int i;
-
+
+ for (i=0;i<200;i++) {
+ in2[i]=EXPAND;
+ ONE[i]=1;
+ }
dpowa(in1,200,in2,out1);
- dpowa(in1,200,0,out2);
- dpowa(in1,200,1,out3);
+ dpowa(in1,200,ZERO,out2);
+ dpowa(in1,200,ONE,out3);
for (i=0;i<200;i++){
assert(( (fabs(out1[i]-res[i]))/(fabs(out1[i])) )<1e-15);
assert(( (fabs(out2[i]-1))/(fabs(out2[i])) )<1e-16);
@@ -481,13 +487,13 @@ void dpowaTest(void) {
void zpowaTest(void) {
double in1R[]=ZSOURCER;
double in1I[]=ZSOURCEI;
- doubleComplex in2=ZEXPAND;
+ doubleComplex in2[200];
double resR[]=ZRESULTR;
double resI[]=ZRESULTI;
doubleComplex *in1,out[200];
int i;
-
+ for (i=0;i<200;i++) in2[i]=ZEXPAND;
in1=DoubleComplexMatrix(in1R,in1I,200);
zpowa(in1,200,in2,out);