 
   if(n_approx > 0)
   {
-    #ifdef INDSOFT
-      float temp = node_eps[tid];
-    #endif
+#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];
-        float4 octopole1 = multipole_data[address + 2];
-      #else
-        float4 monopole  = tex1Dfetch(texMultipole, address);
-        float4 octopole0 = tex1Dfetch(texMultipole, address + 1);
-        float4 octopole1 = tex1Dfetch(texMultipole, address + 2);
-      #endif
+#ifndef TEXTURES
+      float4 monopole  = multipole_data[address    ];
+      float4 octopole0 = multipole_data[address + 1];
+      float4 octopole1 = multipole_data[address + 2];
+#else
+      float4 monopole  = tex1Dfetch(texMultipole, address);
+      float4 octopole0 = tex1Dfetch(texMultipole, address + 1);
+      float4 octopole1 = tex1Dfetch(texMultipole, address + 2);
+#endif
 
       node_mon0[tid] = monopole.w;
       node_mon1[tid] = make_float3(monopole.x,  monopole.y,  monopole.z);
       node_oct0[tid] = make_float3(octopole0.x, octopole0.y, octopole0.z);
       node_oct1[tid] = make_float3(octopole1.x, octopole1.y, octopole1.z);
 
-      #ifdef INDSOFT
-        node_eps[tid]  = octopole0.w;
-      #endif
+#ifdef INDSOFT
+      node_eps[tid]  = octopole0.w;
+#endif
 
     } else {
       //Set non-active memory locations to zero

        if((n_stack1 - c_stack0) >= (LMEM_STACK_SIZE << SHIFT))
        {
          //We overwrote our current stack
	  apprCount = -1; return acc_i;	 
	}
      }
      __syncthreads();


      /***
      **** --> copy nodes1 to nodes0: done by reassigning the pointers
      ***/
      n_nodes0    = n_nodes1;

      n_stack_pre = n_stack0;

      acc_i.y += sh_acc[idx].y;
      acc_i.z += sh_acc[idx].z;
      if (sh_ds2[idx] < ds2_min) {
        ds2_min = sh_ds2[idx];
        ngb     = sh_ngb[idx];
      }
    }
  }
