   int *nodes  = (int*)&direct[   DIM];            //  3*DIM, 13*DIM, 10*DIM
   int *prefix = (int*)&nodes [10*DIM];            // 13*DIM, 15*DIM,  2*DIM
 
+#ifndef OLDPREFIX
+  __shared__ int prefix_shmem[32];
+#endif
+
   float  *node_mon0 = (float* )&nodes    [DIM];   //  4*DIM,  5*DIM,  1*DIM
   float3 *node_mon1 = (float3*)&node_mon0[DIM];   //  5*DIM,  8*DIM,  3*DIM
 #if 0



#define TEXTURES
// #define OLDPREFIX


template<int DIM2, int SHIFT>


#define WARP_SIZE2 5
#define WARP_SIZE  (1 << WARP_SIZE2)

#if 1

