ENSEEIHT/TP-optimisation-numerique-2
1
0
Fork 0
Code Issues Pull requests Packages Projects Releases Wiki Activity

add try/catch

Browse source
This commit is contained in:
relmo 2020-06-16 21:11:17 +02:00
parent 135f3d42dc
commit aa947a8403
2 changed files with 1 additions and 4 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

4
src/Lagrangien_Augmente.jl
View file

@ -18,7 +18,6 @@ Lagrangien_Augmente(algo,fonc,contrainte,gradfonc,hessfonc,grad_contrainte,
* **hessfonc** : la hessienne de la fonction * **hessfonc** : la hessienne de la fonction
* **grad_contrainte** : le gradient de la contrainte * **grad_contrainte** : le gradient de la contrainte
* **hess_contrainte** : la hessienne de la contrainte * **hess_contrainte** : la hessienne de la contrainte
* **phi(x)** : utilisée dans le calcul du gradient du Lagrangien (égale 0 dans le cas des contraintes d'égalités)
* **x0** : la première composante du point de départ du Lagrangien * **x0** : la première composante du point de départ du Lagrangien
* **options** * **options**
1. **epsilon** : utilisé dans les critères d'arrêt 1. **epsilon** : utilisé dans les critères d'arrêt
@ -48,12 +47,11 @@ options = []
contrainte(x) = (x[1]^2) + (x[2]^2) -1.5 contrainte(x) = (x[1]^2) + (x[2]^2) -1.5
grad_contrainte(x) = [2*x[1] ;2*x[2]] grad_contrainte(x) = [2*x[1] ;2*x[2]]
hess_contrainte(x) = [2 0;0 2] hess_contrainte(x) = [2 0;0 2]
phi(x) = 0
output = Lagrangien_Augmente(algo,f,contrainte,gradf,hessf,grad_contrainte,hess_contrainte,phi,x0,options) output = Lagrangien_Augmente(algo,f,contrainte,gradf,hessf,grad_contrainte,hess_contrainte,phi,x0,options)
``` ```
""" """
function Lagrangien_Augmente(algo,fonc::Function,contrainte::Function,gradfonc::Function, function Lagrangien_Augmente(algo,fonc::Function,contrainte::Function,gradfonc::Function,
hessfonc::Function,grad_contrainte::Function,hess_contrainte::Function,phi::Function,x0,options) hessfonc::Function,grad_contrainte::Function,hess_contrainte::Function,x0,options)
if options == [] if options == []
epsilon = 1e-8 epsilon = 1e-8

1
test/runtests.jl
View file

@ -10,7 +10,6 @@ include("../src/Lagrangien_Augmente.jl")
include("../src/Pas_De_Cauchy.jl") include("../src/Pas_De_Cauchy.jl")
include("../src/Regions_De_Confiance.jl") include("../src/Regions_De_Confiance.jl")
# Tester l'algorithme de Newton # Tester l'algorithme de Newton
printstyled("Lancement des tests sur ",bold=true,color=:white) printstyled("Lancement des tests sur ",bold=true,color=:white)
printstyled("l'algorithme de Newton ",bold=true,color=:blue) printstyled("l'algorithme de Newton ",bold=true,color=:blue)