/* * 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 #include #include #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; }