Lapack appelle à inverser une matrice
-
16-09-2019 - |
Question
De ma compréhension, il faut une décomposition / factorisation (LU, QR, Cholesky, etc.), suivi par le calcul inverse de la matrice basée sur la factorisation. Y a-t-il d'autres façons de se déplacer (je suis en train de comprendre si je peux rester avec les 6 fonctions données gratuitement dans la version de tryout CULAtools )? Merci à l'avance pour les commentaires.
La solution
Bien sûr, trouver la Comatrice; c'est un moyen simple d'inverser les petites matricies. La matrice de adjugate est que la transposée de la matrice de co-facteurs, et l'inverse d'une matrice carrée est que le adjugate divisé par le déterminant (scalaire). Consulter ces termes sur Wikipedia si elles ne sont pas familiers.
Si vous travaillez w / grandes matrices, j'acheter le paquet.
Paul
Autres conseils
Les routines LAPACK qui calculent l'inverse de la matrice sont xyyTRI
, où x
indique le type de données ( « S » pour les biens de simple précision, « D » pour les biens à double précision, « C » pour le complexe de simple précision, et « Z » pour complexe double précision) et yy
indique le type de matrice ( « GE » pour le cas général de matrices non symétriques, il y a plus de 20 autres codes à deux lettres pour d'autres types de matrice). Pour les matrices à valeurs réelles, vous devriez habituellement utiliser DGETRI
, et pour les matrices valeurs complexes, vous souhaitez généralement utiliser ZGETRI
.