[Help-glpk] Linear IP solve instability

[Laurence Davies]

Hi, I have a very simple sparse LP which is a L1 norm minimisation. It solves 
fine using the simplex solver but terminates early in the IP solver with a 
numerical instability error.

I've attached the file that causes this.

Unnecessary history: This isn't the first solver I've tried - I've had success 
with CVX-Matlab and Mathematica solvers using various formulations. I tried to 
write my own solver using Dikin's method without success (seems to never 
converge) and I've written my own simplex solver which works fine.

The KKT conditions report high quality except for KKT.PB which reports medium 
quality for the simplex solver and primal infeasible for the IPM.

I've placed nonneg bounds on the residuals u+v and X is bounded by +/- 1e6 for 
the sake of sanity. The latter seemed to be useful in Mathematica. Could anyone 
please make a suggestion of what to try next?

Regards,

Laurence

test_glpk_prob.txt
Description: test_glpk_prob.txt