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


//! Maximum number of threads (width of a dpd_warp)
const unsigned int gpu_aniso_pair_force_max_tpp = 32;

//! 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,
