Structure preserving schemes for Fokker–Planck equations with nonconstant diffusion matrices