Méthode du gradient conjugué

[MacManus]

Bonjour !

Dans l'algorithme du gradient conjugué, il faut au préalable se donner un vecteur de depart , et puis viennent ensuite les itérations .... pour trouver une solution du système linéaire Au=b, avec A symétrique def. pos.

Mais le problème est que je ne sais pas comment choisir ce vecteur de départ...!

Pouvez-vous me dire s'il y a une méthode simple pour le trouver ?

Merci bcp.

[fatal_error]

salut,

ben au pire t'en prends 1 au pif. mettons (1,1,...,1) (j'ai pas osé écrire 0,0,...,0, ca n'aurais pas beaucoup de sens...)
Apres si tu veux qu'il soit déjà dirigé vers ta solution, je sais pas.

[Sylviel]

Par défaut 1 * une norme de b / une norme de A c'est pas mal. (ou une idée de la norme).

[laya]

Si ta matrice est dans et que tu bosses avec la norme 1, tu peux commencer avec le vecteur de : .

[Sylviel]

J'aime pas trop dire (1,1,1,1) parce que si ton b est de l'ordre de 10^4 partout, alors que A a des coeff de l'ordre de l'unité, y'a des chances pour que tu sois loin. En même temps ça n'a pas beaucoup d'importance.

[fatal_error]

uiui, ca suit la logique de ton premier message

[MacManus]

oui fatal :) j'utilisais mon pif aussi.
bon je vais essayer ce que me dit Sylviel.

Merci en tout cas.