Changeset 26 for tests/testeur_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/testeur_lqup.C (modified) (5 diffs)

tests/testeur_lqup.C

@@ -12 +12 @@
 using namespace std;
 //#include "fflas-ffpack/modular-int.h"
-#include "fflas-ffpack/modular-positive.h"
-//#include "fflas-ffpack/modular-balanced.h"
+//#include "fflas-ffpack/modular-positive.h"
+#include "fflas-ffpack/modular-balanced.h"
 #include "timer.h"
 #include "Matio.h"
@@ -82 +82 @@
                         P = new size_t[M];
                         Q = new size_t[N];
-
+                        for (size_t i=0; i<M; ++i) P[i] = 0;
+                        for (size_t i=0; i<N; ++i) Q[i] = 0;
                 }
                 else{
@@ -90 +91 @@
                         P = new size_t[N];
                         Q = new size_t[M];
-
+                        for (size_t i=0; i<N; ++i) P[i] = 0;
+                        for (size_t i=0; i<M; ++i) Q[i] = 0;
                 }
                 
@@ -209 +211 @@
                         FFPACK::applyP( F, FFLAS::FflasLeft, FFLAS::FflasTrans, N,0,R, U, N, Q);
                         for (size_t i=0; i<N; ++i)
-                        F.assign (*(U+i*(N+1)),one);
+                                F.assign (*(U+i*(N+1)),one);
                         if (diag == FFLAS::FflasNonUnit)
                                 for ( size_t i=0; i<R; ++i )
@@ -263 +265 @@
         for (size_t i=0; i<M; ++i)
                 for (size_t j=0; j<N; ++j)
-                        cerr<<i+1<<" "<<j+1<<" "<<((size_t) *(Abis+i*lda+j) )<<endl;
+                        if (!(*(Abis+i*lda+j)))
+                                cerr<<i+1<<" "<<j+1<<" "<<((size_t) *(Abis+i*lda+j) )<<endl;
         cerr<<"0 0 0"<<endl<<endl;