pub fn factorial_size_for_exp() -> usize
Expand description

Return N such that 1/(N-7)! < f64::EPSILON.

This number is needed for the computation of exponential matrix