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
