379 CU_ASSERT_DOUBLE_EQUAL(v_out_top.
lat, g_out_top.
lat, 0.000001);
380 CU_ASSERT_DOUBLE_EQUAL(v_out_top.
lon, g_out_top.
lon, 0.000001);
381 CU_ASSERT_DOUBLE_EQUAL(v_out_bottom.
lat, g_out_bottom.
lat, 0.000001);
382 CU_ASSERT_DOUBLE_EQUAL(v_out_bottom.
lon, g_out_bottom.
lon, 0.000001);
384 gs.
lat = 1.3021240033804449;
385 ge.
lat = 1.3021240033804449;
386 gs.
lon = -1.3387392931438733;
387 ge.
lon = 1.80285336044592;
395 CU_ASSERT_DOUBLE_EQUAL(v_out_top.
lat, g_out_top.
lat, 0.000001);
396 CU_ASSERT_DOUBLE_EQUAL(v_out_top.
lon, g_out_top.
lon, 0.000001);
397 CU_ASSERT_DOUBLE_EQUAL(v_out_bottom.
lat, g_out_bottom.
lat, 0.000001);
398 CU_ASSERT_DOUBLE_EQUAL(v_out_bottom.
lon, g_out_bottom.
lon, 0.000001);
static void point_set(double lon, double lat, GEOGRAPHIC_POINT *p)
int clairaut_geographic(const GEOGRAPHIC_POINT *start, const GEOGRAPHIC_POINT *end, GEOGRAPHIC_POINT *g_top, GEOGRAPHIC_POINT *g_bottom)
Computes the pole of the great circle disk which is the intersection of the great circle with the lin...
int clairaut_cartesian(const POINT3D *start, const POINT3D *end, GEOGRAPHIC_POINT *g_top, GEOGRAPHIC_POINT *g_bottom)
Computes the pole of the great circle disk which is the intersection of the great circle with the lin...
void geog2cart(const GEOGRAPHIC_POINT *g, POINT3D *p)
Convert spherical coordinates to cartesian coordinates on unit sphere.
Point in spherical coordinates on the world.