Re: [Bug-glpk] Assertion at glpipm.c; line 571

[Barry Rountree]

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