Interpolating numerical values (2D numpy array) over a mesh

[dokken]

This is far from a general example representative of your question:

Interpolating numerical values (2D numpy array) over a mesh

as the code involves dolfinx_mpc and solving a problem.
It is unclear why you need this for the problem at hand.

You should also consider the following:

Extract the coordinates of the boundary points and the value of a function on them

Consider the following: from mpi4py import MPI import dolfinx import numpy as np import ufl mesh = dolfinx.mesh.create_unit_square(MPI.COMM_WORLD, 16, 16) v_cg1 = ufl.VectorElement("CG", mesh.ufl_cell(), 1) V = dolfinx.fem.FunctionSpace(mesh, v_cg1) f = dolfinx.fem.Function(V) f.interpolate(lambda x: (0.1*x[1]*abs(x[0]), 0.05*np.sin(2*x[0]))) mesh.topology.create_connectivity(mesh.topology.dim-1, mesh.topology.dim) boundary_facets = dolfinx.mesh.exterior_facet_indices(mesh.topology) boundary_v…

and explain how your relation with the mesh nodes is assumed.