Use std::sqrt() instead of std::pow()

kimwalisch · Mar 11, 2024 · c577322 · c577322
1 parent 9ed8e4c
commit c577322
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/include/imath.hpp b/include/imath.hpp
@@ -150,6 +150,8 @@ inline T iroot(T x)
 
   if (N == 3)
     r = (T) std::cbrt((double) x);
+  else if (N == 4)
+    r = (T) std::sqrt(std::sqrt((double) x));
   else
     r = (T) std::pow((double) x, 1.0 / N);