Module Contents

dtype = 'longdouble'
quadrature(Lmax, dtype=dtype)
Generates the Gauss quadrature grid and weights for spherical harmonics transform.

Returns cos_theta, weights

Will integrate polynomials on (-1,+1) exactly up to degree = 2*Lmax+1.


Lmax (int >=0; spherical-harmonic degree.)

spin2Jacobi(Lmax, m, s, ds=None, dm=None)
harmonics(Lmax, m, s, cos_theta, **kwargs)

Gives spin-wieghted spherical harmonic functions on the Gauss quaduature grid. Returns an array with

shape = ( Lmax - Lmin(m,s) + 1, len(z) )

or (Lmax - Lmin(m,s) + 1,) if z is a single point.

  • Lmax (int >=0; spherical-harmonic degree.)

  • m,s (int) – spherical harmonic parameters.

  • cos_theta (np.ndarray or float.)

  • dtype (output dtype. internal dtype = ‘longdouble’.)

operator(name, dtype=dtype)

Interface to base ShereOperator class.

class SphereOperator(name, radius=1, dtype=dtype)
property radius
property dtype
static identity(dtype=dtype)
static parity(dtype=dtype)
static L(dtype=dtype)
static M(dtype=dtype)
static S(dtype=dtype)
class SphereCodomain(dL=0, dm=0, ds=0, pi=0)