Changeset 51 for tests/test-fgesv.C

Timestamp:

09/26/07 17:13:25 (1 year ago)

Author:

pernet

Message:

Make test-fgesv work. No bugs known.

Files:

• tests/test-fgesv.C (modified) (6 diffs)

tests/test-fgesv.C

@@ -12 +12 @@
 #include <iomanip>
 #include <iostream>
+using namespace std;
+
 #include "fflas-ffpack/modular-balanced.h"
 #include "timer.h"
@@ -18 +20 @@
 
 
-using namespace std;
 
 typedef Modular<double> Field;
@@ -37 +38 @@
         F.init(zero,0.0);
         F.init(one,1.0);
-        Field::Element * A, *B, *B2, *X;
+        Field::Element * A, *B, *B2, *X=NULL;
         A = read_field(F,argv[2],&m,&n);
         B = read_field(F,argv[3],&mb,&nb);
@@ -43 +44 @@
         FFLAS::FFLAS_SIDE side = (atoi(argv[5])) ? FFLAS::FflasRight :  FFLAS::FflasLeft;
                 
-        size_t ldx;
+        size_t ldx=0;
         size_t rhs = (side == FFLAS::FflasLeft) ? nb : mb;
         if (m != n)
@@ -50 +51 @@
                         ldx = nb;
                 } else {
-                        X = new Field::Element[m*mb];
+                        X = new Field::Element[mb*m];
                         ldx = m;
                 }
         
-        if ( ((side == FFLAS::FflasRight) && (m != nb))
-             || ((side == FFLAS::FflasLeft)&&(n != mb)) ) {
+        if ( ((side == FFLAS::FflasRight) && (n != nb))
+             || ((side == FFLAS::FflasLeft)&&(m != mb)) ) {
                 cerr<<"Error in the dimensions of the input matrices"<<endl;
                 exit(-1);
@@ -71 +72 @@
                 else
                         R = FFPACK::fgesv (F, side, m, n, rhs, A, n, X, ldx, B, nb, &info);
+                if (info > 0){
+                        std::cerr<<"System is inconsistent"<<std::endl;
+                        exit(-1);
+                }
+                        
                 t.stop();
                 time+=t.usertime();