> 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
