Re: [firedrake] Mapping nodal values from parent to extruded mesh
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
Thanks David, This is great, I'll try doing the reverse operation. - Tuomas On 03/05/2015 08:20 AM, David Ham wrote:
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 <mailto: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 <mailto: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
_______________________________________________ firedrake mailing list firedrake@imperial.ac.uk https://mailman.ic.ac.uk/mailman/listinfo/firedrake
Hi all, A couple of questions regarding those extruded mesh -> parent mesh copy operations, A while back I got the reverse 2d->3d copy working with a hard-coded pyop2 kernel (like in mesh extrusion). I'm not familiar with COFFEE syntax, is there documentation/examples somewhere? The 3d->2d copy routine in extrusion_extraction branch checks that the function space is CG. Does this method easily generalize to other spaces? DGxCG prisms seem to work OK. I'm also interested in using DGxDG and RTxCG/DG spaces. Cheers, Tuomas On 03/05/2015 10:24 AM, Tuomas Karna wrote:
Thanks David,
This is great, I'll try doing the reverse operation.
- Tuomas
On 03/05/2015 08:20 AM, David Ham wrote:
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 <mailto: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 <mailto: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
_______________________________________________ firedrake mailing list firedrake@imperial.ac.uk https://mailman.ic.ac.uk/mailman/listinfo/firedrake
participants (2)
-
David Ham
-
Tuomas Karna