                                                          float4* rdata_conjqm, 
                                                          unsigned int n_group_bodies,
                                                          unsigned int n_bodies, 
-                                                         unsigned int local_beg,
                                                          unsigned int nmax,
                                                          gpu_boxsize box, 
                                                          float deltaT)
     {
-    unsigned int group_idx = blockIdx.x * blockDim.x + threadIdx.x + local_beg;
+    unsigned int group_idx = blockIdx.x * blockDim.x + threadIdx.x;
     
     if (group_idx < n_group_bodies)
         {

    \param rdata_torque Body torques
    \param n_group_bodies Number of rigid bodies in my group
    \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 box Box dimensions for periodic boundary condition handling
*/

                                                 float4* rdata_torque, 
                                                 unsigned int n_group_bodies,
                                                 unsigned int n_bodies, 
                                                 unsigned int nmax,
                                                 gpu_boxsize box)
    {
    unsigned int group_idx = blockIdx.x;
    
    __shared__ float4 ex_space, ey_space, ez_space;
    float4 *body_force = sum;
