Floating point 高精度数值计算

Floating point 高精度数值计算,floating-point,numeric,numerical-methods,floating-point-precision,Floating Point,Numeric,Numerical Methods,Floating Point Precision,我想探讨我的算法的几个不同实现的数值精度(使用标准的双精度算法)。不幸的是,在许多情况下,我不知道封闭形式的正确结果,因此我需要找到一种方法,使用一些非常高精度的计算来计算基准测试结果 这是一个有趣的项目,所以我的限制是:没有工具预算,最好是Linux平台。我知道Mathematica提供自动错误跟踪和任意精度算法,但我没有许可证。执行速度不是问题,因为这些高精度计算只用于计算基准 编写这些高精度计算的最佳方法是什么?我正在寻找至少四精度,但最好更高。到目前为止,我唯一的想法就是使用。我发现这项

我想探讨我的算法的几个不同实现的数值精度(使用标准的双精度算法)。不幸的是,在许多情况下,我不知道封闭形式的正确结果,因此我需要找到一种方法,使用一些非常高精度的计算来计算基准测试结果

这是一个有趣的项目,所以我的限制是:没有工具预算,最好是Linux平台。我知道Mathematica提供自动错误跟踪和任意精度算法,但我没有许可证。执行速度不是问题,因为这些高精度计算只用于计算基准

编写这些高精度计算的最佳方法是什么?我正在寻找至少四精度,但最好更高。到目前为止,我唯一的想法就是使用。

我发现这项任务非常适合。优势之一

  • 它非常容易编写,并且具有方便的REPL或可以使用IPython笔记本界面
  • 它带有一个本地任意精度
    BigFloat
    类型,用于包装MPFR库,以及一个
    BigInt
    用于包装GMP库
  • “偷窥引擎盖下”,看看其他方法做什么,看看llvm和本机汇编代码,这是非常重要的

我不是把它作为一个答案发布的,因为它只是一个链接,而是提供了您需要的所有任意精度浮点原语。您可以用这个标签找到很多解决方案