Hi John
Forget about Jacobians, just call elem->volume(). If the element is
> inverted, the volume() will be negative...
>
Calculating the volumes worked. I was doing it before, but I was making a
mistake in the indices somewhere. I fixed it and now I have all positive
Jacobians.
Thanks for your
On Wed, May 3, 2017 at 7:21 AM, Harshad Sahasrabudhe
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 t
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
neg
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.
Thanks! Yes the mesh is only tet4. I went through the LibMesh code for
calcul
On Tue, 2 May 2017, Harshad Sahasrabudhe wrote:
> 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 fo
Hi John,
Thanks for replying.
You're importing an "adaptive" mesh? Does it already have hanging nodes?
I meant the mesh was refined adaptively. It doesn't have hanging nodes.
I would assume the ordering of each Tet's nodes in Sentaurus is the same,
> and you just need to figure out the permuta
On Tue, May 2, 2017 at 12:48 PM, Harshad Sahasrabudhe
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 alread
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. I keep getting the following error due to improperly
ordered nodes in the element:
ERROR: negative Jacobian -2.85931e-06 at
