Python和相关的Legendre多项式

Question

I have been searching for a python implementation of the associated Legendre polynomials quite a long time and have found nothing satisfying me. There is an implementation in scipy.special , but it is not vectorized. I have found a solution to use pygsl interface with gsl library, but I had a hard time to get everything compiled.

Does anyone know better solution to get access to associated Legendre polynomials in efficiently vectorized way, ie Legendre functions has to be applied for multidimensional matrices?

Answer 1

scipy.special.lpmv已向量化。