[Bug-glpk] Assertion failure in glp_simplex GLPK-4.18

[Cameron Kellough]

Hi,

I am using GLPK-4.18 with the new version of octave 2.9.12.  I've pasted 
in the octave program in case it is needed as a reference.
The bug however is that in glplpx01.c the assertion on line 345 fails 
because there is no case in the switch statement for error code 0x4.

Here is the error code clipped from glpk.h:   
#define GLP_EBOUND   0x04  /* invalid bounds */

This error is returned by glp_simplex in routine lpx_simplex.  
Additionally, it lokks like many cases are not handled in this case 
statement.
I'd send a patch but I don't know the error code mappings well enough to 
map between the error code sets.

Thanks for a handy piece of software and hope this helps.

Octave program follows: 


a=[ 100.0, 0.0,   0.0,   0.0,   0.0,   0.0,   0.0,    0.0,   0.0,    
0.0,     0.0,     0.0,      0.0,    0.0,     0.0,     0.0,     0.0,     
0.0,     0.0,     0.0,  0.0
    0.0,  11.0,  25.4,  26.0,  21.3,  11.9,   4.3,   0.1,    0.0,    
0.0,     0.0,     0.0,     0.0,    0.0,     0.0,     0.0,     0.0,     
0.0,     0.0,     0.0,  0.0;
    0.0,   0.0,   0.0,   0.0,   0.0,  25.3,  31.7,   25.6,   4.2,    
2.6,     0.0,     0.0,     0.0,    0.0,     0.0,     0.0,     0.0,     
0.0,     0.0,     0.0,   0.0;
    0.0,   0.0,   0.0,   0.0,   0.0,   0.0,   0.0,    3.6,  25.8,   
41.6,    10.3,     4.3,     1.0,    0.0,     0.0,     0.0,     0.0,     
0.0,     0.0,     0.0,   0.0;
    0.0,   0.0,   0.0,   0.0,   0.0,   0.0,   0.0,    0.0,   0.0,    
1.6,    35.5,    38.7,    18.4,    4.4,     1.0,     0.4,     0.0,     
0.0,     0.0,     0.0,   0.0;
    0.0,   0.0,   0.0,   0.0,   0.0,   0.0,   0.0,    0.0,   0.0,    
0.0,     1.5,    13.5,    22.2,   18.1,    18.8,    15.4,     7.5,     
2.6,     0.4,     0.0,   0.0;
    0.0,   0.0,   0.0,   0.0,   0.0,   0.0,   0.0,    0.0,   0.0,    
0.0,     0.0,     0.0,     0.0,    1.9,    21.0,    36.5,    21.1,     
9.8,     5.5,     2.8,   0.0;]'/100;


c= [1,1,1,1,1,1,1];

b= [.2/3.0, .2/3.0 , .2/3.0 , .2/6.0, .2/6.0, .2/6.0, .2/6.0, .2/6.0, 
.2/6.0, .2/4.0, .2/4.0, .2/4.0, .2/4.0, .2/5.0, .2/5.0, .2/5.0, .2/5.0, 
.2/5.0, .2/3.0, .2/3.0, .2/3.0 ]';

ctype =   "LLLLLLLLLLLLLLLLLLLLF";
vartype = "CCCCCCC";
s = 1;
lb = [0,0,0,0,0,0,0]';
ub = [1,1,1,1,1,1,0]';
        
param.msglev = 3;
param.itlim = 100;

[xopt, fopt, status, extra] = ...
             glpk (c, a, b, lb, ub, ctype, vartype, s, param);

xopt
fopt
status
extra

a*xopt*(5.0/7.0)*.75
sum(a*xopt*(5.0/7.0)*.75)

-- 
Cameron Kellough, Research Engineer
ESD Intelligence & Information Systems
SRI International, Huntsville Field Operations
Tel 256-684-0417 Page 1-877-730-9920