Mathc matrices/a17

Sommaire

---

Total Pivoting

Dans cet algorithme le pivot choisi est la plus grande valeur absolue qui se trouve dans la matrice A. Pour cela on fait des échanges de lignes et de colonnes. Attention on ne touche pas à b.

La fonction GJ_TP_mR(Ab); a été introduite pour traiter les deux types de systèmes que la fonction gj_TP_mR(Ab); est capable de résoudre.

L'algorithme consiste à vérifier si la matrice A n'est pas singulière. Si elle n'est pas singulière on appelle la fonction gj_TP_mR(Ab); Si la matrice n'est pas carré, on prend une sous matrice de A qui a r lignes et r colonnes ou r est le nombre de lignes de A pour calculer le déterminant. (voir la fonction c_Ab_subArxr_mR(Ab,A);

Copier la bibliothèque dans votre répertoire de travail avec les fichiers des parties précédentes :

• v_a.h ..................... Déclaration des fichiers h
  • vgj1tp.h ................... Gauss Jordan 1
  • vgc_coef.h .............. Copier des coéfficients
  • vgc_free.h .............. Introduire les variables libres
  • vgp_free.h .............. Imprimer les variables libres

Exemples d'applications :

• gj_r1.c ........... Résoudre : Ax = b; A[r][c] .... :c = r
• gj_r2.c ........... Résoudre : Ax = b; A[r][c] .... :c > r

• gj_r3.c ........... Les variables libres
• gj_r4.c ........... Vérifier avec octave

---