python - scipy.interp2d [分段错误（核心转储）]

Question

我对函数 scipy.interp2d 有疑问。这个例子很简单：

import matplotlib.pyplot as plt
from scipy.interpolate import interp2d
import numpy as np

n=60   #Critical parameter

x = np.linspace(-10,10,n)
X,Y = np.meshgrid(x,x)
r = X**2+Y**2
plt.imshow(r)
plt.show()

f = interp2d(x,x,r)
new_x = np.linspace(-10,10,10)
plt.imshow(f(new_x,new_x))
plt.show()

如果 n=2,...,59，一切正常，但当 n=60 时，出现分段错误。这是一个 scipy 错误？

另一方面，我确实注意到当函数的值很高并且域很小时，这个函数有很大的问题。我相信这个例程的实现存在一些内部问题。

任何提示？

Scipy 版本：0.11.0

Answer 1

pip uninstall numpy
pip uninstall scipy
pip install numpy --no-cache-dir
pip install scipy --no-cache-dir

这对我有用

崩溃的原因 - Python 共享库：RTLD_GLOBAL 段错误