Hello everyone,

Sorry to revive this thread, but I have a similar question that follows up 
his question. In the case of a distributed grid, there might be 
vertices/DoFs on a corner (or even edges in 3D) of a cell. 

For example, see attached figure (at the trailing edge of the airfoil in 
the square, in the yellow square). Domain 0 is not aware that one of its 
corner is on a boundary. To aggravate the issue, the DoFHandler with FE_Q 
distributes the DoFs such that Domain 0 owns that corner DoF.

It seems like I might have to loop over the neighbors of every cell to make 
sure that account for all the corners/edges DoFs that might be on a 
boundary, but that seems inefficient.

Anyone knows of a better way within deal.II to deal with this?

Best regards,

Doug

On Friday, May 1, 2020 at 7:21:42 PM UTC-4 Wolfgang Bangerth wrote:

> On 5/1/20 5:07 PM, Shahab Golshan wrote:
> > I have also another question on the function 
> */find_cells_adjacent_to_vertex/*.
> > As input argument, it needs the vertex index, not the vertex as a point. 
> How 
> > can I find the index of a vertex (to use in this function) from the 
> vertex 
> > location (point).
>
> cell->vertex_index(v)
>
> is your friend, where v=0...2^d :-)
>
> Best
> W.
>
>
> -- 
> ------------------------------------------------------------------------
> Wolfgang Bangerth email: [email protected]
> 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 [email protected].
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/ef658a63-82c1-46df-8c49-a3ea4a9a063fn%40googlegroups.com.

Reply via email to