Purpose

---

  This function computes the geometric mean of a set of numbers. The
  geometric mean is defined as mean = (x1 * x2 * ... *
  xn)^(1/n). This procedure works with the logarithms of x, which
  prevents overflow/underflow issues when working with larger or
  small numbers.

Inputs

---

  data: A vector of data. Must be non-negative and finite

Keyword Parameters

---

  rms: On output, will contain the "logarithmic rms" of the data --
  that is, 10^(stdev(alog10(data))).

Outputs

---

  The geometric mean of the data

Modification History

---

  July 2010: Written by Chris Beaumont