Algorithms for convolutions and p-values (tail sums) that have guaranteed relative error, even for very small values.
Possible Rust port of the R and Python sisFFT algorithms.