Tabulation of function space V

[dokken]

@prusso630

Is there a way from python “V” or from trial function in python “u” to get output of domain points as they start to get modified?

For instance I might light to see output in terms of node points how or if they get modified from a new test function f_n = u * v? (As far as what the end resulting expressions look like…)

f_n = u *v is just a symbolic expression, no computation or modification of dof coordinates has happened yet.

When you call
dolfinx.fem.form(f_n*ufl.dx) C++ code is generated for assembly of the local element tensor.

Calling assemble_matrix on this object assembles the contributions into a global tensor.

If you want to control the integration points, you can set “quadrature_degree” in ufl.Measure, see for instance:

Numerical values from ufl.SpatialCoordinate

ufl.SpatialCoordinate does not imply the “mesh-nodes” (meaning the vertices for linear meshes). The x=ufl.SpatialCoordinate(mesh) represents any quadrature point used during integration in its physical space. What I would suggest to do, is to interpolate your function (non-analytical) into a suitable function space. Say, if you want to actually use the mesh nodes (vertices), you can create a first order Lagrange function and interpolate into that, i.e. V = dolfinx.fem.FunctionSpace(mesh, ("L…

To inspect these coordinates you should look at Crash when trying to obtain coordinates of quadrature points - #2 by dokken
Mapped quadrature points, weights and solution at these quadrature points - #2 by dokken