Problem with is_affine() method of blending maps
I am summing up the result of a discussion in our developer chat here to document this problem. Thanks to Nils for contributing the main result.
While adding an AffineMap2D
blending map, see [fdc929e9], it turned out that there is a problem with the is_affine()
method. The latter is supposed to be used to indicate that the map does not have to be computed on each quadrature point. However, the implementation of that optimization seems to be buggy.
Thus, AffineMap2D
violates the POLA idea and currently returns False
when is_affine()
is called.