How to get Gaussian integration point directly on a mesh element/cell?

[dokken]

See:

How to determine gauss point coordinates in a mesh

Do you want the physical quadrature points for each cell in your function space? These are automatically generated in ffc in fenics. If you want to access them you could create a quadrature function space and tabulate the dof coordinates: from dolfin import * mesh = UnitSquareMesh(1, 1) el = FiniteElement("Quadrature", mesh.ufl_cell(), degree=2, quad_scheme="default") V = FunctionSpace(mesh, el) print(V.tabulate_dof_coordinates()) returning [[ 0.83333333 0.66666667] [ 0.33333333 0.16666…

@xi_wang

create_quadrature

See: Bitbucket