How to determine gauss point coordinates in a mesh

mesh
4

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.16666667]
 [ 0.83333333  0.16666667]
 [ 0.66666667  0.83333333]
 [ 0.16666667  0.33333333]
 [ 0.16666667  0.83333333]]
1 Like