change pow(x,0.5) to std::sqrt

Signed-off-by: Anatoly <Titov_AI@pnpi.nrcki.ru>

diff --git a/src/spirals.cpp b/src/spirals.cpp
index 8f8bdb3c99be7d356dfd4b81fcf9125057381e31..7d0afdf40b16c63eb702630e0aa9111a20a70aa8 100644 (file)
--- a/src/spirals.cpp
+++ b/src/spirals.cpp
@@ -54,10 +54,10 @@ double Fx (double x0, double y0, double z0, double p1, double p2, double p3, std
     double ret = 0;
     for (int i = 0; i < x.size(); i++) {
         ret +=
-           pow (pow (p2 * (x[i][2] - z0) - p3 * (x[i][1] - y0), 2) +
+           std::sqrt(pow (p2 * (x[i][2] - z0) - p3 * (x[i][1] - y0), 2) +
                 pow (p3 * (x[i][0] - x0) - p1 * (x[i][2] - z0), 2) +
-                pow (p1 * (x[i][1] - y0) - p2 * (x[i][0] - x0), 2), 0.5) /
-           pow (p1 * p1 + p2 * p2 + p3 * p3, 0.5);
+                pow (p1 * (x[i][1] - y0) - p2 * (x[i][0] - x0), 2)) /
+           std::sqrt(p1 * p1 + p2 * p2 + p3 * p3);
     }
     return ret;
 }