NonlinearDiffusionNewtonGalerkin
operator.Bases the operator generation on the new pystencils backend (to be released as version 2.0) and gets rid of our own fork. The new backend is significantly more efficient, less prone to bugs, and easily extensible. Vectorized code is also properly typed, which means that we do not need any stunts in our blending kernels to compute the quadrature points from the loop counters.
Includes some additional quality improvements:
HOGType
class because that was nothing more than boilerplate code.Some (mostly) unused features are removed with the idea that they will be reimplemented more generally in the future:
FUSEDROWS
: To be reimplemented as a special case of general blocking (much like CUBES
).reduce_reuse_dist
and elimtmps_and_reorder
.