Normalized Sinc function. sinc(x) = sin(pi*x)/(pi*x).
x
double sinc(double x) { double z = constants.pi * x; return almostEqualI(z, 0.0, 15) ? 1.0 : sin(z) / z; }