def circumradius(self, o):
mesh = coarsen(o.ufl_domain())
return firedrake.Circumradius(mesh)
def circumradius(self, o):
mesh = o.ufl_domain()
hierarchy, level = utils.get_level(mesh)
new_mesh = hierarchy[level - 1]
return firedrake.Circumradius(new_mesh.ufl_domain())
