21*21的矩阵可以输出好的结果，而40*40的矩阵输出NaN

主题：2121的矩阵可以输出好的结果，而4040的矩阵输出NaN

feiyu861010 [专家分：0] 发布于 2011-05-02 14:51:00

我是做偏微分方程数值解的，以前用matlab编程序，这次改用fortran编，可是低分辨率下可以得到理想的结果，高分辨率就不行了，结果会出现NaN，我看了以前的帖子，有这方面的，不过没有解决，请教各位老师帮我想想办法啊

沙发

yeg001 [专家分：14390] 发布于 2011-05-02 19:07:00

具体要看你代码怎么写了. 如果程序一开始设计成有弹性的话应该不会出现这种情况.

板凳

feiyu861010 [专家分：0] 发布于 2011-05-20 20:43:00

问题已解决，是粗心写错了一个数组

3 楼

feiyu861010 [专家分：0] 发布于 2011-05-20 20:46:00

您好，我有个问题想请教一下，就是我想读取一个dat文件，数据是有3行2列，想读取时存到数组中，可是只能读取第一列，读取第二列就错，
open(3,file = 'bbc.dat')
do j=1,2
  do i=1,3
     read(3,*)b(i,j)
  end do
end do

4 楼

yeg001 [专家分：14390] 发布于 2011-05-20 23:08:00

每次读入缓冲是一行的. 每一次read完之后默认会抛掉这一行的缓冲, 下一个read就去读下一行了.
你还是一行就读入两个数据吧.
do i=1,3
read(3,*)b(i,:)
end do

您尚未登录，请登录后再回复。点此登录或注册