Not yet implemented:

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