Function knot 

fn knot(taus: &[f64], i: usize) -> f64

Knot i of the sequence, without materialising it.

The evaluation path runs per query, so it indexes the anchor and the free pillars in place; building the vector here would allocate on every discount lookup, which is the cost the coefficient seam exists to avoid.