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

Timestamp:

Feb 7, 2005, 5:42:27 PM (21 years ago)

Author:

ansari

Message:

Ajout methode calcul de la matrice inverse par Lapack en utilisant la resolution de systeme + MAJ doc - Reza 7 Fev 2005

File:

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

trunk/SophyaExt/LinAlg/intflapack.h

@@ -18 +18 @@
   virtual int LeastSquareSolve(TArray<T>& a, TArray<T> & b); 
   virtual int LeastSquareSolveSVD_DC(TMatrix<T>& a,TMatrix<T>& b,TVector<r_8>& s,int_4& rank,r_8 rcond=-1.); 
+
+  // Calcul de la matrice inverse en utilisant la resolution de syst. lineaire
+  virtual int ComputeInverse(TMatrix<T>& a, TMatrix<T>& ainv);
 
   virtual int SVD(TArray<T>& a, TArray<T> & s); 
@@ -66 +69 @@
 inline int LapackLeastSquareSolve(TArray<T>& a, TArray<T> & b)
 { LapackServer<T> lps; return( lps.LeastSquareSolve(a, b) );  }
+
+/*! \ingroup LinAlg
+    \fn LapackInverse(TMatrix<T>&)
+    \brief Computes the inverse matrix using linear system solver LapackServer::LinSolve. 
+*/
+template <class T>
+inline TMatrix<T> LapackInverse(TMatrix<T>& a)
+{ LapackServer<T> lps; TMatrix<T> ainv; lps.ComputeInverse(a, ainv);  return ainv; }
 
 /*! \ingroup LinAlg