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