虽然我不知道有什么适合的工具箱，但 Julia 是 MATLAB 和 Octave 的一个很好的开源替代品，它确实可以很好地解决这个问题。MATLAB/Octave 的线性代数语法几乎直接转移到 Julia，尽管您需要A(i)像A[i]. 但经过快速翻译后，您可以将 JuliaBigFloats用于可变精度数字。这是 MPFR 库的包装器，它非常标准并且可以。基于 Arb 库的 ArbFloats.jl 中给出了比特大小 <500 的更快实现。最后，如果您只想要 128 位数字，DecFP.jl 提供了一些很好的实现。

走这条路的好处是 Julia 的调度系统使这些只是标准数字的替代品，因此您可以在大多数编写良好的库中使用它们。因此，您可以使用这些高速可变精度数字，只需调用通用矩阵乘法、SVD 分解、\微分方程求解器、优化例程等，并期望它们代替Float64s 工作。