def test_angular_array_conversion():
"""
Test AngularConversion of AngularArray with regular values.
Conversion of [q in 1/nm] to [u in 'None'].
"""
spec_pnt_sys = SpectralPoint(
quantity="wavelength", unit="nm", value=525.)
ang_array_user = AngularArray(
quantity="q", unit="1/nm", dataarray=np.linspace(0., 0.024, 11))
ang_conv = AngularConverter(ang_array_user, spec_pnt_sys)
ang_array_sys = ang_conv.u_to_s(ang_array_user)
aa_sys_expected = np.array([0., 0.20053523, 0.40107046, 0.60160568,
0.80214091, 1.00267614, 1.20321137,
1.4037466, 1.60428183, 1.80481705,
2.00535228])
assert_allclose(
ang_array_sys.data, aa_sys_expected, rtol=2e-07, atol=1e-16)
assert_equal(ang_array_sys.quantity, "u")
assert_equal(ang_array_sys.unit, "None")
# reverse the conversion
ang_array_us2 = ang_conv.s_to_u(ang_array_sys)
assert_allclose(
ang_array_us2.data, ang_array_user.data, rtol=2e-07, atol=1e-16)
assert_equal(ang_array_us2.quantity, "q")
assert_equal(ang_array_us2.unit, "1/nm")
def test_angular_range_conversion():
"""
Test AngularConversion of AngularRange with regular values.
Conversion of [q in 1/nm] to [u in 'None']
"""
spec_pnt_sys = SpectralPoint(
quantity="wavelength", unit="nm", value=525.)
ang_range_user = AngularRange(
quantity="q", unit="1/nm", begin=0., end=0.012,
sampling_points=100, scale="linear")
ang_conv = AngularConverter(ang_range_user, spec_pnt_sys)
ang_range_sys = ang_conv.u_to_s(ang_range_user)
assert_equal(ang_range_sys.begin, 0.)
assert_allclose(ang_range_sys.end, 1.00267614, rtol=2e-07, atol=1e-16)
assert_equal(ang_range_sys.quantity, "u")
assert_equal(ang_range_sys.unit, "None")
assert_equal(ang_range_sys.sampling_points, 100)
assert_equal(ang_range_sys.scale, "linear")
# reverse the conversion
ang_range_us2 = ang_conv.s_to_u(ang_range_sys)
assert_allclose(
ang_range_user.begin, ang_range_us2.begin, rtol=2e-07, atol=1e-16)
assert_allclose(
ang_range_user.end, ang_range_us2.end, rtol=2e-07, atol=1e-16)
assert_equal(ang_range_us2.quantity, "q")
assert_equal(ang_range_us2.unit, "1/nm")
