                         grid.x = 65535;
                         }
 
-                    shared_bytes += sizeof(Scalar)*block_size;
-
                     gpu_compute_pair_aniso_forces_kernel<evaluator, 1, 0>
                         <<<grid, block_size, shared_bytes>>>(pair_args.d_force,
                                                       pair_args.d_torque,


//! 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;
                        }

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