I found the following code to get quadrature points and weights in the standard simplex. But I want to get the quadrature points on mesh elements/cells directly. Is any way in Fenics to do this? What are the choices for scheme (what are the candidate input strings)?

Thanks! This code gives a list of integration points, but how to get the weights? Can I get the points and weights arranged by elements/cells? What is the default integration scheme?

Thanks! I have read the doc and want to use this function to get integration points and weights. The following code gives some error. How to construct FIAT cell with mesh?

from FIAT import *
cell = Cell(mesh, 0)
print(create_quadrature(cell,2))

AttributeError: ‘dolfin.cpp.mesh.Cell’ object has no attribute ‘get_shape’

Yes the weights do not change, but to align it with the nodes correctly need the information of integration scheme, the distribution, and the order of the integration nodes. If I align them manually, I am afraid of using the wrong information and get bugs that are hard to realize.

Isn’t there a way to get both the integration points and weights of a mesh cell directly? Dolfin or dolfinx will be OK for me. I just installed both versions.