f2py：在与python交互时在fortran中指定实际精度吗？

f2py：在与python交互时在fortran中指定实际精度吗？

在您的第一个示例中，我不知道您为什么说f2py似乎不接受双精度，

test.test.r_dp

而
双精度却如此。一个显示带有小数点的值且没有显式dtype的numpy数组是一个双精度数组。

第二个示例显示了F2PY用处理类型定义的局限性

kind=<kind>

。请参阅常见问题解答：http :
//cens.ioc.ee/projects/f2py2e/FAQ.html#q-what-if-fortran-90-pre-uses-type-
spec-kind-kind

要查看发生了什么，请运行

f2py test.f90 -m test

。我得到这个：

Reading fortran pres...    Reading file 'test.f90' (format:free)Post-processing...    Block: test Block: test     Block: input_sp     Block: input_dpPost-processing (stage 2)...    Block: test        Block: unknown_interface Block: test     Block: input_sp     Block: input_dpBuilding modules...    Building module "test"...        Constructing F90 module support for "test"...          Variables: r_dp sp r_sp dp Constructing wrapper function "test.input_sp"...getctype: "real(kind=sp)" is mapped to C "float" (to override define dict(real = dict(sp="<C typespec>")) in /Users/warren/tmp/.f2py_f2cmap file).getctype: "real(kind=sp)" is mapped to C "float" (to override define dict(real = dict(sp="<C typespec>")) in /Users/warren/tmp/.f2py_f2cmap file).getctype: "real(kind=sp)" is mapped to C "float" (to override define dict(real = dict(sp="<C typespec>")) in /Users/warren/tmp/.f2py_f2cmap file).   input_sp(val) Constructing wrapper function "test.input_dp"...getctype: "real(kind=dp)" is mapped to C "float" (to override define dict(real = dict(dp="<C typespec>")) in /Users/warren/tmp/.f2py_f2cmap file).getctype: "real(kind=dp)" is mapped to C "float" (to override define dict(real = dict(dp="<C typespec>")) in /Users/warren/tmp/.f2py_f2cmap file).getctype: "real(kind=dp)" is mapped to C "float" (to override define dict(real = dict(dp="<C typespec>")) in /Users/warren/tmp/.f2py_f2cmap file).   input_dp(val)    Wrote C/API module "test" to file "./testmodule.c"    Fortran 90 wrappers are saved to "./test-f2pywrappers2.f90"

请注意，它会将“ real（kind = sp）”和“ real（kind = dp）”都映射到C“ float”，这是单精度的。

有几种方法可以解决此问题。

方法1

将类型声明分别更改为“ real（kind = 4）”和“ real（kind = 8）”（或“ real * 4”和“ real * 8”）。

当然，这违背了使用的目的

selected_real_kind

，对于某些编译器，对于单精度和双精度而言，4和8不是正确的KIND值。在这种情况下，使用gfortran

sp

为4和

dp

8，因此可以正常工作。

方法二

告诉f2py如何处理这些声明。在f2py常见问题解答中对此进行了说明，这是上面显示的f2py输出中的“ getctype：…”消息中建议的方法。

在这种情况下，您将创建一个名为

.f2py_f2cmap

（在运行f2py的目录中）的文件，其中包含以下行

dict(real=dict(sp='float', dp='double'))

然后f2py将对那些

real(sp)

和

real(dp)

声明做正确的事情。

方法3

它还可以重新排列代码：

module typesimplicit noneinteger, parameter :: sp = selected_real_kind(6,37) ! single precisioninteger, parameter :: dp = selected_real_kind(15,307) ! double precisionreal(sp) :: r_sp = 1.0real(dp) :: r_dp = 1.0_dpend modulemodule inputcontainssubroutine input_sp(val)  use types  real(sp), intent(in) :: val  real(sp) :: x  x = val  write(*,*) xend subroutinesubroutine input_dp(val)  use types  real(dp), intent(in) :: val  real(dp) :: x  x = val  write(*,*) dp, val, xend subroutineend module