Re: [Libmesh-devel] What's going on with reduced_basis_ex4?

2014-01-03 Thread Roy Stogner
On Fri, 3 Jan 2014, Roy Stogner wrote: > I guess the thing to do is fix our METIS interface to handle > disconnected meshes (or possibly just remove an overzealous > assertion), Fixing that assertion fixes the problems I was seeing with SerialMesh, but seems to have uncovered a real regression w

[Libmesh-devel] What's going on with reduced_basis_ex4?

2014-01-03 Thread Roy Stogner
I should have been keeping a closer eye on the buildbot reports over the holidays, I guess. reduced_basis_ex4 is currently dying in online mode when it tries to get Metis to partition eim_data/interpolation_points_mesh.xda, which is a mesh of 21 QUAD4 elements connecting 84 nodes. Yes, read that

Re: [Libmesh-devel] Issue with "duplicate" nodes and EquationSystems.write()

2014-01-03 Thread Derek Gaston
Right - we don't do m->n restart so that works for us. Also - recall that I don't like this reordering business because the un-reordering destroys the original ordering... so that's another reason not to reorder... Unfortunately though it turns out that something about the reordering might be

Re: [Libmesh-devel] Issue with "duplicate" nodes and EquationSystems.write()

2014-01-03 Thread Kirk, Benjamin (JSC-EG311)
Short term yes, but this will break m->n restart I think. I recall something similar in the mesh a few yearns back when two disjoint nodes created an identical index. The fix there was to handle identical keys for separate nodes. If I can extend that to this case all will be good. On Jan 3, 2

Re: [Libmesh-devel] Issue with "duplicate" nodes and EquationSystems.write()

2014-01-03 Thread Derek Gaston
If I just comment out the renumbering in EquationSystems::write() then everything looks fine. Would anyone object to me adding a flag to EquationSystems::write() allowing you to turn off the renumbering? Derek On Fri, Jan 3, 2014 at 10:39 AM, Derek Gaston wrote: > Ok - MeshCommunication::assi

Re: [Libmesh-devel] Issue with "duplicate" nodes and EquationSystems.write()

2014-01-03 Thread Derek Gaston
Ok - MeshCommunication::assign_global_indices() is beyond my comprehension... BUT there are a couple of "#if 0" error checks in there that if un-#if'd they do trigger these things with my test case: Error: nodes with duplicate Hilbert keys! node 7, (x,y,z)=( 0.5, 0.5,0) has Hilber

Re: [Libmesh-devel] Issue with "duplicate" nodes and EquationSystems.write()

2014-01-03 Thread Roy Stogner
On Thu, 2 Jan 2014, Derek Gaston wrote: One of my users brought an issue to my attention: libMesh is not properly handling multiple nodes occupying the same position when writing out xda/xdr files for the solution vector. I am attaching a mesh and main.C that show the issue.  If you look at th