Hello,

I want to define my mesh to be a circle. In Fenics I know there is a command called CircleMesh that makes this very easy.

Q1) Is there a simple way of doing this in Firedrake?

I looked around and see how circlemanifoldmesh can do this with help for ExtrudedMesh. I found some code that can make a ring and I tried running it myself and defining a function on that mesh. Unfortunately, I get an error, see below.

Q2) If one needs to use CircleManifoldMesh is there something else one needs to do to use a circular geometry as your mesh?

Below are the lines that I tried and the error
>>> from firedrake import *
>>> m    = CircleManifoldMesh(20, radius=2)
>>> mesh = ExtrudedMesh(m, 5, extrusion_type='radial')
>>> x    = SpatialCoordinate(mesh)
>>> V    = FunctionSpace(mesh,'CG',1)
>>> U    = Function(V).interpolate(-(pow(x[0]-0.5,2) + pow(x[1]-0.5,2)))
>>> plot(U);
Traceback (most recent call last):
File "<stdin>", line 1, in <module>
File "/home/fpoulin/software/firedrake/src/firedrake/firedrake/plot.py", line 308, in plot
    **kwargs)
File "/home/fpoulin/software/firedrake/src/firedrake/firedrake/plot.py", line 502, in two_dimension_plot
    num_sample_points)
File "/home/fpoulin/software/firedrake/src/firedrake/firedrake/plot.py", line 457, in _two_dimension_triangle_func_val
    raise RuntimeError("Unsupported Functionality")
RuntimeError: Unsupported Functionality

------------------
Francis Poulin
Associate Professor
Department of Applied Mathematics
University of Waterloo

email:           fpoulin@uwaterloo.ca
Web:            https://uwaterloo.ca/poulin-research-group/
Telephone: +1 519 888 4567 x32637