def test_prolate30_direct(testcase):
"""Test direct method on the prolate 30 ellipsoid"""
options = dict(frame=PROLATE30, degrees=True)
(lat1, lon1, azi1, lat2, lon2, azi2, s12) = testcase[:7]
point1 = GeoPoint(lat1, lon1, **options)
point2, az_b = point1.displace(s12, azi1, long_unroll=True, degrees=True)
lat_b, lon_b = point2.latitude_deg, point2.longitude_deg
assert lat_b == approx(lat2, abs=1e-13)
assert lon_b == approx(lon2, abs=1e-13)
assert az_b == approx(azi2, abs=1e-13)
