如何使外部功能与F2PY一起使用？

Question

我试图用f2py编译一块旧的Fortran代码，以便可以在Python中调用。 但是，具有外部功能的零件是行不通的。 这是MWE，首先我有test.f：

      function f(x)
      implicit double precision (a-z)

      f = x * x

      return
      end function f

      subroutine gauss(fun)
      implicit double precision (a-h, j-z)
      ! external fun

      x = 1.5
      write(*,*) fun(x)

      return
      end subroutine gauss

后来用makefile：

f2py -c --quiet --fcompiler=gnu95 \
        --f90flags=“-Wtabs” \
        -m test \
        test.f

最后是从python调用的：

import test
f = lambda x: x
test.gauss(test.f)

给出错误typeerror： test.gauss（）1st参数（fun）不能转换为double。

在 second 尝试中，我在subroutine gauss中删除行外部有趣，并在编译过程中获取以下错误消息：

/tmp/tmpet9sk3e9/src.linux-x86_64-3.7/testmodule.c: In function ‘cb_fun_in_gauss__user__routines’:
/tmp/tmpet9sk3e9/src.linux-x86_64-3.7/testmodule.c:313:8: error: variable or field ‘return_value’ declared void

我用完了想法，任何帮助将不胜感激！

Answer 1

我（部分）弄清楚了：f2py需要您明确指定功能的使用。
以下是：

      function f(x)
      implicit double precision (a-z)

      f = x * x

      return
      end function f

      subroutine gauss(fun)
      implicit double precision (a-h, j-z)
      !! <-------- updates starts HERE
      external fun
!f2py real*8 y
!f2py y = fun(y)
      double precision fun
      !! <-------- updates ends HERE

      x = 1.5
      write(*,*) fun(x)

      return
      end subroutine gauss

进行编译

main:
    f2py -c --quiet --fcompiler=gnu95 \
        --f90flags=“-Wtabs” \
        -m test \
        test.f

然后使用您可以测试

import test
f = lambda x: x
test.gauss(test.f)
test.gauss(f)

，并看到外部功能适用于Fortran和Python函数。

---

两个旁注：

1. 尽管文档> ）有趣的是必要的，我发现代码可以在没有它的情况下工作。
2. 您可以使用相同的语法指定多个外部函数，甚至可以为所有这些语法使用相同的虚拟输入变量y为所有这些变量（虽然不是一个好练习）。