So, apperantly mshr is on its way out. So, how to generate a surface mesh of triangles in either pymesh or gmsh and more specifically, how do I get it into fenics?

And furhtermore, what is wrong with this code, why does it not work?

N=6

sphere=Sphere(Point(0,0,0,1)

mesh=generate_mesh(sphere,N)

surface=BoundaryMesh(mesh,“exterior”)

V=FunctionSpace(surface,‘P’,1)

Y=Expression('3*sqrt(10)*(x[0]*x[0]-x[1] x[1])(7*[x[2]

*[x[2]-1)’,degree=2)*

u_n=project(Y,V)

u=TrialFunction(V)

v=TrialFunction(V)

a=inner(grad(u),grad(v))dx

u_n=project(Y,V)

u=TrialFunction(V)

v=TrialFunction(V)

a=inner(grad(u),grad(v))

*dx*

u=Function(V)

L=u_nvu=Function(V)

L=u_n

solve(a==L,u)

Is there something wrong with how I define forms on my surface?