<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><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 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 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 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 href="mailto:Dune@dune-project.org">Dune@dune-project.org</a><br>http://lists.dune-project.org/mailman/listinfo/dune<br></blockquote></div><br></body></html>