def refine(hspace: HierarchicalSpace,
marked_entities,
type='cells') -> HierarchicalSpace:
"""
Refines the hierarchical space based on the type of marked entities.
:param hspace:
:param marked_entities:
:param type:
:return:
"""
logging.info("""
Refining hierarchical space
===========================
space.nlevels = {}
hmesh.nlevels = {}
""".format(hspace.nlevels, hspace.mesh.nlevels))
marked_cells = marked_entities
new_cells = hspace.mesh.refine(marked_cells)
marked_functions = hspace.functions_to_deactivate_from_cells(marked_cells)
hspace.refine(marked_functions, new_cells)
logging.info("""
After refinement, the space has
===============================
space.nfuncs = {}
space.nelems = {}
""".format(hspace.nfuncs, hspace.mesh.nel))
return hspace
def refine(hspace: HierarchicalSpace,
marked_entities,
type='cells') -> HierarchicalSpace:
"""
Refines the hierarchical space based on the type of marked entities.
:param hspace: hierarchical space
:param marked_entities: indices of the marked entities, in this case, cells
:param type: defaults to cells, not implemented 'function'-refinement yet.
:return: a refined hierarchical space
"""
logging.info("""
Refining hierarchical space
===========================
space.nlevels = {}
hmesh.nlevels = {}
""".format(hspace.nlevels, hspace.mesh.nlevels))
marked_cells = marked_entities
new_cells = hspace.mesh.refine(marked_cells)
marked_functions = hspace.functions_to_deactivate_from_cells(marked_cells)
hspace.refine(marked_functions, new_cells)
logging.info("""
After refinement, the space has
===============================
space.nfuncs = {}
space.nelems = {}
""".format(hspace.nfuncs, hspace.mesh.nel))
return hspace