///|
/// Calculate the reciprocal of square root of the sum of squares of any number of coordinates.
pub fn rnorm(vec : Array[Double]) -> Double {
  1.0 / norm(vec)
}