Re: [Bug-glpk] Error using pyomo with glpk 4.60

[Andrew Makhorin]

On Thu, 2016-06-16 at 21:24 +0300, Andrew Makhorin wrote:
> On Thu, 2016-06-16 at 14:28 +0000, Alaggia Rodino, Sebastian wrote:
> > Hello,
> > 
> >  
> > 
> > I’m using pyomo 3.4.4 in windows 7 64bits with glpk 4.60. While trying
> > to solve an example I get the following errors:
> > 
> >  
> > 
> > D:\ampl\ej_desp_el>pyomo solve --solver=glpk diet.py diet.dat
> > --stream-solver -c
> > 
> >  
> > 
> > [    0.00] Setting up Pyomo environment
> > 
> > [    0.00] Applying Pyomo preprocessing actions
> > 
> > [    0.00] Creating model
> > 
> > [    0.02] Applying solver
> > 
> > GLPSOL: GLPK LP/MIP Solver, v4.60
> > 
> > Parameter(s) specified in the command line:
> > 
> > --write C:\Users\c288093\AppData\Local\Temp\tmpc39o0gw9.glpk.raw
> > --wglp C:\User
> > 
> > s\c288093\AppData\Local\Temp\tmpb9ljgvf8.glpk.glp
> > 
> > --cpxlp C:\Users\c288093\AppData\Local\Temp\tmp0mgqyiyv.pyomo.lp
> > 
> > Reading problem data from 'C:\Users\c288093\AppData\Local\Temp
> > \tmp0mgqyiyv.pyomo.lp'...
> > 
> > 9 rows, 9 columns, 63 non-zeros
> > 
> > 114 lines were read
> > 
> > Writing problem data to 'C:\Users\c288093\AppData\Local\Temp
> > \tmpb9ljgvf8.glpk.glp'...
> > 
> > 109 lines were written
> > 
> > GLPK Simplex Optimizer, v4.60
> > 
> > 9 rows, 9 columns, 63 non-zeros
> > 
> > Preprocessing...
> > 
> > 8 rows, 8 columns, 62 non-zeros
> > 
> > Scaling...
> > 
> > A: min|aij| = 8.000e+000  max|aij| = 7.000e+001  ratio = 8.750e+000
> > 
> > GM: min|aij| = 5.017e-001  max|aij| = 1.993e+000  ratio = 3.972e+000
> > 
> > EQ: min|aij| = 2.549e-001  max|aij| = 1.000e+000  ratio = 3.922e+000
> > 
> > Constructing initial basis...
> > 
> > Size of triangular part is 8
> > 
> >       0: obj =  0.000000000e+000 inf =  2.845e+003 (4)
> > 
> >       3: obj =  1.185625000e+002 inf =  0.000e+000 (0)
> > 
> > *     6: obj =  8.820000000e+001 inf =  1.123e-013 (0)
> > 
> > OPTIMAL LP SOLUTION FOUND
> > 
> > Time used:   0.0 secs
> > 
> > Memory used: 0.0 Mb (42663 bytes)
> > 
> > Writing basic solution to 'C:\Users\c288093\AppData\Local\Temp
> > \tmpc39o0gw9.glpk.raw'...
> > 
> > 27 lines were written
> > 
> > invalid literal for int() with base 10: 'c'
> > 
> > [    0.05] Pyomo Finished
> > 
> > ERROR: Unexpected exception (PermissionError) while running model
> > diet.py:
> > 
> >        [WinError 32] El proceso no tiene acceso al archivo porque est
> > siendo utilizado por otro proceso: 'C:\\Users\\c288093\\AppData\\Local
> > \\Temp\\tmpc39o0gw9.glpk.raw'
> > 
> >  
> > 
> > Traceback (most recent call last):
> > 
> >   File "c:\python34\lib\runpy.py", line 170, in _run_module_as_main
> > 
> >     "__main__", mod_spec)
> > 
> >   File "c:\python34\lib\runpy.py", line 85, in _run_code
> > 
> >     exec(code, run_globals)
> > 
> >   File "C:\Python34\Scripts\pyomo.exe\__main__.py", line 9, in
> > <module>
> > 
> >     sys.exit(main())
> > 
> >   File "c:\python34\lib\site-packages\pyomo\scripting\pyomo_main.py",
> > line 71, in main
> > 
> >     retval = ret[0].func(*ret)
> > 
> >   File "c:\python34\lib\site-packages\pyomo\scripting\plugins
> > \solve.py", line 214, in solve_exec
> > 
> >     name='pyomo solve')
> > 
> >   File "c:\python34\lib\site-packages\pyutilib\workflow\task.py", line
> > 109, in __call__
> > 
> >     self.execute()
> > 
> >   File "c:\python34\lib\site-packages\pyomo\util\_task.py", line 160,
> > in execute
> > 
> >     retval = self._fn(**self._kwds)
> > 
> >   File "c:\python34\lib\site-packages\pyomo\scripting\util.py", line
> > 950, in run_command
> > 
> >     TempfileManager.pop()
> > 
> >   File "c:\python34\lib\site-packages\pyutilib\component\config
> > \tempfiles.py", line 177, in pop
> > 
> >     os.remove(filename)
> > 
> >  
> > 
> > D:\ampl\ej_desp_el>
> > 
> >  
> > 
> > I’ve compile glpk with gcc 4.9.2 (mingw) and 5.3.0 (Cygwin). Also I’ve
> > tried with glpk 4.59, 4.58 and got the same errors. The error
> > disappear in version 4.57. This is happening with other users
> > (https://groups.google.com/forum/#!topic/pyomo-forum/AbpBqHJzdqo and
> > https://github.com/FRED-2/OptiType/issues/28 ).
> > 
> > Pyomo is an algebraic modeling language which supports multiple
> > solvers. You can download pyomo 3.4.4 from
> > http://www.pyomo.org/installation/ 
> > 
> >  
> > 
> > Thanks,
> > 
> > Sebastin
> > 
> > 
> 
> The error occurs in pyomo, not in glpk. Please report it to the pyomo
> developers/maintainers or to an appropriate mailing list. Thanks.
> 
> 
> Andrew Makhorin

Just for your information:

In glpk 4.58 (released on Feb 2016) the solution format was changed.

[from ChangeLog]
        The following source modules were changed due to changes in the
        data formats:
        * src/api/rdsol.c
        * src/api/wrsol.c
        * src/api/rdipt.c
        * src/api/wript.c
        * src/api/rdmip.c
        * src/api/wrmip.c

Now these routines use a more consistent DIMACS-like format rather than
a "raw" format used before.

Please ask the pyomo developers/maintainers to make corresponding
changes in the routines that read glpk solution files. Thanks.


Andrew Makhorin