                                                                      n_bodies, 
                                                                      local_beg,
                                                                      nmax,
-                                                                     block_size,
                                                                      box, 
                                                                      deltaT);
         }

    \param n_bodies Total number of rigid bodies
    \param local_beg Starting body index in this card
    \param nmax Maximum number of particles in a rigid body
    \param block_size Block size
    \param box Box dimensions for periodic boundary condition handling
    \param deltaT Time step
*/

                                                        unsigned int n_bodies, 
                                                        unsigned int local_beg,
                                                        unsigned int nmax,
                                                        gpu_boxsize box,
                                                        float deltaT)
    {
