geocentric_radius(): save a sqrt()

rouault · Mar 16, 2024 · 2d22cf6 · 2d22cf6
1 parent 36da998
commit 2d22cf6
Showing 1 changed file with 3 additions and 4 deletions.
diff --git a/src/conversions/cart.cpp b/src/conversions/cart.cpp
@@ -130,10 +130,9 @@ static double geocentric_radius(double a, double b_div_a, double cosphi,
     const double b_div_a_squared = b_div_a * b_div_a;
     const double b_div_a_squared_mul_sinphi_squared =
         b_div_a_squared * sinphi_squared;
-    return a *
-           sqrt(cosphi_squared +
-                b_div_a_squared * b_div_a_squared_mul_sinphi_squared) /
-           sqrt(cosphi_squared + b_div_a_squared_mul_sinphi_squared);
+    return a * sqrt((cosphi_squared +
+                     b_div_a_squared * b_div_a_squared_mul_sinphi_squared) /
+                    (cosphi_squared + b_div_a_squared_mul_sinphi_squared));
 }
 
 /*********************************************************************/