On 25-Aug-2021, at 7:16 PM, Lawrence Mitchell <wence@gmx.li> wrote:

Dear Praveen,

On 25 Aug 2021, at 14:44, Praveen C <praveen@gmx.net> wrote:

Dear all

I am a new user of firedrake but old user of fenics.

I need your help to compute the mesh diameter, i.e.,

max diameter of all cells in the mesh

or some representative mesh size to use in a mesh convergence study.

If you do:

f = interpolate(CellDiameter(mesh), FunctionSpace(mesh, "DG", 0))

with f.dat.vec_ro as v:
  fh = v.max()

will give you it, I think.

Lawrence