重温：Fortran 程序中的堆栈溢出

Question

在 Compac Fortran 中执行以下代码时出现堆栈溢出。 特别针对该行：

DIMENSION IZ(NX,NY),VX(NX-1,NY-1),VY(NX-1,NY-1),VZ(NX-1,NY-1)

我正在做的是从参数文件中读取 NX,NY,DX,DY 。

有什么建议么？

       PARAMETER(NGMAX=30000)
       INTEGER NX,NY,DX,DY
       OPEN(6,FILE='MGSTAOriggroup15.asc',STATUS='old')
       OPEN(7,FILE='Gravity.asc',STATUS='old')
       OPEN(8,FILE='Gravity200.nor',STATUS='old')
       OPEN(9,FILE='linana.para',STATUS='old')
       OPEN(10,FILE='MGSTAOriggroup15coord.dat',STATUS='unknown')
       OPEN(12,FILE='MGSTAOriggroup15ncoord.dat',STATUS='unknown')
       OPEN(11,FILE='linana.fsn',STATUS='unknown')
       READ(9,*) NX,NY,DX,DY
       CALL ANALYSIS(NX,NY)
       Close(6)
       Close(7)
       Close(8)
       Close(9)
       Close(10)
       Close(11)
       Close(12)
C
       STOP 
     END
C 
      SUBROUTINE ANALYSIS(NX,NY)
      INTEGER NX,NY,DX,DY
      COMMON/COM1/ DX,DY
      PARAMETER(NGMAX=30000)
      DIMENSION KO(NGMAX)
      DIMENSION XLS(NGMAX,100),XLE(NGMAX,100),
     &          YLS(NGMAX,100),YLE(NGMAX,100)
      DIMENSION IZ(NX,NY),VX(NX-1,NY-1),VY(NX-1,NY-1),VZ(NX-1,NY-1)
      """"""""""""""""""""""""""""""""""""""""""""
      To the end

注意：从 (NGMAX=30000) 更改为 (NGMAX=30) 没有帮助。

是的，我已将 NGMAX 的所有实例更改为 30。

下面列出了新版本。

---

    INTEGER NX,NY,DX,DY
    PARAMETER(NGMAX=30)
      COMMON /CM1/ FX(NGMAX),FY(NGMAX),FZ(NGMAX),FR(NGMAX),IL(NGMAX)
      COMMON /CM2/ JST(NGMAX),KST(NGMAX),JDP(NGMAX),KDP(NGMAX)
      COMMON /CM3/ XC(NGMAX),YC(NGMAX),ZC(NGMAX),KTYP(NGMAX)
      COMMON /CM4/ A(4,4),B(4),U(4)
    COMMON /CM5/ KO(NGMAX)
      COMMON /CM6/ XLS(NGMAX,10),XLE(NGMAX,10),
     &          YLS(NGMAX,10),YLE(NGMAX,10)
      OPEN(6,FILE='MGSTAOriggroup15.asc',STATUS='old')
      OPEN(7,FILE='Gravity.asc',STATUS='old')
      OPEN(8,FILE='Gravity200.nor',STATUS='old')
      OPEN(9,FILE='linana.para',STATUS='old')
      OPEN(10,FILE='MGSTAOriggroup15coord.dat',STATUS='unknown')
      OPEN(12,FILE='MGSTAOriggroup15ncoord.dat',STATUS='unknown')
      OPEN(11,FILE='linana.fsn',STATUS='unknown')
      READ(9,*) NX,NY,DX,DY
    CALL ANALYSIS(NX,NY)
      Close(6)
      Close(7)
      Close(8)
      Close(9)
      Close(10)
      Close(11)
      Close(12)
C
       STOP 
     END
