Make numerical stability fix actually be present

git-svn-id: http://svn.osgeo.org/postgis/trunk@8666 b70326c6-7e19-0410-871a-916f4a2858ee

diff --git a/liblwgeom/lwspheroid.c b/liblwgeom/lwspheroid.c
index 5f11555094cf4d4d15323258a0a583fa08fcb124..32c0ac047765fbf26a235f3602e5d125febc3380 100644 (file)
--- a/liblwgeom/lwspheroid.c
+++ b/liblwgeom/lwspheroid.c
@@ -183,7 +183,6 @@ double spheroid_direction(const GEOGRAPHIC_POINT *r, const GEOGRAPHIC_POINT *s,
                else
                        alpha = asin(sin_alpha);
 
-               alpha = asin(sin_alpha);
                cos_alphasq = POW2(cos(alpha));
                cos2_sigma_m = cos(sigma) - (2.0 * sin_u1 * sin_u2 / cos_alphasq);