Hi,
In this plasticity example I want to plot the von Mises stress and plastic strain, but when interpolating I get the error (using dolfinx 0.7.3):
Cannot get interpolation points - no Basix element available. Maybe this is a mixed element?
Here is a MWE:
import dolfinx
import ufl
from mpi4py import MPI
def f(x):
return x[0]
mesh = dolfinx.mesh.create_unit_square(MPI.COMM_WORLD, 10, 10)
el = ufl.VectorElement("Quadrature", mesh.ufl_cell(), degree=1, dim=1, quad_scheme='default')
V = dolfinx.fem.functionspace(mesh, el)
u = dolfinx.fem.Function(V)
u.interpolate(f)
This is a known issue, and should be implemented in the upcoming 0.8 release, but I would like to be able to plot my values before that. Is there a (temporary) different way to still interpolate these values?
Thank you!