Batchifying the calculation of the commutator of the position operator and the DFT+U(+V) term. This allows to compute the current fully batchified.