<br><font size=2 face="sans-serif">Hello,</font>
<br>
<br><font size=2 face="sans-serif">I have an optimization problem with
4 parameters having to be optimized. I suppose the hessian of Lagrangian
is dense that's why I set nele_hess = 4 * (4 + 1) / 2 = 10. Morevoer, I
have 2 equality constraints, that's why I set nele_jac = 4 * 2 = 8</font>
<br>
<br><font size=2 face="sans-serif">If I use my own hessian fonction, called
"eval_h", I obtain a great result. Unfortunately, If I use the
option "hessian_approximation" = "limite-memory", after
several iterations, I have an error saying "the memory can't be read".</font>
<br>
<br><font size=2 face="sans-serif">My program crashes inside the Ipopt
fonction </font><font size=2>UpdateHessian</font><font size=2 face="sans-serif">()</font>
<br>
<br><font size=2 face="sans-serif">My code is the following (I'm using
the C interface) :</font>
<br>
<br><font size=2 color=#008000>// create the IpoptProblem</font>
<br><font size=2>nlp = CreateIpoptProblem(n, x_L, x_U, m, g_L, g_U, nele_jac,
nele_hess,index_style, eval_f, eval_g, eval_grad_f, eval_jac_g, eval_h);</font>
<br>
<br><font size=2>AddIpoptStrOption(nlp,</font><font size=2 color=#a11f12>"hessian_approximation"</font><font size=2>,
"limited-memory");</font>
<br>
<br><font size=2>status = IpoptSolve(nlp, x, g, &obj, mult_g,mult_x_L,
mult_x_U, user_data);</font>
<br>
<br><font size=2>Do you know where could be the problem ? Have you got
any ideas ?</font>
<br>
<br><font size=2>Thank you very much.</font>
<br>
<br><font size=2>Best regards,</font>
<br>
<br><font size=2>Pierre-Lin Pommier</font>
<br>
<br><font size=2 face="sans-serif">_______________________________________________________<br>
Pierre-Lin POMMIER<br>
Manufacture Française des Pneumatiques Michelin <br>
CTE/DTO/SIM/ET CER Ladoux - Bâtiment F32 – 2ème étage
<br>
23 place des Carmes-Dechaux 63040 Clermont-Ferrand Cedex 9 <br>
Tel int : 67 194
Fax int
: 68 541 <br>
Tel ext : (+33) (0)4 73 10 71 94 Fax ext : (+33) (0)4
73 10 85 41 <br>
mail : pierre-lin.pommier@fr.michelin.com<br>
<br>
Confidentialité D3<br>
_______________________________________________________</font>