summaryrefslogtreecommitdiff
path: root/src/matrixOperations/eyes/testMatrixEyes.c
diff options
context:
space:
mode:
Diffstat (limited to 'src/matrixOperations/eyes/testMatrixEyes.c')
-rw-r--r--src/matrixOperations/eyes/testMatrixEyes.c169
1 files changed, 0 insertions, 169 deletions
diff --git a/src/matrixOperations/eyes/testMatrixEyes.c b/src/matrixOperations/eyes/testMatrixEyes.c
deleted file mode 100644
index a4e5711b..00000000
--- a/src/matrixOperations/eyes/testMatrixEyes.c
+++ /dev/null
@@ -1,169 +0,0 @@
-/*
- * Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
- * Copyright (C) 2008-2008 - INRIA - Allan SIMON
- *
- * 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 <stdio.h>
-#include "matrixEyes.h"
-#include <math.h>
-
-#define LEADDIM 10
-
-static void deyesaTest ( void )
-{
- int i = 0 ;
-
- double result = 0 ;
-
- double in [LEADDIM*LEADDIM] ;
-
- deyesa ( in , LEADDIM , LEADDIM ) ;
-
- for ( i = 0 ; i < LEADDIM*LEADDIM ; i++ )
- {
- if ( i%(LEADDIM+1) == 0 )
- result = 1 ;
- else
- result = 0 ;
-
- printf ( "\t\t %d in : %e\tresult : %e\tassert : %e \n" , i, in[i] , result , fabs( in[i] - result) / fabs( in[i]) ) ;
-
- if ( in[i] < 1e-14 && result < 1e-14 )
- assert(1);
- else
- assert ( fabs ( in[i] - result) / fabs( in[i]) < 1e-14 ) ;
- }
-}
-
-
-static void seyesaTest ( void )
-{
- int i = 0 ;
-
- float result = 0 ;
-
- float in[LEADDIM*LEADDIM] ;
-
- seyesa ( in , LEADDIM , LEADDIM ) ;
-
- for ( i = 0 ; i < LEADDIM*LEADDIM ; i++)
- {
- if ( i%(LEADDIM+1) == 0 )
- result = 1.0f ;
- else
- result = 0 ;
-
- printf ( "\t\t %d in : %e\tresult : %e\tassert : %e \n" , i, in[i] , result , fabs( in[i] - result) / fabs( in[i]) ) ;
- if ( in[i] < 1e-6 && result < 1e-6 )
- assert(1);
- else
- assert ( fabs ( in[i] - result) / fabs( in[i]) < 1e-6 ) ;
- }
-}
-
-
-static void zeyesaTest ( void )
-{
- int i = 0 ;
-
- doubleComplex result = DoubleComplex ( 0 , 0) ;
-
- doubleComplex in[LEADDIM*LEADDIM] ;
-
- zeyesa ( in , LEADDIM , LEADDIM ) ;
-
- for ( i = 0 ; i < LEADDIM*LEADDIM ; i++ )
- {
- if ( i%(LEADDIM+1) == 0 )
- result = DoubleComplex ( 1 , 0 ) ;
- else
- result = DoubleComplex ( 0 , 0 ) ;
-
- printf ( "\t\t %d in : %e\tresult : %e \n" , i, zreals( in[i]) , zreals ( result ) ) ;
-
- if ( zreals( in[i]) < 1e-14 && zreals( result) < 1e-14 )
- assert(1);
- else
- assert ( fabs ( zreals(in[i]) - zreals(result)) / fabs( zreals(in[i])) < 1e-14 ) ;
-
- if ( zimags ( in[i]) < 1e-14 )
- assert (1);
- else
- assert (0);
-
- }
-}
-
-
-
-
-
-static void ceyesaTest ( void )
-{
- int i = 0 ;
-
- floatComplex result = FloatComplex ( 0 , 0 ) ;
-
- floatComplex in[LEADDIM*LEADDIM] ;
-
- ceyesa ( in , LEADDIM , LEADDIM ) ;
-
- for ( i = 0 ; i < LEADDIM*LEADDIM ; i++)
- {
- if ( i%(LEADDIM+1) == 0 )
- result = FloatComplex ( 1.0f , 0 ) ;
- else
- result = FloatComplex ( 0 , 0 ) ;
-
- printf ( "\t\t %d in : %e\tresult : %e \n" , i, creals( in[i]) , creals ( result ) ) ;
- if ( creals( in[i]) < 1e-6 && creals( result) < 1e-6 )
- assert(1);
- else
- assert ( fabs ( creals(in[i]) - creals(result)) / fabs( creals(in[i])) < 1e-6 ) ;
-
- if ( cimags ( in[i]) < 1e-6 )
- assert (1);
- else
- assert (0);
-
- }
-}
-
-
-
-static int testExponential(void) {
-
- printf("\n>>>> Matrix Exponential Tests\n");
- printf("\t>>>> Matrix Double Realt Tests\n");
- deyesaTest();
-
- printf("\n\n\t>>>> Matrix Float Realt Tests\n");
- seyesaTest();
-
- printf("\n\n\n");
- printf("\t>>>> Matrix Float Complex Tests\n");
- ceyesaTest();
-
- printf("\n\n\n");
- printf("\t>>>> Matrix Double Complex Tests\n");
- zeyesaTest();
-
- return 0;
-}
-
-
-
-int main(void) {
- assert(testExponential() == 0);
- return 0;
-}
-
-