提问者:小点点

在pyomo中优化Fortran函数


我希望使用Pyomo优化Fortran函数。目标函数和约束都是用Fortran语言编写的。根据这里给出的答案,我们可以使用ExternalFunctionexpression对象。但即使对于最简单的函数,我也无法得到结果。下面给出了一个可复制的示例,其中包括Fortran函数、python(python 2.7.12)脚本、为优化而执行的命令和错误。

Fortran函数文件(funcs.f)-

  SUBROUTINE OBJ1(ARG,OBJ)
  DOUBLE PRECISION, INTENT(IN) :: ARG(2)
  DOUBLE PRECISION, INTENT(OUT) :: OBJ
  OBJ = ARG(1)+ARG(2)
  RETURN
  END SUBROUTINE

Python脚本(pytest.py)-

import funcs
from pyomo.environ import *
from pyomo.opt import *
from pyomo.core import *

m = ConcreteModel()
m.a = Var(RangeSet(1,2),within=NonNegativeReals,bounds=(0,10))
m.f = ExternalFunction(library='funcs.so',function='OBJ1')
expr = m.f(m.a)
m.obj = Objective(expr=expr,sense=minimize)

opt = SolverFactory('ipopt')
results = opt.solve(m,tee=True)
print(results)

在终端执行的命令-

>> f2py -c -m funcs funcs.f
>> python pytest.py

错误-

File "/usr/local/lib/python2.7/dist-packages/pyomo/core/base/external.py", line 160, in load_library
FUNCADD(('funcadd_ASL', self._so))(byref(AE))
AttributeError: /home/utkarsh/Desktop/python/modules/blackboxOptimization/funcs.so: undefined symbol: funcadd_ASL

我只给出了我认为相关的一小部分错误。

有鉴于此,我对以下问题有一个初步的看法-

>

对于我的整个项目,我必须用Fortran本身给出约束条件。约束子例程将返回一个实数,该实数将使用pyomo进行限定。如何对这些类型的约束进行建模?

我假设Pyomo将这个黑盒作为非线性优化。因此,我使用的是ipoptsolver。这个假设正确吗?

软件包的版本-

>

  • 俾莫-5.5。1(VOTD)(Linux 4.4.0-127-generic上的CPython 2.7.12)

    ipopt-ipopt 3.12。8.

    f2py-与numpy 1.16一起安装。2.

    谢谢你的帮助!


  • 共1个答案

    匿名用户

    如果你没有绑定Pyomo,你可以使用优秀的Pygmo包,它为不同类型的问题提供解决方案,包括黑盒解决方案。

    下面是一个关于如何在连续约束单目标测试问题上使用它的小示例:

    import pygmo as pg
    import pandas as pd
    
    
    class Rosenbrock():
        """Rosenbrock function constrained to a disk.
    
        See: https://en.wikipedia.org/wiki/Test_functions_for_optimization
        """
    
        def fitness(self, x):
            """Evaluate fitness.
    
            Instead of the Rosenbrock function you could call your Fortran
            code here e.g. by using F2PY: https://www.numfys.net/howto/F2PY/
            """
            obj = (1-x[0])**2+100*(x[1]-x[0]**2)**2
            ineq = x[0]**2+x[1]**2-2
            return [obj, ineq]
    
        def get_bounds(self):
            """Return boundaries."""
            return ([-1.5]*2, [1.5]*2)
    
        def get_nic(self):
            """Determine number of inequalities."""
            return 1
    
    
    # set up and solve problem
    pro = pg.problem(Rosenbrock())
    pop = pg.population(pro, size=200)
    # see: https://github.com/esa/pagmo2/blob/master/include/pagmo/algorithms/
    algo = pg.algorithm(pg.ihs(gen=10000))
    algo.set_verbosity(100)
    pop = algo.evolve(pop)
    
    # extract solutions
    fits = pd.DataFrame(pop.get_f())
    vectors = pd.DataFrame(pop.get_x())
    best_idx = pop.best_idx()
    best_vector = vectors.loc[best_idx].to_frame().T
    best_fitness = fits.loc[best_idx].to_frame().T
    
    print(best_vector)
    print(best_fitness)
    

    然后,您只需在fitness函数中处理Fortran代码的“接口”。

    希望这有帮助!