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];