rought draft of new version to call Cactus local interp

-- still needs more work


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@408 f88db872-0e4f-0410-b76b-b9085cfa78c5

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_;
 	};