enovol operator in CarpetLib is broken

The implementation of eno_interpolation_type = averages, which is not used by default, is broken and cannot be used since it does not allow restriction. Even when adding restriction the computation of the required source region is done incorrectly and leads to segfaults since data outside of the source region is accessed. The check for regbbox + stencil width to be included in srcbbox in the prolongation operator source file seems to be broken as well since it does not take the difference between coarse and fine grid spacing into account.

Keyword:

Comments (2)