Top panel: dispersion relation of the numerical scheme. Red points are the measured oscillation frequencies normalised by as a function of kl_H with l_H = η_H/c_A. Solid black lines represent the theoretical dispersion relation, while blue points show the relative error of each measure. Bottom panel: damping rate of whistler waves as a function of resolution, compared with Δx and Δx² scaling.