[Ipopt] Precompiled IPOPT MATLAB Mex File Crash

Ian Washington washinid at mcmaster.ca
Mon Oct 29 15:34:33 EDT 2012


I had a similar problem and you probably need to set the matlab 
recognised environment variables

LAPACK_VERSION
BLAS_VERSION

to the location of blas/lapack, which you linked ipopt against. Read 
Tony Kelman's suggestions.

Ian.



On 10/28/2012 09:54 PM, Anil Rao wrote:
> 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 <mailto:anilvrao at gmail.com> (primary); 
> anilvrao at ufl.edu <mailto: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.
>
>
>
> _______________________________________________
> 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/20121029/3115a4a7/attachment-0001.html>


More information about the Ipopt mailing list