Changeset 26 for tests/test-lqup.C

Timestamp:

07/06/07 18:58:19 (2 years ago)

Author:

pernet

Message:

New row echelon form and reduced row echelon form

Files:

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

tests/test-lqup.C

@@ -20 +20 @@
 #include "Matio.h"
 #include "timer.h"
-#include "fflas-ffpack/modular-balanced.h"
+//#include "fflas-ffpack/modular-balanced.h"
 #include "fflas-ffpack/modular-positive.h"
 #include "fflas-ffpack/ffpack.h"
@@ -68 +68 @@
                         A = read_field(F,argv[2],&m,&n);
                 }
-                for (j=0;j<n;j++)
+                for (j=0;j<maxP;j++)
                         P[j]=0;
-                for (j=0;j<m;j++)
+                for (j=0;j<maxQ;j++)
                         Q[j]=0;
                 tim.clear();      
@@ -79 +79 @@
                 timc+=tim;
         }
-        write_field (F,cerr<<"Result = "<<endl, A, m,n,n);
+        //write_field (F,cerr<<"Result = "<<endl, A, m,n,n);
 
         cerr<<"P = [";
@@ -116 +116 @@
                 }
                 
-                write_field(F,cerr<<"L = "<<endl,L,m,m,m);
-                write_field(F,cerr<<"U = "<<endl,U,m,n,n);
+                // 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);
                 for ( int i=0; i<m; ++i )
                         F.assign(*(L+i*(m+1)), one);
 
-                write_field(F,cerr<<"L = "<<endl,L,m,m,m);
-                write_field(F,cerr<<"U = "<<endl,U,m,n,n);
+//              write_field(F,cerr<<"L = "<<endl,L,m,m,m);
+//              write_field(F,cerr<<"U = "<<endl,U,m,n,n);
                 if (diag == FFLAS::FflasNonUnit)
                         for ( int i=0; i<R; ++i )
@@ -134 +134 @@
                         }
                 }
-                write_field(F,cerr<<"L = "<<endl,L,m,m,m);
-                write_field(F,cerr<<"U = "<<endl,U,m,n,n);
+//              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, U, n, P);
@@ -167 +167 @@
                                 F.assign( *(U+i+j*n), zero);
                 }
-                write_field(F,cerr<<"L = "<<endl,L,m,n,n);
-                write_field(F,cerr<<"U = "<<endl,U,n,n,n);
+//              write_field(F,cerr<<"L = "<<endl,L,m,n,n);
+//              write_field(F,cerr<<"U = "<<endl,U,n,n,n);
 
                 FFPACK::applyP( F, FFLAS::FflasLeft, FFLAS::FflasTrans, n,0,R, U, n, Q);
@@ -185 +185 @@
                         }
                 }
-                write_field(F,cerr<<"L = "<<endl,L,m,n,n);
-                write_field(F,cerr<<"U = "<<endl,U,n,n,n);
+//              write_field(F,cerr<<"L = "<<endl,L,m,n,n);
+//              write_field(F,cerr<<"U = "<<endl,U,n,n,n);
 
                 FFPACK::applyP (F, FFLAS::FflasLeft, FFLAS::FflasTrans, n,0,R, L, n, P);
@@ -199 +199 @@
                                 fail=true;
         
-        write_field(F,cerr<<"X = "<<endl,X,m,n,n);
-        write_field(F,cerr<<"B = "<<endl,B,m,n,n);
+//      write_field(F,cerr<<"X = "<<endl,X,m,n,n);
+//      write_field(F,cerr<<"B = "<<endl,B,m,n,n);
         delete[] B;
         if (fail)