Calculate the evaluation and jacobians of geometry mappings of cells
using sum factorization. By avoiding all geometry mapping of intersections
we can avoid the face-twist problem and just make the grid consistens (not
yet done).

TODO:
- Make grid consistent
- Vectorization
- 3D

Not yet implemented:

- Vectorization of quadrature loop
- Anisotropic amount of quadrature points per direction

Note: This is still structured for now

Rename the function and remove the backend switch. Reasoning:

- This function was only called once without the backend switching
  mechanism.
- The function takes different arguments than the pdelab
  pymbolic_quadrature_position.
- It returns only one component of the quadrature point.

Probably not the solution we want to have in the long run but right
now it is not possible to have a backend switch for
default/blockstructured/sumfact. This way the old code should behave
the same way as before and it is possible to switch between those
three things.

Adjust time stepping parameters of navierstokes 3d test

See merge request dominic/dune-perftool!253

Merge branch 'test/adjust-time-step-in-test' of ssh://gitlab.dune-project.org:22022/dominic/dune-perftool into test/adjust-time-step-in-test

Revert "Switch to repositories on submodules which comply with our filters in CI network configuration"

This reverts commit 92bd8051fd7fe925914f2c57c3ee6972c1a9fe73.

This reverts commit 32e8611140c7af63f26497b88c9b1626262e4d4c.

Revert "Switch to repositories on submodules which comply with our filters in CI network configuration"

This reverts commit 92bd8051fd7fe925914f2c57c3ee6972c1a9fe73.

This reverts commit 32e8611140c7af63f26497b88c9b1626262e4d4c.