Fixed StokesStudyForLeszek so that now my bad numbers match his bad numbers.
Added "enhanceFluxes" boolean argument to calls to VGP Navier-Stokes formulation constructors.  False is the right thing, at least on the present implementation.  May axe the argument in the end...
Added generic implementation of regularRefinementPattern that takes a cell topology key as argument.
Added commented-out lines for printing local stiffness matrices (for debugging, and comparing stiffness matrices with Leszek).
Added option for including an algebraic pressure equation.
Changes to BasisReconciliation from yesterday (omitted in earlier commit).