Initialization of FEEvaluation like FEValues

A new way to initialize an FEEvaluation object (fast tensor product evaluation routines) has been introduced. It takes the same arguments as the initialization method of FEValues, namely a Mapping, a FiniteElement, a one-dimensional Quadrature<1> object, and update flags. This enables using FEEvaluation on the fly just as FEValues. Since FEEvaluation treats objects of one base element only, this class can select which components to read.

The concept is explained in the documentation of the FEEvaluation class (around line 1376 in fe_evaluation.h) and an example for the Laplace matrix is given. However, it also works for the Stokes operator as shown by the test program assemble_matrix_02.
I see. I wasn't able to read this from the text above the code snippet, but I also don't quite know the background of it all so it may have been clear to someone who knows it better than I do.