How to use same vector-valued mathematical expressions with both ufl and numpy?

[bleyerj]

Maybe this ?

Evaluating Finite Element Function

The points have the wrong shape (single point, the function expects multiple points): from mpi4py import MPI import numpy as np from dolfinx import mesh from dolfinx import fem from dolfinx.fem import FunctionSpace from dolfinx import default_scalar_type import ufl from dolfinx.fem.petsc import LinearProblem from dolfinx import geometry domain = mesh.create_unit_square(MPI.COMM_WORLD, 8, 8, mesh.CellType.quadrilateral) V = FunctionSpace(domain, ("Lagrange", 1)) uD = fem.Function(V)…