def test_qd2geo_invalid_lat():
A = Apex(date=2000, refh=300)
with pytest.raises(ValueError):
A.qd2geo(91, 0, 0, precision=1e-2)
with pytest.raises(ValueError):
A.qd2geo(-91, 0, 0, precision=1e-2)
A.qd2geo(90, 0, 0, precision=1e-2)
A.qd2geo(-90, 0, 0, precision=1e-2)
assert_allclose(A.qd2geo(90+1e-5, 0, 0, 1e-2), A.qd2geo(90, 0, 0, 1e-2), rtol=0, atol=1e-8)
def test_qd2geo_invalid_lat():
apex_out = Apex(date=2000, refh=300)
with pytest.raises(ValueError):
apex_out.qd2geo(91, 0, 0, precision=1e-2)
with pytest.raises(ValueError):
apex_out.qd2geo(-91, 0, 0, precision=1e-2)
apex_out.qd2geo(90, 0, 0, precision=1e-2)
apex_out.qd2geo(-90, 0, 0, precision=1e-2)
assert_allclose(apex_out.qd2geo(90 + 1e-5, 0, 0, 1e-2),
apex_out.qd2geo(90, 0, 0, 1e-2),
rtol=0,
atol=1e-8)
def test_qd2geo():
A = Apex(date=2000, refh=300)
lat, lon, error = A.qd2geo(60, 15, 100, precision=1e-2)
assert_allclose((lat, lon, error), A._qd2geo(60, 15, 100, 1e-2))
assert type(lat) != np.ndarray
assert type(lon) != np.ndarray
assert type(error) != np.ndarray
def test_qd2geo():
apex_out = Apex(date=2000, refh=300)
lat, lon, error = apex_out.qd2geo(60, 15, 100, precision=1e-2)
assert_allclose((lat, lon, error), apex_out._qd2geo(60, 15, 100, 1e-2))
assert type(lat) != np.ndarray
assert type(lon) != np.ndarray
assert type(error) != np.ndarray
def test_apex2geo():
A = Apex(date=2000, refh=300)
lat, lon, error = A.apex2geo(60, 15, 100, precision=1e-2)
assert_allclose((lat, lon, error),
A.qd2geo(*A.apex2qd(60, 15, 100),
height=100,
precision=1e-2))
assert type(lat) != np.ndarray
assert type(lon) != np.ndarray
assert type(error) != np.ndarray
def test_basevectors_apex_scalar_qd():
A = Apex(date=2000, refh=300)
f1, f2, f3, g1, g2, g3, d1, d2, d3, e1, e2, e3 = A.basevectors_apex(60, 15, 100, coords='qd', precision=1e-2)
glat, glon, _ = A.qd2geo(60, 15, 100, precision=1e-2)
_, _, _, _, f1_, f2_, _, d1_, d2_, d3_, _, e1_, e2_, e3_ = A._geo2apexall(glat, glon, 100)
assert_allclose(f1, f1_)
assert_allclose(f2, f2_)
assert_allclose(d1, d1_)
assert_allclose(d2, d2_)
assert_allclose(d3, d3_)
assert_allclose(e1, e1_)
assert_allclose(e2, e2_)
assert_allclose(e3, e3_)
def test_basevectors_apex_scalar_qd():
apex_out = Apex(date=2000, refh=300)
(f1, f2, f3, g1, g2, g3, d1, d2, d3, e1, e2,
e3) = apex_out.basevectors_apex(60, 15, 100, coords='qd', precision=1e-2)
glat, glon, _ = apex_out.qd2geo(60, 15, 100, precision=1e-2)
(_, _, _, _, f1_, f2_, _, d1_, d2_, d3_, _, e1_, e2_,
e3_) = apex_out._geo2apexall(glat, glon, 100)
assert_allclose(f1, f1_)
assert_allclose(f2, f2_)
assert_allclose(d1, d1_)
assert_allclose(d2, d2_)
assert_allclose(d3, d3_)
assert_allclose(e1, e1_)
assert_allclose(e2, e2_)
assert_allclose(e3, e3_)
def test_basevectors_qd_scalar_qd():
A = Apex(date=2000, refh=300)
glat, glon, _ = A.qd2geo(60, 15, 100, precision=1e-2)
assert_allclose(A.basevectors_qd(60, 15, 100, coords='qd', precision=1e-2), A._basevec(glat, glon, 100))
def test_qd2geo_vectorization():
A = Apex(date=2000, refh=300)
assert A.qd2geo([60, 60], 15, 100)[0].shape == (2,)
assert A.qd2geo(60, [15, 15], 100)[0].shape == (2,)
assert A.qd2geo(60, 15, [100, 100])[0].shape == (2,)
def test_convert_qd2geo():
A = Apex(date=2000, refh=300)
assert_allclose(A.convert(60, 15, 'qd', 'geo', height=100, precision=1e-2), A.qd2geo(60, 15, 100, precision=1e-2)[:-1])
def test_qd2geo_vectorization():
apex_out = Apex(date=2000, refh=300)
assert apex_out.qd2geo([60, 60], 15, 100)[0].shape == (2, )
assert apex_out.qd2geo(60, [15, 15], 100)[0].shape == (2, )
assert apex_out.qd2geo(60, 15, [100, 100])[0].shape == (2, )
def test_convert_qd2geo():
apex_out = Apex(date=2000, refh=300)
assert_allclose(
apex_out.convert(60, 15, 'qd', 'geo', height=100, precision=1e-2),
apex_out.qd2geo(60, 15, 100, precision=1e-2)[:-1])
def test_basevectors_qd_scalar_qd():
apex_out = Apex(date=2000, refh=300)
glat, glon, _ = apex_out.qd2geo(60, 15, 100, precision=1e-2)
assert_allclose(
apex_out.basevectors_qd(60, 15, 100, coords='qd', precision=1e-2),
apex_out._basevec(glat, glon, 100))
