Merge topic 'BSplineTransformIndexDecrement'

f27b680 BUG: Adjust index decrement in BSplineTransform
Kitware · May 25, 2016 · 09c443d · 09c443d
2 parents f70b7cb + f27b680
commit 09c443d
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/Modules/Core/Transform/include/itkBSplineTransform.hxx b/Modules/Core/Transform/include/itkBSplineTransform.hxx
@@ -474,11 +474,11 @@ BSplineTransform<TParametersValueType, NDimensions, VSplineOrder>
   bool inside = true;
   for( unsigned int j = 0; j < SpaceDimension; j++ )
     {
-    ScalarType maxLimit = static_cast<ScalarType>( gridSize[j] ) - 0.5
+    const ScalarType maxLimit = static_cast<ScalarType>( gridSize[j] ) - 0.5
       * static_cast<ScalarType>( SplineOrder - 1 ) - 1.0;
-    if(Math::AlmostEquals( index[j], maxLimit ))
+    if( Math::FloatAlmostEqual( index[j], maxLimit, 4 ) )
       {
-      index[j] = Math::FloatAddULP( index[j], -6 );
+      index[j] = Math::FloatAddULP( maxLimit, -6 );
       }
     else if( index[j] >= maxLimit )
       {