                                                                      nmax,
                                                                      box, 
                                                                      deltaT);
+                                                        
         }
     
     return cudaSuccess;

        dim3 particle_threads(block_size, 1, 1);
        
        gpu_nve_rigid_step_one_particle_sliding_kernel<<< particle_grid, particle_threads >>>(pdata.pos, 
                                                                     pdata.vel, 
                                                                     pdata.image,
                                                                     rigid_data.particle_oldpos,
                                                                     rigid_data.particle_oldvel,
                                                                     rigid_data.virial,

extern "C" __global__ void gpu_nve_rigid_step_one_particle_sliding_kernel(float4* pdata_pos,
                                                        float4* pdata_vel,
                                                        int4* pdata_image,
                                                        float* d_particle_mass,
                                                        float4* d_net_force,
                                                        unsigned int *d_rigid_group,
                                                        float4* d_rigid_orientation,
                                                        float4* d_rigid_com,
                                                        float4* d_rigid_vel,
                                                        float4* d_rigid_angvel,
                                                        int* d_rigid_imagex,
                                                        int* d_rigid_imagey,
                                                        int* d_rigid_imagez,
                                                        unsigned int* d_rigid_particle_idx,
                                                        float4* d_rigid_particle_dis,
                                                        float4* rdata_oldpos,
                                                        float4* rdata_oldvel,
                                                        float *d_virial,
