[Ipopt] Precompiled IPOPT MATLAB Mex File Crash

Anil Rao anilvrao at gmail.com
Sun Oct 28 21:54:37 EDT 2012


I am running a problem using the precompiled MATLAB Mac OS-X mex file that
is compiled with the trunk version of IPOPT and the linear solver MUMPS.
In the process of solving my problem, IPOPT crashes MATLAB, giving me the
MATLAB dialog box that says, "MATLAB has encountered an internal problem
and needs to close."  The iteration output from IPOPT is shown below.  The
final objective function shown in the iteration output is the correct one
for this problem, but you will see that the following message is generated
just prior to the crash of MATLAB:  "Restoration phase is called at almost
feasible point, but acceptable point from iteration 92 could be restored."
As it turns out, I get a MATLAB crash every time I see this last message
REGARDLESS of whether I am using a mex file I have compiled myself or if I
am using the precompiled mex file.  I am wondering if there is a bug
somewhere in either the IPOPT source code or the Matlab mex interface.  It
seems odd to me that IPOPT would hum along just fine and would crash MATLAB
when it enters the aforementioned restoration phase.

Regards,

Anil Rao

-- 
Anil V. Rao
Associate Professor
Department of Mechanical and Aerospace Engineering
University of Florida
Gainesville, FL 32611-6250
Tel:  (352) 392-5523 (office); (352) 672-1529 (cell)
E-mail:  anilvrao at gmail.com (primary); anilvrao at ufl.edu (secondary)


******************************************************************************
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
******************************************************************************

NOTE: You are using Ipopt by default with the MUMPS linear solver.
      Other linear solvers might be more efficient (see Ipopt
documentation).


This is Ipopt version trunk, running with linear solver mumps.

Number of nonzeros in equality constraint Jacobian...:     1402
Number of nonzeros in inequality constraint Jacobian.:        1
Number of nonzeros in Lagrangian Hessian.............:      783

Total number of variables............................:      198
                     variables with only lower bounds:        0
                variables with lower and upper bounds:      198
                     variables with only upper bounds:        0
Total number of equality constraints.................:      160
Total number of inequality constraints...............:        1
        inequality constraints with only lower bounds:        0
   inequality constraints with lower and upper bounds:        1
        inequality constraints with only upper bounds:        0

   0  1.0699999e+02 7.18e-01 8.76e-01   0.0 0.00e+00    -  0.00e+00
0.00e+00   0
   1  1.1255327e+02 5.73e-01 3.64e+02   0.5 4.40e+00    -  9.25e-01
1.80e-01f  1
   2  1.1380862e+02 2.71e-01 5.62e+02   0.8 6.82e-01   2.0 1.00e+00
5.15e-01h  1
   3  1.3499130e+02 1.82e-01 4.76e+02  -5.3 6.28e+00    -  2.37e-01
4.82e-01h  1
   4  2.4094031e+02 9.45e-01 1.78e+03   0.9 1.51e+01    -  5.12e-01
1.00e+00f  1
   5  2.4101437e+02 4.12e-01 1.13e+03   0.7 4.06e-01   2.4 1.00e+00
9.60e-01h  1
   6  2.4174477e+02 2.84e-01 7.95e+02   0.5 4.45e-01   1.9 4.50e-01
3.78e-01f  2
   7  2.4215610e+02 2.36e-01 5.62e+03   1.3 4.33e-01   3.3 9.60e-01
6.36e-01h  1
   8  2.4821337e+02 2.14e-01 4.34e+03   1.0 6.56e+00    -  2.11e-01
1.32e-01h  1
   9  2.6251810e+02 1.57e-01 2.99e+03  -0.7 8.17e+00    -  3.75e-01
2.50e-01h  3
iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du
alpha_pr  ls
  10  2.6782805e+02 1.53e-01 2.73e+03   0.8 1.09e+01    -  2.52e-01
6.97e-02f  1
  11  3.1435221e+02 9.50e-01 1.69e+03   0.5 1.18e+01    -  2.95e-01
5.65e-01f  1
  12  3.5382332e+02 1.31e-01 1.66e+03   0.8 5.64e+00    -  9.09e-01
