def test_ellipsoid():
"""Test Ellipsoid initialiser and verify computed attributes."""
wgs84_polar_radius = 6356752.314245 * u.m
wgs84_ellipsoid1 = Ellipsoid(a=WGS84_EQUATORIAL_RADIUS, f=WGS84_FLATTENING)
wgs84_ellipsoid2 = Ellipsoid(a=WGS84_EQUATORIAL_RADIUS, b=wgs84_polar_radius)
assert qisclose(wgs84_ellipsoid1.b, wgs84_ellipsoid2.b)
assert qisclose(wgs84_ellipsoid1.f, wgs84_ellipsoid2.f)
with pytest.raises(TypeError):
Ellipsoid(a=WGS84_EQUATORIAL_RADIUS)
with pytest.raises(TypeError):
Ellipsoid(a=WGS84_EQUATORIAL_RADIUS, b=wgs84_polar_radius,
f=WGS84_FLATTENING)
def test_ellipsoid():
"""Test Ellipsoid initialiser and verify computed attributes."""
wgs84_polar_radius = 6356752.314245 * u.m
wgs84_ellipsoid1 = Ellipsoid(a=WGS84_EQUATORIAL_RADIUS, f=WGS84_FLATTENING)
wgs84_ellipsoid2 = Ellipsoid(a=WGS84_EQUATORIAL_RADIUS,
b=wgs84_polar_radius)
assert qisclose(wgs84_ellipsoid1.b, wgs84_ellipsoid2.b)
assert qisclose(wgs84_ellipsoid1.f, wgs84_ellipsoid2.f)
with pytest.raises(TypeError):
Ellipsoid(a=WGS84_EQUATORIAL_RADIUS)
with pytest.raises(TypeError):
Ellipsoid(a=WGS84_EQUATORIAL_RADIUS,
b=wgs84_polar_radius,
f=WGS84_FLATTENING)
def test_celestial_body():
"Test CelestialBody initialiser and verify computed attributes"
a = CelestialBody(name='earth', ellipsoid=wgs84, mu=wgs84.mu)
b = CelestialBody.from_reference_ellipsoid(name='earth', ellipsoid=wgs84)
assert a.mu == b.mu == wgs84.mu
ellipsoid = Ellipsoid(a=1 * u.m, b=1 * u.m)
c = CelestialBody(name='earth', ellipsoid=ellipsoid,
mu=GEOCENTRIC_GRAVITATIONAL_CONSTANT)
assert qisclose(c.mass, EARTH_MASS)
def test_celestial_body():
"Test CelestialBody initialiser and verify computed attributes"
a = CelestialBody(name='earth', ellipsoid=wgs84, mu=wgs84.mu)
b = CelestialBody.from_reference_ellipsoid(name='earth', ellipsoid=wgs84)
assert a.mu == b.mu == wgs84.mu
ellipsoid = Ellipsoid(a=1 * u.m, b=1 * u.m)
c = CelestialBody(name='earth',
ellipsoid=ellipsoid,
mu=GEOCENTRIC_GRAVITATIONAL_CONSTANT)
assert qisclose(c.mass, EARTH_MASS)
