High order semi-Lagrangian methods for the BGK equation