Hello,
in the literature, most adaptive mesh refinements bisect marked triangle and so we get two new triangles (and then, due to the closure, a few more). FEniCS divides marked triangles into four triangles (see the example below). Can I change/ modify the refinement algorithm such that it works like the standard newest vertex bisection?
from dolfin import *
import matplotlib.pyplot as plt
mesh = UnitTriangleMesh.create()
plt.figure(1)
plot(mesh)
cell_markers = MeshFunction(“bool”, mesh, mesh.topology().dim())
cell_markers.set_all(True)
mesh = refine(mesh,cell_markers)
plt.figure(2)
plot(mesh)
plt.show()