Local Refinement of SubDomain

[dokken]

First please use 3x` encapsulation to format the code, ie

```python
#add code here
```

Secondly, consider

Adapt for mesh and MeshFunctions

The issue is that when refining the mesh, new facets are introduced into the mesh (that had no previous parent facet, as they are internal to a cell, these facets get large numbers as values. I would like to note that your example is a very good one (good way of presenting a minimal working example), but I have simplified it slightly for easier visual interpretation: from dolfin import * parameters["refinement_algorithm"] = "plaza_with_parent_facets" class Top(SubDomain): def inside(self,x,…