【求助】在子程序定义数组并赋值出错

苍山负雪

主程序定义了 PARAMETER NZ=100
经过 PROGRAM MAIN 传到 SUB1, 再传到SUB2
  
在SUB2中定义数组 A(4*NZ,4*NZ)
DIMENSION A(4*NZ,4*NZ)
DOUBLE PRECISION A
并赋值为0
DO I=1, 4*NZ
      DO J=1,4*NZ
            A(I,J)=0.0
      END DO
END DO
发现赋值过程中A的某一个元素就会突然溢出数值为XX.XXE+203
  
这怎么回事？有没有别的办法来定义并赋值

edodo

苍山负雪 wrote:
主程序定义了 PARAMETER NZ=100
经过 PROGRAM MAIN 传到 SUB1, 再传到SUB2
  
在SUB2中定义数组 A(4*NZ,4*NZ)
DIMENSION A(4*NZ,4*NZ)
DOUBLE PRECISION A
并赋值为0
DO I=1, 4*NZ
      DO J=1,4*NZ
            A(I,J)=0.0
      END DO
END DO
发现赋值过程中A的某一个元素就会突然溢出数值为XX.XXE+203
  
这怎么回事？有没有别的办法来定义并赋值

  
老哥试试这样
将你原来的
在SUB2中定义数组 A(4*NZ,4*NZ)
改为
在SUB2中定义数组 A(400,400)
我觉得可能是4*nz误差造成的

chenyq

与这位兄台见解相同！不过我觉得，如果前面NZ没有定义为整型（Common)，NZ应该不再是整数了。你不妨在循环之前将NZ输出来，看看是否为400。

苍山负雪

输出NZ=100 没有错误
按照edodo的做法也不行
其实在SUB2中还定义了数组B(4*NZ),给B赋值就一点问题都没有 奇怪了:~)
另外必须定义COMMON吗？
NZ作为参数，经常要改变的，直接把数组大小定为400也不合适吧

edodo

哈哈！
那就是这样的问题了
因为开的数组太大400*400=160000
是不是会有问题！：）
老哥
你可以试着改改
拆分；）
有必要非要做成400*400么？
要不来个400*200和400*201
哈哈！看看吧！

苍山负雪

晕倒，我要疯了
试着把NZ改成20，也就是80*80的矩阵，居然也通不过？
见鬼的什么意思萨。

chenyq

见兄台如此见解，本人在机子上尝试了一下，按照给出的程序编译了一番，发现没有什么问题，只是开头与你的有一点不同：PARAMETER (NZ=400)，因为在编译时不加括号会出错。现在我也想不通了！！！

edodo

哈哈！
有鬼了！
大家快来抓鬼呀！
苍山负雪  兄的机器里面有鬼了！

苍山负雪

俺超级faint，费了半天劲，突然发现有一步子程序里把A定义成一维的了。手误手误~~
  
编译的时候也没有提示出来，白白耗费了俺一天时间

edodo

哈哈哈！！
  
老哥那我们找乐！
哈哈！
开个玩笑

苍山负雪

:I 还是谢谢 edodo兄和chenyq兄

edodo

客气

niuton2044

有VC中应该如何定义呢？

jing-qian

1# 苍山负雪
苍山负雪:
I have a same problem as you had. Do you mind taking a look?

I have a very simple program like the following. When I run this code under Visual Studio 2008 with Intel Fortran 11, it reports error of stack overflow.  If I define array m as "integer m(512, 512)" but not "integer m(k(1),k(2))", the code works well. My system is Windows XP 32, and RAM is 1GB. It should not be the problem of memory. In addition, I had tried to change k(1) and k(2) to 500 but not 512 in main program. And the code works well too. I am so confused!


program test

IMPLICIT REAL(A-H, O-Z)

integer k(2)
common /index/ k

k(1) = 512
k(2) = 512

call dary

end program test

subroutine dary

IMPLICIT REAL(A-H, O-Z)
integer k(2)
common /index/ k

integer m(k(1),k(2))

m(1,1) = 2
write(*,*) m(1,1)

end subroutine dary

jing-qian

I already fixed the problem. Thanks!