[Dune] GmshReader in parallel applications

[Benedikt Oswald]

Good morning Alastair,

in fact, we are currently making a parallel gmsh reader operational,
and not only this but a Dune compatible tetrahedral grid manager
for linear and curvilinear higher order elements. If you are interested
in testing it, you are welcome to get back to us. Greetings, Benedikt


On 10/03/15 22:36, Jö Fahlke wrote:
> Am Tue, 10. Mar 2015, 16:57:29 +0000 schrieb Radcliffe, Alastair:
>> Hi fellow Duners,
>>
>> Has anyone out there ever successfully used the Dune::GmshReader in
>> *parallel* executions ?
>>
>> I'm using release 2.3 of all the core modules, and the following
>> code snippet works great in serial (./myprog):
>>
>>   Dune::GridPtr<GridType> gridPtr( Dune::GmshReader<GridType>::read(
>> gridfile, true, false ) );
> Uh, I'm surprised that works.  You're read() returns *GridType (a C-Pointer to
> GridType).  GridPtr is the DGF-Parser...  (Yes, I know, the class name could
> use improvements.)  I'd simply put the result of read() into a
> shared_ptr<GridType> or unique_ptr<GridType>.
>
> Anyway, I asked about this (and more) last month, you can read the thread
> here: http://lists.dune-project.org/pipermail/dune/2015-February/013406.html
>
> Basically, you have to use the methods on the gmsh-reader that take a grid
> factory to fill.  Create the grid factory on all ranks, use the gmsh-reader to
> fill the factory on rank 0 only, then let the gridfactory create the grid on
> all ranks, and finally call loadBalance() on the new grid.  Examples of this
> are in dune-grid/src/gmsh-to-alu (for the old alugrid 1.52).
>
> This is true if you are using UG or the old alugrid 1.52.  For the new
> alugrid, please read the thread above and the alugrid paper.
>
> Regards,
> Jö.
>
>> but hangs in parallel runs (mpiexec -np 2 myprog) ...
>>
>> Any suggestions would be most welcome.
>
>
> _______________________________________________
> Dune mailing list
> Dune at dune-project.org
> http://lists.dune-project.org/mailman/listinfo/dune

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.dune-project.org/pipermail/dune/attachments/20150311/3da32c6c/attachment.htm>