NBody solver for GPU using cell lists