I had tacitly assumed it would create the internal numbering from an arbitrary numbering. I now reordered it and it works fine! Thanks a lot for your help, the issue is resolved.
Christian On Mon, Sep 27, 2010 at 2:36 PM, Wolfgang Bangerth <[email protected]> wrote: > >> The grid it creates looks reasonable and has the outer shape it is >> supposed to have. However, if I run this code, which shows a line for >> each cell with boundary indicators of all faces, >> I mostly get lines like the following: >> >> 0 0 255 255 255 255 >> >> In particular, if I do a grep -v 0, I realize that there is no element >> without a 0 assigned face. In other words, according to this >> triangulation all cells have a face on the outer boundary, which is >> clearly not the way it should be. > > deal.II tries to match all faces of cells with faces of other cells. If it > can't find a matching face, then it marks it as external. My best guess is > that you got the numbering of vertices in each cell wrong (e.g. in counter > clockwise sense, rather than lexicographic, if this were in 2d). You can check > this by just creating a single cell and plotting it in gnuplot -- if the > numbering is wrong you'll see that it is distorted. > > Best > W. > > ------------------------------------------------------------------------- > Wolfgang Bangerth email: [email protected] > www: http://www.math.tamu.edu/~bangerth/ > _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
