IVF中实参虚参类型不一致
使用的一个软件包，是linux下写的，使用了不同数据类型的实参虚参的传递，但是IVF不支持，应该怎么办？跟大家请教一下，谢谢暂时通过取消接口检查解决的问题，但是结果跟Linux下运行出来的结果不太一样。这个是软件说明的文档。
结果不一致，很难说哪里的问题。你最好把linux下的makefile拿出来看看，采用相同的编译参数。 li913 发表于 2018-9-26 13:05
结果不一致，很难说哪里的问题。你最好把linux下的makefile拿出来看看，采用相同的编译参数。 ...
好的，我再看看，主要是不太清楚怎么允许这种数据传递，同时不出错 PROGRAM a
REAL :: x
CALL b(x)
CONTAINS
SUBROUTINE b(i)
INTEGER ::ｉ
END SUBROUTINE b
END PROGRAM a
a 與 b 之間為 explicit interface, (type inconsistent) -> (compile-time fatal error, all compilers)
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, ... chiangtp 发表于 2018-9-26 16:54
PROGRAM a
REAL :: x
CALL b(x)
就是计算的结果有时候不太一样，但有时候又是一样的。
谢谢版主！ 怎的(不太一样)？
1.23456789, 1.23456987,for 4-byte REAL, 這算正常
否則, 你 call FFTPack 的 "codes" 出錯的機會極大 chiangtp 发表于 2018-9-26 19:47
怎的(不太一样)？
1.23456789, 1.23456987,for 4-byte REAL, 這算正常
我只是单纯看的结果，那个软件是进行矩阵分解C=BP的，分解出来的矩阵B是C的某些列，IVF和Linux的对应的列数是不一致的。
虽然结果不太一样，但感觉应该都是对的 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, 都算正常 chiangtp 发表于 2018-9-26 20:04
non-zero (A, B) 的不一樣
相對誤差 abs(A-B)/abs(A or B)
也就是可能跟编译器的精度有关？ 什么叫使用不同数据类型的实参虚参的传递啊……
页:
[1]
2