我不是编程新手,但我是Python和数值计算的新手,所以这可能是一个新手问题。
我使用GMPY2来处理非常大的数字的计算,但是我想使用MatPlotLib来绘制结果。我目前有
import gmpy2
from gmpy2 import mpz, mpq, mpfr, mpc
import numpy as np
import matplotlib.pyplot as plt
gmpy2.get_context().precision = 100
def hop(x, y, n):
if n > 0:
return gmpy2.exp(hop(gmpy2.log(x), gmpy2.log(y), n-1))
elif n == 0:
return gmpy2.add(x, y)
else:
raise ValueError('Negative value of n in hop method.')
t = np.arange(-5.0, 5.0, 0.01)
plt.plot(t, hop(t, t, 0))我假设没有从mpfr到numpy的本机类型的隐式转换。在np.arange中,我尝试设置dtype=mpfr,在plot方法中,我尝试编写hop(mpfr(t), mpfr(t), 0),但当然没有成功。
https://stackoverflow.com/questions/38118252
复制相似问题