在互联网上已经花了几天的时间寻找类似的内容,但我仍然找不到任何描述这个问题的东西。即使阅读了《Writing R Extensions》(否则相当值得推荐),也没有提供解决方案。因此,这是我最紧急的问题:
是否可能通过.Fortran(...)调用将函数(为简单起见,假设一个简单的R函数-实际上,问题甚至更加复杂)作为函数/子例程参数传递给Fortran-如果是这样,如何传递?
我编写了两个简单的函数来测试这个问题,首先是Fortran子例程(专门用于使用我最初打算传递的函数,因此界面中有一些奇怪的维度)。
调用
是否可能通过.Fortran(...)调用将函数(为简单起见,假设一个简单的R函数-实际上,问题甚至更加复杂)作为函数/子例程参数传递给Fortran-如果是这样,如何传递?
我编写了两个简单的函数来测试这个问题,首先是Fortran子例程(专门用于使用我最初打算传递的函数,因此界面中有一些奇怪的维度)。
subroutine foo(o, x)
implicit none
interface
subroutine o(t, y, param, f)
double precision, intent(in) :: t
double precision, dimension(58), intent(in) :: y, param
double precision, dimension(22), intent(out) :: f
end subroutine
end interface
double precision, dimension(22), intent(out) :: x
double precision, dimension(58) :: yt, paramt
integer :: i
do i = 1, 58
yt(i) = rand(0)
paramt(i) = rand(1)
end do
call o(dble(4.2), yt, paramt, x)
end subroutine
还需要一个简单的R函数来传递给上面的函数:
asdf <- function(a, s, d, f){x <- c(a, s, d, f)}
调用
.Fortran("foo", asdf, vector(mode="numeric", length=22))
会导致Error: invalid mode (closure) to pass to Fortran (arg 1)
的错误,而将"asdf"(作为字符串)传递会导致段错误,因为该参数显然不符合预期类型(即函数)。
顺便说一下,我并不指望代码能做出任何有意义的事情(这将是另一个函数的任务),我主要想知道是否可能从R中传递函数(或函数指针),或者我最好立即放弃这种方法,寻找可能奏效的其他方法。
提前感谢您的帮助,
迪安