On Tue, May 2, 2017 at 12:48 PM, Harshad Sahasrabudhe <hsaha...@purdue.edu> wrote:
> Hi all, > > I'm trying to convert an adaptive tetrahedral mesh generated by Sentaurus > from their grd format to Gmsh's msh format for importing in my program > which uses LibMesh. You're importing an "adaptive" mesh? Does it already have hanging nodes? > I keep getting the following error due to improperly > ordered nodes in the element: > > ERROR: negative Jacobian -2.85931e-06 at point (x,y,z)=( 1.54167, 1, > 0.5059) in element 69 > > When I manually reorder the nodes for the above element, the error is fixed > for that element. Is there a simple way to correct the node ordering so > that I can use it in my conversion code? I tried computing the volume using > triple vector product, but I still get this error for some elements. I'm > using first order Lagrange basis functions. > I would assume the ordering of each Tet's nodes in Sentaurus is the same, and you just need to figure out the permutation of the node ids that gives you the libmesh/Gmsh numbering (they are the same). I would just write a simple reader for the Sentaurus file format and skip the Gmsh conversion step as I don't see the need for it. -- 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 Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
