Discussion :

[ClaireA]

Bonsoir tout le monde

pour mon TP de calcul scientifique, je dois écrire un programme (en fortran) qui factorise une matrice A tridiagonale en LU. J'ai donc voulu stocker chacune des trois diagonales dans un vecteur (respectivement up, low et diag). Par contre, au moment de la compilation j'ai des erreurs du type:
error: rank mismatch in array reference
lorsque j'avais écrit par exemple
do i=1,n
diag(i)=A(i,i)
end do.

Je débute en fortran, pourriez vous me donner une piste pour comprendre mon erreur ?
Peut être que j'ai mal utilisé allocate.
Merci beaucoup !

[FR119492]

Salut!
Pour pouvoir te répondre, il nous faudrait voir comment tu as déclaré tes variables.
Jean-Marc Blanc

[Bloodagain]

Comment tu as déclaré "A" dans ton programme?

Sinon tu peux demander au gens de la classe, ils ne mordent pas ...

[genteur slayer]

attention toute fois: après factorisation LU, une matrice tridiagonale ne reste pas nécéssairement tridiagonale... une méthode classique consiste à déclarer A comme pleine et de la modifiée en LU: L est une matrice triangulaire inférieur avec des 1 sur la diag (on stoke pa les "1" on se souvien seullement qu'ils sont là) et on stoque le tout dans la partie triangulaire inférieur de A. U7 est une matrice triangulaire suppérieur (avec valeur sur la diag) que l'on stoke dans la partie triangulaire suppérieur de A (diag incluse cette fois)

[FR119492]

Salut!

je dois écrire un programme (en fortran) qui factorise une matrice A tridiagonale en LU

1. Pourquoi? Que veux-tu faire avec tes matrices L et U?
2. Tes matrices L et U sont-elles définies par la relation L+U=A ou par L*U=A ?
3. Ta matrice A est-elle symétrique, et, si oui, est-elle définie positive?
4. Stocker une matrice tridiagonale dans un tableau à 2 indices est un gaspillage pur et simple de place.

Quand tu auras répondu à ces questions, on pourra aller un peu plus loin.
Jean-Marc Blanc

[genteur slayer]

1. Pourquoi? Que veux-tu faire avec tes matrices L et U?
2. Tes matrices L et U sont-elles définies par la relation L+U=A ou par L*U=A ?
3. Ta matrice A est-elle symétrique, et, si oui, est-elle définie positive?
4. Stocker une matrice tridiagonale dans un tableau à 2 indices est un gaspillage pur et simple de place.

c'est marrant: c'est toujours les mêmes questions (pertinentes) qui reviennent au fil des topics...