On 11/29/23 09:23, Lance Zhang wrote:

//***********************************************************************************************************************************//
         if(cell->at_boundary()){
           continue;
           }//to check whether the cell exists

I don't think this is what you want. Here, you are only asking whether the cell has *any* faces that are at the boundary. But if you want to know whether a cell has a neighbor across a certain face, you need to ask whether the cell is at the boundary *at that particular face*. Daniel's code does that.

Best
 W.

--
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/88a5f345-a088-473e-a893-1ed6736440de%40colostate.edu.

Reply via email to