def METHOD_NAME(self, universe):
    protein = universe.select_atoms('protein')
    total = np.sum([protein.center_of_geometry() for ts in
                    universe.trajectory])
    assert_almost_equal(total, self.ref_sum_centre_of_geometry, self.prec,
                        err_msg="sum of centers of geometry over the "
                                "trajectory do not match")