从R调用Fortran：“大小合理的逻辑数组”

Question

这是我在R和Fortran中遇到的一个有趣的问题。

考虑以下Fortran子例程：

! foo.f90
subroutine fun(n, x)
  implicit none
  integer, intent(in) :: n
  logical, intent(in) :: x(n)
  integer :: i

  write(*,*) pack([(i, i=1, n)], x)

end subroutine fun

用R CMD SHLIB foo.f90R 编译并在R中调用时

dyn.load("foo.so")
x <- c(TRUE ,FALSE)
f <- .Fortran("fun", 2L, as.logical(x))

我收到错误消息：

Fortran runtime error: Funny sized logical array

从Fortran程序调用相同的子例程时，

! test.f90
program test
  implicit none
  logical :: x(2)

  x = [.true., .false.]
  call fun(2, x)

contains

  include "foo.f90"

end program test

没有错误发生：

$ gfortran foo.f90 test.f90 -o test
$ ./test 
           1

潜在的罪魁祸首：我在此示例中使用的Fortran编译器。我使用了GNU Fortran 9.2.1。我注意到gfortran 8.3.0存在相同的问题，但是早期版本7.4.0、6.5.0和5.5.0没有错误。但是，为什么test.f90具有相同子例程的Fortran程序在这里可以工作？（即使当相同的标志使用R CMD SHLIB用于编译foo.f90和test.f90：-fno-optimize-sibling-calls -fpic -ggdb -O0 -Wall -Wextra -pedantic -fimplicit-none -fcheck=all -fbacktrace）。

刚开始时看起来很有趣，但是经过一段时间解决这个问题后，这变得令人沮丧...

在此先感谢您的帮助。

更多细节：

> sessionInfo()
R version 3.6.1 (2019-07-05)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu 18.04.3 LTS

Matrix products: default
BLAS:   /usr/lib/x86_64-linux-gnu/blas/libblas.so.3.7.1
LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.7.1

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
 [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8    
 [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8   
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C            
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

loaded via a namespace (and not attached):
[1] compiler_3.6.1