<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Hi Eike,<br>
the refinement rule that you want is HEX_QUADSECT_0.<br>
There was a bit of implementation missing in UG itself. Please<br>
try the attached patch -- it works for me.<br>
best,<br>
Oliver<br>
<br>
Am 05.11.2012 20:24, schrieb Eike Mueller:
<blockquote
cite="mid:8DE564E1-7302-43E7-AC07-D781B0BF9B6A@bath.ac.uk"
type="cite">
<div>Hi Oliver,</div>
<div><br>
</div>
<div>no, it either crashed or did nothing. I attach the updated
code.</div>
<div><br>
</div>
<div>Thanks, Eike</div>
<div><br>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<br>
<div>
<div>On 5 Nov 2012, at 19:22, Oliver Sander wrote:</div>
<br class="Apple-interchange-newline">
<blockquote type="cite">
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
<div bgcolor="#FFFFFF" text="#000000"> Hi Eike,<br>
as mentioned, this is a hardly tested feature. Can you send
me an updated test<br>
program? Did you find any rule (besides RED) that does
anything?<br>
best,<br>
Oliver<br>
<br>
Am 05.11.2012 20:14, schrieb Eike Mueller:
<blockquote
cite="mid:5599230B-C696-4044-AD8A-3E021C776F63@bath.ac.uk"
type="cite">
<div>Hi Oliver,</div>
<div><br>
</div>
<div>thanks a lot, I think that's what I was looking for.
I've now tried the following RefinementRules, but with
little success:</div>
<div><br>
</div>
<div>RED: I get the previous refinement in all three
directions</div>
<div>BLUE, COARSE: Nothing happens, i.e. elements do not
get refined</div>
<div><br>
</div>
<div>For</div>
<div><br>
</div>
<div>
<div>HEX_BISECT_0_1,</div>
<div>HEX_BISECT_0_2,</div>
<div>HEX_BISECT_0_3,</div>
<div>HEX_TRISECT_0,</div>
<div>HEX_TRISECT_5,</div>
<div>HEX_QUADSECT_0,</div>
<div>HEX_QUADSECT_1,</div>
<div>HEX_QUADSECT_2</div>
</div>
<div><br>
</div>
<div>I always get the following runtime error:</div>
<div><br>
</div>
<div>ERROR in Patterns2Rules: no mapping for HEXAHEDRON
and this pattern!</div>
<div>uggrid: rm.c:2838: UG::INT
UG::D3::Patterns2Rules(UG::D3::ELEMENT*, UG::INT):
Assertion `0' failed.</div>
<div><br>
</div>
<div>so I must still be missing something. I loop over the
grid like this:</div>
<div><br>
</div>
<div>
<div>for (ElementLeafIterator it = ibegin;it!=iend;++it)
{</div>
<div> grid.mark(*it,UG::D3::PRISM_QUADSECT,0);</div>
<div>}</div>
</div>
<div><br>
</div>
<div>Eike</div>
<div><br>
</div>
<div>
<div>On 5 Nov 2012, at 17:01, Oliver Sander wrote:</div>
<br class="Apple-interchange-newline">
<blockquote type="cite">
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
<div bgcolor="#FFFFFF" text="#000000"> <br>
<br>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
Hi Eike,<br>
<br>
Am 05.11.2012 17:45, schrieb Eike Mueller:
<blockquote
cite="mid:89792302-805B-46AF-BCCC-6094120908EA@bath.ac.uk"
type="cite">
<pre wrap="">Dear dune-list,
(1) How can I ensure that the grid is not refined in the radial direction, i.e. I always only end up with one layer of cells in the radial direction? I guess I have to run over the grid and mark cells for refinement, but is this documented anywhere in more detail?</pre>
</blockquote>
UGGrid has a special non-interface method for
anisotropic refinement. Have a look<br>
at the doxygen documentation of the UGGrid class.
There is a method<br>
<br>
bool <a moz-do-not-send="true" class="el"
href="imap://sander@indy5.igpm.rwth-aachen.de:143/fetch%3EUID%3E.INBOX.Sent%3E209?part=1.2.2&filename=a00865.html">mark</a>
(const typename Traits::template <a
moz-do-not-send="true" class="el"
href="imap://sander@indy5.igpm.rwth-aachen.de:143/fetch%3EUID%3E.INBOX.Sent%3E209?part=1.2.3&filename=a00218.html">Codim</a><
0 >::<a moz-do-not-send="true" class="el"
href="imap://sander@indy5.igpm.rwth-aachen.de:143/fetch%3EUID%3E.INBOX.Sent%3E209?part=1.2.4&filename=a00362.html">Entity</a>
&e, typename UG_NS< dim >::RefinementRule
rule, int side=0)<br>
<br>
which marks the entity 'e' for refinement, but with
a special refinement rule.<br>
The list of rules is UG-specific and has not been
properly documented (Freiwillige vor!).<br>
You find all possible enum values in the UG header
ug/gm.h at line 317ff<br>
They are reasonably self-explanatory, but you need
to do some experimenting.<br>
<br>
<blockquote
cite="mid:89792302-805B-46AF-BCCC-6094120908EA@bath.ac.uk"
type="cite">
<pre wrap="">(2) The parallel load balancing does not seem to work. If I run on 24 processes and I call globalRefine(), followed by loadBalance(), some processes end up with no part of the grid. I have checked that I use the parallel version of UG grid.
I attach the code I used to generate the grid, as well as the .vtu files mentioned above.</pre>
</blockquote>
No clue, I'd have to try this out myself.<br>
cheers,<br>
Oliver<br>
<br>
<blockquote
cite="mid:89792302-805B-46AF-BCCC-6094120908EA@bath.ac.uk"
type="cite">
<pre wrap="">Thanks a lot for any help,
Eike
PS: Alternatively, I would also be equaully happy if I can create a PARALLEL 2d surface grid for the surface of the sphere (again, either cubed sphere or a icosahedral grid). However, as far as I can see it, there is no parallel implementation for this, we've tried 2d ALUGRID, parallelised with the metagrid approach, but this does not scale beyond a few hundred processors as the entire grid is stored on each processor and my code runs out of memory.
</pre>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
Dune mailing list
<a moz-do-not-send="true" class="moz-txt-link-abbreviated" href="mailto:Dune@dune-project.org">Dune@dune-project.org</a>
<a moz-do-not-send="true" 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>
<br>
</div>
_______________________________________________<br>
Dune mailing list<br>
<a moz-do-not-send="true"
href="mailto:Dune@dune-project.org">Dune@dune-project.org</a><br>
<a moz-do-not-send="true"
class="moz-txt-link-freetext"
href="http://lists.dune-project.org/mailman/listinfo/dune">http://lists.dune-project.org/mailman/listinfo/dune</a><br>
</blockquote>
</div>
<br>
</blockquote>
<br>
</div>
</blockquote>
</div>
<br>
</blockquote>
<br>
</body>
</html>