diff options
author | torset | 2009-02-27 14:43:06 +0000 |
---|---|---|
committer | torset | 2009-02-27 14:43:06 +0000 |
commit | 032f682c8ee956db344b1aae267a28e344ace94b (patch) | |
tree | 0b91956992576bdb9fc3a3cff3927dc46ad41086 /src/signalProcessing/hilbert/testHilbert.c | |
parent | 8e207d069eb29be177718a7d2e760434d1e40e1f (diff) | |
download | scilab2c-032f682c8ee956db344b1aae267a28e344ace94b.tar.gz scilab2c-032f682c8ee956db344b1aae267a28e344ace94b.tar.bz2 scilab2c-032f682c8ee956db344b1aae267a28e344ace94b.zip |
change precision for tests : put it to 3e-16 for double and 3e-6 for float
Diffstat (limited to 'src/signalProcessing/hilbert/testHilbert.c')
-rw-r--r-- | src/signalProcessing/hilbert/testHilbert.c | 52 |
1 files changed, 26 insertions, 26 deletions
diff --git a/src/signalProcessing/hilbert/testHilbert.c b/src/signalProcessing/hilbert/testHilbert.c index 8a41760d..4743d8f4 100644 --- a/src/signalProcessing/hilbert/testHilbert.c +++ b/src/signalProcessing/hilbert/testHilbert.c @@ -390,9 +390,9 @@ static void dhilbertaTest(void){ /* Matrice 3*3 */ printf(">>> Matrice 3*3 <<<\n"); for (i=0;i<9;i++){ - if (zreals(out3[i])>1e-16) assert( (fabs(zreals(out3[i])-resR3[i]) / fabs(zreals(out3[i])) )<1e-15); + if (zreals(out3[i])>1e-16) assert( (fabs(zreals(out3[i])-resR3[i]) / fabs(zreals(out3[i])) )<3e-16); else assert(1); - if (zimags(out3[i])>1e-16) assert( (fabs(zimags(out3[i])+resI3[i]) / fabs(zimags(out3[i])) )<1e-14); + if (zimags(out3[i])>1e-16) assert( (fabs(zimags(out3[i])+resI3[i]) / fabs(zimags(out3[i])) )<3e-16); else assert(1); } @@ -408,18 +408,18 @@ static void dhilbertaTest(void){ /* Matrice 5*5 */ printf(">>> Matrice 5*5 <<<\n"); for (i=0;i<25;i++){ - if (zreals(out5[i])>1e-16) assert( (fabs(zreals(out5[i])-resR5[i]) / fabs(zreals(out5[i])) )<1e-14); + if (zreals(out5[i])>1e-16) assert( (fabs(zreals(out5[i])-resR5[i]) / fabs(zreals(out5[i])) )<3e-16); else assert(1); - if (zimags(out5[i])>1e-16) assert( (fabs(zimags(out5[i])+resI5[i]) / fabs(zimags(out5[i])) )<1e-13); + if (zimags(out5[i])>1e-16) assert( (fabs(zimags(out5[i])+resI5[i]) / fabs(zimags(out5[i])) )<3e-16); else assert(1); } /* Matrice 7*7 */ printf(">>> Matrice 7*7 <<<\n"); for (i=0;i<49;i++){ - if (zreals(out7[i])>1e-16) assert( (fabs(zreals(out7[i])-resR7[i]) / fabs(zreals(out7[i])) )<1e-13); + if (zreals(out7[i])>1e-16) assert( (fabs(zreals(out7[i])-resR7[i]) / fabs(zreals(out7[i])) )<3e-16); else assert(1); - if (zimags(out7[i])>1e-16) assert( (fabs(zimags(out7[i])+resI7[i]) / fabs(zimags(out7[i])) )<1e-13); + if (zimags(out7[i])>1e-16) assert( (fabs(zimags(out7[i])+resI7[i]) / fabs(zimags(out7[i])) )<3e-16); else assert(1); } @@ -435,7 +435,7 @@ static void dhilbertaTest(void){ /* Vecteur 1*3 */ printf(">>> Matrice 1*3 <<<\n"); for (i=0;i<3;i++){ - if (zreals(out13[i])>1e-16) assert( (fabs(zreals(out13[i])-resR13[i]) / fabs(zreals(out13[i])) )<1e-15); + if (zreals(out13[i])>1e-16) assert( (fabs(zreals(out13[i])-resR13[i]) / fabs(zreals(out13[i])) )<3e-16); else assert(1); if (zimags(out13[i])>1e-16) assert( (fabs(zimags(out13[i])+resI13[i]) / fabs(zimags(out13[i])) )<1e-16); else assert(1); @@ -453,25 +453,25 @@ static void dhilbertaTest(void){ /* Vecteur 1*5 */ printf(">>> Matrice 1*5 <<<\n"); for (i=0;i<5;i++){ - if (zreals(out15[i])>1e-16) assert( (fabs(zreals(out15[i])-resR15[i]) / fabs(zreals(out15[i])) )<1e-15); + if (zreals(out15[i])>1e-16) assert( (fabs(zreals(out15[i])-resR15[i]) / fabs(zreals(out15[i])) )<3e-16); else assert(1); - if (zimags(out15[i])>1e-16) assert( (fabs(zimags(out15[i])+resI15[i]) / fabs(zimags(out15[i])) )<1e-15); + if (zimags(out15[i])>1e-16) assert( (fabs(zimags(out15[i])+resI15[i]) / fabs(zimags(out15[i])) )<3e-16); else assert(1); } /* Vecteur 1*7 */ printf(">>> Matrice 1*7 <<<\n"); for (i=0;i<7;i++){ - if (zreals(out17[i])>1e-16) assert( (fabs(zreals(out17[i])-resR17[i]) / fabs(zreals(out17[i])) )<1e-15); + if (zreals(out17[i])>1e-16) assert( (fabs(zreals(out17[i])-resR17[i]) / fabs(zreals(out17[i])) )<3e-16); else assert(1); - if (zimags(out17[i])>1e-16) assert( (fabs(zimags(out17[i])+resI17[i]) / fabs(zimags(out17[i])) )<1e-15); + if (zimags(out17[i])>1e-16) assert( (fabs(zimags(out17[i])+resI17[i]) / fabs(zimags(out17[i])) )<3e-16); else assert(1); } /* Matrices 3*1 */ printf(">>> Matrice 3*1 <<<\n"); for (i=0;i<3;i++){ - if (zreals(out31[i])>1e-16) assert( (fabs(zreals(out31[i])-resR31[i]) / fabs(zreals(out31[i])) )<1e-15); + if (zreals(out31[i])>1e-16) assert( (fabs(zreals(out31[i])-resR31[i]) / fabs(zreals(out31[i])) )<3e-16); else assert(1); if (zimags(out31[i])>1e-16) assert( (fabs(zimags(out31[i])+resI31[i]) / fabs(zimags(out31[i])) )<1e-16); else assert(1); @@ -480,7 +480,7 @@ static void dhilbertaTest(void){ /* Matrices 3*2 */ printf(">>> Matrice 3*2 <<<\n"); for (i=0;i<6;i++){ - if (zreals(out32[i])>1e-16) assert( (fabs(zreals(out32[i])-resR32[i]) / fabs(zreals(out32[i])) )<1e-15); + if (zreals(out32[i])>1e-16) assert( (fabs(zreals(out32[i])-resR32[i]) / fabs(zreals(out32[i])) )<3e-16); else assert(1); if (zimags(out32[i])>1e-16) assert( (fabs(zimags(out32[i])+resI32[i]) / fabs(zimags(out32[i])) )<1e-16); else assert(1); @@ -492,7 +492,7 @@ static void dhilbertaTest(void){ for (i=0;i<12;i++){ if (zreals(out34[i])>1e-16) assert( (fabs(zreals(out34[i])-resR34[i]) / fabs(zreals(out34[i])) )<1e-16); else assert(1); - if (zimags(out34[i])>1e-16) assert( (fabs(zimags(out34[i])+resI34[i]) / fabs(zimags(out34[i])) )<1e-15); + if (zimags(out34[i])>1e-16) assert( (fabs(zimags(out34[i])+resI34[i]) / fabs(zimags(out34[i])) )<3e-16); else assert(1); } @@ -500,9 +500,9 @@ static void dhilbertaTest(void){ /* Matrices 3*5 */ printf(">>> Matrice 3*5 <<<\n"); for (i=0;i<15;i++){ - if (zreals(out35[i])>1e-16) assert( (fabs(zreals(out35[i])-resR35[i]) / fabs(zreals(out35[i])) )<1e-15); + if (zreals(out35[i])>1e-16) assert( (fabs(zreals(out35[i])-resR35[i]) / fabs(zreals(out35[i])) )<3e-16); else assert(1); - if (zimags(out35[i])>1e-16) assert( (fabs(zimags(out35[i])+resI35[i]) / fabs(zimags(out35[i])) )<1e-14); + if (zimags(out35[i])>1e-16) assert( (fabs(zimags(out35[i])+resI35[i]) / fabs(zimags(out35[i])) )<3e-16); else assert(1); } @@ -510,9 +510,9 @@ static void dhilbertaTest(void){ /* Matrices 3*7 */ printf(">>> Matrice 3*7 <<<\n"); for (i=0;i<21;i++){ - if (zreals(out37[i])>1e-16) assert( (fabs(zreals(out37[i])-resR37[i]) / fabs(zreals(out37[i])) )<1e-13); + if (zreals(out37[i])>1e-16) assert( (fabs(zreals(out37[i])-resR37[i]) / fabs(zreals(out37[i])) )<3e-16); else assert(1); - if (zimags(out37[i])>1e-16) assert( (fabs(zimags(out37[i])+resI37[i]) / fabs(zimags(out37[i])) )<1e-14); + if (zimags(out37[i])>1e-16) assert( (fabs(zimags(out37[i])+resI37[i]) / fabs(zimags(out37[i])) )<3e-16); else assert(1); } @@ -674,7 +674,7 @@ static void shilbertaTest(void){ for (i=0;i<9;i++){ if (creals(out3[i])>1e-16) assert( (fabs(creals(out3[i])-resR3[i]) / fabs(creals(out3[i])) )<1e-6); else assert(1); - if (cimags(out3[i])>1e-16) assert( (fabs(cimags(out3[i])+resI3[i]) / fabs(cimags(out3[i])) )<1e-5); + if (cimags(out3[i])>1e-16) assert( (fabs(cimags(out3[i])+resI3[i]) / fabs(cimags(out3[i])) )<3e-6); else assert(1); } @@ -692,9 +692,9 @@ static void shilbertaTest(void){ printf(">>> Matrice 5*5 <<<\n"); shilberta(in5, 5, 5, out5); for (i=0;i<25;i++){ - if (creals(out5[i])>1e-16) assert( (fabs(creals(out5[i])-resR5[i]) / fabs(creals(out5[i])) )<1e-5); + if (creals(out5[i])>1e-16) assert( (fabs(creals(out5[i])-resR5[i]) / fabs(creals(out5[i])) )<3e-6); else assert(1); - if (cimags(out5[i])>1e-16) assert( (fabs(cimags(out5[i])+resI5[i]) / fabs(cimags(out5[i])) )<1e-5); + if (cimags(out5[i])>1e-16) assert( (fabs(cimags(out5[i])+resI5[i]) / fabs(cimags(out5[i])) )<3e-6); else assert(1); } @@ -704,7 +704,7 @@ static void shilbertaTest(void){ for (i=0;i<49;i++){ if (creals(out7[i])>1e-16) assert( (fabs(creals(out7[i])-resR7[i]) / fabs(creals(out7[i])) )<1e-6); else assert(1); - if (cimags(out7[i])>1e-16) assert( (fabs(cimags(out7[i])+resI7[i]) / fabs(cimags(out7[i])) )<1e-5); + if (cimags(out7[i])>1e-16) assert( (fabs(cimags(out7[i])+resI7[i]) / fabs(cimags(out7[i])) )<3e-6); else assert(1); } @@ -732,7 +732,7 @@ static void shilbertaTest(void){ printf(">>> Matrice 1*4 <<<\n"); shilberta(in14, 1, 4, out14); for (i=0;i<4;i++){ - if (creals(out14[i])>1e-16) assert( (fabs(creals(out14[i])-resR14[i]) / fabs(creals(out14[i])) )<1e-5); + if (creals(out14[i])>1e-16) assert( (fabs(creals(out14[i])-resR14[i]) / fabs(creals(out14[i])) )<3e-6); else assert(1); if (cimags(out14[i])>1e-16) assert( (fabs(cimags(out14[i])+resI14[i]) / fabs(cimags(out14[i])) )<1e-6); else assert(1); @@ -796,7 +796,7 @@ static void shilbertaTest(void){ for (i=0;i<15;i++){ if (creals(out35[i])>1e-16) assert( (fabs(creals(out35[i])-resR35[i]) / fabs(creals(out35[i])) )<1e-6); else assert(1); - if (cimags(out35[i])>1e-16) assert( (fabs(cimags(out35[i])+resI35[i]) / fabs(cimags(out35[i])) )<1e-5); + if (cimags(out35[i])>1e-16) assert( (fabs(cimags(out35[i])+resI35[i]) / fabs(cimags(out35[i])) )<3e-6); else assert(1); } @@ -805,9 +805,9 @@ static void shilbertaTest(void){ printf(">>> Matrice 3*7 <<<\n"); shilberta(in37, 3, 7, out37); for (i=0;i<21;i++){ - if (creals(out37[i])>1e-16) assert( (fabs(creals(out37[i])-resR37[i]) / fabs(creals(out37[i])) )<1e-5); + if (creals(out37[i])>1e-16) assert( (fabs(creals(out37[i])-resR37[i]) / fabs(creals(out37[i])) )<3e-6); else assert(1); - if (cimags(out37[i])>1e-16) assert( (fabs(cimags(out37[i])+resI37[i]) / fabs(cimags(out37[i])) )<1e-5); + if (cimags(out37[i])>1e-16) assert( (fabs(cimags(out37[i])+resI37[i]) / fabs(cimags(out37[i])) )<3e-6); else assert(1); } |