 //! The developer has chosen not to document this variable
 __device__ __constant__ float GPU_rho_coeff[CONSTANT_SIZE];
 
+//! Coefficients of a power expansion of sin(x)/x
+const Scalar sinc_coeff[] = {Scalar(1.0), Scalar(-1.0/6.0), Scalar(1.0/120.0),
+                        Scalar(-1.0/5040.0),Scalar(1.0/362880.0),
+                        Scalar(-1.0/39916800.0)};
+
+//! Coefficients of a power expansion of sin(x)/x
+__device__ __constant__ float gpu_sinc_coeff[6];
+
 /*! \file PPPMForceGPU.cu
   \brief Defines GPU kernel code for calculating the Fourier space forces for the Coulomb interaction. Used by PPPMForceComputeGPU.
 */


//! Texture for reading charge parameters
texture<float, 1, cudaReadModeElementType> pdata_charge_tex;

//! Implements workaround atomic float addition on sm_1x hardware
__device__ inline void atomicFloatAdd(float* address, float value)

        kvec_array[idx] = k;

        float sqk = dot(k,k);
        // omit DC term
        if(idx == 0) {
            vg[0+6*idx] = 0.0f;
            vg[1+6*idx] = 0.0f;
            vg[2+6*idx] = 0.0f;
