[Commits] [svn:einsteintoolkit] NewRad/trunk/src/ (Rev. 16)

[schnetter at cct.lsu.edu]

User: eschnett
Date: 2012/07/03 08:29 PM

Modified:
 /trunk/src/
  extrap.cc, newrad.cc

Log:
 Introduce cctk_ash, retire cctk_lssh
 
 Introduce cctk_ash, describing the process-local array shape that has
 been allocated. This may be larger than cctk_lsh, the process-local
 shape that should be used.
 
 Retire cctk_lssh and related infrastructure to handle staggered grid
 functions.

File Changes:

Directory: /trunk/src/
======================

File [modified]: extrap.cc
Delta lines: +2 -2
===================================================================
--- trunk/src/extrap.cc	2012-06-17 21:43:18 UTC (rev 15)
+++ trunk/src/extrap.cc	2012-07-04 01:29:43 UTC (rev 16)
@@ -94,7 +94,7 @@
 {
   int imin[3], imax[3], is_symbnd[6], is_physbnd[6], is_ipbnd[6];
   GenericFD_GetBoundaryInfo
-    (cctkGH, cctkGH->cctk_lsh, cctkGH->cctk_lssh, cctkGH->cctk_bbox,
+    (cctkGH, cctkGH->cctk_ash, cctkGH->cctk_lsh, cctkGH->cctk_bbox,
      cctkGH->cctk_nghostzones, 
      imin, imax, is_symbnd, is_physbnd, is_ipbnd);
   
@@ -142,7 +142,7 @@
                 break;
               case +1:
                 bmin[d] = imax[d];
-                bmax[d] = cctkGH->cctk_lssh[CCTK_LSSH_IDX(0,d)];
+                bmax[d] = cctkGH->cctk_lsh[d];
                 have_bnd = true;
                 all_physbnd = all_physbnd and is_physbnd[2*d+1];
                 any_physbnd = any_physbnd or is_physbnd[2*d+1];

File [modified]: newrad.cc
Delta lines: +2 -2
===================================================================
--- trunk/src/newrad.cc	2012-06-17 21:43:18 UTC (rev 15)
+++ trunk/src/newrad.cc	2012-07-04 01:29:43 UTC (rev 16)
@@ -282,7 +282,7 @@
 {
   int imin[3], imax[3], is_symbnd[6], is_physbnd[6], is_ipbnd[6];
   GenericFD_GetBoundaryInfo
-    (cctkGH, cctkGH->cctk_lsh, cctkGH->cctk_lssh, cctkGH->cctk_bbox,
+    (cctkGH, cctkGH->cctk_ash, cctkGH->cctk_lsh, cctkGH->cctk_bbox,
      cctkGH->cctk_nghostzones, 
      imin, imax, is_symbnd, is_physbnd, is_ipbnd);
   
@@ -330,7 +330,7 @@
                 break;
               case +1:
                 bmin[d] = imax[d];
-                bmax[d] = cctkGH->cctk_lssh[CCTK_LSSH_IDX(0,d)];
+                bmax[d] = cctkGH->cctk_lsh[d];
                 have_bnd = true;
                 all_physbnd = all_physbnd and is_physbnd[2*d+1];
                 any_physbnd = any_physbnd or is_physbnd[2*d+1];