rfa 0.5.9

#[inline]
pub fn pm(p: &[f64;3])->f64{
    ( p[0]*p[0] + p[1]*p[1] + p[2]*p[2] ).sqrt()
}