Emergency fix to Vincenty's formula at antipodal points

Wetitpig · Sep 5, 2020 · d895e73 · d895e73
1 parent 5d1d33d
commit d895e73
Showing 1 changed file with 0 additions and 4 deletions.
diff --git a/src/vincenty.c b/src/vincenty.c
@@ -190,10 +190,6 @@ void vincenty_inverse(struct Coordinates *location, struct Coordinates *location
 	}
 
 	a = helmertA(calp);
-
-	ssig = hypotl(cu2 * sinl(lambda), cu1 * su2 - su1 * cu2 * cosl(lambda));
-	csig = su1 * su2 + cu1 * cu2 * cosl(lambda);
-	sig = atan2_modified(ssig, csig);
 	b = helmertB(calp, sig);
 
 	s = RAD_MIN * a * (sig + b);