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位下编译 链接也都没什么错误，就是在运行到这里时候报错，我想，应该不是我的数组分配问题吧。

出错的地方就在call dpbtrf('U',NEQUA,MBAND-1,FACT,MBAND,info)，其错误就是上面所说。32位的时候，计算没问题，64为，运行到这里就不行了，过不去了，这个是MKL库函数，看不到错误信息