[Ipopt] Compiling ipopt with ma77
Mehmet Ersin YUMER
meyumer at gmail.com
Wed Jan 23 09:43:51 EST 2013
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> 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>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/34724279/attachment.html>
More information about the Ipopt
mailing list