summaryrefslogtreecommitdiff
path: root/src/signalProcessing/ifft
diff options
context:
space:
mode:
Diffstat (limited to 'src/signalProcessing/ifft')
-rw-r--r--src/signalProcessing/ifft/testDoubleIfft.c64
-rw-r--r--src/signalProcessing/ifft/testFloatIfft.c68
-rw-r--r--src/signalProcessing/ifft/testMatIfft.c36
3 files changed, 82 insertions, 86 deletions
diff --git a/src/signalProcessing/ifft/testDoubleIfft.c b/src/signalProcessing/ifft/testDoubleIfft.c
index bf188e1e..315e6f58 100644
--- a/src/signalProcessing/ifft/testDoubleIfft.c
+++ b/src/signalProcessing/ifft/testDoubleIfft.c
@@ -386,13 +386,13 @@ static void zifftmaTest2 (void )
if ( zreals(out[i]) < 1e-14 && zreals (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zreals(out[i]) - zreals (Result[i]) ) / fabs (zreals (out[i])) < 1e-12 );
+ assert ( fabs( zreals(out[i]) - zreals (Result[i]) ) / fabs (zreals (out[i])) < 3e-16 );
if ( zimags(out[i]) < 1e-14 && zimags (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zimags(out[i]) - zimags (Result[i]) ) / fabs (zimags (out[i])) < 1e-12 ) ;
+ assert ( fabs( zimags(out[i]) - zimags (Result[i]) ) / fabs (zimags (out[i])) < 3e-16 ) ;
}
@@ -441,13 +441,13 @@ static void zifftmaTest3 (void )
if ( zreals(out[i]) < 1e-14 && zreals (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zreals(out[i]) - zreals (Result[i]) ) / fabs (zreals (out[i])) < 1e-12 );
+ assert ( fabs( zreals(out[i]) - zreals (Result[i]) ) / fabs (zreals (out[i])) < 3e-16 );
if ( zimags(out[i]) < 1e-14 && zimags (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zimags(out[i]) - zimags (Result[i]) ) / fabs (zimags (out[i])) < 1e-12 ) ;
+ assert ( fabs( zimags(out[i]) - zimags (Result[i]) ) / fabs (zimags (out[i])) < 3e-16 ) ;
}
@@ -491,13 +491,13 @@ static void zifftmaTest4 (void )
if ( zreals(out[i]) < 1e-14 && zreals (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zreals(out[i]) - zreals (Result[i]) ) / fabs (zreals (out[i])) < 1e-12 );
+ assert ( fabs( zreals(out[i]) - zreals (Result[i]) ) / fabs (zreals (out[i])) < 3e-16 );
if ( zimags(out[i]) < 1e-14 && zimags (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zimags(out[i]) - zimags (Result[i]) ) / fabs (zimags (out[i])) < 1e-12 ) ;
+ assert ( fabs( zimags(out[i]) - zimags (Result[i]) ) / fabs (zimags (out[i])) < 3e-16 ) ;
*/
}
@@ -548,13 +548,13 @@ static void zifftmaTest5 (void )
if ( zreals(out[i]) < 1e-14 && zreals (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zreals(out[i]) - zreals (Result[i]) ) / fabs (zreals (out[i])) < 1e-12 );
+ assert ( fabs( zreals(out[i]) - zreals (Result[i]) ) / fabs (zreals (out[i])) < 3e-16 );
if ( zimags(out[i]) < 1e-14 && zimags (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zimags(out[i]) - zimags (Result[i]) ) / fabs (zimags (out[i])) < 1e-12 ) ;
+ assert ( fabs( zimags(out[i]) - zimags (Result[i]) ) / fabs (zimags (out[i])) < 3e-16 ) ;
}
@@ -605,13 +605,13 @@ static void zifftmaTest6 (void )
if ( zreals(out[i]) < 1e-14 && zreals (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zreals(out[i]) - zreals (Result[i]) ) / fabs (zreals (out[i])) < 1e-12 );
+ assert ( fabs( zreals(out[i]) - zreals (Result[i]) ) / fabs (zreals (out[i])) < 3e-16 );
if ( zimags(out[i]) < 1e-14 && zimags (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zimags(out[i]) - zimags (Result[i]) ) / fabs (zimags (out[i])) < 1e-12 ) ;
+ assert ( fabs( zimags(out[i]) - zimags (Result[i]) ) / fabs (zimags (out[i])) < 3e-16 ) ;
}
@@ -662,13 +662,13 @@ static void zifftmaTest7 (void )
if ( zreals(out[i]) < 1e-14 && zreals (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zreals(out[i]) - zreals (Result[i]) ) / fabs (zreals (out[i])) < 1e-12 );
+ assert ( fabs( zreals(out[i]) - zreals (Result[i]) ) / fabs (zreals (out[i])) < 3e-16 );
if ( zimags(out[i]) < 1e-14 && zimags (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zimags(out[i]) - zimags (Result[i]) ) / fabs (zimags (out[i])) < 1e-12 ) ;
+ assert ( fabs( zimags(out[i]) - zimags (Result[i]) ) / fabs (zimags (out[i])) < 3e-16 ) ;
}
@@ -711,13 +711,13 @@ static void zifftmaTest8 (void )
if ( zreals(out[i]) < 1e-14 && zreals (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zreals(out[i]) - zreals (Result[i]) ) / fabs (zreals (out[i])) < 1e-12 );
+ assert ( fabs( zreals(out[i]) - zreals (Result[i]) ) / fabs (zreals (out[i])) < 3e-16 );
if ( zimags(out[i]) < 1e-14 && zimags (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zimags(out[i]) - zimags (Result[i]) ) / fabs (zimags (out[i])) < 1e-12 ) ;
+ assert ( fabs( zimags(out[i]) - zimags (Result[i]) ) / fabs (zimags (out[i])) < 3e-16 ) ;
*/
}
@@ -761,13 +761,13 @@ static void zifftmaTest9 (void )
if ( zreals(out[i]) < 1e-14 && zreals (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zreals(out[i]) - zreals (Result[i]) ) / fabs (zreals (out[i])) < 1e-12 );
+ assert ( fabs( zreals(out[i]) - zreals (Result[i]) ) / fabs (zreals (out[i])) < 3e-16 );
if ( zimags(out[i]) < 1e-14 && zimags (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zimags(out[i]) - zimags (Result[i]) ) / fabs (zimags (out[i])) < 1e-12 ) ;
+ assert ( fabs( zimags(out[i]) - zimags (Result[i]) ) / fabs (zimags (out[i])) < 3e-16 ) ;
}
@@ -820,13 +820,13 @@ static void zifftmaTest10 (void )
if ( zreals(out[i]) < 1e-14 && tRealResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zreals(out[i]) - tRealResult[i] ) / fabs (zreals (out[i])) < 1e-12 );
+ assert ( fabs( zreals(out[i]) - tRealResult[i] ) / fabs (zreals (out[i])) < 3e-16 );
if ( zimags(out[i]) < 1e-14 && tImagResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zimags(out[i]) - tImagResult[i] ) / fabs (zimags (out[i])) < 1e-12 ) ;
+ assert ( fabs( zimags(out[i]) - tImagResult[i] ) / fabs (zimags (out[i])) < 3e-16 ) ;
}
@@ -877,13 +877,13 @@ static void zifftmaTest11 (void )
if ( zreals(out[i]) < 1e-14 && tRealResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zreals(out[i]) - tRealResult[i] ) / fabs (zreals (out[i])) < 1e-12 );
+ assert ( fabs( zreals(out[i]) - tRealResult[i] ) / fabs (zreals (out[i])) < 3e-16 );
if ( zimags(out[i]) < 1e-14 && tImagResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zimags(out[i]) - tImagResult[i] ) / fabs (zimags (out[i])) < 1e-12 ) ;
+ assert ( fabs( zimags(out[i]) - tImagResult[i] ) / fabs (zimags (out[i])) < 3e-16 ) ;
}
@@ -935,13 +935,13 @@ static void zifftmaTest12 (void )
if ( zreals(out[i]) < 1e-14 && tRealResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zreals(out[i]) - tRealResult[i] ) / fabs (zreals (out[i])) < 1e-12 );
+ assert ( fabs( zreals(out[i]) - tRealResult[i] ) / fabs (zreals (out[i])) < 3e-16 );
if ( zimags(out[i]) < 1e-14 && tImagResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zimags(out[i]) - tImagResult[i] ) / fabs (zimags (out[i])) < 1e-12 ) ;
+ assert ( fabs( zimags(out[i]) - tImagResult[i] ) / fabs (zimags (out[i])) < 3e-16 ) ;
}
@@ -993,13 +993,13 @@ static void zifftmaTest13 (void )
if ( zreals(out[i]) < 1e-14 && tRealResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zreals(out[i]) - tRealResult[i] ) / fabs (zreals (out[i])) < 1e-12 );
+ assert ( fabs( zreals(out[i]) - tRealResult[i] ) / fabs (zreals (out[i])) < 3e-16 );
if ( zimags(out[i]) < 1e-14 && tImagResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zimags(out[i]) - tImagResult[i] ) / fabs (zimags (out[i])) < 1e-12 ) ;
+ assert ( fabs( zimags(out[i]) - tImagResult[i] ) / fabs (zimags (out[i])) < 3e-16 ) ;
}
@@ -1051,13 +1051,13 @@ static void zifftmaTest14 (void )
if ( zreals(out[i]) < 1e-14 && tRealResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zreals(out[i]) - tRealResult[i] ) / fabs (zreals (out[i])) < 1e-12 );
+ assert ( fabs( zreals(out[i]) - tRealResult[i] ) / fabs (zreals (out[i])) < 3e-16 );
if ( zimags(out[i]) < 1e-14 && tImagResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zimags(out[i]) - tImagResult[i] ) / fabs (zimags (out[i])) < 1e-12 ) ;
+ assert ( fabs( zimags(out[i]) - tImagResult[i] ) / fabs (zimags (out[i])) < 3e-16 ) ;
}
@@ -1110,13 +1110,13 @@ static void zifftmaTest15 (void )
if ( zreals(out[i]) < 1e-14 && tRealResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zreals(out[i]) - tRealResult[i] ) / fabs (zreals (out[i])) < 1e-12 );
+ assert ( fabs( zreals(out[i]) - tRealResult[i] ) / fabs (zreals (out[i])) < 3e-16 );
if ( zimags(out[i]) < 1e-14 && tImagResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zimags(out[i]) - tImagResult[i] ) / fabs (zimags (out[i])) < 1e-12 ) ;
+ assert ( fabs( zimags(out[i]) - tImagResult[i] ) / fabs (zimags (out[i])) < 3e-16 ) ;
}
@@ -1169,13 +1169,13 @@ static void zifftmaTest16 (void )
if ( zreals(out[i]) < 1e-14 && tRealResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zreals(out[i]) - tRealResult[i] ) / fabs (zreals (out[i])) < 1e-12 );
+ assert ( fabs( zreals(out[i]) - tRealResult[i] ) / fabs (zreals (out[i])) < 3e-16 );
if ( zimags(out[i]) < 1e-14 && tImagResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zimags(out[i]) - tImagResult[i] ) / fabs (zimags (out[i])) < 1e-12 ) ;
+ assert ( fabs( zimags(out[i]) - tImagResult[i] ) / fabs (zimags (out[i])) < 3e-16 ) ;
}
@@ -1228,13 +1228,13 @@ static void zifftmaTest32 (void )
if ( zreals(out[i]) < 1e-14 && tRealResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zreals(out[i]) - tRealResult[i] ) / fabs (zreals (out[i])) < 1e-12 );
+ assert ( fabs( zreals(out[i]) - tRealResult[i] ) / fabs (zreals (out[i])) < 3e-16 );
if ( zimags(out[i]) < 1e-14 && tImagResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( zimags(out[i]) - tImagResult[i] ) / fabs (zimags (out[i])) < 1e-12 ) ;
+ assert ( fabs( zimags(out[i]) - tImagResult[i] ) / fabs (zimags (out[i])) < 3e-16 ) ;
}
diff --git a/src/signalProcessing/ifft/testFloatIfft.c b/src/signalProcessing/ifft/testFloatIfft.c
index 4dc0f442..e00d5432 100644
--- a/src/signalProcessing/ifft/testFloatIfft.c
+++ b/src/signalProcessing/ifft/testFloatIfft.c
@@ -383,13 +383,13 @@ static void cifftmaTest2 (void )
if ( creals(out[i]) < 1e-14 && creals (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( creals(out[i]) - creals (Result[i]) ) / fabs (creals (out[i])) < 1e-4 );
+ assert ( fabs( creals(out[i]) - creals (Result[i]) ) / fabs (creals (out[i])) < 3e-6 );
if ( cimags(out[i]) < 1e-14 && cimags (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( cimags(out[i]) - cimags (Result[i]) ) / fabs (cimags (out[i])) < 1e-4 ) ;
+ assert ( fabs( cimags(out[i]) - cimags (Result[i]) ) / fabs (cimags (out[i])) < 3e-6 ) ;
}
@@ -434,18 +434,16 @@ static void cifftmaTest3 (void )
cimags (Result[i]),
fabs( creals(out[i]) - creals (Result[i]) ) / fabs (creals (out[i])) ,
fabs( cimags(out[i]) - cimags (Result[i]) ) / fabs (cimags (out[i])));
-/*
if ( creals(out[i]) < 1e-14 && creals (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( creals(out[i]) - creals (Result[i]) ) / fabs (creals (out[i])) < 1e-4 );
+ assert ( fabs( creals(out[i]) - creals (Result[i]) ) / fabs (creals (out[i])) < 3e-6 );
if ( cimags(out[i]) < 1e-14 && cimags (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( cimags(out[i]) - cimags (Result[i]) ) / fabs (cimags (out[i])) < 1e-4 ) ;
-*/
+ assert ( fabs( cimags(out[i]) - cimags (Result[i]) ) / fabs (cimags (out[i])) < 3e-6 ) ;
}
free(out);
free(in);
@@ -482,18 +480,16 @@ static void cifftmaTest4 (void )
i ,creals(out[i]) , cimags(out[i]), creals (Result[i]) , cimags (Result[i]),
fabs( creals(out[i]) - creals (Result[i]) ) / fabs (creals (out[i])) ,
fabs( cimags(out[i]) - cimags (Result[i]) ) / fabs (cimags (out[i])));
-/*
if ( creals(out[i]) < 1e-14 && creals (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( creals(out[i]) - creals (Result[i]) ) / fabs (creals (out[i])) < 1e-4 );
+ assert ( fabs( creals(out[i]) - creals (Result[i]) ) / fabs (creals (out[i])) < 3e-6 );
if ( cimags(out[i]) < 1e-14 && cimags (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( cimags(out[i]) - cimags (Result[i]) ) / fabs (cimags (out[i])) < 1e-4 ) ;
-*/
+ assert ( fabs( cimags(out[i]) - cimags (Result[i]) ) / fabs (cimags (out[i])) < 3e-6 ) ;
}
free(out);
free(in);
@@ -541,13 +537,13 @@ static void cifftmaTest5 (void )
if ( creals(out[i]) < 1e-14 && creals (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( creals(out[i]) - creals (Result[i]) ) / fabs (creals (out[i])) < 1e-4 );
+ assert ( fabs( creals(out[i]) - creals (Result[i]) ) / fabs (creals (out[i])) < 3e-6 );
if ( cimags(out[i]) < 1e-14 && cimags (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( cimags(out[i]) - cimags (Result[i]) ) / fabs (cimags (out[i])) < 1e-4 ) ;
+ assert ( fabs( cimags(out[i]) - cimags (Result[i]) ) / fabs (cimags (out[i])) < 3e-6 ) ;
}
free(out);
@@ -596,13 +592,13 @@ static void cifftmaTest6 (void )
if ( creals(out[i]) < 1e-14 && creals (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( creals(out[i]) - creals (Result[i]) ) / fabs (creals (out[i])) < 1e-4 );
+ assert ( fabs( creals(out[i]) - creals (Result[i]) ) / fabs (creals (out[i])) < 3e-6 );
if ( cimags(out[i]) < 1e-14 && cimags (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( cimags(out[i]) - cimags (Result[i]) ) / fabs (cimags (out[i])) < 1e-4 ) ;
+ assert ( fabs( cimags(out[i]) - cimags (Result[i]) ) / fabs (cimags (out[i])) < 3e-6 ) ;
}
free(out);
@@ -651,13 +647,13 @@ static void cifftmaTest7 (void )
if ( creals(out[i]) < 1e-14 && creals (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( creals(out[i]) - creals (Result[i]) ) / fabs (creals (out[i])) < 1e-4 );
+ assert ( fabs( creals(out[i]) - creals (Result[i]) ) / fabs (creals (out[i])) < 3e-6 );
if ( cimags(out[i]) < 1e-14 && cimags (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( cimags(out[i]) - cimags (Result[i]) ) / fabs (cimags (out[i])) < 1e-4 ) ;
+ assert ( fabs( cimags(out[i]) - cimags (Result[i]) ) / fabs (cimags (out[i])) < 3e-6 ) ;
}
free(out);
@@ -698,13 +694,13 @@ static void cifftmaTest8 (void )
if ( creals(out[i]) < 1e-14 && creals (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( creals(out[i]) - creals (Result[i]) ) / fabs (creals (out[i])) < 1e-4 );
+ assert ( fabs( creals(out[i]) - creals (Result[i]) ) / fabs (creals (out[i])) < 3e-6 );
if ( cimags(out[i]) < 1e-14 && cimags (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( cimags(out[i]) - cimags (Result[i]) ) / fabs (cimags (out[i])) < 1e-4 ) ;
+ assert ( fabs( cimags(out[i]) - cimags (Result[i]) ) / fabs (cimags (out[i])) < 3e-6 ) ;
}
free(out);
@@ -746,13 +742,13 @@ static void cifftmaTest9 (void )
if ( creals(out[i]) < 1e-14 && creals (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( creals(out[i]) - creals (Result[i]) ) / fabs (creals (out[i])) < 1e-4 );
+ assert ( fabs( creals(out[i]) - creals (Result[i]) ) / fabs (creals (out[i])) < 3e-6 );
if ( cimags(out[i]) < 1e-14 && cimags (Result[i]) < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( cimags(out[i]) - cimags (Result[i]) ) / fabs (cimags (out[i])) < 1e-4 ) ;
+ assert ( fabs( cimags(out[i]) - cimags (Result[i]) ) / fabs (cimags (out[i])) < 3e-6 ) ;
}
free(out);
@@ -803,13 +799,13 @@ static void cifftmaTest10 (void )
if ( creals(out[i]) < 1e-14 && tRealResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( creals(out[i]) - tRealResult[i] ) / fabs (creals (out[i])) < 1e-4 );
+ assert ( fabs( creals(out[i]) - tRealResult[i] ) / fabs (creals (out[i])) < 3e-6 );
if ( cimags(out[i]) < 1e-14 && tImagResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( cimags(out[i]) - tImagResult[i] ) / fabs (cimags (out[i])) < 1e-4 ) ;
+ assert ( fabs( cimags(out[i]) - tImagResult[i] ) / fabs (cimags (out[i])) < 3e-6 ) ;
}
free(out);
@@ -858,13 +854,13 @@ static void cifftmaTest11 (void )
if ( creals(out[i]) < 1e-14 && tRealResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( creals(out[i]) - tRealResult[i] ) / fabs (creals (out[i])) < 1e-4 );
+ assert ( fabs( creals(out[i]) - tRealResult[i] ) / fabs (creals (out[i])) < 3e-6 );
if ( cimags(out[i]) < 1e-14 && tImagResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( cimags(out[i]) - tImagResult[i] ) / fabs (cimags (out[i])) < 1e-4 ) ;
+ assert ( fabs( cimags(out[i]) - tImagResult[i] ) / fabs (cimags (out[i])) < 3e-6 ) ;
}
free(out);
@@ -915,13 +911,13 @@ static void cifftmaTest12 (void )
if ( creals(out[i]) < 1e-14 && tRealResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( creals(out[i]) - tRealResult[i] ) / fabs (creals (out[i])) < 1e-4 );
+ assert ( fabs( creals(out[i]) - tRealResult[i] ) / fabs (creals (out[i])) < 3e-6 );
if ( cimags(out[i]) < 1e-14 && tImagResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( cimags(out[i]) - tImagResult[i] ) / fabs (cimags (out[i])) < 1e-4 ) ;
+ assert ( fabs( cimags(out[i]) - tImagResult[i] ) / fabs (cimags (out[i])) < 3e-6 ) ;
}
free(out);
@@ -972,13 +968,13 @@ static void cifftmaTest13 (void )
if ( creals(out[i]) < 1e-14 && tRealResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( creals(out[i]) - tRealResult[i] ) / fabs (creals (out[i])) < 1e-4 );
+ assert ( fabs( creals(out[i]) - tRealResult[i] ) / fabs (creals (out[i])) < 3e-6 );
if ( cimags(out[i]) < 1e-14 && tImagResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( cimags(out[i]) - tImagResult[i] ) / fabs (cimags (out[i])) < 1e-4 ) ;
+ assert ( fabs( cimags(out[i]) - tImagResult[i] ) / fabs (cimags (out[i])) < 3e-6 ) ;
}
free(out);
@@ -1029,13 +1025,13 @@ static void cifftmaTest14 (void )
if ( creals(out[i]) < 1e-14 && tRealResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( creals(out[i]) - tRealResult[i] ) / fabs (creals (out[i])) < 1e-4 );
+ assert ( fabs( creals(out[i]) - tRealResult[i] ) / fabs (creals (out[i])) < 3e-6 );
if ( cimags(out[i]) < 1e-14 && tImagResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( cimags(out[i]) - tImagResult[i] ) / fabs (cimags (out[i])) < 1e-4 ) ;
+ assert ( fabs( cimags(out[i]) - tImagResult[i] ) / fabs (cimags (out[i])) < 3e-6 ) ;
}
free(out);
@@ -1087,13 +1083,13 @@ static void cifftmaTest15 (void )
if ( creals(out[i]) < 1e-14 && tRealResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( creals(out[i]) - tRealResult[i] ) / fabs (creals (out[i])) < 1e-4 );
+ assert ( fabs( creals(out[i]) - tRealResult[i] ) / fabs (creals (out[i])) < 3e-6 );
if ( cimags(out[i]) < 1e-14 && tImagResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( cimags(out[i]) - tImagResult[i] ) / fabs (cimags (out[i])) < 1e-4 ) ;
+ assert ( fabs( cimags(out[i]) - tImagResult[i] ) / fabs (cimags (out[i])) < 3e-6 ) ;
}
free(out);
@@ -1145,13 +1141,13 @@ static void cifftmaTest16 (void )
if ( creals(out[i]) < 1e-14 && tRealResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( creals(out[i]) - tRealResult[i] ) / fabs (creals (out[i])) < 1e-4 );
+ assert ( fabs( creals(out[i]) - tRealResult[i] ) / fabs (creals (out[i])) < 3e-6 );
if ( cimags(out[i]) < 1e-14 && tImagResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( cimags(out[i]) - tImagResult[i] ) / fabs (cimags (out[i])) < 1e-4 ) ;
+ assert ( fabs( cimags(out[i]) - tImagResult[i] ) / fabs (cimags (out[i])) < 3e-6 ) ;
}
free(out);
@@ -1203,13 +1199,13 @@ static void cifftmaTest32 (void )
if ( creals(out[i]) < 1e-14 && tRealResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( creals(out[i]) - tRealResult[i] ) / fabs (creals (out[i])) < 1e-4 );
+ assert ( fabs( creals(out[i]) - tRealResult[i] ) / fabs (creals (out[i])) < 3e-6 );
if ( cimags(out[i]) < 1e-14 && tImagResult[i] < 1e-18 )
assert ( 1 ) ;
else
- assert ( fabs( cimags(out[i]) - tImagResult[i] ) / fabs (cimags (out[i])) < 1e-4 ) ;
+ assert ( fabs( cimags(out[i]) - tImagResult[i] ) / fabs (cimags (out[i])) < 3e-6 ) ;
}
free(out);
diff --git a/src/signalProcessing/ifft/testMatIfft.c b/src/signalProcessing/ifft/testMatIfft.c
index 88de15e5..33aee9be 100644
--- a/src/signalProcessing/ifft/testMatIfft.c
+++ b/src/signalProcessing/ifft/testMatIfft.c
@@ -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-14) assert( (fabs(out1[i]-resR1[i]) / fabs(out1[i])) < 1e-14 );
+ if (out1[i]>1e-16) assert( (fabs(out1[i]-resR1[i]) / fabs(out1[i])) < 3e-16 );
else assert(1);
}
printf(" >>> Matrice 2*6 <<< \n");
difftma(in2, 2, 6, out2);
for (i=0;i<12;i++){
- if (out2[i]>1e-14) assert( (fabs(out2[i]-resR2[i]) / fabs(out2[i])) < 1e-14 );
+ if (out2[i]>1e-16) assert( (fabs(out2[i]-resR2[i]) / fabs(out2[i])) < 3e-16 );
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-14) assert( (fabs(out3[i]-resR3[i]) / fabs(out3[i])) < 1e-14 );
+ if (out3[i]>1e-16) assert( (fabs(out3[i]-resR3[i]) / fabs(out3[i])) < 3e-16 );
else assert(1);
}
printf(" >>> Matrice 4*3 <<< \n");
difftma(in4, 4, 3, out4);
for (i=0;i<12;i++){
- if (out4[i]>1e-14) assert( (fabs(out4[i]-resR4[i]) / fabs(out4[i])) < 1e-14 );
+ if (out4[i]>1e-16) assert( (fabs(out4[i]-resR4[i]) / fabs(out4[i])) < 3e-16 );
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])) < 1e-14 );
+ if (out6[i]>1e-16) assert( (fabs(out6[i]-resR6[i]) / fabs(out6[i])) < 3e-16 );
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])) < 1e-16 );
+ if (out9[i]>1e-16) assert( (fabs(out9[i]-resR9[i]) / fabs(out9[i])) < 3e-16 );
else assert(1);
}
@@ -227,18 +227,18 @@ 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-14) assert( (fabs(zreals(out1[i])-resR1[i]) / fabs(zreals(out1[i]))) < 1e-14 );
+ if (zreals(out1[i])>1e-16) assert( (fabs(zreals(out1[i])-resR1[i]) / fabs(zreals(out1[i]))) < 3e-16 );
else assert(1);
- if (zimags(out1[i])>1e-14) assert( (fabs(zimags(out1[i])+resI1[i]) / fabs(zimags(out1[i]))) < 1e-14 );
+ if (zimags(out1[i])>1e-16) assert( (fabs(zimags(out1[i])+resI1[i]) / fabs(zimags(out1[i]))) < 3e-16 );
else assert(1);
}
printf(" >>> Matrice 2*6 <<< \n");
zifftma(in2, 2, 6, out2);
for (i=0;i<12;i++){
- if (zreals(out2[i])>1e-14) assert( (fabs(zreals(out2[i])-resR2[i]) / fabs(zreals(out2[i]))) < 1e-14 );
+ if (zreals(out2[i])>1e-16) assert( (fabs(zreals(out2[i])-resR2[i]) / fabs(zreals(out2[i]))) < 3e-16 );
else assert(1);
- if (zimags(out2[i])>1e-14) assert( (fabs(zimags(out2[i])+resI2[i]) / fabs(zimags(out2[i]))) < 1e-13 );
+ if (zimags(out2[i])>1e-16) assert( (fabs(zimags(out2[i])+resI2[i]) / fabs(zimags(out2[i]))) < 3e-13 );
else assert(1);
}
@@ -246,37 +246,37 @@ 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-14) assert( (fabs(zreals(out3[i])-resR3[i]) / fabs(zreals(out3[i]))) < 1e-14 );
+ 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-14) 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);
}
printf(" >>> Matrice 4*3 <<< \n");
zifftma(in4, 4, 3, out4);
for (i=0;i<12;i++){
- if (zreals(out4[i])>1e-14) assert( (fabs(zreals(out4[i])-resR4[i]) / fabs(zreals(out4[i]))) < 1e-14 );
+ if (zreals(out4[i])>1e-16) assert( (fabs(zreals(out4[i])-resR4[i]) / fabs(zreals(out4[i]))) < 3e-16 );
else assert(1);
- if (zimags(out4[i])>1e-14) assert( (fabs(zimags(out4[i])+resI4[i]) / fabs(zimags(out4[i]))) < 1e-14 );
+ if (zimags(out4[i])>1e-16) assert( (fabs(zimags(out4[i])+resI4[i]) / fabs(zimags(out4[i]))) < 3e-16 );
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]))) < 1e-14 );
+ if (zreals(out6[i])>1e-16) assert( (fabs(zreals(out6[i])-resR6[i]) / fabs(zreals(out6[i]))) < 3e-16 );
else assert(1);
- if (zimags(out6[i])>1e-16) assert( (fabs(zimags(out6[i])+resI6[i]) / fabs(zimags(out6[i]))) < 1e-14 );
+ if (zimags(out6[i])>1e-16) assert( (fabs(zimags(out6[i])+resI6[i]) / fabs(zimags(out6[i]))) < 3e-16 );
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]))) < 1e-16 );
+ 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-15) assert( (fabs(zimags(out9[i])-resI9[i]) / fabs(zimags(out9[i]))) < 1e-15 );
+ if (zimags(out9[i])>1e-16) assert( (fabs(zimags(out9[i])-resI9[i]) / fabs(zimags(out9[i]))) < 3e-16 );
else assert(1);
}
}
generated by cgit (git 2.34.1) at 2025-03-26 07:14:07 +0000