diff options
-rw-r--r-- | redist/dclhelpers.c | 4 | ||||
-rw-r--r-- | redist/test_dcl.c | 39 |
2 files changed, 41 insertions, 2 deletions
diff --git a/redist/dclhelpers.c b/redist/dclhelpers.c index c3725f8..fb6aba6 100644 --- a/redist/dclhelpers.c +++ b/redist/dclhelpers.c @@ -5,9 +5,9 @@ #include "dclapack.h" #include <alloca.h> -void dclPrint( const DCL_FLOAT * A, int Ac, int n, int m ) +void dclPrint( const DCL_FLOAT * PMATRIX, int PMATRIXc, int n, int m ) { - PRINT( A, n, m ); + PRINT( PMATRIX, n, m ); } void dclIdentity( DCL_FLOAT * I, int Ic, int n ) diff --git a/redist/test_dcl.c b/redist/test_dcl.c new file mode 100644 index 0000000..adea7b5 --- /dev/null +++ b/redist/test_dcl.c @@ -0,0 +1,39 @@ +#include "dclhelpers.h" +#include <stdint.h> +#include <stdio.h> + + +int main() +{ + FLT A[2][4] = { { 0, 1, 2, 3 }, { 4, 5, 6, 7} }; + FLT B[4][2]; + dclPrint( A[0], 4, 2, 4 ); + dclTransp( B[0], 2, A[0], 4, 2, 4 ); + dclPrint( B[0], 2, 4, 2 ); + + int i; + for( i = 0; i < 8; i++ ) + { + printf( "%f\n", ((float*)(B[0]))[i] ); + } + + FLT M[3][3] = { + { .32, 1, 0 }, + { 0, 1, 2 }, + { 1, 0, 1 } }; + FLT Mo[3][3]; + dclInv( Mo[0], 3, M[0], 3, 3 ); + dclPrint( Mo[0], 3, 3, 3 ); + + FLT MM[3][3]; + dclMul( MM[0], 3, M[0], 3, Mo[0], 3, 3, 3, 3 ); + + printf( "The following should be an identity matrix\n" ); + dclPrint( MM[0], 3, 3, 3 ); + +//void dclTransp( DCL_FLOAT * R, int Rc, const DCL_FLOAT * A, int Ac, int n, int m ); + +// dclIdentity( A[0], MATx, 5 ); +// dclPrint( A[0], MATx, MATx, MATy ); +} + |