Changeset 14 for tests/test-lqup.C

Timestamp:

04/09/07 12:02:27 (2 years ago)

Author:

pernet

Message:

 

Files:

• tests/test-lqup.C (modified) (5 diffs)

tests/test-lqup.C

@@ -20 +20 @@
 #include "Matio.h"
 #include "timer.h"
-#include "fflas-ffpack/modular-balanced.h"
+#include "fflas-ffpack/modular-positive.h"
 #include "fflas-ffpack/ffpack.h"
 
@@ -36 +36 @@
                 exit(-1);
         }
-        Field F(atoi(argv[1]));
+        Field F((unsigned long)atoi(argv[1]));
         Field::Element * A;
         
@@ -63 +63 @@
                 tim.clear();      
                 tim.start();    
-                R = FFPACK::LUdivine (F, diag, m, n, A, n, P, Q,
-                                              FFPACK::FfpackLQUP, cutoff);
+                R = FFPACK::LUdivine_gauss (F, diag, m, n, A, n, P, Q,
+                                            FFPACK::FfpackLQUP);
                 tim.stop();
                 timc+=tim;
         }
         //write_field (F,cerr<<"Result = "<<endl, A, m,n,n);
-        
+
+//      cerr<<"P = [";
+//      for (size_t i=0; i<n; ++i)
+//              cerr<<P[i]<<" ";
+//      cerr<<"]"<<endl;
+//      cerr<<"Q = [";
+//      for (size_t i=0; i<m; ++i)
+//              cerr<<Q[i]<<" ";
+//      cerr<<"]"<<endl;
 #if DEBUG
         Field::Element * L = new Field::Element[m*m];
@@ -84 +92 @@
                         F.assign (*(U + i*n + j), *(A+ i*n+j));
         }
+        for (size_t i=R;i<m; ++i)
+                for (size_t j=0; j<n; ++j)
+                        F.assign(*(U+i*n+j), zero);
         for ( int i=0; i<m; ++i ){
                 int j=0;
-                for (; j< ((i<n)?i:n) ; ++j )
+                for (; j< ((i<R)?i:R) ; ++j )
                         F.assign( *(L + i*m+j), *(A+i*n+j));
                 for (; j<m; ++j )
@@ -92 +103 @@
         }
 
-        FFPACK::applyP( F, FFLAS::FflasRight, FFLAS::FflasNoTrans, m,0,m, L, m, Q);
+//      write_field(F,cerr<<"L = "<<endl,L,m,m,m);
+//      write_field(F,cerr<<"U = "<<endl,U,m,n,n);
+        FFPACK::applyP( F, FFLAS::FflasRight, FFLAS::FflasNoTrans, m,0,R, L, m, Q);
         if (diag == FFLAS::FflasNonUnit)
                 for ( int i=0; i<m; ++i )