Hi Roy, Is the mesh entirely Tet4, then? If so then there's only two > equivalence classes of orderings; you could try one, test the > jacobian at any point, and if it's negative then swap two nodes and > you're done.
I tried using the face normals defined by Sentaurus, and I still have the negative Jacobian problem. How can I calculate the Jacobian in LibMesh without getting an error if it's negative? Do I need to create the full mesh or can I just create an element and compute it's Jacobian using FEMap? Thanks! Harshad ------------------------------------------------------------------------------ Check out the vibrant tech community on one of the world's most engaging tech sites, Slashdot.org! http://sdm.link/slashdot _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
