Hi all!

I'm hitting the following issue with libMesh: We have a mesh with
let's say 2 blocks. They are sharing the same side, however each one
has its own set of nodes and these nodes share the same position in
the space - i.e. they have the same coordinates. When we try to save
and XDR file for restart, libMesh executes
MeshTools::Private::globally_renumber_nodes_and_elements(mesh);
(equation_systems_io.C around line 450). This renumbering will trigger
that duplicate index error coming from Hilbert indexing. I suspect
that the renumbering is in place because of adaptivity (?).  I tried
to comment that call out and it "fixed" our problem (no adaptivity in
our run). However, I do not believe that is the proper fix.

Any ideas on how to fix this issue? I'd volunteer to prepare the patch...

Thanks,
--
David Andrs

------------------------------------------------------------------------------
Live Security Virtual Conference
Exclusive live event will cover all the ways today's security and 
threat landscape has changed and how IT managers can respond. Discussions 
will include endpoint security, mobile security and the latest in malware 
threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/
_______________________________________________
Libmesh-devel mailing list
Libmesh-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-devel

Reply via email to