沙发
板凳
SUBROUTINE DEM_KBB_gaojie(NEQUA,KFIXA,AKb,MBAND
. ,FACT,AKib,AKbi,AKbb)
C***********************************************
C
USE CONSTANT
IMPLICIT DOUBLE PRECISION(A-H,O-Z)
C
DIMENSION AKb(KFIXA,KFIXA)
DIMENSION AKbb(KFIXA,KFIXA),AKib(NEQUA,KFIXA),AKbi(KFIXA,NEQUA)
. ,FACT(MBAND,NEQUA)
C
call dpbtrf('U',NEQUA,MBAND-1,FACT,MBAND,info)
call dpbtrs('U',NEQUA,MBAND-1,KFIXA,FACT,MBAND,AKib,NEQUA,info)
call MULTI(AKbi,AKib,KFIXA,NEQUA,KFIXA,AKb)
C
AKb=AKbb-AKb
C
RETURN
END
上面就是所调用的子程序,在32位系统下,运行没有问题的,在64位下编译 链接也都没什么错误,就是在运行到这里时候报错,我想,应该不是我的数组分配问题吧。
3 楼
4 楼
出错的地方就在call dpbtrf('U',NEQUA,MBAND-1,FACT,MBAND,info),其错误就是上面所说。32位的时候,计算没问题,64为,运行到这里就不行了,过不去了,这个是MKL库函数,看不到错误信息
您尚未登录,请登录后再回复。点此登录或注册