Hi Martin,
On 03/06/2014 02:08 AM, Martin Luethi wrote: > Hi David > > David Knezevic <dkneze...@seas.harvard.edu> writes: >> I've found a case where ExodusII_IO::write_discontinuous_exodusII has a >> weird bug. The example is here (on dropbox): > visit shows the same problem (missing elements in exo) for > discontinuous_plot.exo > > But the elements are in the Exodus File. You can inspect it with ncdump > (it's just a NetCDF file after all), and there are enough elements and > nodes within the file (same number as gmv). > > Exodus > num_nodes = 116240 ; > num_elem = 14530 ; > num_el_blk = 2 ; > num_side_sets = 3 ; > num_el_in_blk1 = 4280 ; > num_nod_per_el1 = 8 ; > num_el_in_blk2 = 10250 ; > num_nod_per_el2 = 8 ; > > GMV > nodes 116240 > cells 14530 > > So in Exouds the elements are in two blocks (not sure why). > > The elemnt conneciivities look the same in both cases. But maybe > elem_num_map in Exodus gets confused somehow, there are lots of zeroes > at the end of this block. The array serves as translation between > internally stored element numbers, and those shown in a viewer. > > Do you have the Exodus file format documentation? It's expleined there. Thanks, those comments are very helpful. (I'm cc'ing the list in case anyone else is interested.) Note that test_mesh.exo has two element blocks, so that might be what is triggering the bug... When I use a single-element-block mesh (e.g. via build_cube) in the example, it works fine. I'll look into it more. Thanks, David ------------------------------------------------------------------------------ Subversion Kills Productivity. Get off Subversion & Make the Move to Perforce. With Perforce, you get hassle-free workflows. Merge that actually works. Faster operations. Version large binaries. Built-in WAN optimization and the freedom to use Git, Perforce or both. Make the move to Perforce. http://pubads.g.doubleclick.net/gampad/clk?id=122218951&iu=/4140/ostg.clktrk _______________________________________________ Libmesh-devel mailing list Libmesh-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-devel
