From a2efc211eb5fa79ce3c6e666e83c3e65afd22e46 Mon Sep 17 00:00:00 2001 From: Charles Karney Date: Sun, 22 Sep 2019 16:44:57 -0400 Subject: Update to version 1.50 of the geodesic library. * Allow arbitrarily complex polygons in geod_polygon_*. In the case of self-intersecting polygons the area is accumulated "algebraically", e.g., the areas of the 2 loops in a figure-8 polygon will partially cancel. * Simplify code by using C99 functions remainder and remquo. * More test coverage. Fixes to associated files: * src/pipeline.cpp invoke geod_init with f = es / (1 + sqrt(1 - es)) instead of (the less accurate) f = 1 - sqrt(1 - es) * src/apps/geod_set.cpp remove "#undef f" (a dangling relic?). --- src/apps/geod_set.cpp | 1 - 1 file changed, 1 deletion(-) (limited to 'src/apps/geod_set.cpp') diff --git a/src/apps/geod_set.cpp b/src/apps/geod_set.cpp index cd6b5018..ed7edeb9 100644 --- a/src/apps/geod_set.cpp +++ b/src/apps/geod_set.cpp @@ -46,7 +46,6 @@ geod_set(int argc, char **argv) { /* check if line or arc mode */ if (pj_param(nullptr,start, "tlat_1").i) { double del_S; -#undef f phi1 = pj_param(nullptr,start, "rlat_1").f; lam1 = pj_param(nullptr,start, "rlon_1").f; if (pj_param(nullptr,start, "tlat_2").i) { -- cgit v1.2.3