numpy任意精度线性代数

Question

我有一个numpy的2d数组[中/大型-
例如500x500]。我想找到它的元素智能指数的特征值。问题在于某些值是非常负的（-800，-1000等），并且它们的指数下溢（意味着它们非常接近零，因此numpy将其视为零）。无论如何在numpy中使用任意精度？

我做梦的方式：

import numpy as np

np.set_precision('arbitrary') # <--- Missing part
a = np.array([[-800.21,-600.00],[-600.00,-1000.48]])
ex = np.exp(a)  ## Currently warns about underflow
eigvals, eigvecs = np.linalg.eig(ex)

我已经搜索了gmpy和mpmath的解决方案，但无济于事。任何想法都将受到欢迎。

Answer 1

SymPy可以计算任意精度：

from sympy import exp, N, S
from sympy.matrices import Matrix

data = [[S("-800.21"),S("-600.00")],[S("-600.00"),S("-1000.48")]]
m = Matrix(data)
ex = m.applyfunc(exp).applyfunc(lambda x:N(x, 100))
vecs = ex.eigenvects()
print vecs[0][0] # eigen value
print vecs[1][0] # eigen value
print vecs[0][2] # eigen vect
print vecs[1][2] # eigen vect

输出：

-2.650396553004310816338679447269582701529092549943247237903254759946483528035516341807463648841185335e-261
2.650396553004310816338679447269582701529092549943247237903254759946483528035516341807466621962539464e-261
[[-0.9999999999999999999999999999999999999999999999999999999999999999999999999999999999999994391176386872]
[                                                                                                      1]]
[[1.000000000000000000000000000000000000000000000000000000000000000000000000000000000000000560882361313]
[                                                                                                    1]]

您可以将N（x，100）中的100更改为其他精度，但是，当我尝试1000时，本征vect的计算失败。