I'm getting a negative jacobian error for a mesh created in CUBIT. CUBIT
says that the smallest jacobian in the mesh is around 1.e-4, but when I
load the mesh into libMesh and call fe->reinit(elem) on each element, I get
the error below.
The mesh looks fine to me (it's 1.3MB so I can't attach it here, but I'd be
happy to email it to anyone who is interested). Does anyone have any ideas
why this might be happening?
Thanks,
David
0: libMesh::print_trace(std::ostream&)
1: libMesh::MacroFunctions::report_error(char const*, int, char const*,
char const*)
2: libMesh::FEMap::compute_single_point_map(unsigned int,
std::vector<double, std::allocator<double> > const&, libMesh::Elem const*,
unsigned int, std::vector<libMesh::Node*, std::allocator<libMesh::Node*> >
const&)
3: libMesh::FEMap::compute_map(unsigned int, std::vector<double,
std::allocator<double> > const&, libMesh::Elem const*)
4: libMesh::FE<3u, (libMesh::FEFamily)0>::reinit(libMesh::Elem const*,
std::vector<libMesh::Point, std::allocator<libMesh::Point> > const*,
std::vector<double, std::allocator<double> > const*)
5: ./check_for_neg_jacobian-opt_real() [0x415d7e]
6: __libc_start_main
7: ./check_for_neg_jacobian-opt_real() [0x416496]
[0] /home/dknez/software/libmesh-src/src/fe/fe_map.C, line 618, compiled
May 5 2015 at 08:08:28
ERROR: negative Jacobian: -1.79627e-07 in element 7051
------------------------------------------------------------------------------
One dashboard for servers and applications across Physical-Virtual-Cloud
Widest out-of-the-box monitoring support with 50+ applications
Performance metrics, stats and reports that give you Actionable Insights
Deep dive visibility with transaction tracing using APM Insight.
http://ad.doubleclick.net/ddm/clk/290420510;117567292;y
_______________________________________________
Libmesh-devel mailing list
Libmesh-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-devel
