<p dir="ltr">I don't think it is the reason. The problem is well behaved and convex. All other linear solvers except ma77 are able to solve it correctly.</p>
<div class="gmail_quote">On Jan 23, 2013 3:13 AM, "Mostafa.B" <<a href="mailto:bakhtvar@gmail.com">bakhtvar@gmail.com</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">I had a similar problem with other linear sparse solvers, ipopt failing after only one iteration, however it turned out that the problems was in my optimization model, talking more specifically about my case, it was a division by zero resulting from bad placing of variables, have a check on your model anyways.</div>
<div class="gmail_extra"><br><br><div class="gmail_quote">On Tue, Jan 22, 2013 at 7:25 PM, Mehmet Ersin YUMER <span dir="ltr"><<a href="mailto:meyumer@gmail.com" target="_blank">meyumer@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Stefan,<div><br></div><div>Thanks for your reply. I am still having problems but the problem changed a little. I am able to successfully build ipopt trunk version with ma77 but cannot solve the optimization problem with it, here is a breakdown:</div>
<div><br></div><div>- I checked out the trunk version with svn, and placed a coinhsl directory into $[ipopt_dir]/ThirdParty/HSL/ whose contents are the contents of coinhsl-2012.12.21.tar.gz.</div><div><br></div>
<div>- run ./get.Metis and ./get.Mumps from their directories.</div><div><br></div><div>- Create build directory in $[ipopt_dir].</div><div><br></div><div>- Run ../configure, make, make test, make install from within $[ipopt_dir]/build/</div>
<div><br></div><div>Using the newly created libraries and include files in the build directory, I am able to solve a test problem succesfully with ma27, ma57, and mumps. However, when I try solving the problem with ma77, I get the output below, and Ipopt terminates. As can be understood from above, I did not use any flags while configuring Ipopt. Let me know if you can guess what the problem might be with ma77.</div>
<div><br></div><div>Thanks.</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div>
<p><b>******************************************************************************</b></p>
<p><b>This program contains Ipopt, a library for large-scale nonlinear optimization.</b></p>
<p><b> Ipopt is released as open source code under the Eclipse Public License (EPL).</b></p>
<p><b> For more information visit <a href="http://projects.coin-or.org/Ipopt" target="_blank">http://projects.coin-or.org/Ipopt</a></b></p>
<p><b>******************************************************************************</b></p>
<p><b></b><br></p>
<p><b>This is Ipopt version trunk, running with linear solver ma77.</b></p>
<p><b></b><br></p>
<p><b>Number of nonzeros in equality constraint Jacobian...: 6008</b></p>
<p><b>Number of nonzeros in inequality constraint Jacobian.: 0</b></p>
<p><b>Number of nonzeros in Lagrangian Hessian.............: 158992</b></p>
<p><b></b><br></p>
<p><b>Total number of variables............................: 6008</b></p>
<p><b> variables with only lower bounds: 0</b></p>
<p><b> variables with lower and upper bounds: 6008</b></p>
<p><b> variables with only upper bounds: 0</b></p>
<p><b>Total number of equality constraints.................: 759</b></p>
<p><b>Total number of inequality constraints...............: 0</b></p>
<p><b> inequality constraints with only lower bounds: 0</b></p>
<p><b> inequality constraints with lower and upper bounds: 0</b></p>
<p><b> inequality constraints with only upper bounds: 0</b></p>
<p><b></b><br></p>
<p><b>iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls</b></p>
<p><b> 0 4.0613987e+02 9.58e-03 5.36e+01 0.0 0.00e+00 - 0.00e+00 0.00e+00 0</b></p>
<p><b>WARNING: Problem in step computation; switching to emergency mode.</b></p>
<p><b> 1r 4.0613987e+02 9.58e-03 9.01e+02 0.0 0.00e+00 - 0.00e+00 0.00e+00R 1</b></p>
<p><b>WARNING: Problem in step computation; switching to emergency mode.</b></p>
<p><b>Restoration phase is called at point that is almost feasible,</b></p>
<p><b> with constraint violation 0.000000e+00. Abort.</b></p>
<p><b>Restoration phase in the restoration phase failed.</b></p>
<p><b></b><br></p>
<p><b>Number of Iterations....: 1</b></p>
<p><b></b><br></p>
<p><b> (scaled) (unscaled)</b></p>
<p><b>Objective...............: 4.0613987440593672e+02 4.0613987440593672e+02</b></p>
<p><b>Dual infeasibility......: 5.3643302289191944e+01 5.3643302289191944e+01</b></p>
<p><b>Constraint violation....: 9.5803704550541102e-03 9.5803704550541102e-03</b></p>
<p><b>Complementarity.........: 9.9000001999999998e-01 9.9000001999999998e-01</b></p>
<p><b>Overall NLP error.......: 5.3643302289191944e+01 5.3643302289191944e+01</b></p>
<p><b></b><br></p>
<p><b></b><br></p>
<p><b>Number of objective function evaluations = 2</b></p>
<p><b>Number of objective gradient evaluations = 2</b></p>
<p><b>Number of equality constraint evaluations = 2</b></p>
<p><b>Number of inequality constraint evaluations = 0</b></p>
<p><b>Number of equality constraint Jacobian evaluations = 2</b></p>
<p><b>Number of inequality constraint Jacobian evaluations = 0</b></p>
<p><b>Number of Lagrangian Hessian evaluations = 2</b></p>
<p><b>Total CPU secs in IPOPT (w/o function evaluations) = 0.144</b></p>
<p><b>Total CPU secs in NLP function evaluations = 0.059</b></p>
<p><b></b><br></p>
<p><b>EXIT: Restoration Failed!</b></p>
<p><b></b><br></p>
<p><b></b><br></p>
<p><b>*** The problem FAILED!</b></p></div><div class="gmail_extra"><div><br><br><div class="gmail_quote">On Mon, Jan 21, 2013 at 9:14 PM, Stefan Vigerske <span dir="ltr"><<a href="mailto:stefan@math.hu-berlin.de" target="_blank">stefan@math.hu-berlin.de</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">Hi,<br>
<br>
Ipopt 3.10 does not support MA77 yet.<br>
Ipopt/trunk has been setup to work with MA77, you could try that out.<br>
<br>
Stefan<div><div><br>
<br>
On 01/22/2013 01:38 AM, Mehmet Ersin YUMER wrote:<br>
</div></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div><div>
Hi,<br>
<br>
Up to now I have compiled and used ipopt with ma27 & ma19. Because of an<br>
increase in problem size, I now need to use it with ma77. However, I could<br>
not build ipopt with it. It seems like it is not able to locate the file in<br>
$[ipopt_dir]/Thirdparty/hsl/<br>
<br>
Should the file name be different than: ma77ad.f90 ?<br>
<br>
Source files I am using for ipopt are from version 3.10.0.<br>
<br>
Thanks.<br>
<br>
<br>
<br></div></div>
______________________________<u></u>_________________<br>
Ipopt mailing list<br>
<a href="mailto:Ipopt@list.coin-or.org" target="_blank">Ipopt@list.coin-or.org</a><br>
<a href="http://list.coin-or.org/mailman/listinfo/ipopt" target="_blank">http://list.coin-or.org/<u></u>mailman/listinfo/ipopt</a><br>
<br>
</blockquote>
<br>
</blockquote></div><br><br clear="all"><div><br></div></div><span><font color="#888888">-- <br>Mehmet Ersin Yumer
</font></span></div></div>
<br>_______________________________________________<br>
Ipopt mailing list<br>
<a href="mailto:Ipopt@list.coin-or.org" target="_blank">Ipopt@list.coin-or.org</a><br>
<a href="http://list.coin-or.org/mailman/listinfo/ipopt" target="_blank">http://list.coin-or.org/mailman/listinfo/ipopt</a><br>
<br></blockquote></div><br></div>
</blockquote></div>