My code is configured with —enable-parmesh. Please find the code here: http://paste.ofcode.org/XqU3mKRdgGnQQZRveVwKRr and the l_problem.xda file: http://justpaste.it/qe3i. The libMesh version I’m using is older than the Hilbert keys update. I will update it and see what happens.
Thanks, Miguel A. Salazar de Troya Graduate Scholar, Lawrence Livermore National Laboratory Ph.D. Candidate, University of Illinois at Urbana-Champaign T1878 Rm: 112 Ext: 2-6411 L-code: 288 From: John Peterson <jwpeter...@gmail.com<mailto:jwpeter...@gmail.com>> Date: Tuesday, January 12, 2016 at 10:27 AM To: Miguel Salazar <salazardet...@llnl.gov<mailto:salazardet...@llnl.gov>> Cc: "libmesh-users@lists.sourceforge.net<mailto:libmesh-users@lists.sourceforge.net>" <libmesh-users@lists.sourceforge.net<mailto:libmesh-users@lists.sourceforge.net>> Subject: Re: [Libmesh-users] ParallelMesh, SerialMesh and simulation restarts On Tue, Jan 12, 2016 at 11:17 AM, Salazar De Troya, Miguel <salazardet...@llnl.gov<mailto: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
