Changeset 1343 in Sophya for trunk/SophyaProg/Tests/lpk.cc

Timestamp:

Nov 24, 2000, 10:48:18 AM (25 years ago)

Author:

ansari

Message:

Ajout test LapackServer::SVD() ds lpk.cc - Reza 24/11/2000

File:

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

trunk/SophyaProg/Tests/lpk.cc

@@ -11 +11 @@
 
 
-int lpk_tarr(int n);
-int lpk_tmtx(int n);
+int lpk_rzt_arr(int n);
+int lpk_linsolve_mtx(int n);
+int lpk_svd_mtx(int l, int c, int wsf, bool covu=true);
 
 static double TOLERANCE = 1.e-6;
@@ -25 +26 @@
 
   if (narg < 2) {
-    cout << " lpk - LinAlg/LapackServer test - Usage lpk t/s [size=5] [prtlev=0] [nprtmax=100]\n" 
-         << " t:  lpk_tarr() rztest_lapack with TArray<r_4> \n" 
-         << " s:  lpk_tmtx() LapackServer with TMatrix<r_8> \n" << endl;
+    cout << " lpk - LinAlg/LapackServer test \n" 
+         << " Usage: lpk linsolve/svd/rzt [sizeL,C=5,5] [prtlev=0] \n" 
+         << "             [nprtmax=100] [WorkSpaceSizeFactor=2] \n" 
+         << "   linsolve:  lpk_linsolve_mtx() LapackServer::LinSolve with TMatrix<r_8> \n" 
+         << "   svd:  lpk_svd_mtx() LapackServer::SVD(a,s,u,vt) with TMatrix<r_8> \n" 
+         << "   svds:  lpk_svd_mtx() LapackServer::SVD(a,s); with TMatrix<r_8> \n" 
+         << "   rzt:  lpk_rzt_tarr() rztest_lapack with TArray<r_4> \n" << endl;
     exit(0);
   }
-  int n = 5;
-  int opt = 's';
-  opt = *(arg[1]);
-  if (narg > 2) n = atoi(arg[2]);
+  int l,c;
+  l = c = 5;
+  int wsf = 2;
+  string opt = arg[1];
+  if (narg > 2) sscanf(arg[2], "%d,%d", &l, &c); 
   int nprt = 100;
   int prtlev = 1;
   if (narg > 3) prtlev = atoi(arg[3]);
   if (narg > 4) nprt = atoi(arg[4]);
+  if (narg > 5) wsf = atoi(arg[5]);
 
   int rc = 0;
   BaseArray::SetMaxPrint(nprt, prtlev);
   try {
-    if (opt == 's') rc = lpk_tmtx(n);
-    else rc = lpk_tarr(n);
+    if (opt == "linsolve") rc = lpk_linsolve_mtx(l);
+    else if (opt == "svd") rc = lpk_svd_mtx(l,c,wsf,true);
+    else if (opt == "svds") rc = lpk_svd_mtx(l,c,wsf,false);
+    else if (opt == "rzt") rc = lpk_rzt_arr(l);
+    else { cout << " Unknown option " << opt << " ! " << endl; rc = 66; }
   }
   catch (PThrowable exc) {
@@ -55 +65 @@
   
   PrtTim(" End of lpk LinAlg/Lapack test ");
-  cout << " ---------------  END of Programme -------------- " << endl;
+  cout << " ---------------  END of Programme -------- (Rc= " 
+       << rc << ") --- " << endl;
   return(rc);
 }
 
-
-int lpk_tmtx(int n)
+// -----------------------------------------------------------------------
+/* Nouvelle-Fonction */
+int lpk_linsolve_mtx(int n)
 {
-  int i,j,k;
+  int i,j;
   BaseArray::SetDefaultMemoryMapping(BaseArray::FortranMemoryMapping);
+  cout << " lpk_linsolve_mtx() - Test of LapackServer::LinSolve()  " << endl;
   Matrix a(n,n);
   for(i=0; i<n; i++)
     for(j=0; j<n; j++)  a(j,i) = GauRnd(0., 1.);
-  cout << " ------------ Matrix A = \n " << a << "\n" << endl; 
   
   Vector x(n), b;
@@ -75 +87 @@
   b = a*x;
 
-  //  cout << ":::::::: rztest_lapack - Size=" << n << "  ::::::::: " << endl;
-  cout << " ------- lpk_tmtx() LapackServerTest Using TMatrix<r_8> -------- " << endl; 
+  cout << " ---- lpk_tmtx() LapackServer::LinSolve Test Using TMatrix<r_8> ----- " << endl; 
   cout << " ------------ Matrix A = \n " << a << "\n" << endl; 
   cout << " ------------ Matrix X = \n " << x << "\n" << endl; 
@@ -88 +99 @@
   cout << " ------------ Result B(=X ?) = \n " << b << "\n" << endl; 
   Vector diff = b-x;
+  PrtTim(" End of Compute(diff)");
   cout << " ------------ Vector diff B-X = \n " << diff << "\n" << endl; 
   double min,max;
@@ -102 +114 @@
 }
 
-int  lpk_tarr(int n)
+// -----------------------------------------------------------------------
+/* Nouvelle-Fonction */
+int lpk_svd_mtx(int m, int n, int wsf, bool covu)
 {
-  int i,j,k;
+  BaseArray::SetDefaultMemoryMapping(BaseArray::FortranMemoryMapping);
+  cout << " lpk_svd_mtx() - Test of LapackServer::SVD " << endl;
+  
+  Matrix a(m , n), aa;
+  a = RandomSequence(RandomSequence::Gaussian, 0., 4.);
+  aa = a;
+  Vector s;
+  Matrix u, vt;
+  cout << " ---- lpk_svd_tmtx() LapackServer::SVD Test Using TMatrix<r_8> ---- " << endl; 
+  cout << " ------------ Matrix A = \n " << a << "\n" << endl; 
+
+  cout << "\n   Calling LapackSVD(a,s,u,vt) .... " << endl;
+  PrtTim(" Calling LapackSVD() ");
+  LapackServer<r_8> lpks;
+  lpks.SetWorkSpaceSizeFactor(wsf);
+  if (covu) lpks.SVD(aa, s, u, vt);
+  else lpks.SVD(aa, s);
+  PrtTim(" End LapackSVD() ");
+
+  cout << " ------------ Result S  = \n " << s << "\n" << endl; 
+  if (!covu) return(0);
+  cout << " ------------ Result U  = \n " << u << "\n" << endl; 
+  cout << " ------------ Result VT = \n " << vt << "\n" << endl; 
+  Matrix sm(m,n);
+  int minmn = (m<n) ? m : n ;
+  for(int k=0; k< minmn ; k++) sm(k,k) = s(k);
+  Matrix diff = u*(sm*vt) - a;
+  PrtTim(" End of Compute(diff)");
+  cout << " ------------ Matrix diff U*S*Vt - A = \n " << diff << "\n" << endl; 
+  double min,max;
+  diff.MinMax(min, max);
+  cout << " Min/Max difference Matrix (?=0) , Min= " << min 
+       << " Max= " << max << endl;
+  if ((fabs(min) > TOLERANCE) || (fabs(max) > TOLERANCE)) {
+    cout << " !!! Difference exceeding tolerance (=" << TOLERANCE << ") !!!" 
+         << endl;
+    return(99);
+  } 
+  return(0);
+
+}
+
+int  lpk_rzt_arr(int n)
+{
+  int i,j;
   TArray<r_4> a(n,n);
   for(i=0; i<n; i++)