     gpu_nve_rigid_step_two_body_kernel<<< body_grid, body_threads >>>(rigid_data.vel, 
                                                                       rigid_data.angmom, 
                                                                       rigid_data.angvel,
+                                                                      rigid_data.orientation,
                                                                       rigid_data.conjqm,
+                                                                      rigid_data.body_mass,
+                                                                      rigid_data.moment_inertia,
+                                                                      rigid_data.force,
+                                                                      rigid_data.torque,
+                                                                      rigid_data.body_indices,
                                                                       n_group_bodies,
                                                                       n_bodies, 
                                                                       nmax, 
                                                                       box, 
                                                                       deltaT);
+    
     block_size = 192;
     dim3 particle_grid(group_size/block_size+1, 1, 1);
     dim3 particle_threads(block_size, 1, 1);

                                                           rigid_data.body_imagey, 
                                                           rigid_data.body_imagez,
                                                           rigid_data.conjqm,
                                                           n_group_bodies, 
                                                           n_bodies, 
                                                           box,


//! Flag for invalid particle index, identical to the sentinel value NO_INDEX in RigidData.h
#define INVALID_INDEX 0xffffffff 

/*! Takes the first half-step forward for rigid bodies in the velocity-verlet NVE integration
    \param rdata_com Body center of mass
