Fortran Coder

标题: Bad repeat count in item 1 of list input [打印本页]
作者: BNUPHYKLIANG    时间: 2018-12-4 13:49
标题: Bad repeat count in item 1 of list input
大家好！我在运行一个程序的时候出现了错误： 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位，这个错误会不会与这个有关呢？



作者: li913    时间: 2018-12-4 16:08
最好截个图。既然已经运行了，基本上和编译链接没有关系，不要胡猜。
作者: Jackdaw    时间: 2018-12-4 16:11
本帖最后由 Jackdaw 于 2018-12-4 16:13 编辑

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


Screenshot_2018-12-04_16-10-20.png (146.33 KB, 下载次数: 408)

Screenshot_2018-12-04_16-10-20.png
作者: BNUPHYKLIANG    时间: 2018-12-5 15:53
好的谢谢，我晚上看看
作者: BNUPHYKLIANG    时间: 2018-12-5 16:10
本帖最后由 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      "/)




作者: BNUPHYKLIANG    时间: 2018-12-5 16:19
本帖最后由 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 (55.17 KB, 下载次数: 386)

捕获.PNG
作者: vvt    时间: 2018-12-5 21:44
你要看 PARAM_TYPE 的值。
你的代码里我没有看到自动判断参数类型的代码，可能是写死的。（而不是动态判断的）
作者: BNUPHYKLIANG    时间: 2018-12-7 18:59
好的谢谢，最后我发现问题出在编译器上，gfortran的编译规则应该更严格，而源程序用的时ifort今天我安装了ifort就解决了，代码中的原始数据是字符串，只不过作者在此处用的是复数，所以用了read。问题在于他直接read（字符串，*）复数   这个在ifort里可以，在gfortran里不行（可能因为gfortran认为复数应该有两个部分）
作者: vvt    时间: 2018-12-8 11:02
同一个问题，往往有多种解答。

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

不过，我认为问题出在 PARAM_TYPE(i) 的值。 不应该是 TYPE_COMPLEX
作者: BNUPHYKLIANG    时间: 2018-12-13 14:41
嗯，主要是源程序作者这样设计的，我想的是有可能经过迭代以后那个值是复数，只是输入的初始值是实数
作者: vvt    时间: 2018-12-14 08:59
变量的类型与其定义有关，计算过程不会改变。
（实数是特殊的复数）



欢迎光临 Fortran Coder (http://bbs.fcode.cn/) Powered by Discuz! X3.2