1.00e+00f  1
  13  3.5917148e+02 9.98e-03 9.11e+01   0.0 7.64e-01    -  9.39e-01
1.00e+00f  1
  14  3.3459479e+02 1.25e-01 1.98e+01  -0.1 3.51e+00    -  9.54e-01
1.00e+00f  1
  15  3.1112934e+02 8.73e-02 5.45e+00  -0.5 3.35e+00    -  1.00e+00
1.00e+00f  1
  16  3.1609087e+02 4.78e-02 5.92e+00  -1.1 7.09e-01    -  9.89e-01
1.00e+00h  1
  17  3.1949829e+02 1.49e-02 9.08e-01  -1.4 4.87e-01    -  9.96e-01
1.00e+00h  1
  18  3.2007105e+02 4.47e-03 1.96e+00  -2.6 1.00e-01    -  1.00e+00
8.16e-01h  1
  19  3.2038388e+02 1.64e-03 1.17e-01  -3.2 4.47e-02    -  1.00e+00
1.00e+00h  1
iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du
alpha_pr  ls
  20  3.2038331e+02 2.72e-04 2.33e+00  -4.7 9.19e-03    -  1.00e+00
8.41e-01h  1
  21  3.2038674e+02 3.00e-06 2.27e+00  -6.3 1.52e-03    -  1.00e+00
1.00e+00h  1
  22  3.2038464e+02 7.84e-07 2.21e+00  -7.9 8.58e-04    -  1.00e+00
1.00e+00h  1
  23  3.2038462e+02 2.39e-07 1.25e-01  -9.4 1.57e-05    -  1.00e+00
6.96e-01h  1
  24  3.2038463e+02 7.35e-11 1.63e-03 -10.1 7.36e-07    -  1.00e+00
1.00e+00h  1
  25  3.2038463e+02 1.64e-15 7.09e-06 -11.0 1.80e-09    -  1.00e+00
1.00e+00h  1
  26  3.2038463e+02 1.30e-15 5.71e-06 -11.0 1.40e-10    -  1.00e+00
1.00e+00h  1
  27  3.2038463e+02 1.57e-15 5.36e-06 -11.0 2.25e-10    -  1.00e+00
1.00e+00h  1
  28  3.2038463e+02 1.39e-15 5.36e-06 -11.0 2.85e-10    -  1.00e+00
1.00e+00H  1
  29  3.2038463e+02 9.42e-16 4.11e-06 -11.0 1.18e-10    -  1.00e+00
1.00e+00h  1
iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du
alpha_pr  ls
  30  3.2038463e+02 1.97e-15 4.46e-06 -11.0 1.24e-10    -  1.00e+00
1.00e+00h  1
  31  3.2038463e+02 2.03e-15 1.02e-05 -11.0 6.31e-11    -  1.00e+00
1.00e+00h  1
  32  3.2038463e+02 1.56e-15 8.29e-06 -11.0 1.90e-10    -  1.00e+00
5.00e-01h  2
  33  3.2038463e+02 1.18e-15 1.05e-05 -11.0 1.81e-10    -  1.00e+00
1.00e+00h  1
  34  3.2038463e+02 1.63e-15 1.96e-06 -11.0 1.69e-10    -  1.00e+00
3.12e-02h  6
  35  3.2038463e+02 1.31e-15 5.36e-06 -11.0 1.07e-10    -  1.00e+00
1.00e+00h  1
  36  3.2038463e+02 1.37e-15 6.41e-06 -11.0 4.15e-10    -  1.00e+00
1.00e+00H  1
  37  3.2038463e+02 1.20e-15 1.02e-05 -11.0 2.10e-10    -  1.00e+00
2.50e-01h  3
  38  3.2038463e+02 1.47e-15 5.58e-06 -11.0 3.11e-10    -  1.00e+00
1.00e+00H  1
  39  3.2038463e+02 1.60e-15 8.31e-06 -11.0 2.92e-10    -  1.00e+00
6.10e-05h 15
iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du
alpha_pr  ls
  40  3.2038463e+02 1.56e-15 4.49e-06 -11.0 2.12e-10    -  1.00e+00
3.81e-06h 19
  41  3.2038463e+02 1.60e-15 8.93e-06 -11.0 3.12e-10    -  1.00e+00
