11 # pragma warning (disable: 4996)
22 double a = 6378137, f = 1/298.257223563;
23 double lat1, lon1, azi1, lat2, lon2, azi2, s12;
27 while (scanf(
"%lf %lf %lf %lf", &lat1, &lon1, &azi1, &s12) == 4) {
28 geod_direct(&g, lat1, lon1, azi1, s12, &lat2, &lon2, &azi2);
29 printf(
"%.15f %.15f %.15f\n", lat2, lon2, azi2);
API for the geodesic routines in C.
void GEOD_DLL geod_init(struct geod_geodesic *g, double a, double f)
void GEOD_DLL geod_direct(const struct geod_geodesic *g, double lat1, double lon1, double azi1, double s12, double *plat2, double *plon2, double *pazi2)