[Ipopt] Compiling ipopt with ma77

[Mostafa.B]

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.


On Tue, Jan 22, 2013 at 7:25 PM, Mehmet Ersin YUMER <meyumer at gmail.com>wrote:

> Hi Stefan,
>
> 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:
>
> - 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.
>
> - run ./get.Metis and ./get.Mumps  from their directories.
>
> - Create build directory in $[ipopt_dir].
>
> - Run ../configure,  make,  make test,  make install   from within
> $[ipopt_dir]/build/
>
> 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.
>
> Thanks.
>
>
>
>
>
> *
> ******************************************************************************
> *
>
> *This program contains Ipopt, a library for large-scale nonlinear
> optimization.*
>
> * Ipopt is released as open source code under the Eclipse Public License
> (EPL).*
>
> *         For more information visit http://projects.coin-or.org/Ipopt*
>
> *
> ******************************************************************************
> *
>
> **
>
> *This is Ipopt version trunk, running with linear solver ma77.*
>
> **
>
> *Number of nonzeros in equality constraint Jacobian...:     6008*
>
> *Number of nonzeros in inequality constraint Jacobian.:        0*
>
> *Number of nonzeros in Lagrangian Hessian.............:   158992*
>
> **
>
> *Total number of variables............................:     6008*
>
> *                     variables with only lower bounds:        0*
>
> *                variables with lower and upper bounds:     6008*
>
> *                     variables with only upper bounds:        0*
>
> *Total number of equality constraints.................:      759*
>
> *Total number of inequality constraints...............:        0*
>
> *        inequality constraints with only lower bounds:        0*
>
> *   inequality constraints with lower and upper bounds:        0*
>
> *        inequality constraints with only upper bounds:        0*
>
> **
>
> *iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du
> alpha_pr  ls*
>
> *   0  4.0613987e+02 9.58e-03 5.36e+01   0.0 0.00e+00    -  0.00e+00
> 0.00e+00   0*
>
> *WARNING: Problem in step computation; switching to emergency mode.*
>
> *   1r 4.0613987e+02 9.58e-03 9.01e+02   0.0 0.00e+00    -  0.00e+00
> 0.00e+00R  1*
>
> *WARNING: Problem in step computation; switching to emergency mode.*
>
> *Restoration phase is called at point that is almost feasible,*
>
> *  with constraint violation 0.000000e+00. Abort.*
>
> *Restoration phase in the restoration phase failed.*
>
> **
>
> *Number of Iterations....: 1*
>
> **
>
> *                                   (scaled)                 (unscaled)*
>
> *Objective...............:   4.0613987440593672e+02
> 4.0613987440593672e+02*
>
> *Dual infeasibility......:   5.3643302289191944e+01
> 5.3643302289191944e+01*
>
> *Constraint violation....:   9.5803704550541102e-03
> 9.5803704550541102e-03*
>
> *Complementarity.........:   9.9000001999999998e-01
> 9.9000001999999998e-01*
>
> *Overall NLP error.......:   5.3643302289191944e+01
> 5.3643302289191944e+01*
>
> **
>
> **
>
> *Number of objective function evaluations             = 2*
>
> *Number of objective gradient evaluations             = 2*
>
> *Number of equality constraint evaluations            = 2*
>
> *Number of inequality constraint evaluations          = 0*
>
> *Number of equality constraint Jacobian evaluations   = 2*
>
> *Number of inequality constraint Jacobian evaluations = 0*
>
> *Number of Lagrangian Hessian evaluations             = 2*
>
> *Total CPU secs in IPOPT (w/o function evaluations)   =      0.144*
>
> *Total CPU secs in NLP function evaluations           =      0.059*
>
> **
>
> *EXIT: Restoration Failed!*
>
> **
>
> **
>
> **** The problem FAILED!*
>
>
> On Mon, Jan 21, 2013 at 9:14 PM, Stefan Vigerske <stefan at math.hu-berlin.de
> > wrote:
>
>> Hi,
>>
>> Ipopt 3.10 does not support MA77 yet.
>> Ipopt/trunk has been setup to work with MA77, you could try that out.
>>
>> Stefan
>>
>>
>> On 01/22/2013 01:38 AM, Mehmet Ersin YUMER wrote:
>>
>>> Hi,
>>>
>>> Up to now I have compiled and used ipopt with ma27 & ma19. Because of an
>>> increase in problem size, I now need to use it with ma77. However, I
>>> could
>>> not build ipopt with it. It seems like it is not able to locate the file
>>> in
>>> $[ipopt_dir]/Thirdparty/hsl/
>>>
>>> Should the file name be different than:   ma77ad.f90    ?
>>>
>>> Source files I am using for ipopt are from version 3.10.0.
>>>
>>> Thanks.
>>>
>>>
>>>
>>> ______________________________**_________________
>>> Ipopt mailing list
>>> Ipopt at list.coin-or.org
>>> http://list.coin-or.org/**mailman/listinfo/ipopt<http://list.coin-or.org/mailman/listinfo/ipopt>
>>>
>>>
>>
>
>
> --
> Mehmet Ersin Yumer
>
> _______________________________________________
> Ipopt mailing list
> Ipopt at list.coin-or.org
> http://list.coin-or.org/mailman/listinfo/ipopt
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://list.coin-or.org/pipermail/ipopt/attachments/20130123/ddfda7f7/attachment-0001.html>