factorisation lu

[invite3569df15]

salut

je recherche du code (n'importe langage: c, c++, java...) pour l'algorithme de factorisation LU avec la notation compacte (autre nom connu: méthode de crout)

merci
-----

[martini_bird]

Salut,

tu ne peux pas obtenir ta factorisation simplement avec le pivot de Gauss?
Il y a ce fil qui pourrait peut-être t'intéresser?

Désolé de ne pas t'avoir aidé.

[martini_bird]

Sinon,

google m'a donné ceci (voir bas de page pour un algo en TP).

Cordialement.

[invite88e71a19]

La methode de Crout et la methode de Doolittle sont presque la meme (on pense dans l'une avec le ligne et avec l'autre aux colonnes mais c'est symmetrique).

Si on ecrit



Pour le Doolittle voilà l'algo




for i=1:n
for j=2:i
A(i,j-1)=A(i,j-1)/A(j-1,j-1)
for k=1:j-1
A(i,j)=A(i,j)-A(i,k)*A(k,j)
end
end
for j=i+1:n
for k=1:i-1
A(I,j)=A(i,j)-A(i,k)*A(k,j)
end
end
end

Crout c'est analogue:



avec u_{kk}=1

[A voir en vidéo sur Futura]

[Toufou]

Salut,

Un autre lien avec source en C/C++
http://perso.wanadoo.fr/jean-pierre...._matrices.html