Berk, forgive for delaying my response (I was out of my lab...)<br><br>Just a silly question:<br><br>What do you exactly mean by "Ghost cell"? This silly question is justified by the fact that I thought that ghost cells only make sense when we have subdomains _overlapped_. In these cases, the number of overlapped cell layers should (or could) be considered as "ghost level". <br>
<br>In my case, the subdomains are not overlapped at all. They are just partitioned by Metis and the solver can deal with the computations without the need to use overlapped elements. Thus, if "ghost cells" can be simply defined by "cells somehow touching the parallel interface (even those not overlapped)" I think I could try tagging my parallel interface cells as "ghost level 0" to see what happens...<br>
<br>Renato.<br><br><div class="gmail_quote">On Sun, Apr 5, 2009 at 11:53 PM, Chris Kees <span dir="ltr"><<a href="mailto:christopher.e.kees@usace.army.mil">christopher.e.kees@usace.army.mil</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Berk,<br>
<br>
Thanks for correcting me. I'll try adding the vtkGhostLevels to the XDMF and see how that goes.<br><font color="#888888">
<br>
Chris</font><div><div></div><div class="h5"><br>
On Apr 4, 2009, at 2:50 PM, Berk Geveci wrote:<br>
<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
The right way to deal with this situation is to mark the ghost cells<br>
as ghost. If you create a cell array called vtkGhostLevels and assign<br>
0 to inside cells and 1 to ghost cells, you should not need<br>
MergeBlocks or CleanToGrid. Note that this array has to be of type<br>
unsigned char. There are actual benefits to keeping ghost levels since<br>
some algorithms will produce better results.<br>
<br>
-berk<br>
<br>
On Thu, Apr 2, 2009 at 6:37 PM, Chris Kees<br>
<<a href="mailto:christopher.e.kees@usace.army.mil" target="_blank">christopher.e.kees@usace.army.mil</a>> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
I turned off the overlapping domain decomposition (ghost cells) for a simple<br>
problem and the sequence<br>
<br>
MergeBlocks->CleantoGrid->ExtractSurface->Clip<br>
<br>
shows just the physical boundary of the problem (clipped open so you can see<br>
inside). Also volume visualization and streamline calculation works with no<br>
processor boundary artifacts.<br>
<br>