Evaluate solution at grid points of structured mesh and reshape to tensor

[dokken]

See for instance:

Reordering nodal values in accordance with meshgrid created using numpy

Consider the following minimal example: import numpy as np import matplotlib.pyplot as plt from mpi4py import MPI import dolfinx import ufl # define mesh and function space N = 10 degree = 1 msh = dolfinx.mesh.create_unit_square(MPI.COMM_WORLD,N,N) V = dolfinx.fem.functionspace(msh, ("Lagrange", degree)) x = ufl.SpatialCoordinate(msh) uf = dolfinx.fem.Function(V) uf.interpolate(lambda x : np.sin(2*np.pi*x[0]) + 3*x[1]) # points = msh.geometry.x # I can use either this or tabulate_dof_coordina…

DOLFINx in general does not consider a mesh structured. All meshes are unstructured, and data is reordered with respect to data-locality.