def test__nondim_resist_spitzer(Z):
"""test _nondim_resist_spitzer function"""
alpha = _nondim_resist_spitzer(Z, 'par')
if Z == 1:
alpha_check = 0.5064
rtol = 1e-3
elif Z == 2 or Z == 4 or Z == np.inf:
alpha_check = _nondim_resist_braginskii(0, Z, 'par')
rtol = 2e-2
elif Z == 16:
alpha_check = _nondim_resist_ji_held(0, Z, 'par')
rtol = 2e-2
assert np.isclose(alpha, alpha_check, rtol=rtol)
def test_unmagnetized(self, Z):
"""confirm perp -> par as B -> 0"""
alpha_hat_par = _nondim_resist_braginskii(self.small_hall, Z, 'par')
alpha_hat_perp = _nondim_resist_braginskii(self.small_hall, Z, 'perp')
assert np.isclose(alpha_hat_par, alpha_hat_perp, rtol=1e-3)
def test_cross_vs_ji_held(self, Z):
"""cross should roughly agree with ji-held at hall 0.1"""
alpha_hat_cross_brag = _nondim_resist_braginskii(0.1, Z, 'cross')
alpha_hat_cross_jh = _nondim_resist_ji_held(0.1, Z, 'cross')
assert np.isclose(alpha_hat_cross_brag, alpha_hat_cross_jh, rtol=5e-2)
def test_known_values_par(self, Z, field_orientation, expected):
"""check some known values"""
beta_hat = _nondim_resist_braginskii(self.big_hall, Z,
field_orientation)
assert np.isclose(beta_hat, expected, atol=1e-2)