Fix a bug in the calculation of the cutoff for the 3D-0D FFT in non-cubic cells.
Description
Fix a bug in the calculation of the cutoff for the 3D-0D FFT in non-cubic cells.
This allows to make the FFT method the default for computing the Coulomb integrals used in DFT+U and DFT+U+V. The FFT can also now be used in domain parallelization. This makes the calculation of the Coulomb integrals much faster.
News snippet
Coulomb integrals for DFT+U and DFT+U+V are now computed using FFTs by default.
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.