Struct rgsl::types::series_acceleration::LevinUWorkspace

pub struct LevinUWorkspace { /* private fields */ }

Workspace for Levin U Transform with error estimation

Implementations

impl LevinUWorkspace

pub fn new(n: usize) -> Option<LevinUWorkspace>

This function allocates a workspace for a Levin u-transform of n terms. The size of the workspace is O(2n^2 + 3n).

pub fn accel(
    &mut self,
    array: &[f64],
    sum_accel: &mut f64,
    abserr: &mut f64
) -> Value

This function takes the terms of a series in array of size array_size and computes the extrapolated limit of the series using a Levin u-transform. Additional working space must be provided in w. The extrapolated sum is stored in sum_accel, with an estimate of the absolute error stored in abserr. The actual term-by-term sum is returned in w->sum_plain. The algorithm calculates the truncation error (the difference between two successive extrapolations) and round-off error (propagated from the individual terms) to choose an optimal number of terms for the extrapolation. All the terms of the series passed in through array should be non-zero.

pub fn sum_plain(&self) -> f64

pub fn terms_used(&self) -> usize

pub fn size(&self) -> usize

Trait Implementations

impl Drop for LevinUWorkspace

fn drop(&mut self)

Executes the destructor for this type.