Add support for bubble enhanced P2 element in 2D
The conforming Crouzeix-Raviart element P_2^+-P_1^\text{disc}
for triangles in 2D employs for the velocity component continuous quadratic Lagrange elements extended by a cubic bubble function. This MR adds to HOG support for this P_2^+
part.