Revisiting the multigrid routines
Description
Revisiting the multigrid routines to have an interpolation scheme that is the trilinear interpolation without applying boundary conditions. This makes the code faster.
The interpolation and restriction now work for dimension 1, 2, and 3.
Simple tests on bulk Si show that the new codes converges much faster for small spacings and can lead to a reduction of the run time up to 30%.
News snippet
Improving the multigrid preconditioner.
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.
Edited by Nicolas Tancogne-Dejean