summaryrefslogtreecommitdiff
path: root/2.3-1/src/c/matrixOperations/eye/testMatrixEye.c
diff options
context:
space:
mode:
Diffstat (limited to '2.3-1/src/c/matrixOperations/eye/testMatrixEye.c')
-rw-r--r--2.3-1/src/c/matrixOperations/eye/testMatrixEye.c348
1 files changed, 348 insertions, 0 deletions
diff --git a/2.3-1/src/c/matrixOperations/eye/testMatrixEye.c b/2.3-1/src/c/matrixOperations/eye/testMatrixEye.c
new file mode 100644
index 00000000..a3682743
--- /dev/null
+++ b/2.3-1/src/c/matrixOperations/eye/testMatrixEye.c
@@ -0,0 +1,348 @@
+/*
+ * 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;
+}
+
+
generated by cgit (git 2.34.1) at 2024-12-20 11:29:56 +0000