<br><font size=2 face="sans-serif">Hi Kranthi,</font>
<br>
<br><font size=2 face="sans-serif">If the program segfaults at that point,
you probably gave invalid information about the structure of the Jacobian
(AVAR or NNZA) to Ipopt. &nbsp;Please check that you populate AVAR and
ACON and NNZA correctly.</font>
<br>
<br><font size=2 face="sans-serif">As for LIWMA and LA, those are related
to the work space that MA27 requires for performing and storing the factorization
of the iteration matrix. &nbsp;Since the fill-in of the factorization is
not known before it is attempted, we take a guess of how much it could
be, reserve that memory, try the factorization, and if that turns out not
to be enough, MA27 complaints and we refactorize the matrix after increasing
the size of the work space - this is when that message is printed. &nbsp;If
you really want to know more, please consult the documentation of MA27
(which is available at Harwell, where you also downloaded MA27). &nbsp;Usually,
seeing this message means that the linear system is ill-conditioned, but
in your case that might be because some information about the Jacobian
structure is incorrect (see the seg fault).</font>
<br>
<br><font size=2 face="sans-serif">I hope this helps,</font>
<br>
<br><font size=2 face="sans-serif">Andreas</font>
<br>
<br>
<br>
<br>
<table width=100%>
<tr valign=top>
<td>
<td><font size=1 face="sans-serif"><b>Kranthi K Gade &lt;kranthi@cs.nyu.edu&gt;</b></font>
<br><font size=1 face="sans-serif">Sent by: coin-ipopt-bounces@list.coin-or.org</font>
<p><font size=1 face="sans-serif">12/04/2004 11:34 PM</font>
<td><font size=1 face="Arial">&nbsp; &nbsp; &nbsp; &nbsp; </font>
<br><font size=1 face="sans-serif">&nbsp; &nbsp; &nbsp; &nbsp; To:
&nbsp; &nbsp; &nbsp; &nbsp;coin-ipopt@list.coin-or.org</font>
<br><font size=1 face="sans-serif">&nbsp; &nbsp; &nbsp; &nbsp; cc:
&nbsp; &nbsp; &nbsp; &nbsp;</font>
<br><font size=1 face="sans-serif">&nbsp; &nbsp; &nbsp; &nbsp; Subject:
&nbsp; &nbsp; &nbsp; &nbsp;[Coin-ipopt] LIWMA/LA problem?</font></table>
<br>
<br>
<br><font size=2><tt>Hello all,<br>
<br>
I have been trying to use IPOPT to optimize a problem with around 2000
<br>
variables and 2 constraints. &nbsp;The program gives seg fault on iteration
1.<br>
<br>
Here's the output:<br>
<br>
******************************************************************************<br>
This program contains IPOPT, a program for large-scale nonlinear <br>
optimization.<br>
 &nbsp; &nbsp;IPOPT is released as open source under the Common Public
License (CPL).<br>
 &nbsp; &nbsp;For more information visit www.coin-or.org/Ipopt<br>
******************************************************************************<br>
<br>
Number of variables &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; : &nbsp; &nbsp;
2105<br>
 &nbsp; &nbsp;of which are fixed &nbsp; &nbsp; &nbsp; &nbsp; : &nbsp; &nbsp;
