fortran***在“%a.out'；中出错：free（）：下一个大小无效（正常）：0x0e18b000***

我对fortran是新手。这是我的代码：

  subroutine fft_forward(Q,ier,fftable,ffwork,nfftable,nffwork,     &
 &                       nfft1,nfft2,nfft3,nfftdim1,nfftdim2,       &
 &                       nfftdim3)
  implicit none 

  integer nfft1,nfft2,nfft3,nfftdim1,nfftdim2,nfftdim3
  integer nfftable,nffwork,ier
  double precision,dimension(nfftdim1,nfftdim2,nfftdim3):: Q
  double precision,dimension(nfftable):: fftable
  double precision,dimension(nffwork):: ffwork 

  call pubz3d(Q,ier,nfft1,nfft2,nfft3,nfftdim1,nfftdim2,nfftdim3,   &
 &            fftable,nfftable,ffwork,nffwork)
  write(6,*) 'bye'
  return
  end

我在pubz3d子例程的不同行中放置了一些write语句，它们工作正常（它在子例程的最后一行打印“hi”）。但子例程不返回任何值，甚至打印“bye”

  subroutine pubz3d(Q,ier,nfft1,nfft2,nfft3,nfftdim1,nfftdim2,      &
 &                  nfftdim3,wsave,lensav,work,lenwrk)
  implicit none
  double precision,dimension(nfftdim1,nfftdim2,nfftdim3)::Q
  double precision,dimension(lensav)::wsave      
  double precision,dimension(nfft1)::r
  double precision,dimension(lenwrk)::work
  integer nfft1,nfft2,nfft3,nfftdim1,nfftdim2,nfftdim3,n
  integer lensav,lenwrk,ier
  integer i,j,k

  n = nfftdim1 * nfftdim2 *nfftdim3

  do  k = 1, nfft3
   do  j = 1, nfft2
    do  i = 1,nfft1
      r(i) = Q(i,j,k)
    continue
    call dfft1f(n,1,r,n + 1,wsave,lensav,work,lenwrk,ier)
    do  i = 1,nfft1
     Q(i,j,k) = r(i)
    continue
   continue
  continue

  do  k = 1,nfft3
   do  i = 1,nfft1
    do  j = 1,nfft2
      r(j) = Q(i,j,k)
    continue
    call dfft1f(n,1,r,n + 1,wsave,lensav,work,lenwrk,ier)
    do  j = 1,nfft2
      Q(i,j,k) = r(j)
    continue
   continue
  continue

  do  i = 1, nfft1
   do  j = 1, nfft2
    do  k = 1,nfft3
      r(k) = Q(i,j,k)
    continue
    call dfft1f(n,1,r,n + 1,wsave,lensav,work,lenwrk,ier)
    do  k = 1,nfft3
      Q(i,j,k) = r(k)
    continue
   continue
  continue
  write(6,*) 'hi'
  return
  end

跑步是我得到的

hi
*** Error in `./a.out': free(): invalid next size (normal): 0x0e5a9000 ***
Program received signal SIGABRT: Process abort signal.
Backtrace for this error:
#0  0xB760B163
#1  0xB760B800
#2  0xB77103FF
#3  0xB7710424
#4  0xB7416606
#5  0xB7419A32
#6  0xB7450E52
#7  0xB745B339
#8  0xB745BFAC
#9  0x804BD3B in pubz3d_
#10  0x804BFE4 in fft_forward_
#11  0x8055FF4 in do_pmesh_kspace_
#12  0x8054F2D in MAIN__ at main.f:?
Aborted

---

谢谢

您确定

nfft1>nfft2

和

nfft1>nfft3

？辅助变量

r

的大小为

nfft1

，但您可以在循环中使用它，直到

nfft2

和

nfft3

。如果其中任何一个大于nfft1，则您将得到一个错误。确保

r

的尺寸足够大，可以容纳所有三种循环

---

如果这不是问题所在，请尝试使用一些数组边界检查进行编译（

gfortran

：

-fcheck=all

，

ifort

：

-检查全部以确保安全）。这可能会在出现错误时引发信息更丰富的错误。
您确定
nfft1>nfft2
和
nfft1>nfft3
？辅助变量
r
的大小为
nfft1
，但您可以在循环中使用它，直到
nfft2
和
nfft3
。如果其中任何一个大于nfft1

，则您将得到一个错误。确保

r

的尺寸足够大，可以容纳所有三种循环

---

如果这不是问题所在，请尝试使用一些数组边界检查进行编译（

gfortran

：

-fcheck=all

，

ifort

：

-检查全部以确保安全）。这可能会在出现问题时抛出一个信息量更大的错误。
谢谢你，安德拉斯，这很有效。这是我第一次体验。@Motaharhkia没问题。如果我解决了您的问题，请考虑将我的答案标记为已接受。谢谢您，安德拉斯。这是我第一次体验。@Motaharhkia没问题。如果我解决了你的问题，请考虑将我的答案标记为已接受。