Optim. LDA+U submesh
For DFT+U, the orbitals defined on the submesh are first copied to a local function on the mesh before computing the dot product, as we do for the pseudopotentials. This makes this part of the code faster, as checked using oct-test, for various spacings (from 0.5a.u. to 0.2 a.u.) and values of AOThreshold (from 0.6 to 0.1).
Small optimization of LDA+U
- I have checked that my code follows the Octopus coding standards
- I have added tests for all the new features added in this request.