Hierarchical N-body simulations with auto-tuning for heterogeneous systems