quadrature rule -> numpy
At the moment the buffer is exported as read write, so that the values can be modified, i.e.,
points, weights = rule.get()
weights *= integrationElement(points)
shouldn't work.
We should also consider directly exporting the points buffer on the rule and return modified weights, i.e.,
weights = rules.weights(element) # include integrationElement
values = function( element(rule) )
should perhaps work. '.get' could be kept to access the original weights.