diff options
Diffstat (limited to 'src/patch/patch_interp.hh')
-rw-r--r-- | src/patch/patch_interp.hh | 20 |
1 files changed, 15 insertions, 5 deletions
diff --git a/src/patch/patch_interp.hh b/src/patch/patch_interp.hh index ba62792..2c92da2 100644 --- a/src/patch/patch_interp.hh +++ b/src/patch/patch_interp.hh @@ -152,6 +152,7 @@ private: // min parindex at each iperp // ... this is a reference to an array passed in to our constructor // ==> we do *not* own this! + // ... index is (iperp) const jtutil::array1d<int>& min_parindex_used_; // number of points to interpolate at each iperp @@ -187,6 +188,8 @@ private: fp gridfn_coord_origin_, gridfn_coord_delta_; // + // We use (the 0-origin) parsub as an interpolator point index + // within each interpolation (i.e. at each value of iperp). // The interpolator accesses the gridfn data with the generic // subscripting expression // data[offset + parsub*stride] @@ -196,13 +199,20 @@ private: // Our strategy will be to leave offset unspecified so it // defaults to 0, specify a common stride for all (gfn,iperp), // and specify the appropriate data pointer for each (gfn,iperp). + // Thus data should point to the start of the gridfn data + // we should use for each iperp // - // --> start of gridfn data we can use for each iperp - // subscripts are (gfn,iperp) + // --> start of gridfn data we should use for each iperp + // subscripts are (iperp,gfn) so [gfn] is contiguous at each iperp + // ... n.b. we do *not* own this data! jtutil::array2d<fp *> gridfn_data_ptrs_; - // stride of par coordinate in gridfn data - // (same for all gfn and iperp) - int gridfn_par_stride_; + // CCTK_VARIABLE_ type codes for the gridfns + // ... all are actually CCTK_REAL, but we need the array for + // the interpolator + // ... datatype is CCTK_INT so we can pass this by reference + // to the interpolator + // ... index is (ghosted_gfn) + jtutil::array1d<CCTK_INT> gridfn_type_codes_; }; |