1 2 3 4
/// p-vector inner (=scalar=dot) product. pub fn pdp(a: &[f64; 3], b: &[f64; 3]) -> f64 { a[0] * b[0] + a[1] * b[1] + a[2] * b[2] }