def METHOD_NAME(self, temp_kev, sigmav):
    temp_k = raw_uc(temp_kev, "keV", "K")
    result = reactivity(temp_k, reaction="D-D1", method="Bosch-Hale")
    np.testing.assert_allclose(result, sigmav, rtol=0.003, atol=0)