     /*!
      * \param q The current thread data.
      * \param primitive Index of the intersected primitive.
+     * \param image the current image
      * \returns True If the volumes still overlap after refinement.
      *
      * HOOMD's neighbor lists require additional filtering. This first ensures

     * enabled, the cutoff radius is adjusted based on the diameters of the
     * particles.
     */
    DEVICE bool refine(const ThreadData& q, const int primitive) const
        {
        bool exclude = (q.idx == primitive);


            rc2 *= rc2;

            // compute distance and wrap back into box
            const Scalar3 dr = r - q.position - image;
            const Scalar drsq = dot(dr,dr);

            // exclude if outside the sphere
