pg_utils.numerics.special.eval_jacobi_recur_mpmath

pg_utils.numerics.special.eval_jacobi_recur_mpmath(Nmax: int, alpha: mpf, beta: mpf, z: ndarray, dps: int = 33) → ndarray[source]

Evaluate Jacobi polynomials with recurrence relation up to a degree, to (arbitrary) multi-precision.

This functions generates values for Jacobi polynomials from degree 0 up to a specified degree, using recurrence relations.

Parameters:

Nmax (int) – maximum degree, required to be >= 1
alpha (mpmath.mpf) – alpha index for Jacobi polynomials
beta (mpmath.mpf) – beta index for Jacobi polynomials
z (np.ndarray) – 1-D array of grid points where the Jacobi polynomials are to be evaluated; assumed to be within interval [-1, +1]
dps (int) – number of decimal places for calculation

Returns:

Array with shape (Nmax + 1, z.size), values for Jacobi polynomials at grid points specified in z.

Note

The input parameters alpha and beta as well as z need to at least match the precision of the desired output, otherwise the multi-precision evaluation is meaningless.