<html>
<head>
<meta content="text/html; charset=windows-1252"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Good morning Alastair,<br>
<br>
in fact, we are currently making a parallel gmsh reader operational,<br>
and not only this but a Dune compatible tetrahedral grid manager<br>
for linear and curvilinear higher order elements. If you are
interested<br>
in testing it, you are welcome to get back to us. Greetings,
Benedikt<br>
<br>
<br>
<div class="moz-cite-prefix">On 10/03/15 22:36, Jö Fahlke wrote:<br>
</div>
<blockquote cite="mid:20150310213600.GA3333@paranoia" type="cite">
<pre wrap="">Am Tue, 10. Mar 2015, 16:57:29 +0000 schrieb Radcliffe, Alastair:
</pre>
<blockquote type="cite">
<pre wrap="">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 ) );
</pre>
</blockquote>
<pre wrap="">
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: <a class="moz-txt-link-freetext" href="http://lists.dune-project.org/pipermail/dune/2015-February/013406.html">http://lists.dune-project.org/pipermail/dune/2015-February/013406.html</a>
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ö.
</pre>
<blockquote type="cite">
<pre wrap="">but hangs in parallel runs (mpiexec -np 2 myprog) ...
Any suggestions would be most welcome.
</pre>
</blockquote>
<pre wrap="">
</pre>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
Dune mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Dune@dune-project.org">Dune@dune-project.org</a>
<a class="moz-txt-link-freetext" href="http://lists.dune-project.org/mailman/listinfo/dune">http://lists.dune-project.org/mailman/listinfo/dune</a>
</pre>
</blockquote>
<br>
</body>
</html>