Bugfix DFT+U, domain parallelization, and MKL
Description
Do not call gemv is there is zero point in the submesh, as this causes MKL to return an error about an invalid argument for LDA.
News snippet
Minor fix for DFT+U with domain parallelization for large systems.
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.