Dear David,
Dear libMesh developers,
I have some problems with the new version
of RBEIMConstruction/RBEIMEvaluation AND ParallelMesh.
I guess that the problem was introduced with the following commit
https://github.com/libMesh/libmesh/commit/474509c0d10ff0feb7292410b3cfdd91b743310d
In particular, at line 346 of src/reduced_basis/rb_eim_construction.C,
the optimal_elem_id is broadcast among the processors, in order for the
optimal element to be retrieved at line 356. However, when working in
parallel on a ParallelMesh, the element may not be a local one, causing the
assertion at line 318 of src/mesh/parallel_mesh.C to fail. I think that a
copy of the optimal element should be broadcast between all the processors,
possibly using a procedure similar to the one in MeshSerializer (only for
the optimal element, not for the entire mesh), but my expertise with the
ParallelMesh class is not sufficient for me to implement a bugfix on my own.
Do you have any suggestions?
Thanks,
best regards,
Francesco
------------------------------------------------------------------------------
CenturyLink Cloud: The Leader in Enterprise Cloud Services.
Learn Why More Businesses Are Choosing CenturyLink Cloud For
Critical Workloads, Development Environments & Everything In Between.
Get a Quote or Start a Free Trial Today.
http://pubads.g.doubleclick.net/gampad/clk?id=119420431&iu=/4140/ostg.clktrk
_______________________________________________
Libmesh-devel mailing list
Libmesh-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-devel
