<html>
  <head>
    <meta content="text/html; charset=ISO-8859-1"
      http-equiv="Content-Type">
  </head>
  <body bgcolor="#FFFFFF" text="#000000">
    Hi Andrea,<br>
    there is no way to change the connectivity of any Dune grid
    implementation that<br>
    I know.  Not even a dirty trick.  You need to kill your grid object
    within the time loop<br>
    and create a new one using the grid factory.  That by itself is
    reasonably straightforward.<br>
    The problem of course is how to transfer your data from the old to
    the new grid.<br>
    But the solution to that is problem-specific.  <br>
    good luck,<br>
    Oliver<br>
    <br>
    Am 19.02.2013 19:01, schrieb Sacconi, Andrea:
    <blockquote
      cite="mid:EC7DA3C0917A914BAFB620BB60189EBD1361D774@icexch-m2.ic.ac.uk"
      type="cite">
      <meta http-equiv="Content-Type" content="text/html;
        charset=ISO-8859-1">
      <style id="owaParaStyle" type="text/css">
<!--
p
        {margin-top:0;
        margin-bottom:0}
-->
P {margin-top:0;margin-bottom:0;}</style>
      <div style="direction: ltr;font-family: Tahoma;color:
        #000000;font-size: 10pt;">
        <div style="direction:ltr; font-family:Tahoma; color:#000000;
          font-size:10pt">Hi everyone,<br>
          <br>
          I'd like to ask you a question about the grid.<br>
          <br>
          Suppose you are solving a problem on a given initial grid, and
          after a while some topological changes take place
          (coalescence, breakups, ...). For example, from one connected
          domain you get two, smaller disconnected ones.<br>
          Suppose in addition that you have a software that removes
          tangles, returning the actual disconnected domain without
          self-intersections.<br>
          <br>
          I know that a grid can be neither created nor modified; so
          basically you can only have a reference to a grid, initialised
          by dereferencing a pointer to the the same grid type.<br>
          So, if all the objects (function spaces, grid functions, ...)
          I have in my algorithm are defined over the initial one, is
          there a tricky way to continue my computations on the new
          disconnected domain?<br>
          <br>
          I mean, the new tangle-free grid is returned inside a time
          loop.<br>
          I have initialised a new pointer to the same grid type by
          using a grid factory (since the new disconnected domain is
          returned through its connectivity, I insert vertices and
          elements by hand). However, I can't assign it to the reference
          to the old grid ... Please tell me that there is a (dirty)
          workaround to continue the computations inside the time loop
          ........... otherwise I wouldn't know how to proceed without
          breaking the loop ........<br>
          <br>
          Cheers,<br>
          <div>Andrea<br>
            <div style="font-family:Tahoma; font-size:13px"><br>
              __________________________________________________________<br>
              <br>
              Andrea Sacconi<br>
              PhD student, Applied Mathematics<br>
              AMMP Section, Department of Mathematics, Imperial College
              London,<br>
              London SW7 2AZ, UK<br>
              <a class="moz-txt-link-abbreviated" href="mailto:a.sacconi11@imperial.ac.uk">a.sacconi11@imperial.ac.uk</a><br>
            </div>
          </div>
        </div>
      </div>
      <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>