Algorithme de Gauss-Jordan de Wikipédia : un truc qui cloche ?

Alors du coup normalement c’est bon !

1
2
3
4
5
6
7
8

content.indices.foreach(id_line => {
            if(id_line != id_last_pivot) {
              val tmp = mutable_being_inversed_matrix(id_line)(general_id_column)
              content.indices.foreach(id_column => {
                mutable_being_inversed_matrix(id_line)(id_column) -= mutable_being_inversed_matrix(id_last_pivot)(id_column) * tmp
              })
            }
          })

J’ai rajouté une variable $tmp$ pour garder $A[i,j]$ à la même valeur.

Résultats :

Avec une petite matrice : https://jsfiddle.net/0pe4xab6/

Avec la matrice du topic : https://jsfiddle.net/dhgLqe5y/

La matrice de Wikipédia marche !

Petite question. Maintenant que j’ai réussi à "convertir" la matrice originale en une matrice identité de la bonne façon, je voudrais "convertir" la matrice identité en la matrice inversée, résultat concret de l’élimination de Gauss-Jordan.

Je me contente d’appliquer les mêmes opérations que sur la matrice originale ?

Oui.

Mais genre la division qui permet d’obtenir les 1, je ne la fais pas sur ma matrice identité quand même ? Enfin si je ne me trompe pas elle ne sert à rien pour elle.

Ben… Bien sûr que si…

EDIT : je suis en train de me dire que tu n’as pas compris l’algo en fait, il s’agit d’effectuer les mêmes opérations sur ta matrice à inverser (disons $A$) et sur la matrice $I$ puisqu’en fait, l’algo de Gauss-Jordan est équivalent à multiplier par $A^{-1}$. Si tu Multiplies $A$ par $A^{-1}$ tu obtiens $I$ (c’est la partie que tu as déjà codée), et si tu multiplies $I$ par $A^{-1}$, tu obtiens $A^{-1}$. Il faut effectuer exactement les même opérations sur tes deux matrices pour que l’algo fonctionne, sinon ça n’aurait aucune raison de donner quoique ce soit de particulier…

hello Adri1, désolé de répondre avec un petit peu de retard (j’ai édité l’OP et mis en résolu ce topic mais c’est tout ce que j’ai eu le temps de faire rapidement).

Effectivement, j’ai manqué de recul et fait preuve d’étourderie quand j’ai posté mon dernier message (je parlais de la matrice identité alors que dans ma tête je pensais à l’autre matrice bref, grosse grosse confusion sur le moment).

Par contre depuis le début j’ai bien compris l’idée générale du procédé, là où j’ai eu un peu plus de mal c’est vers la dernière ligne de l’algo de wikipedia, mais là encore je l’avais bien comprise bien avant que j’ai créé ce topic.

A vrai dire, mon algo marchait depuis le début : j’avais juste oublié de stocker certaines valeurs de variables dans des variables "temporaires", pour éviter d’utiliser leur nouvelle valeur dans les calculs qui les impliquent dans une même boucle. C’est ce qui faisait foirer mon algo, mais sinon tout aurait dû marcher nickel depuis le début..................... Edit : d’ailleurs ça a bien été le cas quand j’ai corrigé ces oublis.

Question bête mais pourquoi en Scala ?

Salut @ache, c’est juste parce que ça s’inscrit dans un projet qui est réalisé en Scala!

Logique donc.

Je demandais par curiosité simplement.