pub fn fwd97(x: &mut [f64], tmp: &mut Vec<f64>)
One level of forward CDF 9/7 on x (in place: low half first).
x