discontinuity_jumps

sofia_redux.toolkit.splines.spline_utils.discontinuity_jumps(knot_line, n_knot, degree, b_spline)[source]

Calculates the discontinuity jumps.

Calculates the discontinuity jumps of the kth derivative of the b-splines of degree k at the knots k+2 -> n - k - 1. The results are updated in-place in the b_spline array.

Adapted from the Fortran function fpdisc in the fitpack library.

Parameters:
knot_linenumpy.ndarray (float)

The knot locations in a single dimension of shape (>= n_knots,).

n_knotint

The number of knots in the knot line.

degreeint

The degree of the spline.

b_splinenumpy.ndarray (float)

An array of shape (max(knot_estimate), k + 2) where k is the degree of the spline containing the spline coefficients. Values will be updated in-place.

Returns:
None