def test_vs_ji_held_K2(self):
"""confirm agreement with ji-held model when K=2"""
eta_i_hat_brag = _nondim_visc_i_braginskii(self.big_hall)
eta_i_hat_jh = _nondim_visc_i_ji_held(self.big_hall, 1, 0, 100, K=2)
for idx in [0, 1, 2, 3, 4]:
assert np.isclose(eta_i_hat_brag[idx],
eta_i_hat_jh[idx],
rtol=2e-2)
def test__nondim_visc_i_ji_held(hall, Z, mu, theta, index, expected):
"""test _nondim_visc_i_ji_held function"""
kappa_hat = _nondim_visc_i_ji_held(hall, Z, mu, theta)
kappa_check = expected
assert np.isclose(kappa_hat[index], kappa_check, rtol=2e-2)