C 
      SUBROUTINE ANALYSIS(NX,NY)
    INTEGER NX,NY,DX,DY
      COMMON/COM1/ DX,DY
    PARAMETER(NGMAX=30)
      COMMON /CM1/ FX(NGMAX),FY(NGMAX),FZ(NGMAX),FR(NGMAX),IL(NGMAX)
      COMMON /CM2/ JST(NGMAX),KST(NGMAX),JDP(NGMAX),KDP(NGMAX)
      COMMON /CM3/ XC(NGMAX),YC(NGMAX),ZC(NGMAX),KTYP(NGMAX)
      COMMON /CM4/ A(4,4),B(4),U(4)
    COMMON /CM5/ KO(NGMAX)
      COMMON /CM6/ XLS(NGMAX,10),XLE(NGMAX,10),
     &          YLS(NGMAX,10),YLE(NGMAX,10)
      DIMENSION IZ(NX,NY),VX(NX-1,NY-1),VY(NX-1,NY-1),VZ(NX-1,NY-1)

我很抱歉： 我花了一整天的时间编写代码，终于可以工作了。 然而，声明 NX 和 NY 是我们不需要的。 我想从参数文件中读取这些参数。 代码如下：

---

      PARAMETER(NX=322,NY=399,NGMAX=30000)
C I need to change NX and NY to read from the parameter file
      CHARACTER infile1*80,infile2*80,infile3*80,outfile1*80
     &,outfile2*80,outfile3*80
      DIMENSION KO(NGMAX)
      DIMENSION XLS(NGMAX,100),XLE(NGMAX,100),
     &          YLS(NGMAX,100),YLE(NGMAX,100)
      DIMENSION IZ(NX,NY),VX(NX-1,NY-1),VY(NX-1,NY-1),VZ(NX-1,NY-1)
      DIMENSION FX(NGMAX),FY(NGMAX),FZ(NGMAX),FR(NGMAX),IL(NGMAX)
      DIMENSION JST(NGMAX),KST(NGMAX),JDP(NGMAX),KDP(NGMAX)
      DIMENSION XC(NGMAX),YC(NGMAX),ZC(NGMAX),KTYP(NGMAX)
      DIMENSION A(4,4),B(4),U(4)
C
    CALL getenv('INFILE1',infile1)
    CALL getenv('INFILE2',infile2)
    CALL getenv('INFILE3',infile3)
    CALL getenv('OUTFILE1',outfile1)
    CALL getenv('OUTFILE2',outfile2)
    CALL getenv('OUTFILE3',outfile3)
    OPEN(1,FILE='Alaa1.para',STATUS='old')
    READ(1,*)DX,DY,infile1,infile2,infile3,outfile1,outfile2
     &,outfile3
C
       OPEN(6,FILE=infile1,STATUS='old')
       OPEN(7,FILE=infile2,STATUS='old')
       OPEN(8,FILE=infile3,STATUS='old')
       OPEN(10,FILE=outfile1,STATUS='unknown')
       OPEN(12,FILE=outfile2,STATUS='unknown')
       OPEN(11,FILE=outfile3,STATUS='unknown')

       TO THE END

---

Answer 1

听起来子例程ANALYSIS正在根据传递给它的参数（NX和NY）在堆栈上分配空间。 如果这些值太大，则堆栈上可能没有足够的空间来分配数组。

您可以做的是打印出从 linana.para 文件中读取的 NX 和 NY 值，看看它们对您是否有意义应用。

Answer 2

呵呵。 好吧，从技术上讲，这可能属于 ServerFault，因为这是操作系统问题，但我认为那里也没有人会回答 VMS 问题。 对于我和 VMS 来说已经有一段时间了，对于 FORTRAN 来说甚至更长，但我会尝试一下...

好吧，这就是您所得到的：ANALYSIS 子例程有 4 个数组，每个数组有 (30,000*100)=3,000,000 个单元格，每个单元格为 12M单元或每个 8 字节（我认为？） 为 72Mbytes。 再加上整数数组 IZ 和实数数组 VX、VY 和 VZ 占用的空间，我无法估计，因为我不知道 NX 和 NY 的最大值是多少。

