implémentation algo du pivot de Gauss

[221]

bonjour a tous,

j'essaye d'implémenter l'algo d'élimination par la méthode du pivot de gauss ,
j ai un problème avec la partie triangularisation de la matrice de mon programme, le débogueur n'indique aucune erreur mais le programme ne triangularise pas la matrice .

for (k=0; k<ligne; k++)
    {

                pivot=score[k][k];
                if (pivot!= 0){
                    for (i=k+1; i<ligne; i++)
                    {
                     for (j=0; j< colone; j++)
                        {
                        matrice[i][j]=matrice[i][j] - (matrice[k][j]/pivot)*matrice[i][j];
                        printf("%7d", matrice[i][j]);
                        }
                        printf("\n");
                    }
    }

le code est tiré de l algorithme de triangularisation :


merci,
-----

[Jack]

Je n'ai pas vérifié ton algo, mais d'une manière générale, tu peux utiliser un débogueur qui va te permettre d'exécuter ton programme pas-à-pas en visualisant dynamiquement l'état de tes variables (et de ta matrice par conséquent).

Sinon, à l'ancienne, tu appelles une fonction qui t'affiche ta matrice lorsque tu veux connaitre son état.

A+

[221]

y avait une erreur d affectation dans mon programme que j ai corrigé :

for (k=0; k<ligne; k++)
    {

                pivot=matrice[k][k];
                if (pivot!= 0){
                    for (i=k+1; i<ligne; i++)
                    {
                     for (j=0; j< colone; j++)
                        {
                        matrice[i][j]=matrice[i][j] - (matrice[k][j]/pivot)*matrice[i][j];
                        printf("%7d", matrice[i][j]);
                        }
                        printf("\n");
                    }
    }

ca ne triangularise tjr pas la matrice
merci de m éclairer

[221]

merci jack,
si quelqu un a un algo plus probant plus simple a implémenter il est le bienvenue

[A voir en vidéo sur Futura]

[Jack]

Mon dernier cours d'analyse numérique remontant à une trentaine d'années, je ne me sens plus vraiment capable de ce genre de performances

[221]

j comprends c est de l ordre du souvenir lointain x) .

le fait que le débogueur n'indique rien ne m aide pas trop, je vois vraiment pas ou ca déconne
je ne vois rien d illogique dans :
matrice[i][j]=matrice[i][j] - (matrice[k][j]/pivot)*matrice[i][j];

[jojo150393]

for (i=k+1; i<ligne; i++)
                    {
                     for (j=0; j< colone; j++)

C'est jusqu'à n et par n-1 (ici ligne et colonne). Remplace par <=

[jojo150393]

C'est normale que les indices de cette ligne :

matrice[i][j]=matrice[i][j] - (matrice[k][j]/pivot)*matrice[i][j]

ne correspondent pas aux indices de l'algo ?

[Jack]

j comprends c est de l ordre du souvenir lointain x) .

le fait que le débogueur n'indique rien ne m aide pas trop, je vois vraiment pas ou ca déconne
je ne vois rien d illogique dans :
matrice[i][j]=matrice[i][j] - (matrice[k][j]/pivot)*matrice[i][j];

Tu es sur de cette dernière ligne, parce que si on regarde l'algo que tu as donné, il me semble que c'est plutôt:
matrice[i][j]=matrice[i][j] - (matrice[k][j]/pivot)*matrice[i][k];

PS : grillé par jojo. Je n'avais pas vu ta réponse car j'avais du interrompre ma saisie pendant quelques minutes

[221]

merci jojo150393,
j ai pas vraiment suivi l algo question indices
enfaîte dans la ligne :

matrice[i][j]=matrice[i][j] - (matrice[k][j]/pivot)*matrice[i][j]

-matrice[k][j] est l élément j eme de ma linge K a savoir la ligne du pivot actuel ,pour chaque ligne on a un pivot donc k varie de 0 jusqu au nbr de ligne .
-le pivot de chaque ligne est l element matrice[k][k] qui varie aussi de 0 jusqu a nbr de ligne.
-matrice [i][j] est l élément j eme de la ligne i=k+1, ligne juste en dessous de la ligne du pivot ,il varie de i=k+1 jusqu a nbr ligne.
en gros j ai ca donne
nouvelle linge en dessous du pivot(éléments de la ligne)= éléments de la ligne en dessous du pivot -(éléments de la lignes du pivot /pivot lui meme)*éléments de la ligne du dessous

j espère que c est lisible

[jojo150393]

Je comprend pas désolé. Il faut plus de clarté ou on pourra pas t'aider.