use sqrt for older gcc

Saurav Agarwal · Saurav Agarwal · commit 97af831c473d · 2025-05-25T00:10:50.000-04:00
diff --git a/cppsrc/core/include/CoverageControl/constants.h b/cppsrc/core/include/CoverageControl/constants.h
@@ -48,9 +48,9 @@ namespace CoverageControl {
 double const kEps = 1e-10;          /*!< Epsilon for double comparison */
 float const kEpsf = 1e-6f;          /*!< Epsilon for float comparison */
 double const kLargeEps = 1e-4;      /*!< Large epsilon for double comparison */
-double const kSqrt2 = std::sqrt(2); /*!< Square root of 2 */
-double const kOneBySqrt2 = 1. / std::sqrt(2); /*!< 1 by square root of 2 */
-float const kOneBySqrt2f = 1.f / sqrtf(2.f);  /*!< 1 by square root of 2 */
+double const kSqrt2 = std::sqrt(2.); /*!< Square root of 2 */
+double const kOneBySqrt2 = 1.0 / kSqrt2; /*!< 1 by square root of 2 */
+float const kOneBySqrt2f = 1.0f / std::sqrt(2.0f);  /*!< 1 by square root of 2 */
 double const kInfD =
     std::numeric_limits<double>::infinity(); /*!< Infinity for double */
 constexpr auto kMaxPrecision{std::numeric_limits<long double>::digits10 +
diff --git a/cppsrc/core/src/world_idf.cpp b/cppsrc/core/src/world_idf.cpp
@@ -119,7 +119,7 @@ void WorldIDF::GenerateMapCuda(float const resolution, float const truncation,
     host_dists[i].sigma_y = static_cast<float>(sigma.y());
     float rho = static_cast<float>(normal_distributions_[i].GetRho());
     host_dists[i].rho = rho;
-    host_dists[i].sqrt_one_minus_rho_squared = std::sqrtf(1. - rho * rho);
+    host_dists[i].sqrt_one_minus_rho_squared = std::sqrt(1. - rho * rho);
     host_dists[i].scale =
         static_cast<float>(normal_distributions_[i].GetScale());
   }

Original file line number	Diff line number	Diff line change
`@@ -119,7 +119,7 @@ void WorldIDF::GenerateMapCuda(float const resolution, float const truncation,`
`119`	`119`	`host_dists[i].sigma_y = static_cast<float>(sigma.y());`
`120`	`120`	`float rho = static_cast<float>(normal_distributions_[i].GetRho());`
`121`	`121`	`host_dists[i].rho = rho;`
`122`		`- host_dists[i].sqrt_one_minus_rho_squared = std::sqrtf(1. - rho * rho);`
	`122`	`+ host_dists[i].sqrt_one_minus_rho_squared = std::sqrt(1. - rho * rho);`
`123`	`123`	`host_dists[i].scale =`
`124`	`124`	`static_cast<float>(normal_distributions_[i].GetScale());`
`125`	`125`	`}`