FEniCS / ffcx / 38575778870

88%

push

github

Fix sum factorisation with several or one point quadrature rules (#902)

- 1D tensor factor tables were named FE_TF0, FE_TF1, ... from zero for
  each quadrature rule, so the tables of integrals with different
  quadrature degrees clashed and one rule read the other's tables, giving
  a wrong element tensor. Include the quadrature rule in their name.
- A piecewise table with tensor factors, such as a coefficient on a one
  point rule, is accessed outside the quadrature loop, where the
  quadrature index can't be split by direction, which failed an
  assertion. Its factors are constant over the points, so read point 0.

Co-authored-by: Claude Opus 5.5 (1M context) <noreply@anthropic.com>

1 of 1 new or added line in 1 file covered. (100.0%)

4779 of 5441 relevant lines covered (87.83%)

0.88 hits per line