<div dir="ltr">can be useful to use some marker or flag ?</div><div class="gmail_extra"><br><br><div class="gmail_quote">2013/2/20 <span dir="ltr"><<a href="mailto:dune-request@dune-project.org" target="_blank">dune-request@dune-project.org</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Send Dune mailing list submissions to<br>
<a href="mailto:dune@dune-project.org">dune@dune-project.org</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="http://lists.dune-project.org/mailman/listinfo/dune" target="_blank">http://lists.dune-project.org/mailman/listinfo/dune</a><br>
or, via email, send a message with subject or body 'help' to<br>
<a href="mailto:dune-request@dune-project.org">dune-request@dune-project.org</a><br>
<br>
You can reach the person managing the list at<br>
<a href="mailto:dune-owner@dune-project.org">dune-owner@dune-project.org</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of Dune digest..."<br>
<br>Today's Topics:<br>
<br>
1. Topological changes in the grid (Sacconi, Andrea)<br>
2. Re: Topological changes in the grid (Oliver Sander)<br>
<br><br>---------- Messaggio inoltrato ----------<br>From: "Sacconi, Andrea" <<a href="mailto:a.sacconi11@imperial.ac.uk">a.sacconi11@imperial.ac.uk</a>><br>To: "<a href="mailto:dune@dune-project.org">dune@dune-project.org</a>" <<a href="mailto:dune@dune-project.org">dune@dune-project.org</a>><br>
Cc: <br>Date: Tue, 19 Feb 2013 18:01:03 +0000<br>Subject: [Dune] Topological changes in the grid<br>
<div>
<div style="direction:ltr;font-size:10pt;font-family:Tahoma">
<div style="direction:ltr;font-size:10pt;font-family:Tahoma">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 href="mailto:a.sacconi11@imperial.ac.uk" target="_blank">a.sacconi11@imperial.ac.uk</a><br>
</div>
</div>
</div>
</div>
</div>
<br><br>---------- Messaggio inoltrato ----------<br>From: Oliver Sander <<a href="mailto:sander@igpm.rwth-aachen.de">sander@igpm.rwth-aachen.de</a>><br>To: <a href="mailto:dune@dune-project.org">dune@dune-project.org</a><br>
Cc: <br>Date: Tue, 19 Feb 2013 19:32:25 +0100<br>Subject: Re: [Dune] Topological changes in the grid<br>
<div 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 type="cite">
<div style="direction:ltr;font-size:10pt;font-family:Tahoma">
<div style="direction:ltr;font-size:10pt;font-family:Tahoma">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 href="mailto:a.sacconi11@imperial.ac.uk" target="_blank">a.sacconi11@imperial.ac.uk</a><br>
</div>
</div>
</div>
</div>
<br>
<fieldset></fieldset>
<br>
<pre>_______________________________________________
Dune mailing list
<a href="mailto:Dune@dune-project.org" target="_blank">Dune@dune-project.org</a>
<a href="http://lists.dune-project.org/mailman/listinfo/dune" target="_blank">http://lists.dune-project.org/mailman/listinfo/dune</a>
</pre>
</blockquote>
<br>
</div>
<br>_______________________________________________<br>
Dune mailing list<br>
<a href="mailto:Dune@dune-project.org">Dune@dune-project.org</a><br>
<a href="http://lists.dune-project.org/mailman/listinfo/dune" target="_blank">http://lists.dune-project.org/mailman/listinfo/dune</a><br>
<br></blockquote></div><br></div>