diff options
Diffstat (limited to 'src/c/signalProcessing/ifft/testMatIfft.c')
| -rw-r--r-- | src/c/signalProcessing/ifft/testMatIfft.c | 47 |
1 files changed, 19 insertions, 28 deletions
diff --git a/src/c/signalProcessing/ifft/testMatIfft.c b/src/c/signalProcessing/ifft/testMatIfft.c index 2d43267f..64048a8d 100644 --- a/src/c/signalProcessing/ifft/testMatIfft.c +++ b/src/c/signalProcessing/ifft/testMatIfft.c @@ -1,7 +1,7 @@ /* * Scilab ( http://www.scilab.org/ ) - This file is part of Scilab - * Copyright (C) 2008-2008 - INRIA - Allan SIMON + * Copyright (C) 2008-2008 - INRIA - Arnaud TORSET * * This file must be used under the terms of the CeCILL. * This source file is licensed as described in the file COPYING, which @@ -130,14 +130,14 @@ static void difftmaTest(void){ printf(" >>> Matrice 1*12 <<< \n"); difftma(in1, 1, 12, out1); for (i=0;i<12;i++){ - if (out1[i]>1e-16) assert( (fabs(out1[i]-resR1[i]) / fabs(out1[i])) < 3e-16 ); + if (out1[i]>1e-16) assert( (fabs(out1[i]-resR1[i]) / fabs(out1[i])) < 3e-14 ); else assert(1); } printf(" >>> Matrice 2*6 <<< \n"); difftma(in2, 2, 6, out2); for (i=0;i<12;i++){ - if (out2[i]>1e-16) assert( (fabs(out2[i]-resR2[i]) / fabs(out2[i])) < 3e-16 ); + if (out2[i]>1e-16) assert( (fabs(out2[i]-resR2[i]) / fabs(out2[i])) < 3e-15 ); else assert(1); } @@ -145,28 +145,28 @@ static void difftmaTest(void){ printf(" >>> Matrice 3*4 <<< \n"); difftma(in3, 3, 4, out3); for (i=0;i<12;i++){ - if (out3[i]>1e-16) assert( (fabs(out3[i]-resR3[i]) / fabs(out3[i])) < 3e-16 ); + if (out3[i]>1e-16) assert( (fabs(out3[i]-resR3[i]) / fabs(out3[i])) < 3e-15 ); else assert(1); } printf(" >>> Matrice 4*3 <<< \n"); difftma(in4, 4, 3, out4); for (i=0;i<12;i++){ - if (out4[i]>1e-16) assert( (fabs(out4[i]-resR4[i]) / fabs(out4[i])) < 3e-16 ); + if (out4[i]>1e-16) assert( (fabs(out4[i]-resR4[i]) / fabs(out4[i])) < 3e-15 ); else assert(1); } printf(" >>> Matrice 6*2 <<< \n"); difftma(in6, 6, 2, out6); for (i=0;i<12;i++){ - if (out6[i]>1e-16) assert( (fabs(out6[i]-resR6[i]) / fabs(out6[i])) < 3e-16 ); + if (out6[i]>1e-16) assert( (fabs(out6[i]-resR6[i]) / fabs(out6[i])) < 3e-15 ); else assert(1); } printf(" >>> Matrice 3*3 <<< \n"); difftma(in9, 3, 3, out9); for (i=0;i<9;i++){ - if (out9[i]>1e-16) assert( (fabs(out9[i]-resR9[i]) / fabs(out9[i])) < 3e-16 ); + if (out9[i]>1e-16) assert( (fabs(out9[i]-resR9[i]) / fabs(out9[i])) < 3e-15 ); else assert(1); } @@ -193,10 +193,9 @@ static void zifftmaTest(void){ double resI4[]=IRESULT4; double resR6[]=RRESULT6; double resI6[]=IRESULT6; - double resR9[]=RRESULT9; - double resI9[]=IRESULT9; + - doubleComplex *in1, *in2, *in3, *in4, *in6, *in9, out1[12], out2[12], out3[12], out4[12], out6[12], out9[12]; + doubleComplex *in1, *in2, *in3, *in4, *in6, *in9, out1[12], out2[12], out3[12], out4[12], out6[12]; in1=(doubleComplex *)malloc((unsigned int)12*sizeof(doubleComplex)); in2=(doubleComplex *)malloc((unsigned int)12*sizeof(doubleComplex)); @@ -227,16 +226,16 @@ static void zifftmaTest(void){ printf(" >>> Matrice 1*12 <<< \n"); zifftma(in1, 1, 12, out1); for (i=0;i<12;i++){ - if (zreals(out1[i])>1e-16) assert( (fabs(zreals(out1[i])-resR1[i]) / fabs(zreals(out1[i]))) < 3e-16 ); + if (zreals(out1[i])>1e-16) assert( (fabs(zreals(out1[i])-resR1[i]) / fabs(zreals(out1[i]))) < 3e-14 ); else assert(1); - if (zimags(out1[i])>1e-16) assert( (fabs(zimags(out1[i])+resI1[i]) / fabs(zimags(out1[i]))) < 3e-16 ); + if (zimags(out1[i])>1e-16) assert( (fabs(zimags(out1[i])+resI1[i]) / fabs(zimags(out1[i]))) < 3e-15 ); else assert(1); } printf(" >>> Matrice 2*6 <<< \n"); zifftma(in2, 2, 6, out2); for (i=0;i<12;i++){ - if (zreals(out2[i])>1e-16) assert( (fabs(zreals(out2[i])-resR2[i]) / fabs(zreals(out2[i]))) < 3e-16 ); + if (zreals(out2[i])>1e-16) assert( (fabs(zreals(out2[i])-resR2[i]) / fabs(zreals(out2[i]))) < 3e-15 ); else assert(1); if (zimags(out2[i])>1e-16) assert( (fabs(zimags(out2[i])+resI2[i]) / fabs(zimags(out2[i]))) < 3e-13 ); else assert(1); @@ -246,39 +245,31 @@ static void zifftmaTest(void){ printf(" >>> Matrice 3*4 <<< \n"); zifftma(in3, 3, 4, out3); for (i=0;i<12;i++){ - if (zreals(out3[i])>1e-16) assert( (fabs(zreals(out3[i])-resR3[i]) / fabs(zreals(out3[i]))) < 3e-16 ); + if (zreals(out3[i])>1e-16) assert( (fabs(zreals(out3[i])-resR3[i]) / fabs(zreals(out3[i]))) < 3e-15 ); else assert(1); - if (zimags(out3[i])>1e-16) assert( (fabs(zimags(out3[i])+resI3[i]) / fabs(zimags(out3[i]))) < 3e-16 ); + if (zimags(out3[i])>1e-16) assert( (fabs(zimags(out3[i])+resI3[i]) / fabs(zimags(out3[i]))) < 3e-15 ); else assert(1); } printf(" >>> Matrice 4*3 <<< \n"); zifftma(in4, 4, 3, out4); for (i=0;i<12;i++){ - if (zreals(out4[i])>1e-16) assert( (fabs(zreals(out4[i])-resR4[i]) / fabs(zreals(out4[i]))) < 3e-16 ); + if (zreals(out4[i])>1e-16) assert( (fabs(zreals(out4[i])-resR4[i]) / fabs(zreals(out4[i]))) < 3e-15 ); else assert(1); - if (zimags(out4[i])>1e-16) assert( (fabs(zimags(out4[i])+resI4[i]) / fabs(zimags(out4[i]))) < 3e-16 ); + if (zimags(out4[i])>1e-16) assert( (fabs(zimags(out4[i])+resI4[i]) / fabs(zimags(out4[i]))) < 3e-15 ); else assert(1); } printf(" >>> Matrice 6*2 <<< \n"); zifftma(in6, 6, 2, out6); for (i=0;i<12;i++){ - if (zreals(out6[i])>1e-16) assert( (fabs(zreals(out6[i])-resR6[i]) / fabs(zreals(out6[i]))) < 3e-16 ); + if (zreals(out6[i])>1e-16) assert( (fabs(zreals(out6[i])-resR6[i]) / fabs(zreals(out6[i]))) < 3e-15 ); else assert(1); - if (zimags(out6[i])>1e-16) assert( (fabs(zimags(out6[i])+resI6[i]) / fabs(zimags(out6[i]))) < 3e-16 ); + if (zimags(out6[i])>1e-16) assert( (fabs(zimags(out6[i])+resI6[i]) / fabs(zimags(out6[i]))) < 3e-15 ); else assert(1); } - printf(" >>> Matrice 3*3 <<< \n"); - zifftma(in9, 3, 3, out9); - for (i=0;i<9;i++){ - if (zreals(out9[i])>1e-16) assert( (fabs(zreals(out9[i])-resR9[i]) / fabs(zreals(out9[i]))) < 3e-16 ); - else assert(1); - - if (zimags(out9[i])>1e-16) assert( (fabs(zimags(out9[i])-resI9[i]) / fabs(zimags(out9[i]))) < 3e-16 ); - else assert(1); - } + } |