[localevaluation.elliptic] enable matrixvalued diffusion
[local{functional,operator] use correct quadrature order
* Joe is to blame if this is wrong