[Dune] Topological changes in the grid

[Sacconi, Andrea]

Hi everyone,

I'd like to ask you a question about the grid.

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.
Suppose in addition that you have a software that removes tangles, returning the actual disconnected domain without self-intersections.

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.
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?

I mean, the new tangle-free grid is returned inside a time loop.
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 ........

Cheers,
Andrea

__________________________________________________________

Andrea Sacconi
PhD student, Applied Mathematics
AMMP Section, Department of Mathematics, Imperial College London,
London SW7 2AZ, UK
a.sacconi11 at imperial.ac.uk
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.dune-project.org/pipermail/dune/attachments/20130219/8426d923/attachment.htm>