aboutsummaryrefslogtreecommitdiff
path: root/attic/dave/dclapack_test.c
diff options
context:
space:
mode:
Diffstat (limited to 'attic/dave/dclapack_test.c')
-rw-r--r--attic/dave/dclapack_test.c34
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;
+}
+