// compute the interpolation indexes
for (unsigned int n = 0; n< ImageDimension; n++) {
if (splineOrder & 1) { // Use this index calculation for odd splineOrder
- indx = (long)vcl_floor((float)x[n]) - splineOrder / 2;
+ indx = (long)std::floor((float)x[n]) - splineOrder / 2;
for (unsigned int k = 0; k <= splineOrder; k++) {
evaluateIndex[n][k] = indx++;
}
} else { // Use this index calculation for even splineOrder
- indx = (long)vcl_floor((float)(x[n] + 0.5)) - splineOrder / 2;
+ indx = (long)std::floor((float)(x[n] + 0.5)) - splineOrder / 2;
for (unsigned int k = 0; k <= splineOrder; k++) {
evaluateIndex[n][k] = indx++;
}