Skip to content

rapport du TP6 - #4

Open
bouayoubegrine wants to merge 5 commits into
dpo-mth8408:mainfrom
bouayoubegrine:TP6
Open

bouayoubegrine wants to merge 5 commits into
dpo-mth8408:mainfrom
bouayoubegrine:TP6

Conversation

@bouayoubegrine

Copy link
Copy Markdown

No description provided.

@github-actions

Copy link
Copy Markdown

Status:
Success--Here-is-the-PDF

Comment thread rapport.qmd
t = 1.0
xt = projection!(xk + t .* p, C)
errorsquare = dot(xt - xk, xt - xk)
while obj(nlp, xt) > obj(nlp, xk) - 1.0e-4 * errorsquare / t # alpha = 1.0e-4

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Le errorsquare dans le membre de droite doit être $\Vert x_k - x(t_k) \Vert^2$ mais, tu ne le mets pas à jour. Mais à la fin, il faut quand-même renvoyer $\Vert x_k - x(1) \Vert^2$.

Comment thread rapport.qmd
p = -grad(nlp,xk)
println("Iter\t\t\t||xk - xk(1)||\t\t\t\t objectif")
println(@sprintf("%d\t\t\t%s\t\t\t\t%.3e",
k, " - ", obj(nlp, xk)))

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Finalement, les colonnes ne sont pas vraiment bien alignées... 😢

Comment thread rapport.qmd
C = SimpleBounds(model3.meta.lvar, model3.meta.uvar)
x3 = projected_grad(model3, C)
stats3 = ipopt(model3, print_level=0)
println("Erreur relatif d'écart entre gradient projeté est Ipopt: ", norm(x3 - stats3.solution)/norm(stats3.solution))

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Ok

Comment thread rapport.qmd
x .= C.c .+ (C.r/nd) .* d
end
return x
end

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

👍

Comment thread rapport.qmd

## Explication

- Cette méthode peut être utilisée pour le calcul du pas dans les méthodes de région de confiance, puisqu’on impose la condition d’Armijo, laquelle assure la décroissance de la fonction objectif — point essentiel.

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Ce n'est pas suffisant.

Comment thread rapport.qmd

- Cette méthode peut être utilisée pour le calcul du pas dans les méthodes de région de confiance, puisqu’on impose la condition d’Armijo, laquelle assure la décroissance de la fonction objectif — point essentiel.

- L’avantage de cette méthode est qu’elle permet de trouver un point stationnaire sous certains hypothèses du sous-problème de région de confiance à une tolérance près . À l’inverse, dans la méthode du gradient conjugué, dès qu’un itéré sort de la boule, on projette l’avant-dernier itéré sur la boule, ce qui n’est ni la solution du sous-problème de région de confiance ni même un point stationnaire.

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

On ne le projette pas ; il est dans la boule.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants