Question

我目前正在将函数作为参数传递。

在下面的程序中，我使用内置函数EXP作为参数积分函数。我的编译器给了我以下错误：

integrate1.f90:22.26:

r = integral(-1.0,1.0,EXP,1000);
                      1
Error: Expected a procedure for argument 'f' at (1)

如果我在变量EXP的声明中取消注释r的用法，我就不会收到此错误。

所以看来，如果我不使用内置函数，我就不能用它作为参数这有点奇怪，导致“内置”的一种暗示函数无论如何都被加载。

如果不明确使用EXP功能，如何防止此错误？我是否需要使用USE语句来加载内置的？如果没有别的方法，我有兴趣知道这是由于Fortran规范还是编译器问题？

我正在使用GNU Fortran (Ubuntu/Linaro 4.7.3-1ubuntu1) 4.7.3。

示例：

MODULE MINTEGRATE
CONTAINS
FUNCTION integral(from,to,f,n)
    INTERFACE
        FUNCTION f(y); REAL, INTENT(IN) :: y; END FUNCTION
    END INTERFACE
    REAL :: from,to,integral,width;
    INTEGER :: n;

    width=ABS(to-from)/n;
    integral = 0.0;
    DO i=0,n
        integral = integral+f(from+width*i)*width;
    END DO
END
END

PROGRAM INTEGRATE
USE MINTEGRATE;

!PROCEDURE(EXP), POINTER :: f => EXP; ! using the variable f below works without error
REAL :: r!=EXP(0.0);
r = integral(-1.0,1.0,EXP,1000);
WRITE(*,*) r;
END

Answer 1

“当一个内部函数作为一个实际参数传递给一个过程时，必须使用它的特定名称，并且在调用时，它的参数必须是标量。并非所有特定的内部函数都可以作为实际参数出现。（更多信息，见 Intrinsic Functions Not Allowed as Actual Arguments）。“来自英特尔Fortran手册。

你很幸运，因为单精度real的泛型exp的具体名称也是exp，但是要小心并传递正确的特定函数，或者编写自己的包装器来调用通用名称。例如，如果您想要默认的实对数，则必须使用alog。

您可以通过以下方式通知编译器您的意思是内在函数：

intrinsic exp

放在主程序的声明中。

是否可以使用USE语句加载内置函数？

1 个答案: