所以我有两个版本的流体动力学代码,它们具有相同的基础物理特性。让我们称它们为代码 A 和 B。但是代码 B 更优化且更面向对象。我试图比较两个代码的数值结果,看看它们是否兼容。我这样做的方式是使用相同的初始参数运行两个代码并计算特定的流体动力学值 X。
X 的数量级为 1E6。为了比较这两个代码,我实际上是在做
在哪里是一些错误。测试是给.
这意味着两个代码的前 8 位数字都吐出相同的 X。
这个错误有多严重?我对非整数使用双精度,这样的大小是否仅仅意味着水电代码中涉及的数千次计算的累积舍入误差?还是我缺少的根本区别?