We run ut of PIDs during a typical run on KNL.

Additionally, uses the right number of subrefined intervals in the blockstructured case

Still not beautiful...

This reverts commit b82e6531.

for consistency.

Otherwise debugging of nonlinear problems becomes hard.

We recently learned that usage of sub_elements in UFL is ambiguous...

otherwise my paraview throws me segfault.

To comply with latest pdelab.

Make anisotropic quadrature order possible for non DG examples.

- Use TensorProductElement in one example.

- We still use the sumfact option to swicth to the sum factorization
  code branch. This way it is still possible to easily swicth between
  sumfact and non sumfact code.

- Only TensorProductElements with the same degree in all directions
  will work. Anisotropie and adaption of quadrature rule will happen
  in the next commits.