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

Timestamp:

Aug 11, 2010, 2:58:45 PM (15 years ago)

Author:

ansari

Message:

correction ds lpk (pour least_square_solve), Reza 11/08/2010

File:

• trunk/SophyaProg/Tests/lpk.cc (modified) (1 diff)

trunk/SophyaProg/Tests/lpk.cc

@@ -208 +208 @@
   PrtTim(" Calling LapackLeastSquareSolve() ");
   LapackLeastSquareSolve(aa, bx);
+  Vector solvx(n);
+  solvx = bx(Range(0,n-1));
   PrtTim(" End LapackLeastSquareSolve() ");
-  bx.Share(bx(Range(0,0,n)));
-  if (prtlev > 0) {
-    cout << " ------------ Result X  = \n " << bx << "\n" << endl; 
-    cout << " ------------ X-X_Real  = \n " << bx-x << "\n" << endl; 
-  }
-  Vector diff = b-a*bx;
-  PrtTim(" End of Compute(diff)");
-  if (prtlev > 0) 
-    cout << " ------------ Matrix diff b-a*x = \n " << diff << "\n" << endl; 
+  //  bx.Share(bx(Range(0,0,n)));
+  if (prtlev > 0) {
+    cout << " ------------ Result X  = \n " << solvx << "\n" << endl; 
+    cout << " ------------ X-X_Real  = \n " << solvx-x << "\n" << endl; 
+  }
+  Vector diff = b-a*solvx;
+  PrtTim(" End of Compute(diff)");
+  if (prtlev > 0) 
+    cout << " ------------ Vector diff b-a*x = \n " << diff << "\n" << endl; 
   double mean,sigma;
   MeanSigma(diff, mean, sigma);