使用Python PYOMO中的GAMS / CPLEX

时间:2018-01-12 08:57:06

标签: pyomo gams-math

我注意到Pyomo 5.3提供了一个GAMS求解器插件。 https://github.com/Pyomo/pyomo/blob/master/pyomo/solvers/plugins/solvers/GAMS.py

这非常令人兴奋,因为我们有一个GAMS / CPLEX许可证,我们可以使用CPLEX作为求解器,但只能通过GAMS。使用新的Pyomo-Gams界面,从我的理解中可以在Pyomo中制定问题,并将其翻译为GAMS并由CPLEX解决。

然而,当我使用shell集成测试时,它非常慢(对于30个小型MIP的解决方案为40秒,而对于glpk / ipopt / cbc为6个)。此外,该插件的文档实际上是不存在的。

但是也许你们中的某个人有使用该界面的经验,可以帮助我

  • pyomo实际上是否将pyomo模型转换为gams代码?如果是,我在哪里可以找到gams-file?
  • 翻译的效率如何?如果我想重复解决一个小型模型,我该怎么办?
  • 使用shell或GAMS Python API有什么区别?
  • 有没有找到这方面文件的地方?

  • 此外,似乎conda仅为Linux / Python 3.6或Windows / Python 2.7提供Pyomo 5.3 https://anaconda.org/conda-forge/pyomo/files?version=5.3,所以我不得不使用pip在我的机器上安装Pyomo 5.3。

提前致谢,Theo

import pyomo.environ as pe

# set up the model
model = pe.ConcreteModel()

model.MaxWeight = pe.Param(initialize=0,mutable=True)
model.Item = ['hammer','wrench','screwdriver','towel']

Weight = {'hammer':5,'wrench':7,'screwdriver':4,'towel':3}
Value = {'hammer':8,'wrench':3,'screwdriver':6,'towel':11}

model.x = pe.Var(model.Item,within=pe.Binary)
model.z = pe.Objective(expr=sum(Value[i] * model.x[i] for i in model.Item),sense=pe.maximize)
model.constraint = pe.Constraint(expr=sum(Weight[i]*model.x[i] for i in model.Item) <= model.MaxWeight)

# time execution
solver_list = ['cbc', 'ipopt', 'gams', 'glpk']

for i, solver_name in enumerate(solver_list):
    solver = pe.SolverFactory(solver_name)
    print(solver_name)
    tic = time.time()
    for MaxWeight_i in range(0,30):
        model.MaxWeight = MaxWeight_i
        result = solver.solve(model)

        soln_items = list()
        for i in model.x:
            if pe.value(model.x[i]) > 0.5:
                soln_items.append(i)
        # print("Maximum Weight =", MaxWeight_i, soln_items)

    print("{:7.2f} s".format(time.time()-tic))
    print(" ")

1 个答案:

答案 0 :(得分:0)

这是相当延迟的,但我可以回答你的一些问题。

首先,刚刚为readthedocs上的GAMS界面创建了一个基本文档页面,您可以在http://pyomo.readthedocs.io/en/latest/library_reference/solvers/gams.html找到它。请注意,此位置可能会发生变化,因为我相信我们很快会重新构建文档树,但您应该能够搜索“gams”以便将来再次找到它。如果您认为自己或其他人希望看到更多文档,请告诉我,因为我很乐意提供任何有用的文章。

至于shell接口和Python API接口之间的区别,确实没有。我认为通过使用API​​会有性能提升,但过去似乎并非如此(事实上,我尝试过的一个模型看到shell界面更快)。如果你尝试两种方式并且不这样做,我也很高兴知道这一点。