Added functions for calculating the geometric mean and standard deviation to the relax library.
These are the geometric_mean() and geometric_std() functions of the lib.statistics module. The
implementation is designed to be fast, using numpy array arithmetic rather than Python loops.
