X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=itk%2FitkBSplineInterpolateImageFunctionWithLUT.txx.original;h=877a78494bf33abf6359e2693fa09f5335dd5945;hb=8d51bd1cb7e9ca416b46677dcd73ba76d107a587;hp=3a8744ff6999a55903fabfa547531ab40d6f39a9;hpb=71511076230b1d0a60a638c2c1b6577a34116e02;p=clitk.git diff --git a/itk/itkBSplineInterpolateImageFunctionWithLUT.txx.original b/itk/itkBSplineInterpolateImageFunctionWithLUT.txx.original index 3a8744f..877a784 100644 --- 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; lm_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; } }