Struct Work 

pub struct Work<T: Lane, const B: usize> {
    pub kt: usize,
    pub c: Vec<Vec<Fl<T, B>>>,
    pub aux: Vec<Vec<Fl<T, B>>>,
    pub rate: Rate,
}

The Taylor coefficients of one step about its start x0 and the equations’ scratch series.

c[f][k] is coefficient k of field f, k = 0..=kt. The engine sets order 0 (the state at x0) before Equations::taylor_coeffs and reads orders 0..=kt after it. aux[i] are NAUX series the recursion may keep (products, right-hand sides), each of the same length; their contents are the equations’.

Fields

kt: usize

the Taylor order: coefficients 0..=kt

c: Vec<Vec<Fl<T, B>>>

aux: Vec<Vec<Fl<T, B>>>

rate: Rate

the growth rate of a mode-matrix evaluation, zero otherwise: the recursion passes each t-derivative of a coefficient through rate.shift (crate::mode)

Implementations

impl<T: Lane, const B: usize> Work<T, B>

pub fn new(nf: usize, naux: usize, kt: usize) -> Self

work for nf fields, naux scratch series and order kt; every series has kt + 2 entries, zero