Classes
class	Case

Variables
list	cases

int	n_interp = 37

	cubspl = CubicSpline(case.x, case.y, bc_type=case.bc)

	x_interp = np.linspace(case.x[0], case.x[-1], n_interp)

	y_interp = cubspl(x_interp)

	dy_interp = cubspl(x_interp, 1)

	d2y_interp = cubspl(x_interp, 2)

str	tag = 1 else 'second_deriv'

Detailed Description

This script generates reference data for the cross-check in cubic_spline_test.cpp.
To add more tests, append a new Case object to the cases list.

Variable Documentation

◆ cases

list gen_ref.cases

Initial value:

=  [
    Case(np.sin, np.logspace(-1, 0, 10), 'not-a-knot',
         'sin_not_a_knot.dat'),
    Case(np.cos, [0, 0.1, 1, 1.5, 3, 2*np.pi], 'periodic',
         'cos_periodic.dat'),
    Case(np.exp, np.logspace(-1, 0, 10), ((1, np.exp(0.1)), (1, np.exp(1))),
         'exp_first_deriv.dat'),
    Case(np.log, np.logspace(0, 1, 10), ((2, -1), (2, -0.01)),
         'log_second_deriv.dat'),
    Case(np.sqrt, np.logspace(0, 1, 10), ((1, 0.5), 'not-a-knot'),
         'sqrt_mix_bc.dat'),
    Case(lambda x: np.ones_like(x), [0, 1], 'periodic',
         'two_points_periodic.dat'),
    Case(np.arccos, [0, 0.5], ((1, -1), (1, -2/np.sqrt(3))),
         'two_points_first_deriv.dat'),
    Case(np.sqrt, [1, 4], ((2, 0.5), (2, 0.25)),
         'two_points_second_deriv.dat'),
    Case(np.sin, [0, 0.3, 1], 'not-a-knot',
         'three_points_not_a_knot.dat'),
]

◆ cubspl

gen_ref.cubspl = CubicSpline(case.x, case.y, bc_type=case.bc)

◆ d2y_interp

gen_ref.d2y_interp = cubspl(x_interp, 2)

◆ dy_interp

gen_ref.dy_interp = cubspl(x_interp, 1)

◆ n_interp

int gen_ref.n_interp = 37

◆ tag

str gen_ref.tag = 1 else 'second_deriv'

◆ x_interp

gen_ref.x_interp = np.linspace(case.x[0], case.x[-1], n_interp)

◆ y_interp

gen_ref.y_interp = cubspl(x_interp)

Classes

Variables