def METHOD_NAME(self):
    # Modified gmpe
    gmpe = ChiouYoungs2014()
    mgmpe = CY14SiteTerm(gmpe_name='ChiouYoungs2014')
    imt = PGA()
    # Compute results
    mean, stds = mgmpe.get_mean_and_stddevs(
        self.ctx, self.ctx, self.ctx, imt, self.stdt)
    # Compute the expected results
    mean_expected, stds_expected = gmpe.get_mean_and_stddevs(
        self.ctx, self.ctx, self.ctx, imt, self.stdt)
    # Test that for reference soil conditions the modified GMPE gives the
    # same results of the original gmpe
    np.testing.assert_almost_equal(mean, mean_expected, decimal=7)
    # Here we use a quite large tolerance since in the site term we take
    # the std from the calculation of motion on reference rock. This
    # does not match the std that the same GMM computes for soft soils
    # with the same remaining conditions
    np.testing.assert_almost_equal(stds, stds_expected, decimal=1)