[Help-glpk] (no subject)

[Epiphanie Nyirabahizi/nyirabahizie/O/VCU]

Hello all,

what has been changed into the glpk package that makes my code erroneous?

The code used to run without any problem, now, if I run it the element of my 
matrix are not recognizable.


library(glpk)
lp <- lpx_create_prob()
lpx_set_prob_name(lp, "Endpoint matter")
lpx_set_obj_dir(lp, LPX_MIN)

lpx_add_rows(lp,4)
lpx_set_row_name(lp, 1, "br")
lpx_set_row_bnds(lp, 1, LPX_UP,0.0,514)
#lpx_set_row_bnds(lp, 1,LPX_FX , .514,.514)

lpx_set_row_name(lp, 2, "bl")
#lpx_set_row_bnds(lp, 2, LPX_UP,0.0,.174)
lpx_set_row_bnds(lp, 2,LPX_FX ,.174,.174)

lpx_set_row_name(lp, 3, "ac")
lpx_set_row_bnds(lp, 3, LPX_UP,0.0,.53)
#lpx_set_row_bnds(lp, 3,LPX_FX , .53,.53)


lpx_set_row_name(lp, 4, "ta")
lpx_set_row_bnds(lp, 4, LPX_UP,0.0, 1.52)
#lpx_set_row_bnds(lp, 4,LPX_FX , 1.52,1.52)

 


lpx_add_cols(lp,2)
lpx_set_col_name(lp, 1, "x1")
lpx_set_col_bnds(lp, 1, LPX_LO, 0.0, 0.0)
#lpx_set_obj_coef(lp, 1, -1)

lpx_set_col_name(lp, 2, "x2")
lpx_set_col_bnds(lp, 2, LPX_LO, 0.0, 0.0)
#lpx_set_obj_coef(lp, 2, -1)

     ia[1] = 1; ja[1] = 1; ar[1]= 5.934e-03; # a[1,1]=-5.934e-03
     ia[2] = 1; ja[2] = 2; ar[2]= 1.098e-01; # a[1,2]=-1.098e-01
     ia[3] = 2; ja[3] = 1; ar[3]= 4.309e-02; # a[2,1]=-4.309e-02
     ia[4] = 2; ja[4] = 2; ar[4]= 7.204e-02; # a[2,2]=-7.204e-02
     ia[5] = 3; ja[5] = 1; ar[5]= 2.227e-03; # a[3,1]=-2.227e-03
     ia[6] = 3; ja[6] = 2; ar[6]= 3.604e-02; # a[3,2]=-3.604e-02
     ia[7] = 4; ja[7] = 1; ar[7]= -2.580e-03; # a[4,1]=2.580e-03
     ia[8] = 4; ja[8] = 2; ar[8]= -7.392e-03; # a[4,2]=7.392e-03

lpx_load_matrix(lp, 8, ia, ja, ar);
lpx_simplex(lp);
LPX_OPT
LPX_FEAS
lpx_get_status(lp)

Hello all,

what has been changed into the glpk package that makes my code erroneous?

The code used to run without any problem, now, if I run it the element of my matrix are not recognizable.

library(glpk)
lp <- lpx_create_prob()
lpx_set_prob_name(lp, "Endpoint matter")
lpx_set_obj_dir(lp, LPX_MIN)

lpx_add_rows(lp,4)
lpx_set_row_name(lp, 1, "br")
lpx_set_row_bnds(lp, 1, LPX_UP,0.0,514)
#lpx_set_row_bnds(lp, 1,LPX_FX , .514,.514)

lpx_set_row_name(lp, 2, "bl")
#lpx_set_row_bnds(lp, 2, LPX_UP,0.0,.174)
lpx_set_row_bnds(lp, 2,LPX_FX ,.174,.174)

lpx_set_row_name(lp, 3, "ac")
lpx_set_row_bnds(lp, 3, LPX_UP,0.0,.53)
#lpx_set_row_bnds(lp, 3,LPX_FX , .53,.53)

lpx_set_row_name(lp, 4, "ta")
lpx_set_row_bnds(lp, 4, LPX_UP,0.0, 1.52)
#lpx_set_row_bnds(lp, 4,LPX_FX , 1.52,1.52)

 

lpx_add_cols(lp,2)
lpx_set_col_name(lp, 1, "x1")
lpx_set_col_bnds(lp, 1, LPX_LO, 0.0, 0.0)
#lpx_set_obj_coef(lp, 1, -1)

lpx_set_col_name(lp, 2, "x2")
lpx_set_col_bnds(lp, 2, LPX_LO, 0.0, 0.0)
#lpx_set_obj_coef(lp, 2, -1)

     ia[1] = 1; ja[1] = 1; ar[1]= 5.934e-03; # a[1,1]=-5.934e-03
     ia[2] = 1; ja[2] = 2; ar[2]= 1.098e-01; # a[1,2]=-1.098e-01
     ia[3] = 2; ja[3] = 1; ar[3]= 4.309e-02; # a[2,1]=-4.309e-02
     ia[4] = 2; ja[4] = 2; ar[4]= 7.204e-02; # a[2,2]=-7.204e-02
     ia[5] = 3; ja[5] = 1; ar[5]= 2.227e-03; # a[3,1]=-2.227e-03
     ia[6] = 3; ja[6] = 2; ar[6]= 3.604e-02; # a[3,2]=-3.604e-02
     ia[7] = 4; ja[7] = 1; ar[7]= -2.580e-03; # a[4,1]=2.580e-03
     ia[8] = 4; ja[8] = 2; ar[8]= -7.392e-03; # a[4,2]=7.392e-03

lpx_load_matrix(lp, 8, ia, ja, ar);
lpx_simplex(lp);
LPX_OPT
LPX_FEAS
lpx_get_status(lp)