diff --git a/GeoLib/IO/XmlIO/Qt/XmlStnInterface.cpp b/GeoLib/IO/XmlIO/Qt/XmlStnInterface.cpp
index 465c4cf1fa0e0700a8bc4b8c365aa4527948f279..df88156fb89f290b38c4995a6d32d0480f1a81c9 100644
--- a/GeoLib/IO/XmlIO/Qt/XmlStnInterface.cpp
+++ b/GeoLib/IO/XmlIO/Qt/XmlStnInterface.cpp
@@ -208,7 +208,9 @@ void XmlStnInterface::readStratigraphy( const QDomNode &stratRoot,
             /* add other horizon features here */
 
             double depth (horizon.attribute("z").toDouble());
-            if (fabs(depth - depth_check) > std::numeric_limits<double>::epsilon()) // skip soil-layer if its thickness is zero
+            if (std::abs(depth - depth_check) >
+                std::numeric_limits<double>::
+                    epsilon())  // skip soil-layer if its thickness is zero
             {
                 borehole->addSoilLayer(horizon.attribute("x").toDouble(),
                                        horizon.attribute("y").toDouble(),
@@ -326,7 +328,7 @@ void XmlStnInterface::writeBoreholeData(QDomDocument &doc,
     boreholeTag.appendChild(stationDepthTag);
     QDomText stationDepthText = doc.createTextNode(QString::number(borehole->getDepth(), 'f'));
     stationDepthTag.appendChild(stationDepthText);
-    if (fabs(borehole->getDate()) > 0)
+    if (std::abs(borehole->getDate()) > 0)
     {
         QDomElement stationDateTag = doc.createElement("bdate");
         boreholeTag.appendChild(stationDateTag);
diff --git a/GeoLib/Raster.cpp b/GeoLib/Raster.cpp
index 7336958d5def40d875c52278044664ae7256d3fd..48623dd40387a6842e21d5b7fd178e290a7af35d 100644
--- a/GeoLib/Raster.cpp
+++ b/GeoLib/Raster.cpp
@@ -88,8 +88,8 @@ double Raster::interpolateValueAtPoint(MathLib::Point3d const& pnt) const
     double const yIdx (std::floor(yPos));    //  so not to over- or underflow.
 
     // weights for bilinear interpolation
-    double const xShift = std::fabs((xPos - xIdx) - 0.5);
-    double const yShift = std::fabs((yPos - yIdx) - 0.5);
+    double const xShift = std::abs((xPos - xIdx) - 0.5);
+    double const yShift = std::abs((yPos - yIdx) - 0.5);
     std::array<double,4> weight = {{ (1-xShift)*(1-yShift), xShift*(1-yShift), xShift*yShift, (1-xShift)*yShift }};
 
     // neighbors to include in interpolation
@@ -119,7 +119,8 @@ double Raster::interpolateValueAtPoint(MathLib::Point3d const& pnt) const
         }
 
         // remove no data values
-        if (std::fabs(pix_val[j] - _header.no_data) < std::numeric_limits<double>::epsilon())
+        if (std::abs(pix_val[j] - _header.no_data) <
+            std::numeric_limits<double>::epsilon())
         {
             weight[j] = 0;
             no_data_count++;