&nbsp; &nbsp;0<br>
Number of constraints &nbsp; &nbsp; &nbsp; &nbsp; : &nbsp; &nbsp; &nbsp;
&nbsp;2<br>
Number of lower bounds &nbsp; &nbsp; &nbsp; &nbsp;: &nbsp; &nbsp; &nbsp;
&nbsp;2<br>
Number of upper bounds &nbsp; &nbsp; &nbsp; &nbsp;: &nbsp; &nbsp; &nbsp;
&nbsp;0<br>
Number of nonzeros in Jacobian: &nbsp; &nbsp; 4210<br>
Number of nonzeros in Hessian : &nbsp; &nbsp; 2105<br>
<br>
ITER &nbsp; &nbsp; ERR &nbsp; &nbsp; &nbsp; MU &nbsp; &nbsp; &nbsp;||C||
&nbsp; ||YPY|| &nbsp;||PZ|| &nbsp; &nbsp; ||D|| &nbsp; ALFA(V) <br>
ALFA(X) &nbsp; &nbsp; NU &nbsp; #LS &nbsp; &nbsp; &nbsp; F &nbsp; &nbsp;
&nbsp;#cor &nbsp; Regu &nbsp; &nbsp;CPU(s)<br>
0 .987E+05p .100E+00 .140E+06 .000E+00 .000E+00 &nbsp;.000E+00 .000E+00
<br>
.000E+00 &nbsp;.000E+00 &nbsp;0 -.57231356E+03 &nbsp;0 .000E+00 .313E-01<br>
<br>
ma27_call: &nbsp;LIWMA increased from 73675 to 1307435<br>
ma27_call: &nbsp;LA increased from 63145 to 1307435<br>
<br>
1 .134E+06p .100E+00 .138E+06 .000E+00 .000E+00 &nbsp;.621E+04 .976E+00d
<br>
.100E+01h .000E+00 &nbsp;1 -.82994516E+03 &nbsp;2 .100E-03 .938E-01<br>
<br>
Can somebody tell me what are these LIWMA and LA? &nbsp;And is there a
<br>
problem if they are increased? &nbsp;Also is there a good source (FAQ list
or <br>
something) to which I can refer to when faced with such problems?<br>
<br>
I am attaching some debugging information just in case you might find it
<br>
useful.<br>
<br>
Thanks a lot.<br>
kranthi<br>
Program received signal SIGSEGV, Segmentation fault.<br>
0x0043b113 in constr_order_a__ (nzorig=0x4fe8e0, norig=0x3d8340, m=0x3d8344,
nind=0x4fe0f0, ivar=0x75edd4, nfix=0x4fe120, nz=0x4fe9b0, a=0x650020, acon=0x740048,
avar=0x744210,<br>
 &nbsp; &nbsp;nzc=0x4fe9a0, liw=0x22f240, iw=0x760eb8, ierr=0x22fee4) at
constr.F:2310<br>
2310 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;k = IW(p_ivar1+AVAR(ln2))<br>
Current language: &nbsp;auto; currently fortran<br>
(gdb) bt<br>
#0 &nbsp;0x0043b113 in constr_order_a__ (nzorig=0x4fe8e0, norig=0x3d8340,
m=0x3d8344, nind=0x4fe0f0, ivar=0x75edd4, nfix=0x4fe120, nz=0x4fe9b0, a=0x650020,
acon=0x740048, avar=0x744210,<br>
 &nbsp; &nbsp;nzc=0x4fe9a0, liw=0x22f240, iw=0x760eb8, ierr=0x22fee4) at
constr.F:2310<br>
#1 &nbsp;0x0043a9b1 in constr_new_a__ (task=0x4f4a70, nzorig=0x4fe8e0,
norig=0x3d8340, xorig=0x3d85f8, m=0x3d8344, nind=0x4fe0f0, ivar=0x75edd4,
nfix=0x4fe120, ifix=0x740048, cscale=0x700f00,<br>
 &nbsp; &nbsp;nz=0x4fe9b0, a=0x650020, acon=0x740048, avar=0x744210, nzc=0x4fe9a0,
liw=0x22f428, iw=0x760eb8, ierr=0x22fee4, ev_f__=0x402d86 &lt;EV_F&gt;,
ev_c__=0x402e30 &lt;EV_C&gt;, ev_g__=0x402ddb &lt;EV_G&gt;,<br>
 &nbsp; &nbsp;ev_a__=0x402e8e &lt;EV_A&gt;, ev_h__=0x402f01 &lt;EV_H&gt;,
ev_hlv__=0x402fa5 &lt;EV_HLV&gt;, ev_hov__=0x402ff0 &lt;EV_HOV&gt;, ev_hcv__=0x403034
&lt;EV_HCV&gt;, dat=0x0, idat=0x3d8340) at constr.F:1818<br>
#2 &nbsp;0x0043a443 in constr_ (task=0x4f06f8, iter=0x4f06f4, n=0x4fe0e0,
nind=0x4fe0f0, m=0x3d8344, ivar=0x75edd4, nfix=0x4fe120, ifix=0x740048,
norig=0x3d8340, xorig=0x3d85f8,<br>
 &nbsp; &nbsp;cscale=0x700f00, vin=0x22f840, vout=0x700f10, ivec1=0x22f8a4,
