 
 //! Maximum number of threads (width of a dpd_warp)
 const unsigned int gpu_aniso_pair_force_max_tpp = 32;
-
-//! CTA reduce
-template<typename T>
-__device__ static T aniso_warp_reduce(unsigned int NT, int tid, T x, volatile T* shared)
-    {
-    shared[tid] = x;
-
-    for (int dest_count = NT/2; dest_count >= 1; dest_count /= 2)
-        {
-        if (tid < dest_count)
-           {
-            x += shared[dest_count + tid];
-            shared[tid] = x;
-            }
-         }
-    T total = shared[0];
-    return total;
-    }
-
 
 //! Wraps arguments to gpu_cgpf
 struct a_pair_args_t

                        grid.x = 65535;
                        }

                    shared_bytes += sizeof(Scalar)*block_size;

                    gpu_compute_pair_aniso_forces_kernel<evaluator, 0, 0>
                        <<<grid, block_size, shared_bytes>>>(pair_args.d_force,
                                                      pair_args.d_torque,

                        grid.y = grid.x/65535 + 1;
                        grid.x = 65535;
                        }

                    gpu_compute_pair_aniso_forces_kernel<evaluator, 0, 1>
                        <<<grid, block_size, shared_bytes>>>(pair_args.d_force,
