Would it be correct to say then that your 2d mesh would be a single face of
the extruded 3d mesh?
that is correct.
If you work with the svn repository rather than the last release, then there
is a function GridTools::extract_boundary_mesh that could give you the 2d
mesh given the 3d mesh, along with a mapping from one to the other. That
might make your life marginally simpler.
Thank you! I will look into this function.
To facilitate the mapping from 3d cells to 2d cell, you may wish to keep a
vector<Triangulation<2>::active_cell_iterator> with all the cells of the 2d
mesh somewhere, and use the user_pointer field of the cell iterators of your
3d mesh to store pointers into this array.
Thank, I didn't know about the user_pointer field. I was actually
planning on keeping a Vector (corresponding to cell-centered DoFs of the
3D grid) that contains the index of the corresponding DoF (or cell,
etc.) on the 2D grid. Aside from the time to generate this Vector, is
there any reason this approach would be ill-advised?
Thanks,
Dan
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
