diff options
author | ultramn <dchapm2@umbc.edu> | 2016-12-16 13:50:55 -0800 |
---|---|---|
committer | ultramn <dchapm2@umbc.edu> | 2016-12-16 13:50:55 -0800 |
commit | 93ad6a810479a5cc09a809f2ea23d549c7cc9c2a (patch) | |
tree | 2df4787ec7386670e96d7c7846506e8a052e290e /dave/dclapack_test.c | |
parent | 4dbacd89da4fe29ec6ffcf0bd6e1fc4d942b4e27 (diff) | |
download | libsurvive-93ad6a810479a5cc09a809f2ea23d549c7cc9c2a.tar.gz libsurvive-93ad6a810479a5cc09a809f2ea23d549c7cc9c2a.tar.bz2 |
Added dclapack.h
Diffstat (limited to 'dave/dclapack_test.c')
-rw-r--r-- | dave/dclapack_test.c | 33 |
1 files changed, 33 insertions, 0 deletions
diff --git a/dave/dclapack_test.c b/dave/dclapack_test.c new file mode 100644 index 0000000..38ae651 --- /dev/null +++ b/dave/dclapack_test.c @@ -0,0 +1,33 @@ +#define FLOAT float +#define ORDER 50 +#include "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 = 12; + 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(Ainv,n,n); + MUL(A,Ainv,Prod,n,n,n); + PRINT(Prod,n,n); + + return 0; +} + |