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, 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.
Thanks
Miguel
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
------------------------------------------------------------------------------
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
