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

Timestamp:

Jul 23, 2004, 12:53:35 PM (21 years ago)

Author:

cmv

Message:

add SVD decomp by Divide and Conquer (cmv 23/07/04)

File:

• trunk/SophyaExt/LinAlg/intflapack.h (modified) (2 diffs)

trunk/SophyaExt/LinAlg/intflapack.h

@@ -20 +20 @@
   virtual int SVD(TArray<T>& a, TArray<T> & s); 
   virtual int SVD(TArray<T>& a, TArray<T> & s, TArray<T> & u, TArray<T> & vt);
+  virtual int SVD_DC(TMatrix<T>& a, TVector<T>& s, TMatrix<T>& u, TMatrix<T>& vt);
  
   virtual int LapackEigenSym(TArray<T>& a, TVector<r_8>& b, bool eigenvector=true);
@@ -83 +84 @@
 
 /*! \ingroup LinAlg
+    \fn LapackSVD_DC(TMatrix<T>&, TVector<T>&, TMatrix<T>&, TMatrix<T>&)
+    \brief SVD decomposition DC using LapackServer. 
+*/
+template <class T>
+inline int LapackSVD_DC(TMatrix<T>& a, TVector<T>& s, TMatrix<T>& u, TMatrix<T>& vt)
+{ LapackServer<T> lps; return( lps.SVD_DC(a, s, u, vt) ); }
+
+
+/*! \ingroup LinAlg
     \fn LapackEigenSym(TArray<T>&, TArray<T> &)
     \brief Compute the eigenvalues and eigenvectors of A (symetric or hermitian).