ivec2=0x22f8a4, lrs=0x4fe020, rs=0x650020, lis=0x4fe028, is=0x740048, lrw=0x22f6f8,
rw=0x700f20, liw=0x22f6fc,<br>
 &nbsp; &nbsp;iw=0x760eb8, ierr=0x22fee4, ev_f__=0x402d86 &lt;EV_F&gt;,
ev_c__=0x402e30 &lt;EV_C&gt;, ev_g__=0x402ddb &lt;EV_G&gt;, ev_a__=0x402e8e
&lt;EV_A&gt;, ev_h__=0x402f01 &lt;EV_H&gt;, ev_hlv__=0x402fa5 &lt;EV_HLV&gt;,<br>
 &nbsp; &nbsp;ev_hov__=0x402ff0 &lt;EV_HOV&gt;, ev_hcv__=0x403034 &lt;EV_HCV&gt;,
dat=0x0, idat=0x3d8340) at constr.F:1502<br>
#3 &nbsp;0x00407cba in get_scale__ (iter=0x22ff64, norig=0x3d8340, xorig=0x3d85f8,
n=0x4fe0e0, x=0x6fcd28, nind=0x4fe0f0, m=0x3d8344, ivar=0x75edd4, nfix=0x4fe120,
ifix=0x740048,<br>
 &nbsp; &nbsp;nlbo=0x3d8348, ilbo=0x3d83a0, bnds_lo__=0x3d83c0, nubo=0x3d834c,
iubo=0xbaadf00d, bnds_uo__=0xbaadf00d, cscale=0x700f00, kconstr=0x4fe020,
lrs=0x4fe040, rs=0x650020, lis=0x4fe050,<br>
 &nbsp; &nbsp;is=0x740048, lrw=0x22fa0c, rw=0x700f10, liw=0x22fa14, iw=0x760eb8,
ierr=0x22fee4, ev_f__=0x402d86 &lt;EV_F&gt;, ev_c__=0x402e30 &lt;EV_C&gt;,
ev_g__=0x402ddb &lt;EV_G&gt;, ev_a__=0x402e8e &lt;EV_A&gt;,<br>
 &nbsp; &nbsp;ev_h__=0x402f01 &lt;EV_H&gt;, ev_hlv__=0x402fa5 &lt;EV_HLV&gt;,
ev_hov__=0x402ff0 &lt;EV_HOV&gt;, ev_hcv__=0x403034 &lt;EV_HCV&gt;, dat=0x0,
idat=0x3d8340) at get_scale.F:362<br>
#4 &nbsp;0x004045d7 in ipoptmain_ (norig=0x3d8340, xorig=0x3d85f8, m=0x3d8344,
nlbo=0x3d8348, ilbo=0x3d83a0, bnds_lo__=0x3d83c0, nubo=0x3d834c, iubo=0xbaadf00d,
bnds_uo__=0xbaadf00d,<br>
 &nbsp; &nbsp;v_lo__=0x3dc7d8, v_uo__=0x3dc800, lam=0x3dc820, c=0x3dc848,
lrw=0x22fce8, rw=0x650020, liw=0x22fce4, iw=0x740048, iter=0x22ff64, ierr=0x22fee4,
ev_f__=0x402d86 &lt;EV_F&gt;,<br>
 &nbsp; &nbsp;ev_c__=0x402e30 &lt;EV_C&gt;, ev_g__=0x402ddb &lt;EV_G&gt;,
ev_a__=0x402e8e &lt;EV_A&gt;, ev_h__=0x402f01 &lt;EV_H&gt;, ev_hlv__=0x402fa5
&lt;EV_HLV&gt;, ev_hov__=0x402ff0 &lt;EV_HOV&gt;, ev_hcv__=0x403034 &lt;EV_HCV&gt;,<br>
 &nbsp; &nbsp;dat=0x0, idat=0x3d8340) at ipopt.F:632<br>
