24小时热门版块排行榜    

查看: 2906  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

todd0226

新虫 (初入文坛)

[求助] fortran 调用子程序时,出现program exception-access violation 已有1人参与

大家好,本人fortran新人,编了一段程序,但是在运行时出现了如题所示的错误!
因为程序太长,而且错综复杂,不方便贴出,敬请见谅!
下面我进行一些相关的描述,请各位指导!
我通过debug调试出问题的出处:错误信息出现在subroutine info_nd(...)这一行
其中:step,nelm, nj,nnode 是在主程序中声明的整数常量
           剩余的变量全是数组,全部是动态变量,而且我全部在主程序中进行了allocate,各个数组的维数以及上下界和子程序中声明的完全一样。
调试时,屏幕输出了“info_nd调试开始”,没有输出“info_nd调试第一点”,而且绿色箭头指向了subroutine info_nd(...),所以我断定问题出处在这一行。
我网上搜了一下,这种问题的可能原因之一是数组越界,但是我感觉我的这种定义数组的形式应该不会造成这个问题(?)
此外,有这样一个现象:nj = 6904,nelm=6804的时候,无论debug还是release都会在subroutine info_nd(...)这一行出现如题目所说的问题;但是当我把nj,nelm调小了之后(<6000),发现debug下可以运行成功,但是release仍然出现上述问题。
所以在此请教大家,恳请予以指导!谢谢!


主程序
program main
integer,parameter:: step=10,nj=6904,nelm=6804,nnode=4
......
......   
do istp = 1,step
      ......
      write(*,*)'info_nd调试开始'
      call info_nd(step,istp,nelm,nj,nnode,maxgint,id_el_nd,id_nd_el,num_nd_el,gint_t,gint,elarea,stgp,svars_gp,sig11_nd,sig22_nd,sig12_nd,mises_nd)
      ......
end do
......
end   !主程序结束

出现问题的子程序
subroutine info_nd(step,istp,nelm,nj,nnode,maxgint,id_el_nd,id_nd_el,num_nd_el,gint_t,gint,elarea,stgp,svars_gp,sig11_nd,sig22_nd,sig12_nd,mises_nd)
implicit none

integer   step,istp,nelm,nj,nnode
integer::   id_el_nd(nelm,nnode),id_nd_el(nj,nelm),num_nd_el(nj),gint_t(nelm+1,step),gint(nelm,step),maxgint(step)
real*8::   svars_gp(nelm,20*24,step),elarea(nelm),stgp(3,24*nelm,step)
real*8::   sig_nd_temp(nnode,3),sig11_nd_temp(nj,nelm),sig22_nd_temp(nj,nelm),sig12_nd_temp(nj,nelm)
real*8::   sig11_nd(nj,step),sig22_nd(nj,step),sig12_nd(nj,step),mises_nd(nj,step),area(nj)
real*8::   NTN(nnode,nnode),invNTN(nnode,nnode)
real*8,allocatable::   sig_gp(:,,N(:,,NT(:,,invN(:,
real*8    sigp,tigp
integer   gint_el,gint_iel
integer   iel,igp,i,j,ii

write(*,*)'info_nd调试第一点'
......
end
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

todd0226

新虫 (初入文坛)

引用回帖:
7楼: Originally posted by pippi6 at 2014-11-04 17:20:43
No idea。已经不是fortran问题了...

好的,多谢
8楼2014-11-04 21:08:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

【答案】应助回帖

感谢参与,应助指数 +1
在你的makefile (或编译 option)里加上 /check:bounds  /traceback
就会在runtime 时间显示在何处越界。
2楼2014-11-03 18:50:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

todd0226

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by pippi6 at 2014-11-03 18:50:05
在你的makefile (或编译 option)里加上 /check:bounds  /traceback
就会在runtime 时间显示在何处越界。

不好意思,对这个我了解的不多。
你说的在Makefile里加上 /check:bounds  /traceback,Makefile是我编写的源程序文件吗?加在哪个位置?
编译option,这是编译器的选项吗?我用的是Intel visual Fortran编译器,没有找到这一项呀?
3楼2014-11-03 21:35:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

todd0226

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by pippi6 at 2014-11-03 18:50:05
在你的makefile (或编译 option)里加上 /check:bounds  /traceback
就会在runtime 时间显示在何处越界。

我在Intel visual Fortran中找到了/check:bound /trackback,然后调试发现没有出现数组越界的问题。
在运行到subroutine info_nd(...)时,
提示信息说:其原因可能是堆被破坏,这说明NASA_mb05_p1.exe中或它               所加载的任何DLL中有Bug。原因也可能是用户在NASA_mb05_p1.exe具有焦点时按下了F12

其中,不是后一个原因,在运行时我并没有碰到键盘,关于第一个原因,能从中得到哪些信息呢?
4楼2014-11-04 10:12:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 极限拉扯,如何抉择 +11 Isak0320 2026-10-12 15/750 2026-10-14 21:55 by 猫鱼不是鱼
[基金申请] 申请科研项目纯靠实力能成功吗? 200+5 David1982 2026-10-13 8/400 2026-10-14 21:38 by lhw_GPA4.0
[硕博家园] 中科院海洋所招收海洋材料相关领域博士研究生及研究助理 +3 sujiner 2026-10-08 4/200 2026-10-14 20:40 by qicaiyun
[论文投稿] 二审4个多月了还是Awaiting Reviewer Scores 10+3 纳斯达克ML 2026-10-11 4/200 2026-10-14 17:30 by babero
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 13/650 2026-10-14 15:12 by zuesli
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见