Changeset 1342 in Sophya for trunk/SophyaExt/LinAlg/intflapack.h

Timestamp:

Nov 24, 2000, 10:47:37 AM (25 years ago)

Author:

ansari

Message:

Ajout de SVD ds LapackServer - Reza 24/11/2000

File:

• trunk/SophyaExt/LinAlg/intflapack.h (modified) (1 diff)

trunk/SophyaExt/LinAlg/intflapack.h

@@ -10 +10 @@
 class LapackServer { 
 public:
-  int LinSolve(TArray<T>& a, TArray<T> & b); 
+  LapackServer();
+  virtual ~LapackServer();
+
+  virtual int LinSolve(TArray<T>& a, TArray<T> & b); 
+  virtual int SVD(TArray<T>& a, TArray<T> & s); 
+  virtual int SVD(TArray<T>& a, TArray<T> & s, TArray<T> & u, TArray<T> & vt); 
+
+  inline void SetWorkSpaceSizeFactor(int f = 2)
+  { wspace_size_factor = (f > 1) ? f : 1; }
+  inline int  GetWorkSpaceSizeFactor() 
+  { return wspace_size_factor; }
+
+private:
+  int SVDDriver(TArray<T>& a, TArray<T> & s, 
+                TArray<T>* up=NULL, TArray<T> * vtp=NULL);
+
+  int wspace_size_factor;
 };
 
 template <class T>
 inline int LapackLinSolve(TArray<T>& a, TArray<T> & b)
-  { LapackServer<T> lps; return( lps.LinSolve(a, b) );  }
+{ LapackServer<T> lps; return( lps.LinSolve(a, b) );  }
+
+template <class T>
+inline int LapackSVD(TArray<T>& a, TArray<T> & s)
+{ LapackServer<T> lps; return( lps.SVD(a, s) ); }
+
+template <class T>
+inline int LapackSVD(TArray<T>& a, TArray<T> & s, TArray<T> & u, TArray<T> & vt)
+{ LapackServer<T> lps; return( lps.SVD(a, s, u, vt) ); }