传递字符串以在 Fortran 子例程中执行

Question

在下面的子程序中，我想传递一个名为str 的字符串变量。如果是'poly'、'gaus'、'slat'，那么它有一个预定义的操作（fval = 参见下面的代码）。我想让用户指定要使用的函数并将其作为字符串变量传递。

那是……

如果str = '3*cos(i*t)'，那么我希望fval 等于3*cos(i*t)。如何让 Fortran 将输入的字符串解释为由 Fortran 执行的命令？

subroutine f(fval, i, t, str)
implicit none
integer, parameter :: ikind = selected_int_kind(8)
integer, parameter :: dbl = selected_real_kind(15,307)

integer(kind = ikind) :: i
real(kind = dbl) :: fval, t
character str*100

if(str .eq. 'poly') then
    fval = t**i
elseif(str .eq. 'slat') then
    fval = exp(-i*t)
elseif(str .eq. 'gaus') then
    fval = exp(-i*t*t)
else
    fval = ???
endif

end subroutine

Answer 1

你不能。不容易。不过，您可以做两件事：

• 首先，您可以使用函数指针。 http://www.macresearch.org/advanced_fortran_90_callbacks_with_the_transfer_function 有一个简单的例子——这个想法是你可以传递在别处定义的函数的名称。这可能足以解决您的问题。

• 您可以调用一个库来解析字符串并计算表达式。我对fortran一无所知，但是您的fortran编译器可能支持调用c例程。对于 gfortran（也许还有其他人），这看起来可以工作 - http://www.gnu.org/s/libmatheval/（http://www.gnu.org/s/libmatheval/manual/libmatheval.html#Fortran-sample-program 的示例 fortran 代码）

Answer 2

其实很简单

subroutine f(fval, i, t, str) ... character(len=*), intent(in) :: str ... end subroutine

诀窍是使用intent(in) 修饰符将虚拟字符串参数定义为“未知长度”。