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