请教python(numpy)与fortran在计算精度方面的差异

lon91ong · 发表于 2014-8-28 09:04:36

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由 lon91ong 于 2014-8-29 10:15 编辑

因为在我自己的程序中发现, 同样都是单精度计算的前提下, python (numpy)计算的结果似乎不如fortran的精度高

我是做有限差分计算的时候发现这个问题的, 相同的网格划分条件下, fortran的计算结果可以看到明显的频散现象
而python的计算几乎没有频散, 就好像图像结果被"柔化"处理过的一样

但事实上是没有任何后续处理的, 两种程序的逻辑是完全相同的

百思不得其解

烦请大侠解惑, 感激万分!

补充两个图:

问题解决了, 原来是Fortran程序有点细节方面的问题, 震源函数中应该用浮点数的地方没加小数点

freekiller · 发表于 2014-8-28 09:27:33

不懂帮顶，羡慕一下会Python的

lon91ong · 发表于 2014-8-28 17:31:22

麻烦英语好的大侠帮忙把问题翻译过来发Stack Overflow, 上次发过的问题被评成负分

burf · 发表于 2014-9-9 13:09:47

精度好像不一样，python默认float64 fortran好像是float32？记不清了……你试试看，设置好经度后没什么区别…我用f2py没发现精度问题…… 但是貌似如果没设置好的话python内部计算都会出问题…之前做一个lorenz63的小试验，有一个矩阵里存在10e-6量级的误差…结果最后结果完全不一样Orz 不过这点误差对于大型系统应该不算什么，应该会自动慢慢消耗掉，但是系统越简单越容易受影响吧

		自动登录	找回密码
密码			立即注册

[经验总结] 请教python(numpy)与fortran在计算精度方面的差异

登录后查看更多精彩内容~