3.81e-06h 19
  42  3.2038463e+02 2.06e-15 3.21e-06 -11.0 2.11e-10    -  1.00e+00
1.00e+00H  1
  43  3.2038463e+02 1.62e-15 3.40e-06 -11.0 1.19e-10    -  1.00e+00
6.25e-02h  5
  44  3.2038463e+02 1.43e-15 2.52e-06 -11.0 2.46e-10    -  1.00e+00
3.12e-02h  6
  45  3.2038463e+02 1.90e-15 5.36e-06 -11.0 1.13e-10    -  1.00e+00
1.00e+00h  1
  46  3.2038463e+02 1.40e-15 6.06e-06 -11.0 2.41e-10    -  1.00e+00
1.00e+00H  1
  47  3.2038463e+02 1.93e-15 1.44e-06 -11.0 1.03e-10    -  1.00e+00
1.00e+00H  1
  48  3.2038463e+02 2.03e-15 3.80e-06 -11.0 1.95e-10    -  1.00e+00
1.00e+00h  1
  49  3.2038463e+02 1.50e-15 7.58e-06 -11.0 1.62e-10    -  1.00e+00
5.00e-01h  2
iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du
alpha_pr  ls
  50  3.2038463e+02 1.99e-15 1.05e-05 -11.0 2.60e-10    -  1.00e+00
2.50e-01h  3
  51  3.2038463e+02 2.45e-15 1.02e-05 -11.0 1.36e-10    -  1.00e+00
1.00e+00H  1
  52  3.2038463e+02 1.57e-15 4.89e-06 -11.0 2.00e-10    -  1.00e+00
5.00e-01h  2
  53  3.2038463e+02 1.65e-15 4.91e-06 -11.0 1.04e-10    -  1.00e+00
1.56e-02h  7
  54  3.2038463e+02 9.54e-16 5.15e-06 -11.0 1.36e-10    -  1.00e+00
1.00e+00H  1
  55  3.2038463e+02 1.68e-15 1.43e-05 -11.0 2.60e-10    -  1.00e+00
1.00e+00h  1
  56  3.2038463e+02 1.49e-15 7.86e-06 -11.0 1.28e-10    -  1.00e+00
5.00e-01h  2
  57  3.2038463e+02 8.75e-16 9.11e-06 -11.0 8.05e-11    -  1.00e+00
5.00e-01h  2
  58  3.2038463e+02 2.86e-15 1.46e-05 -11.0 1.08e-10    -  1.00e+00
1.00e+00H  1
  59  3.2038463e+02 1.79e-15 8.21e-06 -11.0 8.21e-11    -  1.00e+00
1.00e+00h  1
iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du
alpha_pr  ls
  60  3.2038463e+02 1.97e-15 1.32e-05 -11.0 1.49e-10    -  1.00e+00
1.25e-01h  4
  61  3.2038463e+02 1.09e-15 5.36e-06 -11.0 1.61e-10    -  1.00e+00
3.12e-02h  6
  62  3.2038463e+02 1.78e-15 1.01e-05 -11.0 8.98e-11    -  1.00e+00
1.00e+00H  1
  63  3.2038463e+02 2.03e-15 5.91e-06 -11.0 9.23e-11    -  1.00e+00
1.00e+00h  1
  64  3.2038463e+02 1.17e-15 1.04e-05 -11.0 9.45e-11    -  1.00e+00
5.00e-01h  2
  65  3.2038463e+02 2.45e-15 9.37e-06 -11.0 1.22e-10    -  1.00e+00
1.00e+00h  1
  66  3.2038463e+02 1.81e-15 5.17e-06 -11.0 1.33e-10    -  1.00e+00
1.00e+00h  1
  67  3.2038463e+02 1.32e-15 5.59e-06 -11.0 1.73e-10    -  1.00e+00
7.81e-03h  8
  68  3.2038463e+02 1.47e-15 8.57e-06 -11.0 1.79e-10    -  1.00e+00
5.00e-01h  2
  69  3.2038463e+02 1.30e-15 1.09e-05 -11.0 2.43e-10    -  1.00e+00
1.56e-02h  7
iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du
alpha_pr  ls
  70  3.2038463e+02 1.78e-15 3.78e-06 -11.0 2.27e-10    -  1.00e+00
