diff options
Diffstat (limited to 'src/signalProcessing/lev/testDoubleLev.c')
| -rw-r--r-- | src/signalProcessing/lev/testDoubleLev.c | 26 |
1 files changed, 13 insertions, 13 deletions
diff --git a/src/signalProcessing/lev/testDoubleLev.c b/src/signalProcessing/lev/testDoubleLev.c index d9d4b2f5..6967890e 100644 --- a/src/signalProcessing/lev/testDoubleLev.c +++ b/src/signalProcessing/lev/testDoubleLev.c @@ -33,7 +33,7 @@ static void dlevaTest(void) { ; for (i=0;i<3;i++) { - if (out1[i]>1e-16) assert((fabs(out1[i]-result1[i]) / fabs (out1[i]))<1e-15); + if (out1[i]>1e-16) assert((fabs(out1[i]-result1[i]) / fabs (out1[i]))<3e-16); else assert(1); } @@ -75,11 +75,11 @@ static void dlevaTest(void) { else assert(1); } - if (out2>1e-16) assert((fabs(out2-result2) / fabs (out2))<1e-15); + if (out2>1e-16) assert((fabs(out2-result2) / fabs (out2))<3e-16); for (i=0;i<7;i++) { - if (out3[i]>1e-16) assert((fabs(out3[i]-result3[i]) / fabs (out3[i]))<1e-14); + if (out3[i]>1e-16) assert((fabs(out3[i]-result3[i]) / fabs (out3[i]))<3e-16); else assert(1); } @@ -110,19 +110,19 @@ static void zlevaTest(void) { out2 = zleva(in,4,out1,out3); for (i=0;i<3;i++){ - if (zreals(out1[i])>1e-16) assert((fabs(zreals(out1[i])-result1R[i]) / fabs(zreals(out1[i]))) <1e-15); + if (zreals(out1[i])>1e-16) assert((fabs(zreals(out1[i])-result1R[i]) / fabs(zreals(out1[i]))) <3e-16); else assert(1); if (zimags(out1[i])>1e-16) assert((fabs(zimags(out1[i])-result1I[i]) / fabs(zimags(out1[i]))) <1e-16); else assert(1); } - if (zreals(out2)>1e-16) assert((fabs(zreals(out2)-result2R) / fabs(zreals(out2))) <1e-15); + if (zreals(out2)>1e-16) assert((fabs(zreals(out2)-result2R) / fabs(zreals(out2))) <3e-16); else assert(1); - if (zimags(out2)>1e-16) assert((fabs(zimags(out2)-result2I) / fabs(zimags(out2))) <1e-15); + if (zimags(out2)>1e-16) assert((fabs(zimags(out2)-result2I) / fabs(zimags(out2))) <3e-16); else assert(1); for (i=0;i<3;i++){ - if (zreals(out3[i])>1e-16) assert((fabs(zreals(out3[i])-result3R[i]) / fabs(zreals(out3[i]))) <1e-15); + if (zreals(out3[i])>1e-16) assert((fabs(zreals(out3[i])-result3R[i]) / fabs(zreals(out3[i]))) <3e-16); else assert(1); if (zimags(out3[i])>1e-16) assert((fabs(zimags(out3[i])-result3I[i]) / fabs(zimags(out3[i]))) <1e-16); else assert(1); @@ -162,21 +162,21 @@ static void zlevaTest(void) { for (i=0;i<5;i++) printf("out : %1.20f+%1.20f\n",zreals(out1[i]),zimags(out1[i])); for (i=0;i<5;i++){ - if (zreals(out1[i])>1e-16) assert((fabs(zreals(out1[i])-result1R[i]) / fabs(zreals(out1[i]))) <1e-14); + if (zreals(out1[i])>1e-16) assert((fabs(zreals(out1[i])-result1R[i]) / fabs(zreals(out1[i]))) <3e-16); else assert(1); - if (zimags(out1[i])>1e-16) assert((fabs(zimags(out1[i])-result1I[i]) / fabs(zimags(out1[i]))) <1e-14); + if (zimags(out1[i])>1e-16) assert((fabs(zimags(out1[i])-result1I[i]) / fabs(zimags(out1[i]))) <3e-16); else assert(1); } - if (zreals(out2)>1e-16) assert((fabs(zreals(out2)-result2R) / fabs(zreals(out2))) <1e-14); + if (zreals(out2)>1e-16) assert((fabs(zreals(out2)-result2R) / fabs(zreals(out2))) <3e-16); else assert(1); - if (zimags(out2)>1e-16) assert((fabs(zimags(out2)-result2I) / fabs(zimags(out2))) <1e-14); + if (zimags(out2)>1e-16) assert((fabs(zimags(out2)-result2I) / fabs(zimags(out2))) <3e-16); else assert(1); for (i=0;i<5;i++){ - if (zreals(out3[i])>1e-16) assert((fabs(zreals(out3[i])-result3R[i]) / fabs(zreals(out3[i]))) <1e-15); + if (zreals(out3[i])>1e-16) assert((fabs(zreals(out3[i])-result3R[i]) / fabs(zreals(out3[i]))) <3e-16); else assert(1); - if (zimags(out3[i])>1e-16) assert((fabs(zimags(out3[i])-result3I[i]) / fabs(zimags(out3[i]))) <1e-14); + if (zimags(out3[i])>1e-16) assert((fabs(zimags(out3[i])-result3I[i]) / fabs(zimags(out3[i]))) <3e-16); else assert(1); } } |