Scribd
Importer
3 vues1 page

Algorithme du Gradient Conjugué

L'algorithme Sequential_CG est une méthode itérative pour résoudre des systèmes d'équations linéaires de la forme Ax = b. Il commence par initialiser un vecteur x0 et calcule le résidu r0, puis itère jusqu'à ce que la norme du résidu soit inférieure à un seuil ε. À chaque itération, il met à jour les valeurs de x, r et p en fonction des calculs basés sur les résidus et les directions de recherche.

Transféré par

jessica.akono
Copyright
© All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd
3 vues1 page

Algorithme du Gradient Conjugué

L'algorithme Sequential_CG est une méthode itérative pour résoudre des systèmes d'équations linéaires de la forme Ax = b. Il commence par initialiser un vecteur x0 et calcule le résidu r0, puis itère jusqu'à ce que la norme du résidu soit inférieure à un seuil ε. À chaque itération, il met à jour les valeurs de x, r et p en fonction des calculs basés sur les résidus et les directions de recherche.

Transféré par

jessica.akono
Copyright
© All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd

Algorithme: Sequential_CG (A, b, 𝜀)

Début
Initialiser x0 dans IRn;

𝑟₀ = 𝑏 – 𝐴𝑥₀;

𝑝₀ = - 𝑟₀;

k = 0;

Tantque ‖𝑟ₖ‖ > 𝜀 faire


αₖ = (𝑟ₖᵀ𝑟ₖ) / (𝑝ₖᵀ𝐴𝑝ₖ);

𝑥ₖ₊₁ = 𝑥ₖ + αₖ𝑝ₖ;

𝑟ₖ₊₁ = 𝑟ₖ − αₖ𝐴𝑝ₖ;

βₖ = (𝑟ₖ₊₁ᵀ𝑟ₖ₊₁) / (𝑟ₖᵀ𝑟ₖ);

𝑝ₖ₊₁ = - 𝑟ₖ₊₁ + βₖ𝑝ₖ;

𝑘 = 𝑘 + 1;
Fintque
Retourner xk-1

Fin

Vous aimerez peut-être aussi