无论如何，VAX 上至少有 144,000 个页面，Alpha 上至少有 24,000 个页面（加上其他 4 个阵列）。 您的进程地址表至少需要是该值的 2 倍（实际上要多一点），因为用户堆栈只能获得该值的一半。

因此，转到命令提示符并输入“SHOW PROC/MEM”并告诉我它的内容（我想我是对的）。

注意：NX、NY 只是 DIMENSION 的参数，它是占用堆栈的数组。

将数组设为全局可能行不通，因为堆和静态分配空间并不比堆栈大。

好的，我看到您的最大 NX、NY 值为 5,000。 这又是 4*25,000,000 或 100,000,000。 每个 8 字节即 800MB 加上上面的 72MB。 坏消息是，无论 VMS 进程设置如何，VAX 在架构上都限制为 1GB 的用户堆栈，因此我们希望浮点数为 4 字节（而不是 8 字节），或者您使用的是 Alpha。

Answer 3

假设异常发生在您提到的行上，那么您的程序的堆栈对于在其上分配的变量来说不够大。 您需要在编译器中找到该标志来增加程序上的默认堆栈分配。 或者，如果您根本不使用递归，则可以将这些变量设为全局变量，这样它们就不会在堆栈上分配。

请注意您分配的大小：

PARAMETER(NGMAX=30000)
DIMENSION KO(NGMAX)
DIMENSION XLS(NGMAX,100),XLE(NGMAX,100),YLS(NGMAX,100),YLE(NGMAX,100)
DIMENSION IZ(NX,NY),VX(NX-1,NY-1),VY(NX-1,NY-1),VZ(NX-1,NY-1)

您为四个数组 XLS、中的每一个数组分配 30000 * 100 个值（我假设是实数）的空间XLE、YLS 和 YLE。 那是 12,000,000 个值。 假设 32 位值，则内存为 48 MB。 但我依稀记得在某些平台上，实数是80位？ 在这种情况下，内存将是 120 Meg。

您真的需要您的数组在一维长度为 30,000 个元素吗？ 如果没有，减少分配可能会对您有所帮助。

查看您正在使用的 FORTRAN 编译器的编译器选项。 我记得 VAX/VMS FORTRAN 编译器有一个命令行选项可以使所有变量静态化，也就是说不在堆栈上分配。 假设您可以保证不会进行任何递归，这可能是解决您的问题的最快解决方案。

另一个选择是将所有大型数组放入 COMMON 块中，但请注意，公共块中的数组无法动态调整大小。

Answer 4

看起来您的代码是 Fortran77 如果是的话，以下内容应该有意义，如果不是，它可能会给您一些线索。

在 Fortran77 中，语言标准中没有动态内存分配支持，许多编译器确实对此提供了一些支持，但通过使用编译器特定的扩展来实现。 您的代码不会触发类似扩展的任何记忆。

执行您想要的操作的一种方法是生成一个主程序，该程序将声明数组空间，然后将数组作为参数传递给子例程。 该主程序将被编译并与其余代码链接。 每组值都会产生一个新的主程序。

例如 - 对于 NX=2456 & 纽约=789

  PROGRAM MAIN
  INTEGER NX, NY
  DIMENSION IZ(2456,789),VX(2455,788),VY(2455,788),VZ(2455,788)
  NX=2456
  NY=789
  CALL ANALYSIS(NX,NY,IZ,VX,VY,VZ)
  STOP
  END

  SUBROUTINE ANALYSIS(NX,NY,IZ,VX,VY,VZ)
  INTEGER NX,NY,DX,DY
  DIMENSION IZ(NX,NY),VX(NX-1,NY-1),VY(NX-1,NY-1),VZ(NX-1,NY-1)

...
返回
我不记得

了，但可能可以在主程序中使用参数语句，以便您完全参数化上面的主程序。