There's some more strange things going on with the pyramid... In the definition of the geometric element (src/geom/cell_pyramid5.C) the numbering of the base is *clockwise* instead of following the usual counter-clockwise numbering used by e.g. Hexes. The build_side() code in the class is self-consistent with this numbering, but the numbering of the Lagrange basis functions is not, nor is the element defined in reference_elements/3D/one_pyramid.xda (this reference element should be inverted with respect to the definition in the Pyramid5 class, as far as I can tell.) Some questions...
Remind me again what was the reason for numbering it this way? I'm sure there was one, I just can't remember what it was. I'll change the basis function numbering and the reference element definition. Is there any other obvious part of the library to fix up? Thanks, -- John ------------------------------------------------------------------------- This SF.Net email is sponsored by the Moblin Your Move Developer's challenge Build the coolest Linux based applications with Moblin SDK & win great prizes Grand prize is a trip for two to an Open Source event anywhere in the world http://moblin-contest.org/redirect.php?banner_id=100&url=/ _______________________________________________ Libmesh-devel mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-devel
