2

我正在求解一组模拟方程（IMODE = 1，SOLVER = 3）。IPOPT 求解器求解到可接受的水平并退出，但 gekko 为此返回错误并返回我的解决方案。根据 IPOPT 文档，可接受级别的容差为 1.0e-6，这与 gekko 使用的 OTOL 和 RTOL 的默认值（以及我正在使用的值）相同。我能够修改 gekko.py 源代码以让我的答案返回，但这样做我绕过了所有类型的错误。我不希望绕过所有错误，因为它们显然有助于调试其他问题，例如不可行性。是否有我缺少的 m.solve 选项，或者当 IPOPT 解决到可接受的水平时不触发错误的另一种方法？

4

1 回答 1

2

处理求解器错误的一种方法是将求解命令包装在try,except语句中。APPINFO输出可能会指导您遇到的错误类型，并让您对“不可行的解决方案”、“解决到可接受的水平”或其他 IPOPT 错误代码做出不同的响应。

try:
   m.solve(disp=True)
except:
   print('Solver error, looking at APPINFO')
   if m.options.APPINFO==1:
      print('APPINFO=1')
   elif m.options.APPINFO==2:
      print('APPINFO=2')

另一种选择是尝试不同的求解器，例如 APOT 或 BPOPT。

m.options.SOLVER = 1

编辑：APPINFO Gekko 引发求解器异常时，不会更新参数。请尝试以下操作debug=0：

m.solve(disp=True,debug=0)
if m.options.APPINFO!=0:
   print('Solver error, looking at APPINFO')
if m.options.APPINFO==1:
   print('APPINFO=1')
elif m.options.APPINFO==2:
   print('APPINFO=2')

我刚刚更新了 Gekko，以便远程解决方案也将绕过引发的异常并使用APPINFO信息完成对选项文件的处理。本地运行时APPINFO信息在run目录下，用in读入。options.jsonload_JSONgk_post_solve.py

于 2019-09-20T14:10:01.023 回答