How to obtain coordinates of vertices of mesh?

[Alpha-Legion]

I used

VDG = VectorFunctionSpace(mesh, 'DG', 0)
vert_coords = assemble(interpolate(mesh.coordinates, VDG))

But it gives me a function of shape (num_cells, 3).
I expected it to be of shape (num_vertices, 3).
What's wrong?

[connorjward]

A DG0 function has a single DoF in the middle of each cell. A CG1 function would give you the result that you expect. Note that mesh.coordinates is already such a function though so the interpolation isn't needed.

[Alpha-Legion]

Thank you for point this out to me. I misunderstood the definition of the function space DG0.

[connorjward]

defelement (https://defelement.org/elements/discontinuous-lagrange.html) is a great resource for this sort of thing.

[Alpha-Legion]

Great! Thanks for recommending it to me.