On Tue, Jan 12, 2016 at 11:17 AM, Salazar De Troya, Miguel < salazardet...@llnl.gov> wrote:
> Hello > > I’m attaching a simple code that illustrates the problem. I open a mesh > (l_problem.xda) and perform a uniform refinement (if I don’t, there are no > problems), then I create an ExplicitSystem and assign random values to its > solution. I print the solution for verification purposes. I save both the > equation system and the refined mesh. Then I rerun the code reading the > saved data and the saved mesh and print the solution to compare it with the > previous run. Running the code in parallel, if my mesh is SerialMesh, both > solutions coincide, however if my mesh is plain Mesh, By "plain mesh" I guess you mean ParallelMesh and that you configured libmesh with --enable-parmesh? > I get this assertion error: > > > Assertion `cnt < io_buffer.size()' failed. > > cnt = 12755 > > io_buffer.size() = 12755 > > > The program would be run with the following flags > > mpirun –n 2 ./program –uniref 1 –read_sol 0 > output1.txt > > > Rerun: > > mpirun –n 2 ./program –uniref 1 –read_sol 1 > output2.txt > > > In output2.txt, the vector printed after “initial solution” doesn’t > coincide with the vector printed after “final solution” in output1.txt if: > Mesh is used instead of SerialMesh, there’s a uniform refinement, (-uniref > 1), and it is run in parallel. I guess the issue has to do on how > ParallelMesh reads the mesh. > > It's possible there is an issue with reading XDA files with ParallelMesh in parallel... I don't think that feature gets used very often. What version of libmesh are you using? Roy recently fixed some Hilbert keys stuff, not sure if that would be related or not. Also, please send us (after making it as minimal as possible) your test code so we can try to reproduce the error? -- John ------------------------------------------------------------------------------ Site24x7 APM Insight: Get Deep Visibility into Application Performance APM + Mobile APM + RUM: Monitor 3 App instances at just $35/Month Monitor end-to-end web transactions and take corrective actions now Troubleshoot faster and improve end-user experience. Signup Now! http://pubads.g.doubleclick.net/gampad/clk?id=267308311&iu=/4140 _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users