Struct Cubic

pub struct Cubic<T>
where
    T: Float + Debug,
{ /* private fields */ }

Cubic Spline.

Cubic spline with natural boundary conditions. The resulting curve is piecewise cubic on each interval, with matching first and second derivatives at the supplied data-points. The second derivative is chosen to be zero at the first and last point.

Example

let xa = [0.0, 1.0, 2.0];
let ya = [0.0, 2.0, 4.0];
let interp = Cubic::new(&xa, &ya)?;

Reference

Numerical Algorithms with C - Gisela Engeln-Mullges, Frank Uhlig - 1996 - Algorithm 10.1, pg 254

Trait Implementations

impl<T> Interpolation<T> for Cubic<T>

where T: Float + Debug + Scalar + Lapack,

const MIN_SIZE: usize = 3usize

The minimum number of points required by the interpolator. For example, Akima spline interpolation requires a minimum of 5 points.

const NAME: &'static str = "cubic"

The name of the interpolator.

fn new(xa: &[T], ya: &[T]) -> Result<Self, InterpolationError>

where Self: Sized,

Creates a new interpolator.

fn eval( &self, xa: &[T], ya: &[T], x: T, acc: &mut Accelerator, ) -> Result<T, DomainError>

Returns the interpolated value y for a given point x, using the data arrays xa and ya and the Accelerator acc.

fn eval_deriv( &self, xa: &[T], ya: &[T], x: T, acc: &mut Accelerator, ) -> Result<T, DomainError>

Returns the derivative dx/dy of an interpolated function for a given point x, using the data arrays xa and ya and the Accelerator acc.

fn eval_deriv2( &self, xa: &[T], ya: &[T], x: T, acc: &mut Accelerator, ) -> Result<T, DomainError>

Returns the second derivative d²x/dy² of an interpolated function for a given point x, using the data arrays xa and ya and the Accelerator acc.

fn eval_integ( &self, xa: &[T], ya: &[T], a: T, b: T, acc: &mut Accelerator, ) -> Result<T, DomainError>

Returns the numerical integral of an interpolated function over the range [a,b], using the data arrays xa and ya and the Accelerator acc.