Hi Tuomas,

This wasn't there this morning but I've implemented one of the cases (pulling out the top and bottom maps). The result is in the extrusion_extraction branches of both PyOP2 and Firedrake. The 2d->3d operation would be rather similar except that you'd have to interrogate the fiat_element on the extruded Function in order to determine which extruded nodes are "on top of" the 2d nodes you have. Unfortunately I won't have time to do that one soon (huge amounts to do before SIAM CSE next week) but feel free to have a try and complain when it doesn't work!

Cheers,

David

On 4 March 2015 at 22:23, Tuomas Karna <tuomas.karna@gmail.com> wrote:

Hi all,

I'd need to copy nodal values between fields on parent and extruded
meshes. For example, copy 2d->3d (constant over vertical) or 3d->2d
(extract surface/bottom level). The horizontal function space is the
same. Is there an easy way to do this? Is there a map of extruded nodes
somewhere?

Thanks,

Tuomas

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

Dr David Ham
Departments of Mathematics and Computing
Imperial College London

http://www.imperial.ac.uk/people/david.ham