<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Apr 23, 2015 at 2:57 PM, Jason Sarich <span dir="ltr"><<a href="mailto:jason.sarich@gmail.com" target="_blank">jason.sarich@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div dir="ltr">Hi David,<div><br></div><div>The IPM code is a naive implementation and not very battle-tested. There is nothing specifically broken, but the linear solve operator becomes very ill-conditioned and for most cases will require a direct solver (e.g. superlu). This is the major problem, and if you're ok with that then it may work for you.<br></div><div>Unfortunately, we don't have anything else to recommend for generally constrained problems.</div></div></blockquote><div><br></div><div><br></div><div>A quick follow-up question: I notice that the "toy" example I referred to earlier in this thread uses VecCreateSeq and MatCreateSeq, so I just wanted to check if TAO IPM is restricted to sequential data? (My guess is that it's not restricted to sequential data, and that Seq is used only because the problem is so small, but I just wanted to check.)</div><div><br></div><div>Thanks,</div><div>David</div><div><br></div></div></div></div>