***捕获了段错误***地址0x2aaeb4b6f440，当R调用Fortran时导致“内存未映射”

问题描述

我有一个非常具体的错误，因此谷歌搜索没有帮助，很抱歉，我不知道如何为该问题提供一个简单的可生产示例。下面的代码可以在我的本地计算机上完美运行，但是在HPC上会产生此错误：

 *** caught segfault ***
address 0x2ad718ba0440,cause 'memory not mapped'

Traceback:
 1: array(.Fortran("hus_vertical_interpolation",m = as.integer(DIM[1]),n = as.integer(DIM[2]),o = as.integer(DIM[3]),p = as.integer(DIM[4]),req = as.integer(length(req_press_levels)),hus_on_model_level = as.numeric(spec_hum_data[]),pres = as.numeric(req_press_levels),pressure_full_level = as.numeric(pressure[]),hus_on_press_level = as.numeric(output_array[]))$hus_on_press_level,dim = output_DIM)
 2: Specific_humidity_afterburner(spec_hum_file = q_nc.files[x],req_press_levels = required_PLev)

Possible actions:
1: abort (with core dump,if enabled)
2: normal R exit
3: exit R without saving workspace
4: exit R saving workspace
Selection: 

该代码应该：

1. 遍历NetCDF文件的向量，并将文件名spec_hum_file传递给功能Specific_humidity_afterburner。
2. 该函数读取NetCDF文件，提取数据并传递到第一个编译的子例程，进行数学运算并返回值。
3. 获取结果，将其传递给另一个FORTRAN子例程，然后返回第二个结果。
4. 将第二个结果写入新的NetCDF文件。

在步骤3中发生错误。R函数是：

Specific_humidity_afterburner<-function(spec_hum_file,req_press_levels){
  
  require(ff)
  require(ncdf4)
  require(stringi)
  require(DescTools)
  library(stringr)
  library(magrittr)
  #1============================================================================ 
  #Reading data from netCDF file
  #2============================================================================ 
  #Reading other variables
  #3============================================================================
  # First Fortran subroutine 
  #4============================================================================
  
  #load vertical interpolate subroutine for specific humidity
  dyn.load("spec_hum_afterburner/vintp2p_afterburner_hus.so")
  #check
  is.loaded("hus_vertical_interpolation")
  
  DIM<-dim(spec_hum_data)
  
  output_DIM<-c(DIM[1],DIM[2],length(req_press_levels),DIM[4])
  
  output_array<-ff(array(0.00,dim =output_DIM),dim =output_DIM)
  
  result<- array(.Fortran("hus_vertical_interpolation",m=as.integer(DIM[1]),n=as.integer(DIM[2]),o=as.integer(DIM[3]),p=as.integer(DIM[4]),pres=as.numeric(req_press_levels),pressure_full_level=as.numeric(pressure[]),hus_on_model_level=as.numeric(spec_hum_data[]),hus_on_press_level=as.numeric(output_array[]))$hus_on_press_level,dim =output_DIM)
  
  DIMNAMES<-dimnames(spec_hum_data)
  
  DIMNAMES[["lev"]]<-req_press_levels
  
  Specific_humidity<- ff(result,dim = output_DIM,dimnames =DIMNAMES )
  rm(result)
  
  #5============================================================================
  # Writing NetCDF file of the interpolated values 
}

Fortran子例程：

subroutine hus_vertical_interpolation(m,n,o,p,req,pres,&
    pressure_full_level,hus_on_model_level,hus_on_press_level)
implicit none
integer :: m,req
integer :: x,y,s,t,plev
double precision :: pres(req),hus_on_model_level(m,p)
double precision :: pressure_full_level(m,p)
double precision :: delta_hus,delta_p,grad_hus_p,diff_p
double precision,intent(out) :: hus_on_press_level(m,p)
real :: arg = -1.0,NaN
NaN= sqrt(arg)
do plev=1,req
 do t=1,p 
  do x=1,m
   do y=1,n 
    do s=1,o    
!above uppest level
    if(pres(plev) .LT. pressure_full_level(x,1,t)) then 
     hus_on_press_level(x,plev,t) = NaN
    end if 
! in between levels
    if(pres(plev) .GE. pressure_full_level(x,t) .AND. pres(plev) .LE. &
     pressure_full_level(x,s+1,t) ) then 
     delta_hus = hus_on_model_level(x,t) - hus_on_model_level(x,t)
     delta_p = log(pressure_full_level(x,t))&
     - log(pressure_full_level(x,t))
     grad_hus_p = delta_hus /delta_p
     diff_p = log(pres(plev)) - log(pressure_full_level(x,t))
     hus_on_press_level(x,t) = hus_on_model_level(x,t)&
     + grad_hus_p * diff_p 
    end if  
! extrapolation below the ground  
    if(pres(plev) .GT. pressure_full_level(x,t)) then
     hus_on_press_level(x,t)
    end if
    end do 
   end do 
  end do    
 end do
end do   
end subroutine hus_vertical_interpolation

Fortran子例程使用以下命令编译：

gfortran -fPIC -shared -ffree-form vintp2p_afterburner_hus.f90 -o vintp2p_afterburner_hus.so

例如，错误行为是不可预测的，可能发生在循环的索引1、2、8等。我们试图将大数组作为最后一个变量移交给Fortran子例程，它最大程度地减少了错误的发生。 此外，NetCDF文件的大小约为2GB。值得一提的是，这些模块是使用EasyBuild构建的，因此，按照HPC支持团队的说法，不可能发生冲突。据我们所知，我们尝试了许多解决方案，但没有任何进展！

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）