bilinear finite elements
[operators] set default quadrature order for L2ProductFunctionalQ1 to 2
[operators] choose right quadrature in L2ProductQ1
[operators] remove 'order' argument from DiffusionOperatorQ1's __init__ method