Remove vcl_math calls

[clitk.git] / itk / itkBSplineInterpolateImageFunctionWithLUT.txx.original

diff --git a/itk/itkBSplineInterpolateImageFunctionWithLUT.txx.original b/itk/itkBSplineInterpolateImageFunctionWithLUT.txx.original
index 3a8744ff6999a55903fabfa547531ab40d6f39a9..877a78494bf33abf6359e2693fa09f5335dd5945 100644 (file)
--- a/itk/itkBSplineInterpolateImageFunctionWithLUT.txx.original
+++ b/itk/itkBSplineInterpolateImageFunctionWithLUT.txx.original
@@ -259,7 +259,7 @@ GetSampleIndexOfPixelPosition(const ContinuousIndexType & x, IndexType & Evaluat
   // bool mChange = false;
 
   // Compute t1 = distance to floor 
-  TCoefficientType t1 = x[l]- vcl_floor(x[l]);
+  TCoefficientType t1 = x[l]- std::floor(x[l]);
     
   // Compute index in precomputed weights table
   TCoefficientType t2 = mSamplingFactors[l]*t1;
@@ -311,13 +311,13 @@ EvaluateAtContinuousIndex(const ContinuousIndexType & x) const {
   long indx;
   for (unsigned int l=0; l<d; l++)  {
     if (mSplineOrders[l] & 1) {  // Use this index calculation for odd splineOrder (like cubic)
-      indx = (long)vcl_floor(x[l]) - mSplineOrders[l] / 2 ; //this->m_SplineOrder / 2;
+      indx = (long)std::floor(x[l]) - mSplineOrders[l] / 2 ; //this->m_SplineOrder / 2;
       EvaluateIndex[l] = indx;
     }
     else { // Use this index calculation for even splineOrder
       if (mSplineOrders[l] == 0) EvaluateIndex[l] = (long)rint(x[l]);
       else {
-       indx = (long)vcl_floor((x[l]+ 0.5)) - mSplineOrders[l] / 2; //this->m_SplineOrder / 2;
+       indx = (long)std::floor((x[l]+ 0.5)) - mSplineOrders[l] / 2; //this->m_SplineOrder / 2;
        EvaluateIndex[l] = indx;
       }
     }