[Ipopt] Compiling ipopt with ma77

Jonathan Hogg jonathan.hogg at stfc.ac.uk
Wed Jan 23 10:11:42 EST 2013


It could well be a bug in the MA77 interface - it hasn't had that much 
use. Can you supply us with a code that reproduces the issue?

Regards,

Jonathan.

On 23/01/13 14:43, Mehmet Ersin YUMER wrote:
>
> 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.
>
> On Jan 23, 2013 3:13 AM, "Mostafa.B" <bakhtvar at gmail.com 
> <mailto:bakhtvar at gmail.com>> wrote:
>
>     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 <mailto: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 <mailto: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 <mailto:Ipopt at list.coin-or.org>
>                 http://list.coin-or.org/mailman/listinfo/ipopt
>
>
>
>
>
>         -- 
>         Mehmet Ersin Yumer
>
>         _______________________________________________
>         Ipopt mailing list
>         Ipopt at list.coin-or.org <mailto:Ipopt at list.coin-or.org>
>         http://list.coin-or.org/mailman/listinfo/ipopt
>
>
>
>
> _______________________________________________
> Ipopt mailing list
> Ipopt at list.coin-or.org
> http://list.coin-or.org/mailman/listinfo/ipopt


-- 
Scanned by iCritical.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://list.coin-or.org/pipermail/ipopt/attachments/20130123/72c3cb54/attachment-0001.html>


More information about the Ipopt mailing list