def METHOD_NAME(self):
    # radial
    r = 2.0
    R = 1.0
    isotropic = Anisotropy(anisotropy_type="isotropic")
    kwargs_iso = {}
    beta = isotropic.beta_r(r, **kwargs_iso)
    k = isotropic.K(r, R, **kwargs_iso)
    f = isotropic.anisotropy_solution(r, **kwargs_iso)
    assert f == 1
    print(beta, "test")
    const = Anisotropy(anisotropy_type="const")
    kwargs = {"beta": beta}
    k_const = const.K(r, R, **kwargs)
    print(k, k_const)
    npt.assert_almost_equal(k, k_const, decimal=5)