Changeset 1042 in Sophya

Timestamp:

Jun 8, 2000, 4:55:52 PM (25 years ago)

Author:

ansari

Message:

Modif + correction LapackLinSolve, Reza 8/6/2000

Location:

trunk/SophyaExt/LinAlg

Files:

• intflapack.cc (modified) (3 diffs)
• intflapack.h (modified) (1 diff)

trunk/SophyaExt/LinAlg/intflapack.cc

@@ -15 +15 @@
 
 template <class T>
-TArray<T>& LapackServer<T>::LinSolve(TArray<T>& a, TArray<T> & b)
+int LapackServer<T>::LinSolve(TArray<T>& a, TArray<T> & b)
 {
   if ( ( a.NbDimensions() != 2 ) || ( b.NbDimensions() != 2 ) )
@@ -26 +26 @@
   if ( a.Size(rowa) !=  a.Size(cola)) 
     throw(SzMismatchError("LapackServer::LinSolve(a,b) a Not a square Array"));
-  if ( a.Size(rowa) !=  b.Size(rowa)) 
+  if ( a.Size(rowa) !=  b.Size(rowb)) 
     throw(SzMismatchError("LapackServer::LinSolve(a,b) RowSize(a <> b) "));
 
@@ -56 +56 @@
   }
   delete[] ipiv;
-  return(b);
+  return(info);
 }
 

trunk/SophyaExt/LinAlg/intflapack.h

@@ -10 +10 @@
 class LapackServer { 
 public:
-  TArray<T>& LinSolve(TArray<T>& a, TArray<T> & b); 
+  int LinSolve(TArray<T>& a, TArray<T> & b); 
 };
 
 template <class T>
-inline TArray<T>& LapackLinSolve(TArray<T>& a, TArray<T> & b)
+inline int LapackLinSolve(TArray<T>& a, TArray<T> & b)
   { LapackServer<T> lps; return( lps.LinSolve(a, b) );  }