 
   /*********** shared memory distribution **********/
 
-                                                  //  begin,    end,   size
-                                                  // -----------------------
+  //  begin,    end,   size
+  // -----------------------
   int *approx = (int*)&shmem [     0];            //  0*DIM,  2*DIM,  2*DIM
   int *direct = (int*)&approx[ 2*DIM];            //  2*DIM,  3*DIM,  1*DIM
   int *nodes  = (int*)&direct[   DIM];            //  3*DIM, 13*DIM, 10*DIM


  if(n_approx > 0)
  {
#ifdef INDSOFT
    float temp = node_eps[tid];
#endif

    if (tid < n_approx) {
      int address      = (approx[tid] << 1) + approx[tid];
#ifndef TEXTURES
      float4 monopole  = multipole_data[address    ];
      float4 octopole0 = multipole_data[address + 1];

    }
  }
  __syncthreads();
#endif

  return acc_i;
}
