On Fri, May 7, 2010 at 3:03 PM, Karen Lee <[email protected]> wrote: > Dear libmesh users and developers, > > 1) I'm having a problem where I'm reading in the solution obtained from a > previous successful libmesh run to try to get the gradient of the solution > and I'm getting the following error: > > ERROR: negative Jacobian: -2.55565e-10 in element 0 > [0] src/fe/fe_map.C, line 341, compiled Mar 30 2010 at 16:03:26 > terminate called after throwing an instance of 'libMesh::LogicError' > what(): Error in libMesh internal logic > > My first element consists of these points: > > 1 0.0635691 0.0267701 -0.0401152 4.876378533e-11 > 2 0.063976 0.0263403 -0.0393085 4.866489885e-11 > 3 0.063604 0.0261774 -0.0397662 4.875840464e-11 > 4 0.0639932 0.0260264 -0.0402555 4.881947696e-11 > > which to me look distinct enough...
Distinct, yes, but have you checked whether it is inverted? Assuming this is a tet, it will be inverted (and thus have negative jacobians) if you take one of the vertices and 'pull it through' the opposite face. You say the solution came from a previous libmesh run, what mesh format did you write in? -- John ------------------------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
