Re: [Help-glpk] MIP Solvers (i.e CBC, CPLEX, GLPK, GUROBI, LPSOLVE, SCIPC, SCIPL, SCIPS and XPRESS) Benchmark

[Noli Sicad]

Hi Jiadong,

> You can post your question about Cbc to the mailing list of Cbc
>
> http://list.coin-or.org/mailman/listinfo/cbc

Yes, I will post some questions in CBC ML later on.

It seems that the generated FreeMPS and MPS could not be solved as
well in GLPK/Mathrog CLI (i.e. glpsol).

However, nothing wrong with GLPK/MathProg generated FreeMPS and MPS
MIP models. I check the other MIP models in example (e.g. bpp.mod)

glpsol --math bpp.mod --wfreemps bpp_fmps.mps

And this

glpsol --freemps  bpp_fmps.mps

Got the expected results.

However, this problem (below). The objective function = 0

###################
Nolis-MacBook-Pro:Case_Studies nsicad$ glpsol --mps
TimberHarvestModel_0025p_mps.mps
GLPSOL: GLPK LP/MIP Solver, v4.47
Parameter(s) specified in the command line:
 --mps TimberHarvestModel_0025p_mps.mps
Reading problem data from `TimberHarvestModel_0025p_mps.mps'...
Problem: TimberHa
Objective: R0000001
464 rows, 255 columns, 1421 non-zeros
249 integer variables, all of which are binary
1727 records were read
GLPK Integer Optimizer, v4.47
464 rows, 255 columns, 1421 non-zeros
249 integer variables, all of which are binary
Preprocessing...
418 rows, 255 columns, 1181 non-zeros
249 integer variables, all of which are binary
Scaling...
 A: min|aij| =  9.750e-01  max|aij| =  1.022e+04  ratio =  1.048e+04
GM: min|aij| =  4.933e-01  max|aij| =  2.027e+00  ratio =  4.109e+00
EQ: min|aij| =  2.441e-01  max|aij| =  1.000e+00  ratio =  4.097e+00
2N: min|aij| =  1.233e-01  max|aij| =  1.491e+00  ratio =  1.209e+01
Constructing initial basis...
Size of triangular part = 418
Solving LP relaxation...
GLPK Simplex Optimizer, v4.47
418 rows, 255 columns, 1181 non-zeros
*     0: obj =   3.468641975e+03  infeas =  0.000e+00 (0)
*    10: obj =   0.000000000e+00  infeas =  0.000e+00 (0)
OPTIMAL SOLUTION FOUND
Integer optimization begins...
+    10: mip =     not found yet >=              -inf        (1; 0)
+    10: >>>>>   0.000000000e+00 >=   0.000000000e+00   0.0% (1; 0)
+    10: mip =   0.000000000e+00 >=     tree is empty   0.0% (0; 1)
INTEGER OPTIMAL SOLUTION FOUND
Time used:   0.0 secs
Memory used: 0.6 Mb (614630 bytes)
Nolis-MacBook-Pro:Case_Studies nsicad$

############

Here's the real GLPK/MatProg result (below) of
TimberHarvestModel_0025p.mod (below).

The generated CPLEX lp MIP model by glpsol seems fine.

Anybody wants the CPLEX LP model generated by glpsol? You can download
it in this link.

https://gist.github.com/4108883

Thanks.

Noli


######

Nolis-MacBook-Pro:Case_Studies nsicad$ glpsol --math
TimberHarvestModel_0025p.mod --wglp Timber_glp.lp
GLPSOL: GLPK LP/MIP Solver, v4.47
Parameter(s) specified in the command line:
 --math TimberHarvestModel_0025p.mod --wglp Timber_glp.lp
Reading model section from TimberHarvestModel_0025p.mod...
Reading data section from TimberHarvestModel_0025p.mod...
599 lines were read

.

Generating Evenflow_Harvest_Volume_HV_Alpha...
Generating Evenflow_Harvest_Volume_HV_Beta...
Generating Adjacent_Constraint...
Model has been successfully generated
GLPK Integer Optimizer, v4.47
464 rows, 255 columns, 1421 non-zeros
249 integer variables, all of which are binary
Preprocessing...
418 rows, 255 columns, 1181 non-zeros
249 integer variables, all of which are binary
Scaling...
 A: min|aij| =  9.750e-01  max|aij| =  1.022e+04  ratio =  1.048e+04
GM: min|aij| =  4.933e-01  max|aij| =  2.027e+00  ratio =  4.109e+00
EQ: min|aij| =  2.441e-01  max|aij| =  1.000e+00  ratio =  4.097e+00
2N: min|aij| =  1.233e-01  max|aij| =  1.491e+00  ratio =  1.209e+01
Constructing initial basis...
Size of triangular part = 418
Solving LP relaxation...
GLPK Simplex Optimizer, v4.47
418 rows, 255 columns, 1181 non-zeros
*     0: obj =   3.468641975e+03  infeas =  0.000e+00 (0)
*   204: obj =   3.501953549e+05  infeas =  1.776e-15 (0)
OPTIMAL SOLUTION FOUND
Integer optimization begins...
+   204: mip =     not found yet <=              +inf        (1; 0)
+ 17013: >>>>>   3.141913580e+05 <=   3.486538570e+05  11.0% (3503; 74)
+ 46887: >>>>>   3.415597531e+05 <=   3.484158711e+05   2.0% (10089; 341)
+ 52871: >>>>>   3.427400000e+05 <=   3.484158711e+05   1.7% (7912; 7000)
+ 53854: >>>>>   3.441562963e+05 <=   3.484158711e+05   1.2% (7209; 8753)
+106063: mip =   3.441562963e+05 <=   3.481560380e+05   1.2% (15588; 11915)
+128901: >>>>>   3.448318519e+05 <=   3.480337672e+05   0.9% (19388; 12300)
+167570: mip =   3.448318519e+05 <=   3.476660019e+05   0.8% (22566; 17049)
+203556: mip =   3.448318519e+05 <=   3.476660019e+05   0.8% (26936; 18135)
+234790: mip =   3.448318519e+05 <=   3.476660019e+05   0.8% (31096; 18959)
+261625: mip =   3.448318519e+05 <=   3.476642484e+05   0.8% (34977; 19593)
+287744: mip =   3.448318519e+05 <=   3.476615391e+05   0.8% (38722; 20168)
+309101: mip =   3.448318519e+05 <=   3.476264847e+05   0.8% (41288; 20781)
+326960: mip =   3.448318519e+05 <=   3.475600735e+05   0.8% (43029; 21420)
+348232: mip =   3.448318519e+05 <=   3.475600735e+05   0.8% (45593; 21997)
+363004: mip =   3.448318519e+05 <=   3.472124818e+05   0.7% (46683; 22656)
Time used: 60.0 secs.  Memory used: 61.0 Mb.

.

.

Time used: 2220.4 secs.  Memory used: 251.0 Mb.
+2021285: mip =   3.455566667e+05 <=   3.460048582e+05   0.1% (30924; 342242)
+2027054: mip =   3.455566667e+05 <=   3.460048582e+05   0.1% (29983; 344881)
+2033045: mip =   3.455566667e+05 <=   3.460048582e+05   0.1% (29011; 347508)
+2039477: mip =   3.455566667e+05 <=   3.460048582e+05   0.1% (27974; 350209)
+2046348: mip =   3.455566667e+05 <=   3.460048582e+05   0.1% (26964; 353021)
+2053682: mip =   3.455566667e+05 <=   3.460048582e+05   0.1% (25910; 356029)
+2061480: mip =   3.455566667e+05 <=   3.460048582e+05   0.1% (24786; 359225)
+2070126: mip =   3.455566667e+05 <=   3.460048582e+05   0.1% (23608; 362694)
+2079241: mip =   3.455566667e+05 <=   3.460048582e+05   0.1% (22294; 366177)
+2089328: mip =   3.455566667e+05 <=   3.460048582e+05   0.1% (20946; 370068)
+2100109: mip =   3.455566667e+05 <=   3.460001615e+05   0.1% (19351; 375226)
+2110555: mip =   3.455566667e+05 <=   3.459963620e+05   0.1% (17450; 381677)
Time used: 2280.4 secs.  Memory used: 251.0 Mb.
+2120248: mip =   3.455566667e+05 <=   3.459332288e+05   0.1% (15117; 391249)
+2128163: mip =   3.455566667e+05 <=   3.459332288e+05   0.1% (11986; 406653)
+2139212: mip =   3.455566667e+05 <=   3.459332288e+05   0.1% (7739; 427947)
+2156837: mip =   3.455566667e+05 <=     tree is empty   0.0% (0; 490335)
INTEGER OPTIMAL SOLUTION FOUND
Time used:   2297.0 secs
Memory used: 251.1 Mb (263266271 bytes)
#################################


Discounted Net Value = 20864096.52


STAND  HarvestVolume
1   20393.14
2   19963.78
3   19620.31
4   19130.90
5   18798.87
6   18349.14
Writing resultCut...