Changeset 1248 in Sophya

Timestamp:

Oct 19, 2000, 5:05:20 PM (25 years ago)

Author:

ansari

Message:

Amelioration programme test lapack - Reza 19/10/2000

File:

• trunk/SophyaProg/Tests/lpk.cc (modified) (4 diffs)

trunk/SophyaProg/Tests/lpk.cc

@@ -11 +11 @@
 
 
-void lpk_tarr(int n);
-void lpk_tmtx(int n);
+int lpk_tarr(int n);
+int lpk_tmtx(int n);
 
+static double TOLERANCE = 1.e-6;
 
 int main(int narg, char* arg[])
@@ -38 +39 @@
   if (narg > 4) nprt = atoi(arg[4]);
 
+  int rc = 0;
   BaseArray::SetMaxPrint(nprt, prtlev);
   try {
-    if (opt == 's') lpk_tmtx(n);
-    else lpk_tarr(n);
+    if (opt == 's') rc = lpk_tmtx(n);
+    else rc = lpk_tarr(n);
   }
   catch (PThrowable exc) {
     cerr << " catched Exception (lpk.cc) " << exc.Msg() << endl;
+    rc = 77;
   }  
   catch (...) {
     cerr << " catched unknown (...) exception (lpk.cc) " << endl; 
+    rc = 78; 
   } 
   
   PrtTim(" End of lpk LinAlg/Lapack test ");
   cout << " ---------------  END of Programme -------------- " << endl;
+  return(rc);
 }
 
 
-void lpk_tmtx(int n)
+int lpk_tmtx(int n)
 {
   int i,j,k;
@@ -84 +89 @@
   Vector diff = b-x;
   cout << " ------------ Vector diff B-X = \n " << diff << "\n" << endl; 
+  double min,max;
+  diff.MinMax(min, max);
+  cout << " Min/Max difference Vector (?=0) , Min= " << min 
+       << " Max= " << max << endl;
+  if ((fabs(min) > TOLERANCE) || (fabs(max) > TOLERANCE)) {
+    cout << " !!! Difference exceeding tolerance (=" << TOLERANCE << ") !!!" 
+         << endl;
+    return(99);
+  } 
+  return(0);
 
 }
 
-void lpk_tarr(int n)
+int  lpk_tarr(int n)
 {
   int i,j,k;
@@ -114 +129 @@
 
   cout << " ------------ Result B(=X ?) = \n " << b << "\n" << endl; 
+  return(0);
 }