Fortran Coder

查看: 46|回复: 8

[求助] 萌新求救

[复制链接]

8

帖子

2

主题

0

精华

入门

F 币
112 元
贡献
29 点

规矩勋章

发表于 6 天前 | 显示全部楼层 |阅读模式
大家好!我在运行一个程序的时候出现了错误: Fortran runtime error: Bad repeat count in item 1 of list input   (对应的程序是:read(curr%defaultval,*) curr%ival)
我上网找只有一个相关的内容而且那个情况是:后面的变量有多个,但是被read的文件内的值只有一个,所以报错。

而我的curr%ival 只需要一个变量。所以应该不是这个原因。
另外
  这个程序本来是在ia32上用ifort编译运行的,非常成熟。
但是当时在make的时候我只有gfortran ,所以我修改了Makefile 用LIBS = -llapack -lblas 来链接静态库,
一直无法通过,直到我把LIBS改成MKLLIB才make通过。我想运行的也许错误和这个有关?
  我的操作系统是Ubuntu18 64位,这个错误会不会与这个有关呢?



回复

使用道具 举报

282

帖子

1

主题

0

精华

宗师

F 币
1563 元
贡献
1048 点
发表于 6 天前 | 显示全部楼层
最好截个图。既然已经运行了,基本上和编译链接没有关系,不要胡猜。

44

帖子

6

主题

0

精华

专家

超凡脱俗

F 币
317 元
贡献
168 点
发表于 6 天前 | 显示全部楼层
本帖最后由 Jackdaw 于 2018-12-4 16:13 编辑

文件内容和变量类型不对应,检查输入文件吧


Screenshot_2018-12-04_16-10-20.png
天下英雄出我辈,一入江湖岁月催。

鸿图霸业谈笑间,不胜人生一场醉。

8

帖子

2

主题

0

精华

入门

F 币
112 元
贡献
29 点

规矩勋章

 楼主| 发表于 5 天前 | 显示全部楼层
好的谢谢,我晚上看看

8

帖子

2

主题

0

精华

入门

F 币
112 元
贡献
29 点

规矩勋章

 楼主| 发表于 5 天前 | 显示全部楼层
本帖最后由 BNUPHYKLIANG 于 2018-12-5 16:39 编辑
Jackdaw 发表于 2018-12-4 16:11
文件内容和变量类型不对应,检查输入文件吧

感谢!
[Fortran] 纯文本查看 复制代码
do i = 1, N_Param
       curr => cfg%record(i)
       curr%id         = i
       curr%pname      = PARAM_NAME(i)
       curr%ptype      = PARAM_TYPE(i)
       curr%isArray    = PARAM_ARRAY(i)
       curr%defaultval = PARAM_DVAL(i)

       if (curr%ptype .eq. TYPE_REAL) then
          read(curr%defaultval,*) curr%rval
       else if (curr%ptype .eq. TYPE_INTEGER) then
          read(curr%defaultval,*) curr%ival
       else if (curr%ptype .eq. TYPE_COMPLEX) then
          read(curr%defaultval,*) curr%cval     
       end if

这就是出错处的那个程序,提示倒数第二行read(curr%defaultval,*) curr%cval 出错,错误就
Fortran runtime error: Bad repeat count in item 1 of list input

按理说,它提前判断了数据类型,不应该再出错。
curr定义如下:

[Fortran] 纯文本查看 复制代码
 type(Param),pointer    :: curr
    integer :: i


curr=>cfg%record(i)
  对于报错的 read(curr%defaultval,*) curr%cval
  因为curr%defaultval = PARAM_DVAL(i)
中我在想是不是i的值太大导致 PARAM_DVAL(i) 中的数据已经没有了(但是它判断数据是个复数啊), PARAM_DVAL如下:
[Fortran] 纯文本查看 复制代码
character(len=*), parameter :: PARAM_DVAL(N_Param) =  &
       &(/"-1      ", "HSF.in  ", "HSF.out ", "12      ", "0.0     ", &
       &  "0       ",  &
       &  "0.001   ", "0.0     ", "8       ", "0.125   ", "0.001   ", &
       &  "quest   ", "geom.def", "8       ", "0.0     ", "8       ", &
       &  "16      ", "10      ", "50      ", "4       ", "12      ", &
       &  "12      ", "5000    ", "30      ", "0       ", "1000    ", &
       &  "12      ", "4       ", "4       ", "2       ", "0.0001  ", &
       &  "0       ", "1.0     ", "10      "/)




8

帖子

2

主题

0

精华

入门

F 币
112 元
贡献
29 点

规矩勋章

 楼主| 发表于 5 天前 | 显示全部楼层
