如何解决***捕获了段错误***地址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:
该代码应该:
- 遍历NetCDF文件的向量,并将文件名
spec_hum_file
传递给功能Specific_humidity_afterburner
。 - 该函数读取NetCDF文件,提取数据并传递到第一个编译的子例程,进行数学运算并返回值。
- 获取结果,将其传递给另一个FORTRAN子例程,然后返回第二个结果。
- 将第二个结果写入新的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 举报,一经查实,本站将立刻删除。