summaryrefslogtreecommitdiff
path: root/src/elementaryFunctions/int/testDoubleInt.c
diff options
context:
space:
mode:
Diffstat (limited to 'src/elementaryFunctions/int/testDoubleInt.c')
-rw-r--r--src/elementaryFunctions/int/testDoubleInt.c160
1 files changed, 0 insertions, 160 deletions
diff --git a/src/elementaryFunctions/int/testDoubleInt.c b/src/elementaryFunctions/int/testDoubleInt.c
deleted file mode 100644
index f5bf78b8..00000000
--- a/src/elementaryFunctions/int/testDoubleInt.c
+++ /dev/null
@@ -1,160 +0,0 @@
-/*
- * Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
- * Copyright (C) 2008 - INRIA - Arnaud TORSET
- *
- * This file must be used under the terms of the CeCILL.
- * This source file is licensed as described in the file COPYING, which
- * you should have received as part of this distribution. The terms
- * are also available at
- * http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
- *
- */
-
-
-#include <assert.h>
-#include <math.h>
-#include <stdio.h>
-#include "int.h"
-
-static void dintsTest(void) {
- /* FIXME : results could be differents of what wanted because of C number
- where 1 can be 0.99999999999992347
- This example fail :
- for(i=-2;i<-1;i+=.1){
- out = dints(i);
- }
- In the case of -1, -1 wasn't -1 exactly so error
- */
- assert ((fabs( dints(-2) - (-2))/fabs(dints(-2)))<1e-16);
- assert ((fabs( dints(-1.9) - (-1))/fabs(dints(-1.9)))<1e-16);
- assert ((fabs( dints(-1.8) - (-1))/fabs(dints(-1.8)))<1e-16);
- assert ((fabs( dints(-1.7) - (-1))/fabs(dints(-1.7)))<1e-16);
- assert ((fabs( dints(-1.6) - (-1))/fabs(dints(-1.6)))<1e-16);
- assert ((fabs( dints(-1.5) - (-1))/fabs(dints(-1.5)))<1e-16);
- assert ((fabs( dints(-1.4) - (-1))/fabs(dints(-1.4)))<1e-16);
- assert ((fabs( dints(-1.3) - (-1))/fabs(dints(-1.3)))<1e-16);
- assert ((fabs( dints(-1.2) - (-1))/fabs(dints(-1.2)))<1e-16);
- assert ((fabs( dints(-1.1) - (-1))/fabs(dints(-1.1)))<1e-16);
- assert ((fabs( dints(-1.0) - (-1))/fabs(dints(-1.0)))<1e-16);
- assert (fabs( dints(-.9))<1e-16);
- assert ((fabs( dints(1.0) - (1))/fabs(dints(1.0)))<1e-16);
- assert ((fabs( dints(1.1) - (1))/fabs(dints(1.1)))<1e-16);
- assert ((fabs( dints(1.2) - (1))/fabs(dints(1.2)))<1e-16);
- assert ((fabs( dints(1.3) - (1))/fabs(dints(1.3)))<1e-16);
- assert ((fabs( dints(1.4) - (1))/fabs(dints(1.4)))<1e-16);
- assert ((fabs( dints(1.5) - (1))/fabs(dints(1.5)))<1e-16);
- assert ((fabs( dints(1.6) - (1))/fabs(dints(1.6)))<1e-16);
- assert ((fabs( dints(1.7) - (1))/fabs(dints(1.7)))<1e-16);
- assert ((fabs( dints(1.8) - (1))/fabs(dints(1.8)))<1e-16);
- assert ((fabs( dints(1.9) - (1))/fabs(dints(1.9)))<1e-16);
- assert ((fabs( dints(2.0) - (2))/fabs(dints(2.0)))<1e-16);
-
-}
-
-static void zintsTest(void) {
- doubleComplex in, out;
-
-
- /* tests allant de -2 + i a -1 + 2*i, les reels décroissants de 0.1, les imaginaires croissant de 0.1
- + 1 test supplementaire : -0.9 + 0.9*i*/
-
- in=DoubleComplex(-2,1);
- out = zints(in);
- assert ((fabs( zreals(out) - (-2))/fabs(zreals(out)))<1e-16);
- assert ((fabs( zimags(out) - (1))/fabs(zimags(out)))<1e-16);
- in=DoubleComplex(-1.9,1.1);
- out = zints(in);
- assert ((fabs( zreals(out) - (-1))/fabs(zreals(out)))<1e-16);
- assert ((fabs( zimags(out) - (1))/fabs(zimags(out)))<1e-16);
- in=DoubleComplex(-1.8,1.2);
- out = zints(in);
- assert ((fabs( zreals(out) - (-1))/fabs(zreals(out)))<1e-16);
- assert ((fabs( zimags(out) - (1))/fabs(zimags(out)))<1e-16);
- in=DoubleComplex(-1.7,1.3);
- out = zints(in);
- assert ((fabs( zreals(out) - (-1))/fabs(zreals(out)))<1e-16);
- assert ((fabs( zimags(out) - (1))/fabs(zimags(out)))<1e-16);
- in=DoubleComplex(-1.6,1.4);
- out = zints(in);
- assert ((fabs( zreals(out) - (-1))/fabs(zreals(out)))<1e-16);
- assert ((fabs( zimags(out) - (1))/fabs(zimags(out)))<1e-16);
- in=DoubleComplex(-1.5,1.5);
- out = zints(in);
- assert ((fabs( zreals(out) - (-1))/fabs(zreals(out)))<1e-16);
- assert ((fabs( zimags(out) - (1))/fabs(zimags(out)))<1e-16);
- in=DoubleComplex(-1.4,1.6);
- out = zints(in);
- assert ((fabs( zreals(out) - (-1))/fabs(zreals(out)))<1e-16);
- assert ((fabs( zimags(out) - (1))/fabs(zimags(out)))<1e-16);
- in=DoubleComplex(-1.3,1.7);
- out = zints(in);
- assert ((fabs( zreals(out) - (-1))/fabs(zreals(out)))<1e-16);
- assert ((fabs( zimags(out) - (1))/fabs(zimags(out)))<1e-16);
- in=DoubleComplex(-1.2,1.8);
- out = zints(in);
- assert ((fabs( zreals(out) - (-1))/fabs(zreals(out)))<1e-16);
- assert ((fabs( zimags(out) - (1))/fabs(zimags(out)))<1e-16);
- in=DoubleComplex(-1.4,1.9);
- out = zints(in);
- assert ((fabs( zreals(out) - (-1))/fabs(zreals(out)))<1e-16);
- assert ((fabs( zimags(out) - (1))/fabs(zimags(out)))<1e-16);
- in=DoubleComplex(-1.0,2.0);
- out = zints(in);
- assert ((fabs( zreals(out) - (-1))/fabs(zreals(out)))<1e-16);
- assert ((fabs( zimags(out) - (2))/fabs(zimags(out)))<1e-16);
- in=DoubleComplex(-.9,.9);
- out = zints(in);
- assert (fabs( zreals(out))<1e-16);
- assert (fabs( zimags(out))<1e-16);
-}
-
-static void dintaTest(void) {
- int i;
- double in[10]={1.2, 2.0,-2.0,-3.1,-4.5,6.1,4.1,4.9,-0.7, 0.5};
- double res[10]={1,2,-2,-3,-4,6,4,4,0,0};
- double out[10];
-
- dinta(in,10,out);
-
- for (i=0;i<10;i++)
- if (out[i]!=0) assert(fabs(out[i]-res[i])/fabs(out[i]) <1e-16);
- else assert(fabs(out[i]-res[i]) <1e-16);
-
-
-
-}
-
-static void zintaTest(void) {
- int i;
- double inR[10]={1.2, 2.0,-2.0,-3.1,-4.5,6.1,4.1,4.9,-0.7, 0.5};
- double inI[10]={5.4,1.8,1.4,-4.0,-2.7,-0.3,1.5,1.4,7.3,-4.6};
- double resR[10]={1,2,-2,-3,-4,6,4,4,0,0};
- double resI[10]={5,1,1,-4,-2,0,1,1,7,-4};
- doubleComplex out[10],*in;
-
- in=DoubleComplexMatrix(inR,inI,10);
- zinta(in,10,out);
-
- for (i=0;i<10;i++){
- if(zreals(out[i])!=0) assert(fabs(zreals(out[i])-resR[i])/ fabs(zreals(out[i]))<1e-16);
- assert(fabs(zreals(out[i])-resR[i]) <1e-16);
- if(zimags(out[i])!=0) assert(fabs(zimags(out[i])-resI[i])/ fabs(zimags(out[i]))<1e-16);
- assert(fabs(zimags(out[i])-resI[i]) <1e-16);
- }
-}
-
-
-static int testInt() {
- printf("\n>>>> Double Int Tests\n");
- dintsTest();
- zintsTest();
- dintaTest();
- zintaTest();
-
- return 0;
-}
-
-int main() {
- assert(testInt() == 0);
- return 0;
-}