Mesh: QuadTree: Hanging NodesΒΆ

You can give the refine method a function, which is evaluated on every cell of the TreeMesh.

Occasionally it is useful to initially refine to a constant level (e.g. 3 in this 32x32 mesh). This means the function is first evaluated on an 8x8 mesh (2^3).

from SimPEG import Mesh
import matplotlib.pyplot as plt

def run(plotIt=True):
    M = Mesh.TreeMesh([8, 8])

    def refine(cell):
        xyz =
        dist = ((xyz - [0.25, 0.25])**2).sum()**0.5
        if dist < 0.25:
            return 3
        return 2

    if plotIt:
        M.plotGrid(nodes=True, cells=True, facesX=True)
            'Cell Centers',
            'Hanging Nodes',
            'X faces',
            'Hanging X faces'

if __name__ == '__main__':

Total running time of the script: ( 0 minutes 0.224 seconds)

Generated by Sphinx-Gallery