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