summaryrefslogtreecommitdiff
path: root/src/signalProcessing/hilbert/testHilbert.c
diff options
context:
space:
mode:
authortorset2009-02-27 14:43:06 +0000
committertorset2009-02-27 14:43:06 +0000
commit032f682c8ee956db344b1aae267a28e344ace94b (patch)
tree0b91956992576bdb9fc3a3cff3927dc46ad41086 /src/signalProcessing/hilbert/testHilbert.c
parent8e207d069eb29be177718a7d2e760434d1e40e1f (diff)
downloadscilab2c-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.c52
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);
}