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

[schnetter at cct.lsu.edu]

User: eschnett
Date: 2013/01/14 11:53 AM

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

Log:
 Support array padding (take cctk_ash into account)

File Changes:

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

File [modified]: extrap.cc
Delta lines: +7 -3
===================================================================
--- trunk/src/extrap.cc	2013-01-10 23:25:45 UTC (rev 18)
+++ trunk/src/extrap.cc	2013-01-14 17:53:39 UTC (rev 19)
@@ -24,14 +24,18 @@
   int const nj = cctkGH->cctk_lsh[1];
   int const nk = cctkGH->cctk_lsh[2];
   
+  int const ai = cctkGH->cctk_ash[0];
+  int const aj = cctkGH->cctk_ash[1];
+  int const ak = cctkGH->cctk_ash[2];
+  
   int const si = -dir[0];
   int const sj = -dir[1];
   int const sk = -dir[2];
   
   int const di = 1;
-  int const dj = ni;
-  int const dk = ni*nj;
-  int const np = ni*nj*nk;
+  int const dj = ai;
+  int const dk = ai*aj;
+  int const np = ai*aj*ak;
   
   int const dind = si*di + sj*dj + sk*dk;
   

File [modified]: newrad.cc
Delta lines: +6 -2
===================================================================
--- trunk/src/newrad.cc	2013-01-10 23:25:45 UTC (rev 18)
+++ trunk/src/newrad.cc	2013-01-14 17:53:39 UTC (rev 19)
@@ -32,13 +32,17 @@
   int const nj = cctkGH->cctk_lsh[1];
   int const nk = cctkGH->cctk_lsh[2];
   
+  int const ai = cctkGH->cctk_ash[0];
+  int const aj = cctkGH->cctk_ash[1];
+  int const ak = cctkGH->cctk_ash[2];
+  
   int const si = dir[0];
   int const sj = dir[1];
   int const sk = dir[2];
   
   int const di = 1;
-  int const dj = ni;
-  int const dk = ni*nj;
+  int const dj = ai;
+  int const dk = ai*aj;
   
   CCTK_REAL const dx = cctkGH->cctk_delta_space[0] / cctkGH->cctk_levfac[0];
   CCTK_REAL const dy = cctkGH->cctk_delta_space[1] / cctkGH->cctk_levfac[1];