On Mon March 27 2006 04:28, Andrew Makhorin wrote: > > Assertion failed: x[j] > 0.0 && z[j] > 0.0; file glpipm.c; line 571
> Thank you for the bug report. > > Looks like the bug appears due to excessive round-off errors in the > solution of the normal system, because the problem is badly scaled. > > Could you please gzip the problem data causing the bug in mps format > and send it to me? Eh, that was the hint I needed. I hadn't used mps format before, so I decided to have a look at it before I sent it off to you. It had lots of "Inf"s, which tends to happen when I divide by zero (funny, that). I modified my data so that I'm dividing by a sufficiently small non-zero number and the interior-point method found the optimal solution as expected. I'm appending the file in case you want to check for anything else that might be going on. Thanks for the help! Barry > > Andrew Makhorin _______________________________________________ Bug-glpk mailing list [email protected] http://lists.gnu.org/mailman/listinfo/bug-glpk
