4

嗯，这就是我今天的问题...

我正在编写一个模块过程，该过程具有作为参数的函数。这个模块看起来像这样：

module Integ
    implicit none
     <variables declaration>
contains
    function Integral(a,b,f) result(res)
        real, intent(in)     ::a, b
        real                 ::res

        interface
            pure function f(x)
                real, intent(in) :: x
                real             :: f
            endfunction
        endinterface


     <more code of function Integral>

    endfunction Integral

endmodule Integ

嗯，到这里为止，一切都很好。当我尝试将此函数与Fortran 内在函数一起使用时，就会出现此问题。即，在这段代码中：

program main

use Integ

implicit none

real   ::res,a,b

a=3.0; b=4.0

res=Integral(a,b,sin)  !<- This line does not work

!res=Integral(a,b,sen) !<- This line does work

contains
    function sen(x)
        real, intent(in)   :: x
        real               :: sen

        sen=sin(x)
    endfunction

endprogram

第一行不起作用，给出错误消息：

main.f90(17): error #6404: This name does not have a type, and must have an explicit type.   [SIN]
r=Int1DMonteCarlo(0.0,1.0,sin,10000)
--------------------------^

main.f90(17): error #6637: This actual argument must be the name of an external user function or the name of an intrinsic function.   [SIN]
r=Int1DMonteCarlo(0.0,1.0,sin,10000)
--------------------------^

但是第二行（在片段中注释）确实如此。

这些错误对我来说非常令人迷惑，因为sin它是一个 Fortran 内在函数（与错误 nr 2 相矛盾的东西），因此在每个范围内都是显式的（与错误 nr 1 相矛盾的东西）。

显然我做错了什么，但我不知道是什么。

所以我想问：

  • 可以调用具有内部函数作为实际参数的模块过程吗？
  • 除了在过程中声明接口之外，我还失去了什么？

如果您有兴趣，这是模块的完整源代码，这是主要的源代码

对不起，如果我问一个愚蠢的问题。我想我正在按照我现在正在阅读的书籍（Metcalf，Fortran V:II 的数字食谱）告诉我的方式做事。

感谢您的时间！

4

1 回答 1

5

在主程序中使用内在语句来声明实际参数sin是内在的。Fortran 标准中对内在属性的描述中详细说明了这一要求。

着眼于未来，您最好围绕内在函数编写自己的包装函数 - 创建一个简单调用 sin 的函数 mysin。

于 2013-10-08T20:14:33.967 回答