Fortran Coder

查看: 127|回复: 12

[子程序] IVF中实参虚参类型不一致

[复制链接]

16

帖子

4

主题

0

精华

入门

F 币
90 元
贡献
48 点
发表于 2018-9-25 20:06:10 | 显示全部楼层 |阅读模式
使用的一个软件包,是linux下写的,使用了不同数据类型的实参虚参的传递,但是IVF不支持,应该怎么办?跟大家请教一下,谢谢暂时通过取消接口检查解决的问题,但是结果跟Linux下运行出来的结果不太一样。

这个是软件说明的文档。

1.jpg

回复

使用道具 举报

269

帖子

1

主题

0

精华

宗师

F 币
1504 元
贡献
1014 点
发表于 2018-9-26 13:05:30 | 显示全部楼层
结果不一致,很难说哪里的问题。你最好把linux下的makefile拿出来看看,采用相同的编译参数。

16

帖子

4

主题

0

精华

入门

F 币
90 元
贡献
48 点
 楼主| 发表于 2018-9-26 15:51:43 | 显示全部楼层
li913 发表于 2018-9-26 13:05
结果不一致,很难说哪里的问题。你最好把linux下的makefile拿出来看看,采用相同的编译参数。 ...

好的,我再看看,主要是不太清楚怎么允许这种数据传递,同时不出错

93

帖子

3

主题

0

精华

版主

F 币
456 元
贡献
279 点

贡献勋章管理勋章帅哥勋章元老勋章星光勋章规矩勋章

发表于 2018-9-26 16:54:36 | 显示全部楼层
[Fortran] 纯文本查看 复制代码
PROGRAM a
  REAL :: x
  CALL b(x)
CONTAINS
  SUBROUTINE b(i)
      INTEGER ::i
  END SUBROUTINE b
END PROGRAM a

a 與 b 之間為 explicit interface, (type inconsistent) -> (compile-time fatal error, all compilers)

[Fortran] 纯文本查看 复制代码
PROGRAM aa
  REAL ::xx
  CALL bb(xx)
END PROGRAM aa
SUBROUTINE bb(ii)
  INTEGER :: ii
END SUBROUTINE bb

aa 與 bb 之間為 implicit interface, (type inconsistent) -> (compile-time fatal error or warn 是 compiler or compile-time-option dependent)

aa 與 bb 放在兩個不同的file (aa.f90, bb.f90), 分開compile (王不見王, 沒有type inconsistent的問題), 就OK

==========================
IVF结果跟Linux下运行出来的结果不太一样
==========================
怎的(不太一样)?

不一樣的程度遠高於precision的容忍度, 你 call FFTPack 的 "codes" 出錯的機會極大
uninitialised variables, array out-of-bounds, unsaved local variables, ...

16

帖子

4

主题

0

精华

入门

F 币
90 元
贡献
48 点
 楼主| 发表于 2018-9-26 17:06:11 | 显示全部楼层
chiangtp 发表于 2018-9-26 16:54
[mw_shl_code=fortran,true]PROGRAM a
  REAL :: x
  CALL b(x)

就是计算的结果有时候不太一样,但有时候又是一样的。
谢谢版主!

93

帖子

3

主题

0

精华

版主

F 币
456 元
贡献
279 点

贡献勋章管理勋章帅哥勋章元老勋章星光勋章规矩勋章

发表于 2018-9-26 19:47:56 | 显示全部楼层
怎的(不太一样)?

1.23456789, 1.23456987,  for 4-byte REAL, 這算正常

否則, 你 call FFTPack 的 "codes" 出錯的機會極大

16

帖子

4

主题

0

精华

入门

F 币
90 元
贡献
48 点
 楼主| 发表于 2018-9-26 19:55:09 | 显示全部楼层
chiangtp 发表于 2018-9-26 19:47
怎的(不太一样)?

1.23456789, 1.23456987,  for 4-byte REAL, 這算正常

我只是单纯看的结果,那个软件是进行矩阵分解C=BP的,分解出来的矩阵B是C的某些列,IVF和Linux的对应的列数是不一致的。
虽然结果不太一样,但感觉应该都是对的

93

帖子

3

主题

0

精华

版主

F 币
456 元
贡献
279 点

贡献勋章管理勋章帅哥勋章元老勋章星光勋章规矩勋章

发表于 2018-9-26 20:04:27 | 显示全部楼层
non-zero (A, B) 的不一樣

相對誤差 abs(A-B)/abs(A or B)
< 1.0E-6 (for 4-Byte), or,
< 1.0E-15 (for 8-byte)

的 REAL computations, 都算正常

16

帖子

4

主题

0

精华

入门

F 币
90 元
贡献
48 点
 楼主| 发表于 2018-9-26 20:18:43 | 显示全部楼层
chiangtp 发表于 2018-9-26 20:04
non-zero (A, B) 的不一樣

相對誤差 abs(A-B)/abs(A or B)

也就是可能跟编译器的精度有关?

8

帖子

1

主题

0

精华

入门

F 币
48 元
贡献
21 点
发表于 5 天前 | 显示全部楼层
什么叫使用不同数据类型的实参虚参的传递啊……
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

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

GMT+8, 2018-10-24 00:37

Powered by Discuz! X3.2

© 2001-2017 Comsenz Inc.

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