Moving source code

1 files changed, 0 insertions, 348 deletions

diff --git a/src/matrixOperations/eye/testMatrixEye.c b/src/matrixOperations/eye/testMatrixEye.c
deleted file mode 100644
index a3682743..00000000
--- a/src/matrixOperations/eye/testMatrixEye.c
+++ /dev/null
@@ -1,348 +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 <math.h>
-
-#include "eye.h"
-
-#define LEADDIM 10
-#define ROWS    10
-#define COLS    5
-
-#define ROWS2    5
-#define COLS2    10
-
-static void deyeaTest ( void )
-{
-    int i, j = 0 ;
-
-    double result = 0 ;
-    double in [LEADDIM*LEADDIM] ;
-    double in2 [ROWS * COLS];
-    double in3 [ROWS2 * COLS2];
-
-    deyea ( 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]) < 3e-16 ) ;
-        }
-    
-    deyea ( in2, ROWS, COLS );
-    for ( i = 0 ; i <  ROWS ; i++ )
-        {
-            for ( j = 0 ; j < COLS ; ++j )
-            {
-                if ( i == j )
-                    result = 1 ;
-                else
-                    result = 0 ;
-
-                printf ( "\t\t in[%d, %d] : %e\tresult : %e\tassert : %e \n" , i, j,
-                        in[i * ROWS + j] , result ,
-                        fabs( in[i * ROWS + j] - result) / fabs( in[i * ROWS + j]) ) ;
-
-                if ( in[i * ROWS + j] < 1e-14 && result < 1e-14 )
-                    assert(1);
-                else
-                    assert ( fabs ( in[i * ROWS + j] - result) / fabs( in[i * ROWS + j]) < 3e-16 ) ;
-            }
-        }
-
-    deyea ( in3, ROWS2, COLS2);
-    for ( j = 0 ; j < COLS2 ; ++j )
-        {
-            for ( i = 0 ; i <  ROWS2 ; i++ )
-            {
-                if ( i == j )
-                    result = 1 ;
-                else
-                    result = 0 ;
-
-                printf ( "\t\t in[%d, %d] : %e\tresult : %e\tassert : %e \n" , i, j,
-                        in3[j * ROWS2 + i] , result ,
-                        fabs( in3[j * ROWS2 + i] - result) / fabs( in3[j * ROWS2 + i]) ) ;
-
-                if ( in3[j * ROWS2 + i] < 1e-14 && result < 1e-14 )
-                    assert(1);
-                else
-                    assert ( fabs ( in3[j * ROWS2 + i] - result) / fabs( in3[j * ROWS2 + i]) < 3e-16 ) ;
-            }
-        }
-    
-}
-
-
-static void seyeaTest ( void )
-{
-    int i, j = 0 ;
-
-    float result = 0 ;
-
-    float in[LEADDIM*LEADDIM] ;
-    float in2 [ROWS * COLS];
-    float in3 [ROWS2 * COLS2];
-    seyea ( 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 ) ;
-        }
-    
-    seyea ( in2, ROWS, COLS );
-    for ( i = 0 ; i <  ROWS ; i++ )
-        {
-            for ( j = 0 ; j < COLS ; ++j )
-            {
-                if ( i == j )
-                    result = 1 ;
-                else
-                    result = 0 ;
-
-                printf ( "\t\t in[%d, %d] : %e\tresult : %e\tassert : %e \n" , i, j,
-                        in[i * ROWS + j] , result ,
-                        fabs( in[i * ROWS + j] - result) / fabs( in[i * ROWS + j]) ) ;
-
-                if ( in[i * ROWS + j] < 1e-8 && result < 1e-8 )
-                    assert(1);
-                else
-                    assert ( fabs ( in[i * ROWS + j] - result) / fabs( in[i * ROWS + j]) < 1e-8 ) ;
-            }
-        }
-
-    seyea ( in3, ROWS2, COLS2);
-    for ( j = 0 ; j < COLS2 ; ++j )
-        {
-            for ( i = 0 ; i <  ROWS2 ; i++ )
-            {
-                if ( i == j )
-                    result = 1 ;
-                else
-                    result = 0 ;
-
-                printf ( "\t\t in[%d, %d] : %e\tresult : %e\tassert : %e \n" , i, j,
-                        in3[j * ROWS2 + i] , result ,
-                        fabs( in3[j * ROWS2 + i] - result) / fabs( in3[j * ROWS2 + i]) ) ;
-
-                if ( in3[j * ROWS2 + i] < 1e-14 && result < 1e-14 )
-                    assert(1);
-                else
-                    assert ( fabs ( in3[j * ROWS2 + i] - result) / fabs( in3[j * ROWS2 + i]) < 1e-14 ) ;
-            }
-        }
-
-}
-
-
-static void zeyeaTest ( void )
-{
-    int i = 0 ;
-
-    doubleComplex result = DoubleComplex ( 0 , 0) ;
-
-    doubleComplex in[LEADDIM*LEADDIM] ;
-
-    zeyea ( 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])) < 3e-16) ;
-
-        if ( zimags ( in[i]) < 1e-14 )
-            assert (1);
-        else
-            assert (0);
-
-        }
-}
-
-
-
-
-
-static void ceyeaTest ( void )
-{
-    int i = 0 ;
-
-    floatComplex result = FloatComplex ( 0 , 0  ) ;
-
-    floatComplex in[LEADDIM*LEADDIM] ;
-
-    ceyea ( 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);
-
-        }
-}
-
-/*  EYES TEST  */
-
-static void ceyesTest ( void )
-{
-  
-
-    floatComplex result = FloatComplex ( 1.0f , 0  ) ;
-
-    floatComplex in = FloatComplex ( LEADDIM, LEADDIM ) ;
-
-    in = ceyes ( in  ) ;
-
-
-   
-     assert (  (creals(in) - creals(result)) / fabs( creals(in)) < 1e-6 ) ;
-    assert  (  cimags (in) == 0) ;
-}
-
-
-
-
-
-static void zeyesTest ( void )
-{
-
-    doubleComplex result = DoubleComplex ( 1 , 0  ) ;
-
-    doubleComplex in = DoubleComplex ( LEADDIM , LEADDIM ) ;
-
-    in = zeyes ( in  ) ;
-
-
-     assert (  (zreals(in) - zreals(result)) / fabs( zreals(in)) < 1e-16 ) ;
-    assert  (  zimags (in) == 0) ;
-}
-
-
-
-static void deyesTest ( void )
-{
-
-
-    double result =  1  ;
-
-    double in = LEADDIM ;
-
-    in = deyes ( in  ) ;
-
-
-     assert ( (in - result) / in < 1e-16 ) ; ;
-}
-
-
-
-static void seyesTest ( void )
-{
-
-
-    float result =  1.0f  ;
-
-    float in = 154.0f ;
-
-    in = seyes ( in  ) ;
-
-
-     assert ( (in - result) / in < 1e-8 ) ;
-
-}
-
-
-
-static int testEye(void) {
-
-  printf("\n>>>> Matrix Eye Tests\n");
-  printf("\t>>>> Matrix Double Realt Tests\n");
-  deyeaTest();
-
-  printf("\n\n\t>>>> Matrix Float Realt Tests\n");
-  seyeaTest();
-
-  printf("\n\n\n");
-  printf("\t>>>> Matrix Float Complex Tests\n");
-  ceyeaTest();
-
-  printf("\n\n\n");
-  printf("\t>>>> Matrix Double Complex Tests\n");
-  zeyeaTest();
-    
-    
-    
-    
-  printf("\t>>>> Matrix Double Realt Tests\n");
-  deyesTest();
-
-  printf("\n\n\t>>>> Matrix Float Realt Tests\n");
-  seyesTest();
-
-  printf("\n\n\n");
-  printf("\t>>>> Matrix Float Complex Tests\n");
-  ceyesTest();
-
-  printf("\n\n\n");
-  printf("\t>>>> Matrix Double Complex Tests\n");
-  zeyesTest();
-
-  return 0;
-}
-
-
-
-int main(void) {
-  assert(testEye() == 0);
-  return 0;
-}
-
-