def METHOD_NAME(self):
    anisoClass = Anisotropy(anisotropy_type="const")
    kwargs = {"beta": 1.0}
    k = anisoClass.K(self._r_array, self._R_array, **kwargs)
    npt.assert_almost_equal(k[0], 0.61418484930437822, decimal=5)
    kwargs = {"beta": -0.49}
    k = anisoClass.K(self._r_array, self._R_array, **kwargs)
    npt.assert_almost_equal(k[0], 0.7645553632433857, decimal=5)
    anisoClass = Anisotropy(anisotropy_type="Colin")
    kwargs = {"r_ani": 1}
    k = anisoClass.K(self._r_array, self._R_array, **kwargs)
    npt.assert_almost_equal(k[0], 0.91696135187291117, decimal=5)
    k = anisoClass.K(self._r_array, self._R_array - 0.001, **kwargs)
    npt.assert_almost_equal(k[0], 0.91696135187291117, decimal=2)
    k = anisoClass.K(self._r_array, self._R_array + 0.001, **kwargs)
    npt.assert_almost_equal(k[0], 0.91696135187291117, decimal=2)
    anisoClass = Anisotropy(anisotropy_type="radial")
    kwargs = {}
    k = anisoClass.K(self._r_array, self._R_array, **kwargs)
    npt.assert_almost_equal(k[0], 0.61418484930437856, decimal=5)
    anisoClass = Anisotropy(anisotropy_type="isotropic")
    kwargs = {}
    k = anisoClass.K(self._r_array, self._R_array, **kwargs)
    npt.assert_almost_equal(k[0], 0.8660254037844386, decimal=5)
    anisoClass = Anisotropy(anisotropy_type="OM")
    kwargs = {"r_ani": 1}
    k = anisoClass.K(self._r_array, self._R_array, **kwargs)
    npt.assert_almost_equal(k[0], 0.95827704196894481, decimal=5)
    R = 0.1
    k = anisoClass.K(R, R, **kwargs)
    npt.assert_almost_equal(k, 0, decimal=5)