Privately someone pointed out a couple of mistakes in my
suggested fix. The corrected lines read:
if (interpolationType != NO_INTERPOLATION)
{ Ex -= (0.5 / (DBL)Data->Sx); // Result to voxel center
Ey -= (0.5 / (DBL)Data->Sy);
Ez -= (0.5 / (DBL)Data->Sz);
}
Bill P.