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

Reply via email to