草庐IT

Vandermonde

全部标签

python - 计算 Vandermonde 矩阵的有效方法

我正在计算Vandermondematrix对于相当大的一维阵列。做到这一点的自然而干净的方法是使用np.vander().但是,我发现这大约是。比基于列表理解的方法慢2.5倍。In[43]:x=np.arange(5000)In[44]:N=4In[45]:%timeitnp.vander(x,N,increasing=True)155µs±205nsperloop(mean±std.dev.of7runs,10000loopseach)#oneofthelistedapproachesfromthedocumentationIn[46]:%timeitnp.flip(np.colu