Added dclapack.h

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;
+}
+