Hello
Thank you for what you are doing.
I have a recurrent problem when solving mixed-integer LPs.
Attached are an LP file example and its corresponding Output file.
Here is the glpsol's terminal output:
GLPSOL: GLPK LP/MIP Solver, v4.45
Parameter(s) specified in the command line:
--lp example.lp --output example.out
Reading problem data from `example.lp'...
example.lp:1282: warning: missing final end of line
952 rows, 634 columns, 2536 non-zeros
317 integer variables, all of which are binary
1282 lines were read
GLPK Integer Optimizer, v4.45
952 rows, 634 columns, 2536 non-zeros
317 integer variables, all of which are binary
Preprocessing...
952 rows, 634 columns, 2536 non-zeros
317 integer variables, all of which are binary
Scaling...
A: min|aij| = 1.000e+00 max|aij| = 4.134e+04 ratio = 4.134e+04
GM: min|aij| = 1.655e-01 max|aij| = 6.043e+00 ratio = 3.652e+01
EQ: min|aij| = 2.738e-02 max|aij| = 1.000e+00 ratio = 3.652e+01
2N: min|aij| = 3.027e-02 max|aij| = 1.262e+00 ratio = 4.167e+01
Constructing initial basis...
Size of triangular part = 952
Solving LP relaxation...
GLPK Simplex Optimizer, v4.45
952 rows, 634 columns, 2536 non-zeros
0: obj = 6.675300000e+03 infeas = 2.794e+02 (0)
500: obj = 9.252882224e+03 infeas = 1.295e+02 (0)
Error: unable to factorize the basis matrix (2)
Sorry, basis recovery procedure not implemented yet
glp_intopt: cannot solve LP relaxation
Time used: 0.0 secs
Memory used: 1.0 Mb (1026045 bytes)
Writing MIP solution to `example.out'...
example.lp