On Wed, May 3, 2017 at 7:21 AM, Harshad Sahasrabudhe <[email protected]> wrote:
> 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? > Forget about Jacobians, just call elem->volume(). If the element is inverted, the volume() will be negative... -- John ------------------------------------------------------------------------------ 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