1.56e-02h  7
  71  3.2038463e+02 1.81e-15 1.41e-05 -11.0 3.91e-10    -  1.00e+00
1.00e+00h  1
  72  3.2038463e+02 1.58e-15 1.47e-05 -11.0 2.98e-10    -  1.00e+00
6.25e-02h  5
  73  3.2038463e+02 1.44e-15 1.37e-05 -11.0 1.90e-10    -  1.00e+00
1.56e-02h  7
  74  3.2038463e+02 2.20e-15 3.20e-06 -11.0 2.94e-10    -  1.00e+00
3.12e-02h  6
  75  3.2038463e+02 2.03e-15 8.62e-06 -11.0 1.40e-10    -  1.00e+00
1.00e+00h  1
  76  3.2038463e+02 2.57e-15 2.12e-06 -11.0 3.28e-10    -  1.00e+00
1.00e+00h  1
  77  3.2038463e+02 1.54e-15 1.47e-07 -11.0 2.14e-10    -  1.00e+00
6.25e-02h  5
  78  3.2038463e+02 1.79e-15 4.11e-06 -11.0 1.06e-10    -  1.00e+00
1.25e-01h  4
  79  3.2038463e+02 1.79e-15 9.27e-06 -11.0 1.84e-10    -  1.00e+00
3.81e-06h 19
iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du
alpha_pr  ls
  80  3.2038463e+02 1.17e-15 4.49e-06 -11.0 1.04e-10    -  1.00e+00
1.25e-01h  4
  81  3.2038463e+02 1.09e-15 8.56e-06 -11.0 1.75e-10    -  1.00e+00
9.77e-04h 11
  82  3.2038463e+02 1.28e-15 4.14e-06 -11.0 2.37e-10    -  1.00e+00
5.00e-01h  2
  83  3.2038463e+02 1.88e-15 9.35e-06 -11.0 1.58e-10    -  1.00e+00
5.00e-01h  2
  84  3.2038463e+02 1.47e-15 1.05e-05 -11.0 1.97e-10    -  1.00e+00
1.25e-01h  4
  85  3.2038463e+02 1.58e-15 9.26e-06 -11.0 1.55e-10    -  1.00e+00
1.00e+00H  1
  86  3.2038463e+02 1.41e-15 4.81e-06 -11.0 1.42e-10    -  1.00e+00
1.00e+00h  1
  87  3.2038463e+02 1.49e-15 3.55e-06 -11.0 2.13e-10    -  1.00e+00
1.25e-01h  4
  88  3.2038463e+02 1.50e-15 5.70e-06 -11.0 9.55e-11    -  1.00e+00
3.12e-02h  6
  89  3.2038463e+02 1.41e-15 4.51e-06 -11.0 2.62e-10    -  1.00e+00
1.25e-01h  4
iter    objective    inf_pr   inf_du lg(mu)  ||d||  lg(rg) alpha_du
alpha_pr  ls
  90  3.2038463e+02 1.15e-15 1.24e-05 -11.0 2.20e-10    -  1.00e+00
6.25e-02h  5
  91  3.2038463e+02 9.47e-16 8.30e-06 -11.0 2.47e-10    -  1.00e+00
5.00e-01h  2
  92  3.2038463e+02 2.49e-15 6.99e-07 -11.0 1.55e-10    -  1.00e+00
1.00e+00h  1
  93  3.2038463e+02 1.54e-15 1.39e-06 -11.0 3.53e-10    -  1.00e+00
2.50e-01h  3
  94  3.2038463e+02 1.43e-15 6.26e-06 -11.0 1.39e-10    -  1.00e+00
1.00e+00H  1
  95  3.2038463e+02 1.22e-15 8.28e-06 -11.0 5.58e-11    -  1.00e+00
1.00e+00h  1
  96  3.2038463e+02 1.81e-15 8.28e-06 -11.0 2.78e-10    -  1.00e+00
1.00e+00h  1
Restoration phase is called at almost feasible point,
  but acceptable point from iteration 92 could be restored.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://list.coin-or.org/pipermail/ipopt/attachments/20121028/19365b38/attachment.html>


More information about the Ipopt mailing list