[Ipopt-tickets] [Ipopt] #14: Termination code 11 - when using ampl and Pardiso linear solver

[Ipopt]

#14: Termination code 11 -  when using ampl and Pardiso linear solver
------------------------------+---------------------------------------------
Reporter:  rodrigo            |       Owner:  ipopt-team
    Type:  defect             |      Status:  new       
Priority:  high               |   Component:  Ipopt     
 Version:  3.0 - C++ Version  |    Severity:  major     
Keywords:                     |  
------------------------------+---------------------------------------------
 I've been trying to get faster solutions from an dynamic optimization
 problem involving a reactive distillation column for my graduate research
 project. I've obtained some solutions with MA27, but it takes very long so
 I got an academic license for the Pardiso solvers to try it out. I got the
 solver from
 http://www.computational.unibas.ch/cs/scicomp/software/pardiso/.

 I compiled IPOPT with

 {{{
 $ ./configure -with-pardiso="/home/rush/lib/libpardiso_GNU_IA32.so"
 $ make ; make install
 }}}


 I added the license file for pardiso, and added/created the
 PARDISO_LIC_PATH variable. Also, I'm using AMPL, as an interface to IPOPT.

 I'm using linux (CentOS 4.3 www.centos.org) with gcc version 3.4.5
 20051201 (Red Hat 3.4.5-2). With a P4 with Hyper Threading enabled and 2Gb
 of RAM. I'm using kernel version 2.6.9-34.0.1.ELsmp.

 I'm using the following options file:

 {{{
 max_iter 5000
 output_file output.log
 tol 1e-6

 linear_solver pardiso
 mu_strategy adaptive
 linear_system_scaling mc19
 constr_mult_init_max 10000
 }}}


 After a some iterations I get:


 {{{
 Ipopt 3.1.0: halt_on_ampl_error yes
 ******************************************************************************
 This program contains Ipopt, a library for large-scale nonlinear
 optimization.
  Ipopt is released as open source code under the Common Public License
 (CPL).
          For more information visit http://projects.coin-or.org/Ipopt
 ******************************************************************************

 Number of nonzeros in equality constraint Jacobian...:   166317
 Number of nonzeros in inequality constraint Jacobian.:        0
 Number of nonzeros in Lagrangian Hessian.............:   165840

 Total number of variables............................:    19233
                      variables with only lower bounds:        0
                 variables with lower and upper bounds:    14013
                      variables with only upper bounds:        0
 Total number of equality constraints.................:    19113
 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
   80  4.0848058e+05 2.19e-02 4.02e+05  -4.0 1.96e-01   0.0 1.00e+00 1.55e-
 03h  1
   81  3.9552385e+05 2.44e-01 2.62e+05  -3.4 1.68e+00  -4.0 1.00e+00 4.04e-
 01f  1
   82  3.8613839e+05 6.64e-01 4.05e+03  -3.6 8.99e-01  -4.0 1.00e+00
 1.00e+00h  1
   83  3.9630983e+05 4.79e-01 9.46e+04  -2.4 5.99e-01  -4.0 6.99e-01 5.35e-
 01f  1
   84  4.0972022e+05 4.54e-03 2.20e+04  -2.6 1.22e-01   0.0 8.66e-01
 1.00e+00h  1
   85  4.1001553e+05 1.09e-05 2.37e+00  -2.7 1.63e-03   2.0 1.00e+00
 1.00e+00h  1
 error running ipopt:
         termination code 11
 <BREAK>
 }}}


 I think this might be a bug.... so, please help.

 Thanks all,
   Rodrigo Lopez-Negrete

-- 
Ticket URL: <https://projects.coin-or.org/Ipopt/ticket/14>
Ipopt <http://www.coin-or.org>
Ipopt Tickets