Re: [Users] Howto? Outer loop over coarse maps, inner loop over fine maps

Ok. What I'm trying to do is compute the inter-grid truncation error at every point (where it's defined, on the coarse grid):

2h 2h 2h h 2h h h tau := L I u - I L u h h h

The first term on the right hand side is the one for which I need to find the intersection of bounding boxesŠ

Skyping tomorrow is fine.

On 4/15/12 5:54 PM, "Erik Schnetter" eschnetter@perimeterinstitute.ca wrote:

Scott

In Carpet, fine boxes don't necessarily overlap a particular coarse box. Only the set of fine boxes overlaps (is contained in) the set of coarse boxes.

Also, the meaning of cctk_bbox and what Carpet calls a bbox are distinct. cctk_bbox contains booleans indicating whether the current regions contains a piece of the boundary (or whether it is in the interior of the grid), while Carpet's bboxes contain (integer) coordinates describing the rectangles (but in 3D) that make up the refined regions.

I don't quite understand what you are trying to do. Should we Skype tomorrow?

-erik

On Sun, Apr 15, 2012 at 4:47 PM, Scott Hawley scott.hawley@belmont.edu wrote:

...

Erik,

How does one "nest" loops over grids at different levels? In my way of doing things, it was something like...

l = level;

gblm1 = mgh[l-1]; while (gblm1 != NULL) { Bbox coarse_bbox = gblm1->bbox;

   gbl = mgh[l];
   while (gbl != NULL) {
     Bbox fine_bbox = gbl->bbox;

     // On the regions for which the finer grid overlaps the coarse

grid, // do some operations...

    gbl = gbl->next;
  }

gblm1 = gblm1->next; }

I have an idea of something to try, below, but I'm worried that I'm abusing the SWITCH_TO_LEVEL macroŠ

//reflevel = finer gridŠ SWITCH_TO_LEVEL( cctkGH, reflevel - 1 ) { BEGIN_MAP_LOOP (cctkGH, CCTK_GF) { BEGIN_LOCAL_COMPONENT_LOOP (cctkGH, CCTK_GF) { DECLARE_CCTK_ARGUMENTS;

                 SWITCH_TO_LEVEL( cctkGH, reflevel ) {
                       BEGIN_MAP_LOOP (cctkGH, CCTK_GF) {
                              BEGIN_LOCAL_COMPONENT_LOOP (cctkGH,

CCTK_GF) { DECLARE_CCTK_ARGUMENTS;

                                   // where fine grid bbox overlaps

coarse bboxŠ. // do operations on coarse gridŠ SWITCH_TO_LEVEL( cctkGH, reflevel - 1 ) { // do operationsŠ, e.g. call "L" operator } END_SWITCH_TO_LEVEL;

                             } END_LOCAL_COMPONENT_LOOP;
                      } END_MAP_LOOP;
               } END_SWITCH_TO_LEVEL;


   } END_LOCAL_COMPONENT_LOOP;
} END_MAP_LOOP;

} END_SWITCH_TO_LEVEL;

Or?

-- Erik Schnetter eschnetter@perimeterinstitute.ca http://www.perimeterinstitute.ca/personal/eschnetter/ AIM: eschnett247, Skype: eschnett, Google Talk: schnetter@gmail.com