High-order diamond differencing schemes for the Boltzmann Fokker–Planck equation in 1D and 2D Cartesian geometries