To elaborate, Seems like Q elements work according to the following lines:
parametrize_cell_types_tp
@pytest.mark.parametrize("family", ["Q"])
@pytest.mark.parametrize("degree", [2, 3, 4])
def test_P_tp(family, degree, cell_type, datadir):
mesh = get_mesh(cell_type, datadir)
V = FunctionSpace(mesh, (family, degree))
run_scalar_test(mesh, V, degree)