diff options
Diffstat (limited to 'attic/dave/dclapack_test.c')
-rw-r--r-- | attic/dave/dclapack_test.c | 34 |
1 files changed, 34 insertions, 0 deletions
diff --git a/attic/dave/dclapack_test.c b/attic/dave/dclapack_test.c new file mode 100644 index 0000000..3f48b08 --- /dev/null +++ b/attic/dave/dclapack_test.c @@ -0,0 +1,34 @@ +#define FLOAT float +#define ORDER 50 +#include "../redist/dclapack.h" + +#include <stdio.h> +#include <stdlib.h> + +int main() +{ + float A[ORDER][ORDER]; + float Ainv[ORDER][ORDER]; + float Prod[ORDER][ORDER]; + + int i, j, n = 3; + srand(7779); + + for(i=0; i<n; i++) { + for(j=0; j<n; j++) { + A[i][j] = (float)rand() / RAND_MAX; + } + } + + for (i=0; i<10000; i++) { + INV(A,Ainv,n); + } + + PRINT(A,n,n); + PRINT(Ainv,n,n); + MUL(A,Ainv,Prod,n,n,n); + PRINT(Prod,n,n); + + return 0; +} + |