Lawrence, Got it, thanks! Justin On Wed, Jan 11, 2017 at 4:38 PM, Lawrence Mitchell < lawrence.mitchell@imperial.ac.uk> wrote:
Hi Justin,
On 11 Jan 2017, at 21:48, Justin Chang <jychang48@gmail.com> wrote:
Hi all,
Is there a way to pass -dm_refine <number> into Firedrake? That or what's the best way to do a parallel mesh refinement after the initial mesh is created and distributed (I think I recall seeing another mail about this last month).
If you just want the refined distributed mesh, you can use this bit of code:
from firedrake.mg.impl import filter_exterior_facet_labels coarse = Mesh(...) dm = coarse._plex # uniform refinement dm.setRefinementUniform(True) fdm = dm.refine() fdm.removeLabel("interior_facets") fdm.removeLabel("op2_core") fdm.removeLabel("op2_non_core") fdm.removeLabel("op2_exec_halo") fdm.removeLabel("op2_non_exec_halo") filter_exterior_facet_labels(fdm)
fine = Mesh(fdm, dim=coarse.ufl_cell().geometric_dimension(), distribute=False)
We could certainly wrap this up nicely.
Cheers,
Lawrence
_______________________________________________ firedrake mailing list firedrake@imperial.ac.uk https://mailman.ic.ac.uk/mailman/listinfo/firedrake