dedalus.libraries.dedalus_sphere.sphere

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.

Parameters

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.

Parameters

• 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)