-// #include "support_kernels.cu"
+// #include "support_kernels.cu0
 #include <stdio.h>
 #include "../profiling/bonsai_timing.h"
 PROF_MODULE(dev_approximate_gravity);

          nodesM[laneId + i] = 0;                                          //but do not touch those parts which has already been filled

#if 0  /* the following gives different result than then one in else */
       /* the results become the same if I uncomment printf above */
        if (flag)
        {
          nodesM[offset] = child; 

            sh_mass[laneId] = posj.w;
            sh_pos [laneId] = make_float3(posj.x, posj.y, posj.z);

            for (int i = 0; i < WARP_SIZE; i++)
              for (int k = 0; k < NI; k++)
                acc_i[k] = add_acc(acc_i[k], pos_i[k], sh_mass[i], sh_pos[i], eps2);
            direCount += WARP_SIZE*NI;
          }

