Optim. LDA+U submesh
Description
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).
News snippet
Small optimization of LDA+U
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.