提高numpy narrays的精度

Question

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

B是对称负半定矩阵，I是单位矩阵。 问题是，B的值可能太大（>；1e25），当用对角矩阵减去时，有些元素不会，因为有足够的精度。 有没有办法克服这个问题，或者我需要用另一个像mpmath这样的库来重写它？你知道吗

Answer 1

您需要重新实现，但实现可能很简单：更改数值表示形式。因为您没有指定要使用的基类型或实际值范围（1e25与什么？），我不能给出明确的建议。你知道吗

但是，请注意numpy有几种精度更高的类型，包括longdouble（80位）、float96和float128。如果您需要超过128位（38位精度），那么mpmath就是最好的选择。你知道吗