On top of Guido's suggestion you will also have to set the proper direction_flag on each cell, see
http://www.dealii.org/developer/doxygen/deal.II/DEALGlossary.html#GlossDirectionFlag Best -Sebastian On Fri, Jul 29, 2011 at 4:28 PM, Guido Kanschat <[email protected]> wrote: > The best way is changing your coarse mesh so that vertices are in > lexicographical order. Then, deal.II will not have to second-guess your > intentions. > > The orientation of the normal vectors comes only into play in Mapping and is > accessed by FEValues::cell_normal_vector. > > Best, > Guido > > On 07/29/2011 03:36 PM, Daniel Castañon wrote: >> >> Dear all, >> >> I have to solve a PDE on the Torus surface. After some work I was able >> to create the entire Triangulation<2,3> and Boundary<2,3> classes.I >> had to call the create_trianguation_compatibility() function because >> the coarsest mesh (the Triangulation<,>) was created with the right >> hand >> numbering scheme (anticlockwise order) for vertices in cells. >> >> Now, I want to check that the normals of the coarsest mesh are >> effectively pointing outwards from the surface, as I originally took >> care using the right hand convention. I have read in the >> create_trianguation() source code that this function takes care of >> making the numbering of the vertices in such a way that the normal >> field is continuous. So, I think that I just have to pick only one >> cell, take its first and second vertices, take the cross product and >> see if is pointing outwards from the surface. Is this correct? >> >> Thanks. >> >> Daniel. >> >> >> >> >> _______________________________________________ >> dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii > > _______________________________________________ > dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii > _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
