public void GeodSolve71()
{
// Check that DirectLine sets s13.
GeodesicLine line = Geodesic.WGS84.DirectLine(1, 2, 45, 1e7);
GeodesicData dir =
line.Position(0.5 * line.Distance(),
GeodesicMask.STANDARD | GeodesicMask.LONG_UNROLL);
AssertEquals(dir.lat2, 30.92625, 0.5e-5);
AssertEquals(dir.lon2, 37.54640, 0.5e-5);
AssertEquals(dir.azi2, 55.43104, 0.5e-5);
}
public void GeodSolve69()
{
// Check for InverseLine if line is slightly west of S and that s13 is
// correctly set.
GeodesicLine line =
Geodesic.WGS84.InverseLine(-5, -0.000000000000002, -10, 180);
GeodesicData dir =
line.Position(2e7, GeodesicMask.STANDARD | GeodesicMask.LONG_UNROLL);
AssertEquals(dir.lat2, 4.96445, 0.5e-5);
AssertEquals(dir.lon2, -180.00000, 0.5e-5);
AssertEquals(dir.azi2, -0.00000, 0.5e-5);
dir = line.Position(0.5 * line.Distance(),
GeodesicMask.STANDARD | GeodesicMask.LONG_UNROLL);
AssertEquals(dir.lat2, -87.52461, 0.5e-5);
AssertEquals(dir.lon2, -0.00000, 0.5e-5);
AssertEquals(dir.azi2, -180.00000, 0.5e-5);
}