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)