Preconditioners for non-orthogonal cells
Description
In the case of non-orthogonal cells, we use the stencil star general to generate the weights for the preconditioners instead of the stencil star.
The weight of the prefilter preconditioner are now computed in a generic way, which also solve a bug with different spacing, as the previous formula assumed identical spacings along each directions.
News snippet
In the case of non-orthogonal cells, we use the stencil star general to generate the weights for the preconditioners instead of the stencil star.
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 Martin Lueders