[Carpet] missing data for grid component

[Bela Szilagyi]

Giovanni ,

I believe in Carpet you need to give some more space between the different 
refinement level boundaries.  This depends on the type of buffer zones you 
use and the actual width of your buffer size as well.

What you need to keep in mind is that a number of points (buffer size + ghost 
size) at the edge of any refinement level are interpolated onto from the 
coarser level, using an interpolation stencil that is set by the order of 
prolongation in space.  That is, the boundary points of any refinement level 
have to be well inside the next coarser level.  

On a very practical level, this would translate into, say,  outer boundaries 
at 120 (with grid-step 0f 2), the next refinement level boundaries at 60 
(grid-step of 1) and then the finest level at 30, as you have it.  This way 
there are 30 grid-points between the boundary of a refinement level and the 
next coarser level which , I believe, should be more than enough.

Hope this helps.

Bela.



On Tuesday 06 March 2007 12:59, Giovanni Corvino wrote:
> Hi, I have a problem with a compnent of my grid.
>
> I'm running with 2 refinement levels like
>
> [ [ ([-40.0,-40.0,0.0]:[40.0,40.0,18.0]:[1,1,1]) ],
>  [ ([-30.0,-30.0,0.0]:[30,30,8]:[0.5,0.5,0.5]) ] ]"
>
> the outer boundary is at 50 (bitant symmetry).
>
> First of all I get the following warning while generating the initial
> conditions:
>
> INFO (Carpet): Recomposing the grid hierarchy for map 0...
> Not properly nested: rl=1, required distance=4
> Received from level 0:
> bboxset<CCTK_INT4,3>:size=113049,setsize=5,set={([14,14,6]:[30,186,54]:[2,2
>,2]),([32,14,6]:[186,30,54]:[2,2,2]),([170,32,6]:[186,186,54]:[2,2,2]),([32,
>170,6]:[168,186,54]:[2,2,2]),([32,32,38]:[168,168,54]:[2,2,2])} Taboo region
> on level 1:
> bboxset<CCTK_INT4,3>:size=238289,setsize=5,set={([6,6,-2]:[38,194,62]:[2,2,
>2]),([40,6,-2]:[194,38,62]:[2,2,2]),([162,40,-2]:[194,194,62]:[2,2,2]),([40,
>162,-2]:[160,194,62]:[2,2,2]),([40,40,30]:[160,160,62]:[2,2,2])} Sent to
> level 2:
> bboxset<CCTK_INT4,3>:size=35262,setsize=5,set={([36,36,8]:[46,164,32]:[2,2,
>2]),([48,36,8]:[164,46,32]:[2,2,2]),([154,48,8]:[164,164,32]:[2,2,2]),([48,1
>54,8]:[152,164,32]:[2,2,2]),([48,48,22]:[152,152,32]:[2,2,2])} Overlap on
> level 1:
> bboxset<CCTK_INT4,3>:size=13994,setsize=5,set={([36,36,8]:[38,164,32]:[2,2,
>2]),([40,36,8]:[164,38,32]:[2,2,2]),([162,40,8]:[164,164,32]:[2,2,2]),([40,1
>62,8]:[160,164,32]:[2,2,2]),([40,40,30]:[160,160,32]:[2,2,2])}
>
>
> and then wheh I look at the data I see that one or more of the
> components of the grid has not been computed, on all refinement levels.
> Sometimes is only one component, sometimes two. What can it depend on?
> Caould it be something tied to the boudary conditions of the boxes?
>
> Thanks, Giovanni Corvino
>
>
> _______________________________________________
> developers mailing list
> developers at lists.carpetcode.org
> http://lists.carpetcode.org/listinfo/developers

-- 
Bela Szilagyi
----------------------------------------------------
Max-Planck-Institut für Gravitationsphysik
Albert-Einstein-Institut
Tel: +49 331 567 7189
Fax: +49 331 567 7252
----------------------------------------------------