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. <br>
<br>Regards,<br><br>Anil Rao<br><br>-- <br>Anil V. Rao<br>
Associate Professor<br>
Department of Mechanical and Aerospace Engineering<br>
University of Florida<br>
Gainesville, FL 32611-6250<br>Tel: (352) 392-5523 (office); (352) 672-1529 (cell)<br>E-mail: <a href="mailto:anilvrao@gmail.com" target="_blank">anilvrao@gmail.com</a> (primary); <a href="mailto:anilvrao@ufl.edu" target="_blank">anilvrao@ufl.edu</a> (secondary)<br>
<br><br>******************************************************************************<br>This program contains Ipopt, a library for large-scale nonlinear optimization.<br> Ipopt is released as open source code under the Eclipse Public License (EPL).<br>
For more information visit <a href="http://projects.coin-or.org/Ipopt">http://projects.coin-or.org/Ipopt</a><br>******************************************************************************<br><br>NOTE: You are using Ipopt by default with the MUMPS linear solver.<br>
Other linear solvers might be more efficient (see Ipopt documentation).<br><br><br>This is Ipopt version trunk, running with linear solver mumps.<br><br>Number of nonzeros in equality constraint Jacobian...: 1402<br>
Number of nonzeros in inequality constraint Jacobian.: 1<br>Number of nonzeros in Lagrangian Hessian.............: 783<br><br>Total number of variables............................: 198<br> variables with only lower bounds: 0<br>
variables with lower and upper bounds: 198<br> variables with only upper bounds: 0<br>Total number of equality constraints.................: 160<br>Total number of inequality constraints...............: 1<br>
inequality constraints with only lower bounds: 0<br> inequality constraints with lower and upper bounds: 1<br> inequality constraints with only upper bounds: 0<br><br> 0 1.0699999e+02 7.18e-01 8.76e-01 0.0 0.00e+00 - 0.00e+00 0.00e+00 0<br>
1 1.1255327e+02 5.73e-01 3.64e+02 0.5 4.40e+00 - 9.25e-01 1.80e-01f 1<br> 2 1.1380862e+02 2.71e-01 5.62e+02 0.8 6.82e-01 2.0 1.00e+00 5.15e-01h 1<br> 3 1.3499130e+02 1.82e-01 4.76e+02 -5.3 6.28e+00 - 2.37e-01 4.82e-01h 1<br>
4 2.4094031e+02 9.45e-01 1.78e+03 0.9 1.51e+01 - 5.12e-01 1.00e+00f 1<br> 5 2.4101437e+02 4.12e-01 1.13e+03 0.7 4.06e-01 2.4 1.00e+00 9.60e-01h 1<br> 6 2.4174477e+02 2.84e-01 7.95e+02 0.5 4.45e-01 1.9 4.50e-01 3.78e-01f 2<br>
7 2.4215610e+02 2.36e-01 5.62e+03 1.3 4.33e-01 3.3 9.60e-01 6.36e-01h 1<br> 8 2.4821337e+02 2.14e-01 4.34e+03 1.0 6.56e+00 - 2.11e-01 1.32e-01h 1<br> 9 2.6251810e+02 1.57e-01 2.99e+03 -0.7 8.17e+00 - 3.75e-01 2.50e-01h 3<br>
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls<br> 10 2.6782805e+02 1.53e-01 2.73e+03 0.8 1.09e+01 - 2.52e-01 6.97e-02f 1<br> 11 3.1435221e+02 9.50e-01 1.69e+03 0.5 1.18e+01 - 2.95e-01 5.65e-01f 1<br>
12 3.5382332e+02 1.31e-01 1.66e+03 0.8 5.64e+00 - 9.09e-01 1.00e+00f 1<br> 13 3.5917148e+02 9.98e-03 9.11e+01 0.0 7.64e-01 - 9.39e-01 1.00e+00f 1<br> 14 3.3459479e+02 1.25e-01 1.98e+01 -0.1 3.51e+00 - 9.54e-01 1.00e+00f 1<br>
15 3.1112934e+02 8.73e-02 5.45e+00 -0.5 3.35e+00 - 1.00e+00 1.00e+00f 1<br> 16 3.1609087e+02 4.78e-02 5.92e+00 -1.1 7.09e-01 - 9.89e-01 1.00e+00h 1<br> 17 3.1949829e+02 1.49e-02 9.08e-01 -1.4 4.87e-01 - 9.96e-01 1.00e+00h 1<br>
18 3.2007105e+02 4.47e-03 1.96e+00 -2.6 1.00e-01 - 1.00e+00 8.16e-01h 1<br> 19 3.2038388e+02 1.64e-03 1.17e-01 -3.2 4.47e-02 - 1.00e+00 1.00e+00h 1<br>iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls<br>
20 3.2038331e+02 2.72e-04 2.33e+00 -4.7 9.19e-03 - 1.00e+00 8.41e-01h 1<br> 21 3.2038674e+02 3.00e-06 2.27e+00 -6.3 1.52e-03 - 1.00e+00 1.00e+00h 1<br> 22 3.2038464e+02 7.84e-07 2.21e+00 -7.9 8.58e-04 - 1.00e+00 1.00e+00h 1<br>
23 3.2038462e+02 2.39e-07 1.25e-01 -9.4 1.57e-05 - 1.00e+00 6.96e-01h 1<br> 24 3.2038463e+02 7.35e-11 1.63e-03 -10.1 7.36e-07 - 1.00e+00 1.00e+00h 1<br> 25 3.2038463e+02 1.64e-15 7.09e-06 -11.0 1.80e-09 - 1.00e+00 1.00e+00h 1<br>
26 3.2038463e+02 1.30e-15 5.71e-06 -11.0 1.40e-10 - 1.00e+00 1.00e+00h 1<br> 27 3.2038463e+02 1.57e-15 5.36e-06 -11.0 2.25e-10 - 1.00e+00 1.00e+00h 1<br> 28 3.2038463e+02 1.39e-15 5.36e-06 -11.0 2.85e-10 - 1.00e+00 1.00e+00H 1<br>
29 3.2038463e+02 9.42e-16 4.11e-06 -11.0 1.18e-10 - 1.00e+00 1.00e+00h 1<br>iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls<br> 30 3.2038463e+02 1.97e-15 4.46e-06 -11.0 1.24e-10 - 1.00e+00 1.00e+00h 1<br>
31 3.2038463e+02 2.03e-15 1.02e-05 -11.0 6.31e-11 - 1.00e+00 1.00e+00h 1<br> 32 3.2038463e+02 1.56e-15 8.29e-06 -11.0 1.90e-10 - 1.00e+00 5.00e-01h 2<br> 33 3.2038463e+02 1.18e-15 1.05e-05 -11.0 1.81e-10 - 1.00e+00 1.00e+00h 1<br>
34 3.2038463e+02 1.63e-15 1.96e-06 -11.0 1.69e-10 - 1.00e+00 3.12e-02h 6<br> 35 3.2038463e+02 1.31e-15 5.36e-06 -11.0 1.07e-10 - 1.00e+00 1.00e+00h 1<br> 36 3.2038463e+02 1.37e-15 6.41e-06 -11.0 4.15e-10 - 1.00e+00 1.00e+00H 1<br>
37 3.2038463e+02 1.20e-15 1.02e-05 -11.0 2.10e-10 - 1.00e+00 2.50e-01h 3<br> 38 3.2038463e+02 1.47e-15 5.58e-06 -11.0 3.11e-10 - 1.00e+00 1.00e+00H 1<br> 39 3.2038463e+02 1.60e-15 8.31e-06 -11.0 2.92e-10 - 1.00e+00 6.10e-05h 15<br>
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls<br> 40 3.2038463e+02 1.56e-15 4.49e-06 -11.0 2.12e-10 - 1.00e+00 3.81e-06h 19<br> 41 3.2038463e+02 1.60e-15 8.93e-06 -11.0 3.12e-10 - 1.00e+00 3.81e-06h 19<br>
42 3.2038463e+02 2.06e-15 3.21e-06 -11.0 2.11e-10 - 1.00e+00 1.00e+00H 1<br> 43 3.2038463e+02 1.62e-15 3.40e-06 -11.0 1.19e-10 - 1.00e+00 6.25e-02h 5<br> 44 3.2038463e+02 1.43e-15 2.52e-06 -11.0 2.46e-10 - 1.00e+00 3.12e-02h 6<br>
45 3.2038463e+02 1.90e-15 5.36e-06 -11.0 1.13e-10 - 1.00e+00 1.00e+00h 1<br> 46 3.2038463e+02 1.40e-15 6.06e-06 -11.0 2.41e-10 - 1.00e+00 1.00e+00H 1<br> 47 3.2038463e+02 1.93e-15 1.44e-06 -11.0 1.03e-10 - 1.00e+00 1.00e+00H 1<br>
48 3.2038463e+02 2.03e-15 3.80e-06 -11.0 1.95e-10 - 1.00e+00 1.00e+00h 1<br> 49 3.2038463e+02 1.50e-15 7.58e-06 -11.0 1.62e-10 - 1.00e+00 5.00e-01h 2<br>iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls<br>
50 3.2038463e+02 1.99e-15 1.05e-05 -11.0 2.60e-10 - 1.00e+00 2.50e-01h 3<br> 51 3.2038463e+02 2.45e-15 1.02e-05 -11.0 1.36e-10 - 1.00e+00 1.00e+00H 1<br> 52 3.2038463e+02 1.57e-15 4.89e-06 -11.0 2.00e-10 - 1.00e+00 5.00e-01h 2<br>
53 3.2038463e+02 1.65e-15 4.91e-06 -11.0 1.04e-10 - 1.00e+00 1.56e-02h 7<br> 54 3.2038463e+02 9.54e-16 5.15e-06 -11.0 1.36e-10 - 1.00e+00 1.00e+00H 1<br> 55 3.2038463e+02 1.68e-15 1.43e-05 -11.0 2.60e-10 - 1.00e+00 1.00e+00h 1<br>
56 3.2038463e+02 1.49e-15 7.86e-06 -11.0 1.28e-10 - 1.00e+00 5.00e-01h 2<br> 57 3.2038463e+02 8.75e-16 9.11e-06 -11.0 8.05e-11 - 1.00e+00 5.00e-01h 2<br> 58 3.2038463e+02 2.86e-15 1.46e-05 -11.0 1.08e-10 - 1.00e+00 1.00e+00H 1<br>
59 3.2038463e+02 1.79e-15 8.21e-06 -11.0 8.21e-11 - 1.00e+00 1.00e+00h 1<br>iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls<br> 60 3.2038463e+02 1.97e-15 1.32e-05 -11.0 1.49e-10 - 1.00e+00 1.25e-01h 4<br>
61 3.2038463e+02 1.09e-15 5.36e-06 -11.0 1.61e-10 - 1.00e+00 3.12e-02h 6<br> 62 3.2038463e+02 1.78e-15 1.01e-05 -11.0 8.98e-11 - 1.00e+00 1.00e+00H 1<br> 63 3.2038463e+02 2.03e-15 5.91e-06 -11.0 9.23e-11 - 1.00e+00 1.00e+00h 1<br>
64 3.2038463e+02 1.17e-15 1.04e-05 -11.0 9.45e-11 - 1.00e+00 5.00e-01h 2<br> 65 3.2038463e+02 2.45e-15 9.37e-06 -11.0 1.22e-10 - 1.00e+00 1.00e+00h 1<br> 66 3.2038463e+02 1.81e-15 5.17e-06 -11.0 1.33e-10 - 1.00e+00 1.00e+00h 1<br>
67 3.2038463e+02 1.32e-15 5.59e-06 -11.0 1.73e-10 - 1.00e+00 7.81e-03h 8<br> 68 3.2038463e+02 1.47e-15 8.57e-06 -11.0 1.79e-10 - 1.00e+00 5.00e-01h 2<br> 69 3.2038463e+02 1.30e-15 1.09e-05 -11.0 2.43e-10 - 1.00e+00 1.56e-02h 7<br>
iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls<br> 70 3.2038463e+02 1.78e-15 3.78e-06 -11.0 2.27e-10 - 1.00e+00 1.56e-02h 7<br> 71 3.2038463e+02 1.81e-15 1.41e-05 -11.0 3.91e-10 - 1.00e+00 1.00e+00h 1<br>
72 3.2038463e+02 1.58e-15 1.47e-05 -11.0 2.98e-10 - 1.00e+00 6.25e-02h 5<br> 73 3.2038463e+02 1.44e-15 1.37e-05 -11.0 1.90e-10 - 1.00e+00 1.56e-02h 7<br> 74 3.2038463e+02 2.20e-15 3.20e-06 -11.0 2.94e-10 - 1.00e+00 3.12e-02h 6<br>
75 3.2038463e+02 2.03e-15 8.62e-06 -11.0 1.40e-10 - 1.00e+00 1.00e+00h 1<br> 76 3.2038463e+02 2.57e-15 2.12e-06 -11.0 3.28e-10 - 1.00e+00 1.00e+00h 1<br> 77 3.2038463e+02 1.54e-15 1.47e-07 -11.0 2.14e-10 - 1.00e+00 6.25e-02h 5<br>
78 3.2038463e+02 1.79e-15 4.11e-06 -11.0 1.06e-10 - 1.00e+00 1.25e-01h 4<br> 79 3.2038463e+02 1.79e-15 9.27e-06 -11.0 1.84e-10 - 1.00e+00 3.81e-06h 19<br>iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls<br>
80 3.2038463e+02 1.17e-15 4.49e-06 -11.0 1.04e-10 - 1.00e+00 1.25e-01h 4<br> 81 3.2038463e+02 1.09e-15 8.56e-06 -11.0 1.75e-10 - 1.00e+00 9.77e-04h 11<br> 82 3.2038463e+02 1.28e-15 4.14e-06 -11.0 2.37e-10 - 1.00e+00 5.00e-01h 2<br>
83 3.2038463e+02 1.88e-15 9.35e-06 -11.0 1.58e-10 - 1.00e+00 5.00e-01h 2<br> 84 3.2038463e+02 1.47e-15 1.05e-05 -11.0 1.97e-10 - 1.00e+00 1.25e-01h 4<br> 85 3.2038463e+02 1.58e-15 9.26e-06 -11.0 1.55e-10 - 1.00e+00 1.00e+00H 1<br>
86 3.2038463e+02 1.41e-15 4.81e-06 -11.0 1.42e-10 - 1.00e+00 1.00e+00h 1<br> 87 3.2038463e+02 1.49e-15 3.55e-06 -11.0 2.13e-10 - 1.00e+00 1.25e-01h 4<br> 88 3.2038463e+02 1.50e-15 5.70e-06 -11.0 9.55e-11 - 1.00e+00 3.12e-02h 6<br>
89 3.2038463e+02 1.41e-15 4.51e-06 -11.0 2.62e-10 - 1.00e+00 1.25e-01h 4<br>iter objective inf_pr inf_du lg(mu) ||d|| lg(rg) alpha_du alpha_pr ls<br> 90 3.2038463e+02 1.15e-15 1.24e-05 -11.0 2.20e-10 - 1.00e+00 6.25e-02h 5<br>
91 3.2038463e+02 9.47e-16 8.30e-06 -11.0 2.47e-10 - 1.00e+00 5.00e-01h 2<br> 92 3.2038463e+02 2.49e-15 6.99e-07 -11.0 1.55e-10 - 1.00e+00 1.00e+00h 1<br> 93 3.2038463e+02 1.54e-15 1.39e-06 -11.0 3.53e-10 - 1.00e+00 2.50e-01h 3<br>
94 3.2038463e+02 1.43e-15 6.26e-06 -11.0 1.39e-10 - 1.00e+00 1.00e+00H 1<br> 95 3.2038463e+02 1.22e-15 8.28e-06 -11.0 5.58e-11 - 1.00e+00 1.00e+00h 1<br> 96 3.2038463e+02 1.81e-15 8.28e-06 -11.0 2.78e-10 - 1.00e+00 1.00e+00h 1<br>
Restoration phase is called at almost feasible point,<br> but acceptable point from iteration 92 could be restored.<br clear="all"><br>