Dear George,

Well one way to do something like this is to observe that the node of the DG0 space is at the centre of mass. Accordingly if you project mesh.coordinates into a VectorFunctionSpace with element DG 0 then you get a Function whose n-th value is the centre of mass of cell n.

Regards,

David

On Fri, 18 Nov 2016 at 22:21 George Ovchinnikov <lives9@gmail.com> wrote:

Hello,

There is a nice function, which gives me cell by its coordinates
"mesh.locate_cell(x)"

I wonder, does it have some kind of inverse analogue, so I could get
coordinates for a cell (given by id of cell)?

Sincerely,
George

_______________________________________________
firedrake mailing list
firedrake@imperial.ac.uk
https://mailman.ic.ac.uk/mailman/listinfo/firedrake