diff options
author | torset | 2009-02-19 14:12:35 +0000 |
---|---|---|
committer | torset | 2009-02-19 14:12:35 +0000 |
commit | 6bec7e5598ad5cf78e56d82daf5a61f7d87839a8 (patch) | |
tree | d0951a78dd5be3730b9514c5b2d95386a2f133b3 /src/elementaryFunctions/pow/testDoublePow.c | |
parent | 31ebf3501c23dcd9d9d4b2452e91b6da5623060d (diff) | |
download | scilab2c-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.c | 18 |
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); |