     
     unsigned int n_windows = nmax / window_size + 1;
     
+    // slide the window throughout the block
     for (unsigned int start = 0; start < n_windows; start++)
     {
         unsigned int iwindow = threadIdx.x / window_size;

//! Shared memory for body force and torque reduction, required allocation when the kernel is called
extern __shared__ float4 sum[];

//! Takes the 2nd 1/2 step forward in the velocity-verlet NVE integration scheme
/*! \param rdata_force Particle data to step forward in time
    \param rdata_torque List of pointers to forces on each particle
    \param n_bodies Number of forces listed in \a force_data_ptrs

            body_force[ilocal].z = fi.z;
            body_force[ilocal].w = fi.w;
            
            ri.x = pos.x - com.x;
            ri.y = pos.y - com.y;
            ri.z = pos.z - com.z;
