Changeset 21 for tests

Timestamp:

05/01/07 19:23:34 (2 years ago)

Author:

pernet

Message:

Introduction of the new automatic choice of underlying BLAS, for any finite finite field implementation:
float/double is chosen depending on the prime, the dimension, and efficiency considerations.

Files:

• tests/test-fgemv.C (modified) (3 diffs)

tests/test-fgemv.C

@@ -7 +7 @@
 //-------------------------------------------------------------------------
 
-#define DEBUG 0
+#define DEBUG 1
 #define TIME 1
 
@@ -45 +45 @@
         Field::Element * A;
         Field::Element * b;
-        size_t lda;
-        size_t ldb;
         
         b = read_field(F,argv[3],&n,&k);
@@ -67 +65 @@
         }
 
+#if DEBUG
+        Field::Element *d = new Field::Element[n];
+        for (int i=0; i<m; ++i)
+                F.mul (d[i], beta, b[i]);
+        for (int i=0; i<m; ++i)
+                F.mulin (b[i], alpha);
+        for (int i=0; i<m; ++i)
+                for (int j=0; j<n; ++j)
+                        F.axpyin (d[i], *(A+i*m+j), b[j]);
+        bool fail = false;
+        for (int i=0; i<m; ++i)
+                if (!F.areEqual(d[i], c[i]))
+                        fail = true;
+
+        if (fail)
+                cerr<<"FAIL"<<endl;
+        else
+                cerr<<"PASS"<<endl;
+        delete[] d;
+#endif
+        delete[] A;
+        delete[] b;
+        delete[] c;
 #if TIME
         double mflops = (2.0*(m*n/1000000.0)*nbit/tim.usertime());