@staticmethod
def center_of_mass(points):
    """
        Finds the center of mass of the points given.

        To find the "geometric center" of a gamut lets points be only
        the vertices of the gamut.

        Thanks to: http://stackoverflow.com/questions/8917478/center-of-mass-of-a-numpy-array-how-to-make-less-verbose

        Parameters
        ----------
        points : ndarray
            Shape(N, 3), a list of points

        Returns
        -------
        center : ndarray
            Shape(3,), the coordinate of the center of mass.
        """
    cm = points.sum(0) / points.shape[0]
    for i in range(points.shape[0]):
        points[(i), :] -= cm
    return cm
