1 2 3 4
/// Modulus of p-vector. pub fn pm(p: [f64; 3]) -> f64 { (p[0] * p[0] + p[1] * p[1] + p[2] * p[2]).sqrt() }