Description

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

News snippet

Improved support of GPUs for DFT+U calculations.

Checklist

I have checked that my code follows the Octopus coding standards
I have added tests for all the new features added in this request.