[Dune] Assertion error when running Dune greed in parallel

[Guichard, Roland]

Dear Markus,

Thank you for your prompt response.

Sure, you can find attached a stripped out version of the main. I only have inherited that code, not developed it myself.

You’ll see that the same path to the dgf file is passed to all MPI ranks but I don’t understand why initialising the DGFGridFactory is rank dependent.

Besides, any comments/suggestions on the code are very welcome. Part of my job is also to refactor it for a better use. The main issue is that there is such a gap between Dune2.3 and Dune2.6 which makes it uneasy to track (owing that I have never used Dune before :)

Thank you again for your time

Kind regards,

Dr. Roland Guichard
Research Software Engineer
UCL-RITS
Internal Extension: 86947
External Number: 02031086947







On 5 Jul 2019, at 12:00, Markus Blatt <markus at dr-blatt.de<mailto:markus at dr-blatt.de>> wrote:

Hi,

would you mind to share your setup code for the grid?

Markus

On Fri, Jul 05, 2019 at 09:24:33AM +0000, Guichard, Roland wrote:
Dear Dune people,

I am currently working on a project to update a home made Dune module (peits) to the newest Dune modules. Peits was previously running fine with Dune modules (common, fem, grid, alugrid, geometry) v2.3.

Now I am trying to update to the same Dune modules v2.6. I have managed to compile (under GNU compiler) and run serially Peits with Dune2.6 but the issue appears when I run in parallel using MPI. This occurs when the gridPtr is called,

And trying to create the grid from a dgf file, especially the assertion:

https://gitlab.dune-project.org/core/dune-grid/blob/master/dune/grid/io/file/dgfparser/gridptr.hh#L281

Basically, what I have narrowed down is:


 *   The MPI rank 0 runs fine. It fetches all parameters from files, calls the DGFGridFactory, initializes it etc.
 *   The MPI rank 1 comes to the same point but I have checked that the initialisation of the DGFGridFactory is not setting the necessary parameters the same way, I.e. the elParam_ size is set to 0 (whereas ~57000 in MPI rank 0) causing the above assertion to fail.

I wanted to know if there is a dedicated procedure to create a grid in parallel which I am missing at this point. Could you advise on this issue ?

That would be really helpful.

Kind regards,

Dr. Roland Guichard
Research Software Engineer
UCL-RITS
Internal Extension: 86947
External Number: 02031086947








_______________________________________________
Dune mailing list
Dune at lists.dune-project.org<mailto:Dune at lists.dune-project.org>
https://lists.dune-project.org/mailman/listinfo/dune


--
Using PT-Scotch with DUNE https://t.co/hZH8lxf5ZF

Dr. Markus Blatt - HPC-Simulation-Software & Services http://www.dr-blatt.de
Pedettistr. 38, 85072 Eichstätt, Germany
Tel.: +49 (0) 160 97590858

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.dune-project.org/pipermail/dune/attachments/20190705/9c11ee3a/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: peits-main.cc
Type: application/octet-stream
Size: 4546 bytes
Desc: peits-main.cc
URL: <https://lists.dune-project.org/pipermail/dune/attachments/20190705/9c11ee3a/attachment.obj>