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
>