Support #592: Minor refactoring of RefFEltInLocalNumbering. - Bug #589: Clarify which dimension should be used in intermediate matrices of local variational operators. Now the heat may work in 2D. - Bug #586: Correct a bug in grad-grad operator: work matrix allocation assumed a vectorial unknown. - Feature #185: Shape function label is no longer glued to the unknown in general, but to an unknown in a finite element space. So you might have displacement/P1 in one finite element space and displacement/P2 in another, which wasn't possible previously (or you would have had to cheat and name differently the unknown). - Feature #583 - #26: Heavy refactoring of quadrature rules management; now each operator may use its own quadrature rules. - Feature #579 - #588: Introduce the new policy AtQuadraturePoint for parameters; add the possibility to evolve with time (not for the whole Parameter class, just for Parameters that inherit statically from this new policy). - Feature #578: Parameter: add a Write() functionality. Content of this Write() is not entirely satisfactorily though: composite parameters just specify the values of their components, and Lua functions are not written inside (see #581). - Support #576: Few minor corrections suggested by cppcheck.