On 5/9/20 3:55 AM, cd wrote:
To start with I have a Triangulation<dim> volume_mesh and from this, I
use GridGenerator::extract_boundary_mesh to generate Triangulation<dim - 1,
dim> surface_mesh.
Then I solve a std fem elasticity problem and compute volume_solution vector.
Further, I also solve a meshfree problem using only the surface_mesh
information and from this, I have a surface_solution.
Now I need a dof map to compare the two solution vectors. How can I get this
dof map?
The GridGenerator::extract_boundary_mesh() function returns this kind of map
that connects the cells of the volume and the boundary mesh. If you call the
function with two DoFHandler objects (one of which is based on an empty
triangulation), then you will return a map between DoFHandler iterators that
you can at a later point use to call
cell->get_dof_indices()
to see which degrees of freedom are associated with each volume or surface cell.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/e38f32b1-4b71-c90e-f580-82259424f96f%40colostate.edu.