This onleiner makes my richards example faster by 70%.

This fixes a code generation performance bug in sum factorized code.
It also reveals that statistic profiling could be superhelpful to us.

We can disable sum factorization on boundary integrals.
This commit allows us to use a different mixin for those integrals.

The realize_output interface method was not really necessary and just adds a
layer of indiretion.