Method TaylorHoodNavierStokes::jacobian_volume is never called
The class TaylorHoodNavierStokes implements the method jacobian_volume, but the corresponding flag doJacobianVolume is never set. Consequently, the user is tacitly served the numerical approximation of the Jacobian. Is this a bug or a feature?