caching quadratures from numpy arrays
Write a class Quadrature
that takes numpy arrays for points and weights and that can be passed to an Operator
/Scheme
or a Space
(for codegen):
scheme.setQuadratureOrders([0,quadrature],None)
Should be done so that the dune-python
quadrature techniques used for dune-geometry
or QuadPy
can be directly used.