I would like to change the coordinate of mesh based on the marker.
Specifically, I would like to change the coordinates at the boundary of 3D mesh.
Following is the small example of what I want to do.
from fenics import * mesh = UnitSquareMesh(32,32) bdry= MeshFunction("size_t", mesh, mesh.geometry().dim() - 1, mesh.domains()) CompiledSubDomain("x < 0 + 1e-12").mark(bdry, 1)
I know it is possible to change the coordinate with something like
x = mesh.coordinates() x[:,1] = some manipulation
However, is it possible to do the same manipulation based on the boundary marker ??