                                                            rigid_data.conjqm,
                                                            n_group_bodies, 
                                                            n_bodies, 
-                                                           local_beg,
                                                            box, 
                                                            deltaT);
     

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

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