<div dir="ltr"><div><div dir="auto"> The crash doesn't occur with the newer HSL solvers. Here’s what I’ve found so far - I would be glad to hear if I’ve gotten anything wrong, and suggest that a summary such as the one below might be useful if included with the documentation about the solvers. It is hard to do large-scale optimization with the 32-bit limit.<br></div></div><div dir="auto"><br></div><div dir="auto">Internal support for 64-bit integers in the 32-bit interface:</div><div dir="auto">MA27: no</div><div dir="auto">MA57: no, but at some point Mathworks used to use a modified version that did support 64-bit integers before switching to UMFPACK. I don't know whether this version is available somewhere.<br></div><div dir="auto"><div dir="auto">MUMPS: no as of 4.10.0 (the version downloaded by Ipopt), yes as of 5.1.0. I haven't tried to use 5.1.0+ with Ipopt yet.<br></div><div dir="auto"></div>WSMP: no (it failed in the same place as MA27, MA57, and MUMPS 4.10.0)<div dir="auto"></div>MA86: yes (it worked where MA27, MA57, MUMPS 4.10.0, and WSMP failed)</div><div dir="auto">MA97: yes (presumably the same as MA86, though I haven't tried it yet)<br></div>PARDISO: ? (I've had difficulty linking against MKL on my system, so I haven't tried it)<br><div><br></div><div>    Thanks,</div><div>    Sean Hardesty<br></div><div><br><div class="gmail_quote"><div>On Wed, Jun 13, 2018 at 3:22 PM Sean Hardesty <<a href="mailto:sean.hardesty@gmail.com" target="_blank">sean.hardesty@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div><div>I am experiencing a crash during Ipopt's call to the linear solver, using packages provided with Ubuntu 16.04 (there are packages for Ipopt and MUMPS). The same thing happens with source builds of MUMPS, MA27, and MA57.<br><br>For my problem, the total number of variables / constraints and non-zero entries are all on the order of 10^5 to 10^6, so 32-bit integers are perfectly adequate for specifying the constraint Jacobian in triplet format. What I believe is happening is that factoring the system creates so much fill-in that the solver will fail unless it internally uses 64-bit integers for indexing.<br><br>Are any of the solvers supported by Ipopt free of this issue? I am aware from previous comments by Jonathan Hogg that the more modern HSL solvers use 64-bit integers, but are there any problems with the Ipopt interface? MUMPS gives vague instructions for compiling with 64-bit integers:<br><br><a href="http://mumps.enseeiht.fr/index.php?page=faq#3.3" target="_blank">http://mumps.enseeiht.fr/index<wbr>.php?page=faq#3.3</a><br><br>For HSL or MUMPS, do I have to hack metis 4 to use 64-bit integers, like in this post?<br><br><a href="http://matrixprogramming.rudnyi.ru/2008/05/metis-64bit/" target="_blank">http://matrixprogramming.rudny<wbr>i.ru/2008/05/metis-64bit/</a><br><br></div>    Thanks,<br></div>    Sean Hardesty<br></div>
</blockquote></div></div>
</div>