diff options
Diffstat (limited to 'src/matrixOperations/determ')
| -rw-r--r-- | src/matrixOperations/determ/testDoubleDeterm.c | 40 | ||||
| -rw-r--r-- | src/matrixOperations/determ/testFloatDeterm.c | 10 |
2 files changed, 25 insertions, 25 deletions
diff --git a/src/matrixOperations/determ/testDoubleDeterm.c b/src/matrixOperations/determ/testDoubleDeterm.c index 45372375..b8462514 100644 --- a/src/matrixOperations/determ/testDoubleDeterm.c +++ b/src/matrixOperations/determ/testDoubleDeterm.c @@ -378,13 +378,13 @@ static void ddetermaTest(void){ out8=ddeterma(in8,8); out10=ddeterma(in10,10); - assert((fabs(out2-res2)/fabs(out2))<1e-6); - assert((fabs(out3-res3)/fabs(out3))<1e-7); - assert((fabs(out4-res4)/fabs(out4))<1e-6); - assert((fabs(out5-res5)/fabs(out5))<1e-6); - assert((fabs(out6-res6)/fabs(out6))<1e-6); - assert((fabs(out8-res8)/fabs(out8))<1e-6); - assert((fabs(out10-res10)/fabs(out10))<1e-6); + assert((fabs(out2-res2)/fabs(out2))<3e-16); + assert((fabs(out3-res3)/fabs(out3))<3e-16); + assert((fabs(out4-res4)/fabs(out4))<3e-16); + assert((fabs(out5-res5)/fabs(out5))<3e-16); + assert((fabs(out6-res6)/fabs(out6))<3e-16); + assert((fabs(out8-res8)/fabs(out8))<3e-16); + assert((fabs(out10-res10)/fabs(out10))<3e-16); } static void zdetermaTest(void){ @@ -415,26 +415,26 @@ static void zdetermaTest(void){ out10=zdeterma(in10,10); - assert((fabs(zreals(out2)-resR2)/fabs(zreals(out2)))<1e-15); - assert((fabs(zimags(out2)-resI2)/fabs(zimags(out2)))<1e-15); + assert((fabs(zreals(out2)-resR2)/fabs(zreals(out2)))<3e-16); + assert((fabs(zimags(out2)-resI2)/fabs(zimags(out2)))<3e-16); - assert((fabs(zreals(out3)-resR3)/fabs(zreals(out3)))<1e-15); + assert((fabs(zreals(out3)-resR3)/fabs(zreals(out3)))<3e-16); assert((fabs(zimags(out3)-resI3)/fabs(zimags(out3)))<1e-16); - assert((fabs(zreals(out4)-resR4)/fabs(zreals(out4)))<1e-14); - assert((fabs(zimags(out4)-resI4)/fabs(zimags(out4)))<1e-14); + assert((fabs(zreals(out4)-resR4)/fabs(zreals(out4)))<3e-16); + assert((fabs(zimags(out4)-resI4)/fabs(zimags(out4)))<3e-16); - assert((fabs(zreals(out5)-resR5)/fabs(zreals(out5)))<1e-15); - assert((fabs(zimags(out5)-resI5)/fabs(zimags(out5)))<1e-15); + assert((fabs(zreals(out5)-resR5)/fabs(zreals(out5)))<3e-16); + assert((fabs(zimags(out5)-resI5)/fabs(zimags(out5)))<3e-16); - assert((fabs(zreals(out6)-resR6)/fabs(zreals(out6)))<1e-14); - assert((fabs(zimags(out6)-resI6)/fabs(zimags(out6)))<1e-13); + assert((fabs(zreals(out6)-resR6)/fabs(zreals(out6)))<3e-16); + assert((fabs(zimags(out6)-resI6)/fabs(zimags(out6)))<3e-16); - assert((fabs(zreals(out8)-resR8)/fabs(zreals(out8)))<1e-15); - assert((fabs(zimags(out8)-resI8)/fabs(zimags(out8)))<1e-15); + assert((fabs(zreals(out8)-resR8)/fabs(zreals(out8)))<3e-16); + assert((fabs(zimags(out8)-resI8)/fabs(zimags(out8)))<3e-16); - assert((fabs(zreals(out10)-resR10)/fabs(zreals(out10)))<1e-14); - assert((fabs(zimags(out10)-resI10)/fabs(zimags(out10)))<1e-14); + assert((fabs(zreals(out10)-resR10)/fabs(zreals(out10)))<3e-16); + assert((fabs(zimags(out10)-resI10)/fabs(zimags(out10)))<3e-16); } static int determTest(void){ diff --git a/src/matrixOperations/determ/testFloatDeterm.c b/src/matrixOperations/determ/testFloatDeterm.c index 76de97fc..d429d020 100644 --- a/src/matrixOperations/determ/testFloatDeterm.c +++ b/src/matrixOperations/determ/testFloatDeterm.c @@ -381,11 +381,11 @@ static void sdetermaTest(void){ assert((fabs(out2-res2)/fabs(out2))<1e-6); assert((fabs(out3-res3)/fabs(out3))<1e-7); - assert((fabs(out4-res4)/fabs(out4))<1e-5); + assert((fabs(out4-res4)/fabs(out4))<3e-6); assert((fabs(out5-res5)/fabs(out5))<1e-6); assert((fabs(out6-res6)/fabs(out6))<1e-6); assert((fabs(out8-res8)/fabs(out8))<1e-6); - assert((fabs(out10-res10)/fabs(out10))<1e-5); + assert((fabs(out10-res10)/fabs(out10))<3e-6); } @@ -426,16 +426,16 @@ static void cdetermaTest(void){ assert((fabs(creals(out4)-resR4)/fabs(creals(out4)))<1e-6); assert((fabs(cimags(out4)-resI4)/fabs(cimags(out4)))<1e-6); - assert((fabs(creals(out5)-resR5)/fabs(creals(out5)))<1e-5); + assert((fabs(creals(out5)-resR5)/fabs(creals(out5)))<3e-6); assert((fabs(cimags(out5)-resI5)/fabs(cimags(out5)))<1e-6); assert((fabs(creals(out6)-resR6)/fabs(creals(out6)))<1e-6); - assert((fabs(cimags(out6)-resI6)/fabs(cimags(out6)))<1e-5); + assert((fabs(cimags(out6)-resI6)/fabs(cimags(out6)))<3e-6); assert((fabs(creals(out8)-resR8)/fabs(creals(out8)))<1e-6); assert((fabs(cimags(out8)-resI8)/fabs(cimags(out8)))<1e-6); - assert((fabs(creals(out10)-resR10)/fabs(creals(out10)))<1e-5); + assert((fabs(creals(out10)-resR10)/fabs(creals(out10)))<3e-6); assert((fabs(cimags(out10)-resI10)/fabs(cimags(out10)))<1e-6); } |