quadrature points

Question

Hi,

is there a way to extract the quadrature points?
Specifically, is there a way to extract the coordinates for the quadrature points associated with the VectorFunctionSpace like:

Quad = VectorFunctionSpace(mesh, "Quadrature", 2)

Thanks!

Answer 1

Hi, consider

from dolfin import *

mesh = UnitSquareMesh(10, 10)
gdim = mesh.geometry().dim()
V = VectorFunctionSpace(mesh, 'Quadrature', 2)
# Quad points
xq = V.dofmap().tabulate_all_coordinates(mesh).reshape((-1, gdim))
# You might want to remove the duplicates
xq0 = xq[V.sub(0).dofmap().dofs()]  

Comments

Thanks Mirok! That works !!