Skip to content

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

Merge request reports