Fortran DMOD vs Python numpy.remainder [重复]

Question

为什么在 Fortran 中使用“DMOD”和在 Python 中使用“numpy.remainder”会得到不同的结果？

在 Fortran 中：

DOUBLE PRECISION aaa,bbb
aaa = 2.11561
bbb = 1.19841
out = DMOD(aaa, bbb)
print*,out

我得到：0.917199850

在python中：

import numpy as np
print(np.remainder(2.11561,1.19841))

我得到：0.9172000000000002

Answer 1

如果你想在 Fortran 中实现双精度，你也需要这样写常量。你可以这样做，通过使用_dp 后缀：

integer, parameter:: dp=kind(0.d0)
aaa = 2.11561_dp

或者您可以在号码中使用d：

aaa = 2.11561d0

e 用于单精度，d 用于双精度。单精度是默认的……即使变量是双精度的，常量也可能是单精度的！在 C 语言中也是如此。

区别在于Python和C默认使用双精度，你必须要求单精度。