Re: Equality constraint not fulfilled

[Heinrich Schuchardt]

On 3/1/20 7:07 AM, palik imre wrote:
> Hi all,
>
> The following behaviour looks like a bug to me.  But I am new to GLPK, so sorry 
> if I am missing something obvious.
> I am running GLPK on 4.65-1.
> My problem matrix and constraints are (lower & upper bounds after the colon):
>
> 0.25028125163259074     0.0319722332178770316   0.582456482564010392    
> 0.00913492377653629672  0.166416137875431552    0       0       0       0       
> 0       0       0       0       0       0       0       0       : 
> 24.9998000101579372,24.9998000101579372
> -0.00682567831681596404 2.2945894737689132      -0.125954731367528411   
> 0.655596992505403819    -0.0359870661050081253  0       0       0       0       
> 0       0       0       0       0       0       0       0       : 
> -1.79769313486231571e+308,1.58396121568478776
> 0.00682567831681596404  -2.2945894737689132     0.125954731367528411    
> -0.655596992505403819   0.0359870661050081253   0       0       0       0       
> 0       0       0       0       0       0       0       0       : 
> -1.79769313486231571e+308,2.51116750466182204
> 0       0.0103065434721448315   0.166347686907376502    0.0412261738885793261   
> 0.665390747629506007    0.0103065434721448315   0.166347686907376502    0       
> 0       0       0       0       0       0       0       0       0       : 
> 49.9984822684764154,49.9984822684764154
> 0       0.654269049545594727    -0.0405370974913239479  2.61707619818237891     
> -0.162148389965295792   0.654269049545594727    -0.0405370974913239479  0       
> 0       0       0       0       0       0       0       0       0       : 
> -1.79769313486231571e+308,0.5129008702960558
> 0       -0.654269049545594727   0.0405370974913239479   -2.61707619818237891    
> 0.162148389965295792    -0.654269049545594727   0.0405370974913239479   0       
> 0       0       0       0       0       0       0       0       0       : 
> -1.79769313486231571e+308,3.55059259243481229
> 0       0       0       0.011522375091690247    0.166267894225025192    
> 0.0460895003667609882   0.665071576900100769    0.011522375091690247    
> 0.166267894225025192    0       0       0       0       0       0       0       
> 0       : 74.9943203809251315,74.9943203809251315
> 0       0       0       0.647572320594917983    -0.0448768008500337601  
> 2.59028928237967193     -0.17950720340013504    0.647572320594917983    
> -0.0448768008500337601  0       0       0       0       0       0       0       
> 0       : -1.79769313486231571e+308,-1.32410359077575279
> 0       0       0       -0.647572320594917983   0.0448768008500337601   
> -2.59028928237967193    0.17950720340013504     -0.647572320594917983   
> 0.0448768008500337601   0       0       0       0       0       0       0       
> 0       : -1.79769313486231571e+308,5.1984190796412193
> 0       0       0       0       0       0.0126959681730639627   
> 0.166182400301374655    0.0507838726922558506   0.664729601205498621    
> 0.0126959681730639627   0.166182400301374655    0       0       0       0       
> 0       0       : 99.9861467700102793,99.9861467700102793
> 0       0       0       0       0       0.653681253141257912    
> -0.0499398033134639263  2.61472501256503165     -0.199759213253855705   
> 0.653681253141257912    -0.0499398033134639263  0       0       0       0       
> 0       0       : -1.79769313486231571e+308,-3.47173985969024734
> 0       0       0       0       0       -0.653681253141257912   
> 0.0499398033134639263   -2.61472501256503165    0.199759213253855705    
> -0.653681253141257912   0.0499398033134639263   0       0       0       0       
> 0       0       : -1.79769313486231571e+308,7.5370952474605204
> 0       0       0       0       0       0       0       0.0135250388967059626   
> 0.166116980169458778    0.0541001555868238504   0.664467920677835111    
> 0.0135250388967059626   0.166116980169458778    0       0       0       0       
> : 124.977451930488201,124.977451930488201
> 0       0       0       0       0       0       0       0.658242656713506791    
> -0.0535933023008207457  2.63297062685402716     -0.214373209203282983   
> 0.658242656713506791    -0.0535933023008207457  0       0       0       0       
> : -1.79769313486231571e+308,-5.47507908163918344
> 0       0       0       0       0       0       0       -0.658242656713506791   
> 0.0535933023008207457   -2.63297062685402716    0.214373209203282983    
> -0.658242656713506791   0.0535933023008207457   0       0       0       0       
> : -1.79769313486231571e+308,9.68536004575470244
> 0       0       0       0       0       0       0       0       0       
> 0.013978388658102664    0.166079446134375086    0.0559135546324106558   
> 0.664317784537500344    0.013978388658102664    0.166079446134375086    0       
> 0       : 149.971776746037904,149.971776746037904
> 0       0       0       0       0       0       0       0       0       
> 0.658990315976552665    -0.0554651582303146498  2.63596126390621066     
> -0.221860632921258599   0.658990315976552665    -0.0554651582303146498  0       
> 0       : -1.79769313486231571e+308,-6.87337560268436576
> 0       0       0       0       0       0       0       0       0       
> -0.658990315976552665   0.0554651582303146498   -2.63596126390621066    
> 0.221860632921258599    -0.658990315976552665   0.0554651582303146498   0       
> 0       : -1.79769313486231571e+308,11.1108076471177775
> 0       0       0       0       0       0       0       0       0       0       
> 0       0.0140578740265756147   0.166072736942674326    0.0492025590930146445   
> 0.581254579299360086    0.021086811039863422    0.249109105414011489    : 
> 174.971063955562613,174.971063955562613
> 0       0       0       0       0       0       0       0       0       0       
> 0       0.657056999936078912    -0.0556191504001617432  2.2996994997762763      
> -0.194667026400566073   0.985585499904118478    -0.0834287256002426147  : 
> -1.79769313486231571e+308,-7.16309631969488692
> 0       0       0       0       0       0       0       0       0       0       
> 0       -0.657056999936078912   0.0556191504001617432   -2.2996994997762763     
> 0.194667026400566073    -0.985585499904118478   0.0834287256002426147   : 
> -1.79769313486231571e+308,11.3428432442749312
> -1.0011489951401733     0       0       0       0       0       0       0       
> 0       0       0       0       0       0       0       0       0       : 
> -1.79769313486231571e+308,-0
> 0.151777076402054012    0.995155720849257519    -0.0983111960103961735  
> 1.99031144169851504     -0.196622392020792347   0       0       0       0       
> 0       0       0       0       0       0       0       0       : 
> -1.79769313486231571e+308,0
> 0.110233680363780617    -0.999937238194396172   0.0112035562284168166   
> -1.99987447638879234    0.0224071124568336333   0       0       0       0       
> 0       0       0       0       0       0       0       0       : 
> -1.79769313486231571e+308,0
> 0       -1.98887754788353255    0.210632617428514662    0       0       
> 1.98887754788353255     -0.210632617428514662   0       0       0       0       
> 0       0       0       0       0       0       : -1.79769313486231571e+308,0
> 0       1.99966711057415236     -0.0364889968075429449  0       0       
> -1.99966711057415236    0.0364889968075429449   0       0       0       0       
> 0       0       0       0       0       0       : -1.79769313486231571e+308,0
> 0       0       0       -1.98728454107014496    0.225166944340467912    0       
> 0       1.98728454107014496     -0.225166944340467912   0       0       0       
> 0       0       0       0       0       : -1.79769313486231571e+308,0
> 0       0       0       1.99934691568998035     -0.051106855909585816   0       
> 0       -1.99934691568998035    0.051106855909585816    0       0       0       
> 0       0       0       0       0       : -1.79769313486231571e+308,0
> 0       0       0       0       0       -1.98564529350688668    
> 0.239191906999276949    0       0       1.98564529350688668     
> -0.239191906999276949   0       0       0       0       0       0       : 
> -1.79769313486231571e+308,0
> 0       0       0       0       0       1.99893626199613905     
> -0.065221319190149174   0       0       -1.99893626199613905    
> 0.065221319190149174    0       0       0       0       0       0       : 
> -1.79769313486231571e+308,0
> 0       0       0       0       0       0       0       -1.98442703642495855    
> 0.249097043550612074    0       0       1.98442703642495855     
> -0.249097043550612074   0       0       0       0       : 
> -1.79769313486231571e+308,0
> 0       0       0       0       0       0       0       1.9985859302833584      
> -0.0751949418072969805  0       0       -1.9985859302833584     
> 0.0751949418072969805   0       0       0       0       : 
> -1.79769313486231571e+308,0
> 0       0       0       0       0       0       0       0       0       
> -1.9837397586267167     0.25451241628575888     0       0       
> 1.9837397586267167      -0.25451241628575888    0       0       : 
> -1.79769313486231571e+308,0
> 0       0       0       0       0       0       0       0       0       
> 1.99837324861762089     -0.080649607621220587   0       0       
> -1.99837324861762089    0.080649607621220587    0       0       : 
> -1.79769313486231571e+308,0
> 0       0       0       0       0       0       0       0       0       0       
> 0       -1.98361771971734635    0.255461821067954564    -0.991808859858673175   
> 0.127730910533977282    2.97542657957601975     -0.383192731601931846   : 
> -1.79769313486231571e+308,0
> 0       0       0       0       0       0       0       0       0       0       
> 0       1.99833442018210539     -0.0816060360233772097  0.999167210091052693    
> -0.0408030180116886049  -2.99750163027315786    0.122409054035065815    : 
> -1.79769313486231571e+308,0
> 0       0       0       0       0       0       0       0       0       0       
> 0       0       0       0       0       -0.991830567901905269   
> 0.127562238046312809    : -1.79769313486231571e+308,11.029104867750366
> 0       0       0       0       0       0       0       0       0       0       
> 0       0       0       0       0       0.999174134752683019    
> -0.040633095393131792   : -1.79769313486231571e+308,6.41618650102796551
>
> All the columns are free.  The objective function takes all columns with weight 
> 1.  I am trying to solve the objective function.
>
> After running glp_simplex(lp, NULL); I read the column values via 
> glp_get_col_prim().
> The values are the following:
>
> 39.3951087338556718
> 1.65434914606154071
> 10.4852585667029814
> 2.00638515206249668
> 53.8500338174320206
> 5.10799011664556435
> 73.7645217616148869
> 6.63476476684575633
> 100.122219274645659
> 8.77186408234085313
> 124.323763060647408
> 10.8665289705831594
> 150.647055883489259
> 12.8743226374663653
> 170.619978638680635
> 14.5327174097183143
> 199.456348673906149

39.3951087338557 * 0.250281251632591 = 9.85985712211141
1.65434914606154 * 0.031972233217877 = 0.0528932367216753
10.485258566703 * 0.58245648256401 = 6.10720682353598
2.0063851520625 * 0.0091349237765363 = 0.018328175430465
53.850033817432 * 0.166416137875432 = 8.96151465235842

Which adds up to 24.9998000101579

Best regards

Heinrich

> Substituting these values back to the first row I am getting 
> 46.3759814591898873.  Which is quite far from the equality constraint.
>
> Imre Palik