#5 &nbsp;0x00403965 in ipopt_ (n=0x3d8340, x=0x3d85f8, m=0x3d8344, nlb=0x3d8348,
ilb=0x3d83a0, bnds_l__=0x3d83c0, nub=0x3d834c, iub=0xbaadf00d, bnds_u__=0xbaadf00d,
v_l__=0x3dc7d8,<br>
 &nbsp; &nbsp;v_u__=0x3dc800, lam=0x3dc820, c=0x3dc848, lrw=0x22feec, rw=0x0,
liw=0x22fef0, iw=0x0, iter=0x22ff64, ierr=0x22fee4, ev_f__=0x402d86 &lt;EV_F&gt;,
ev_c__=0x402e30 &lt;EV_C&gt;,<br>
 &nbsp; &nbsp;ev_g__=0x402ddb &lt;EV_G&gt;, ev_a__=0x402e8e &lt;EV_A&gt;,
ev_h__=0x402f01 &lt;EV_H&gt;, ev_hlv__=0x402fa5 &lt;EV_HLV&gt;, ev_hov__=0x402ff0
&lt;EV_HOV&gt;, ev_hcv__=0x403034 &lt;EV_HCV&gt;, dat=0x0, idat=0x3d8340,<br>
 &nbsp; &nbsp;nargs=0x3d837c, args=0x3d83e8, cargs=0x3d8450, __g77_length_cargs=40)
at ipopt.F:285<br>
Exception: STATUS_ACCESS_VIOLATION at eip=004321E9<br>
eax=00000000 ebx=10288B60 ecx=10166C50 edx=00000000 esi=101209D0 edi=00000000<br>
ebp=0022E498 esp=0022E430 program=C:\cygwin\bin\gdb.exe, pid 3724, thread
main<br>
cs=001B ds=0023 es=0023 fs=003B gs=0000 ss=0023<br>
Stack trace:<br>
Frame &nbsp; &nbsp; Function &nbsp;Args<br>
0022E498 &nbsp;004321E9 &nbsp;(10288B00, 101276E0, FFFFFFFF, 10131160)<br>
0022E4B8 &nbsp;004322DE &nbsp;(0022E670, 00000006, 0022E5F8, 00422B64)<br>
0022E4C8 &nbsp;00422CA1 &nbsp;(10129DC8, 0022E630, 0022E5D4, 0042EFF9)<br>
0022E5F8 &nbsp;00422B64 &nbsp;(00422C90, 10129DC8, 0022E630, 0022E628)<br>
0022E638 &nbsp;00422CFA &nbsp;(004322A0, 0022E670, 00431D7E, 00000006)<br>
0022E698 &nbsp;00432AEB &nbsp;(101276E0, 00000006, 00000000, 00000001)<br>
0022E768 &nbsp;0043241A &nbsp;(101276E0, 00000006, 00000000, 00000001)<br>
0022E7A8 &nbsp;004338FB &nbsp;(00000000, 00000000, 00000001, 101A13A0)<br>
0022E7E8 &nbsp;00433B7E &nbsp;(00000000, 00000001, 0022E818, 00407B6C)<br>
0022E7F8 &nbsp;00405879 &nbsp;(10114EC8, 00000000, 00000001, 00000001)<br>
0022E818 &nbsp;00407B6C &nbsp;(10114EC8, 00000000, 00000001, 00000000)<br>
0022E858 &nbsp;004230FB &nbsp;(101002AA, 00000001, 0022E888, 005B8594)<br>
0022E888 &nbsp;00484B4B &nbsp;(101002A8, 10191938, 0022E8B0, 005C852C)<br>
0022E8C8 &nbsp;00484FC5 &nbsp;(1019F220, 1019F220, 0022F0A4, 1019F220)<br>
0022E8E8 &nbsp;005C6ED5 &nbsp;(0022E8A8, 101235E8, 0022E908, 004849AD)<br>
0022E8F8 &nbsp;0048438B &nbsp;(00000000, 0000000C, 0022E928, 0049D35E)<br>
End of stack trace (more stack frames may be present)<br>
Segmentation fault (core dumped)_______________________________________________<br>
Coin-ipopt mailing list<br>
Coin-ipopt@list.coin-or.org<br>
http://list.coin-or.org/mailman/listinfo/coin-ipopt<br>
</tt></font>
<br>