                                                                                              rigid_data.torque,
                                                                                              n_group_bodies, 
                                                                                              n_bodies, 
-                                                                                             local_beg,
                                                                                              nmax,
                                                                                              window_size,
                                                                                              box);       

                                                                     rigid_data.virial,
                                                                     n_group_bodies,
                                                                     n_bodies, 
                                                                     local_beg,
                                                                     box, 
                                                                     deltaT);
        }

    \param rdata_conjqm Conjugate quaternion momentum
    \param n_group_bodies Number of rigid bodies in my group
    \param n_bodies Total number of rigid bodies
    \param deltaT Timestep 
    \param box Box dimensions for periodic boundary condition handling
*/
