On Fri, Aug 1, 2014 at 10:50 AM, John Peterson <jwpeter...@gmail.com> wrote: > On Fri, Aug 1, 2014 at 10:28 AM, John Peterson <jwpeter...@gmail.com> wrote: >> >> What I don't understand is why your original exodus file, which also >> has num_dim=3, actually displays correctly in paraview! This might >> take a big longer to track down... > > Here's something... your original file references side "5" of the TRI3 > elements in question: > elem_ss1 = 1, 17, 33 ; > side_ss1 = 5, 5, 5 ; > > > while the libmesh-generated version references side "1" > elem_ss1 = 1, 17, 33 ; > side_ss1 = 1, 1, 1 ; > > So while I don't understand what side "5" even means on a 3-sided > element, clearly it displays properly in paraview! > > Otherwise, all the element numbers in the sets themselves match > between the two files.
Yeah, and if you run in DEBUG mode you'll hit an assert with your input mesh file: Assertion `static_cast<size_t>(i) < side_map_size' failed. static_cast<size_t>(i) = 4 side_map_size = 3 Stack frames: 6 0: 0 libmesh_dbg.0.dylib 0x000000010e190831 libMesh::print_trace(std::ostream&) + 49 1: 1 libmesh_dbg.0.dylib 0x000000010e18c0a7 libMesh::MacroFunctions::report_error(char const*, int, char const*, char const*) + 55 2: 2 libmesh_dbg.0.dylib 0x000000010e566ea4 libMesh::ExodusII_IO_Helper::Conversion::get_side_map(int) const + 164 3: 3 libmesh_dbg.0.dylib 0x000000010e564289 libMesh::ExodusII_IO::read(std::string const&) + 4041 4: 4 exodus_rw-dbg 0x000000010b279065 main + 293 5: 5 libdyld.dylib 0x00007fff85ef05fd start + 1 [0] ./include/libmesh/exodusII_io_helper.h, line 653, compiled Aug 1 2014 at 09:56:52 -- John ------------------------------------------------------------------------------ Want fast and easy access to all the code in your enterprise? Index and search up to 200,000 lines of code with a free copy of Black Duck Code Sight - the same software that powers the world's largest code search on Ohloh, the Black Duck Open Hub! Try it now. http://p.sf.net/sfu/bds _______________________________________________ Libmesh-devel mailing list Libmesh-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-devel
