On 11/16/2017 07:55 PM, 曾元圆 wrote:
> Thank you for your reply, Bruno.
> Now I understand what you mean. But I'm still curious about how the cell # can
> figure out the correct "neighbor" of it. You mentioned that when you ask for
> the neighbor of cell #, you get a cell that has the same level of refinement
> or that is coarser. But in anisotropic refinement we don't identify a coarser
> cell though the "level" it belongs to, which is the case in isotropic refinement.
> I suppose we have stored a "neighbor" of cell # correctly before ask for that.
> But how do we in fact set the link between cell# and it's neighbor, i.e. the
> right whole mother cell? I just want to know how these cells belonging to
> different "levels" are connected. I know perhaps it is not easy to explain in
> a few words, so it's preferable if there are some available thesis or
> documentation related to this.
I don't know whether that's documented in great detail anywhere, but it
happens at the end of refining each cell:
https://github.com/dealii/dealii/blob/master/source/grid/tria.cc#L4493
In other words, each cell's neighbors are the children of its parent's neighbors.
We're always interested in improving documentation, and if you think that
there is a good place to add this sort of information, we'd be quite happy to
include whatever you come up with!
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email:
bang...@colostate.edu
www:
http://www.math.colostate.edu/~bangerth/