python - 提高 numpy narrays 的精度

Question

我已经使用 numpy 编写了这个优化过程，我必须执行以下计算：

B 是对称负半定矩阵，I 是单位矩阵。问题是 B 的值可能太大（> 1e25），当减去对角矩阵时，一些元素不会因为有足够的精度。有什么办法可以克服这个问题，还是我需要使用另一个库（如 mpmath）重写它？

Answer 1

您需要重新实现，但实现可能很简单：更改您的数字表示。由于您没有指定您使用的基本类型或实际值范围（1e25 与什么？），我无法做出明确的建议。

但是，请注意，它numpy有几种精度更高的类型，包括 longdouble（80 位）、float96 和 float128。如果您需要超过 128 位（38 位精度），那么mpmath就是要走的路。