Estimation du conditionnement d'une matrice symétrique indéfinie
L'algorithme qui estime le conditionnement d'une matrice symétrique définie positive est similaire à celui qui calcule le conditionnement d'une matrice quelconque. La seule différence réside dans l'utilisation de la décomposition LDLT au lieu de la décomposition LU. Par conséquent, nous n'en détaillerons pas les principes ici, car ils sont accessibles via le lien ci-dessus.
Ce module comprend deux sous-programmes : `SMatrixLDLTRCond` et `SMatrixRCond`. Le premier estime le conditionnement d'une matrice à partir de sa décomposition LDLT, tandis que le second estime le conditionnement d'une matrice dont la décomposition LDLT n'a pas encore été calculée. Les normes 1 et ∞ des matrices symétriques étant équivalentes, il n'existe pas de sous-programme distinct pour chaque norme.
Cet algorithme est issu de la bibliothèque LAPACK.