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))