Single-cone real-space finite difference schemes for the Dirac von Neumann equation