本帖最后由 BNUPHYKLIANG 于 2018-12-5 17:42 编辑
li913 发表于 2018-12-4 16:08
最好截个图。既然已经运行了,基本上和编译链接没有关系,不要胡猜。

好的,谢谢您!
那和编译器也没关系吗?我看网上有个外国论坛说了什么用ifort 和gFortran 编译运行的结果不一样,并且有一次  类似read(....)x,y,z 出错,ifort read了3个数据,gfortran  read了1个数据截图如下
图中错误代码是:
[Fortran] 纯文本查看 复制代码
do i = 1, N_Param
       curr => cfg%record(i)
       curr%id         = i
       curr%pname      = PARAM_NAME(i)
       curr%ptype      = PARAM_TYPE(i)
       curr%isArray    = PARAM_ARRAY(i)
       curr%defaultval = PARAM_DVAL(i)

       if (curr%ptype .eq. TYPE_REAL) then
          read(curr%defaultval,*) curr%rval
       else if (curr%ptype .eq. TYPE_INTEGER) then
          read(curr%defaultval,*) curr%ival
       else if (curr%ptype .eq. TYPE_COMPLEX) then
          read(curr%defaultval,*) curr%cval
       end if

就是倒数第二行出错,根据curr%defaultval = PARAM_DVAL(i) 我找到了 PARAM_DVAL如下:
[Fortran] 纯文本查看 复制代码
character(len=*), parameter :: PARAM_DVAL(N_Param) =  &
       &(/"-1      ", "HSF.in  ", "HSF.out ", "12      ", "0.0     ", &
       &  "0       ",  &
       &  "0.001   ", "0.0     ", "8       ", "0.125   ", "0.001   ", &
       &  "quest   ", "geom.def", "8       ", "0.0     ", "8       ", &
       &  "16      ", "10      ", "50      ", "4       ", "12      ", &
       &  "12      ", "5000    ", "30      ", "0       ", "1000    ", &
       &  "12      ", "4       ", "4       ", "2       ", "0.0001  ", &
       &  "0       ", "1.0     ", "10      "/)

本来我先想办法看看是第几次循环出错
在循环中加入了print *,i :就在do i=1,N_Param 后面。
[Fortran] 纯文本查看 复制代码
 do i = 1, N_Param
        print *, i
       curr => cfg%record(i)
       curr%id         = i
       curr%pname      = PARAM_NAME(i)
       curr%ptype      = PARAM_TYPE(i)
       curr%isArray    = PARAM_ARRAY(i)
       curr%defaultval = PARAM_DVAL(i)

       if (curr%ptype .eq. TYPE_REAL) then
          read(curr%defaultval,*) curr%rval
       else if (curr%ptype .eq. TYPE_INTEGER) then
          read(curr%defaultval,*) curr%ival
       else if (curr%ptype .eq. TYPE_COMPLEX) then
          read(curr%defaultval,*) curr%cval
       end if



再make并运行,出现了变化,从1显示到第33,也就是第33次循环除了问题,我看第33个数明明是1.0 不应该判断它是复数啊,而在下面出错时却是
else if (curr%ptype .eq. TYPE_COMPLEX) then

          read(curr%defaultval,*) curr%cval


所以是不是在程序运行过程中数据被改变了


捕获.PNG

613

帖子

0

主题

0

精华

大师

F 币
592 元
贡献
367 点

规矩勋章元老勋章新人勋章水王勋章

QQ
发表于 5 天前 | 显示全部楼层
你要看 PARAM_TYPE 的值。
你的代码里我没有看到自动判断参数类型的代码,可能是写死的。(而不是动态判断的)

8

帖子

2

主题

0

精华

入门

F 币
112 元
贡献
29 点

规矩勋章

 楼主| 发表于 3 天前 | 显示全部楼层
好的谢谢,最后我发现问题出在编译器上,gfortran的编译规则应该更严格,而源程序用的时ifort今天我安装了ifort就解决了,代码中的原始数据是字符串,只不过作者在此处用的是复数,所以用了read。问题在于他直接read(字符串,*)复数   这个在ifort里可以,在gfortran里不行(可能因为gfortran认为复数应该有两个部分)

613

帖子

0

主题

0

精华

大师

F 币
592 元
贡献
367 点

规矩勋章元老勋章新人勋章水王勋章

QQ
发表于 前天 11:02 | 显示全部楼层
同一个问题,往往有多种解答。

楼主把它归结于“编译器差异”,倒不是不可以。

不过,我认为问题出在 PARAM_TYPE(i) 的值。 不应该是 TYPE_COMPLEX
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

QQ|捐赠本站|Archiver|关于我们 About Us|群聊|Fcode

GMT+8, 2018-12-10 17:26

Powered by Discuz! X3.2

© 2001-2017 Comsenz Inc.

快速回复 返回顶部 返回列表