 
         int nl_pre = 0;                                                     //Number of leaves that have already been processed
 
-#define NJMAX (DIM*4)
         while (n_bodies > 0) 
         {
           int nb    = min(n_bodies, NJMAX - n_direct);                    //Make sure number of bides to be extracted does not exceed

  return y;
}


#endif

/*********** Forces *************/


          }

        }
        directM[tid] = body_list[tid];
        __syncthreads();
#endif
      } //end lvl
