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

Timestamp:

Jul 21, 2004, 5:52:36 PM (21 years ago)

Author:

cmv

Message:

• Introduction de l'interface Lapack d'inversion des matrices symetriques
• Introduction de l'interface Lapack de recherche de valeurs et vecteurs propres (cas general, symetrique et hermitique)
• Introduction d'un fonction d'interface pour le calculateur de workspace (ilaenv_)
• Commentaires sur les diverses methodes et sur les matrices FORTRAN
• Pour tester cf Tests/tsttminv.cc

(cmv, 21/07/04)

File:

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

trunk/SophyaExt/LinAlg/intflapack.h

@@ -4 +4 @@
 #include "machdefs.h"
 #include "tarray.h"
+#include "tvector.h"
 
 namespace SOPHYA {
@@ -18 +19 @@
 
   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(TArray<T>& a, TArray<T> & s, TArray<T> & u, TArray<T> & vt);
+ 
+  virtual int LapackEigenSym(TArray<T>& a, TVector<r_8>& b, bool eigenvector=true);
+  virtual int LapackEigen(TArray<T>& a, TVector< complex<r_8> >& eval, TMatrix<T>& evec, bool eigenvector);
 
   //! Set the workspace size factor for LAPACK routines
@@ -78 +82 @@
 
 
+/*! \ingroup LinAlg
+    \fn LapackEigenSym(TArray<T>&, TArray<T> &)
+    \brief Compute the eigenvalues and eigenvectors of A (symetric or hermitian). 
+*/
+template <class T>
+inline int LapackEigenSym(TArray<T>& a, TVector<r_8>& b, bool eigenvector=true)
+{ LapackServer<T> lps; return( lps.LapackEigenSym(a,b,eigenvector) );  }
+
+/*! \ingroup LinAlg
+    \fn LapackEigen(TArray<T>&, TArray<T> &)
+    \brief Compute the eigenvalues and (right) eigenvectors of A (general square matrix). 
+*/
+template <class T>
+inline int LapackEigen(TArray<T>& a, TVector< complex<r_8> >& eval, TMatrix<T>& evec, bool eigenvector=true)
+{ LapackServer<T> lps; return( lps.LapackEigen(a,eval,evec,eigenvector) );  }
+
 } // Fin du namespace