Fix a bug where N_overlap_points was being interpreted with different

semantics (a factor of 2) in different places:
N_overlap_points is now cleanly defined in patch.hh (previous checkin)
N_extend_points is now used for how much to extend a single side


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

1 files changed, 14 insertions, 6 deletions

diff --git a/src/patch/patch_system.hh b/src/patch/patch_system.hh
index 5531925..3044f3e 100644
--- a/src/patch/patch_system.hh
+++ b/src/patch/patch_system.hh
@@ -59,7 +59,7 @@ struct	patch_info
 	// ... result refers to internal static buffer;
 	//     the usual caveats about lifetimes/overwriting apply
 	const grid_arrays::grid_array_pars& grid_array_pars
-		(int N_ghost_points, int N_overlap_points, fp delta_drho_dsigma)
+		(int N_ghost_points, int N_extend_points, fp delta_drho_dsigma)
 		const;
 
 	// compute and return reference to  struct grid::grid_pars
@@ -67,7 +67,7 @@ struct	patch_info
 	// ... result refers to internal static buffer;
 	//     the usual caveats about lifetimes/overwriting apply
 	const grid::grid_pars& grid_pars
-		(int N_overlap_points, fp delta_drho_dsigma)
+		(int N_extend_points, fp delta_drho_dsigma)
 		const;
 	};
 
@@ -275,6 +275,10 @@ public:
 	// from patch to patch.  But in practice we'd probably never
 	// use that generality...
 	//
+	// The present implementation also restricts N_overlap_points
+	// to be even.  N.b. the semantics of N_overlap points are
+	// as documented in patch.hh for patch::edge_adjacent_to_patch().
+	//
 	patch_system(fp origin_x_in, fp origin_y_in, fp origin_z_in,
 		     enum patch_system_type type_in,
 		     int N_ghost_points, int N_overlap_points,
@@ -290,19 +294,23 @@ private:
 
 	// create/interlink patches, ghost zones, and frontiers
 	void setup_full_sphere_patch_system
-		(int N_ghost_points, int N_overlap_points, fp delta_drho_dsigma,
+		(int N_ghost_points, int N_extend_points,
+		 int N_overlap_points, fp delta_drho_dsigma,
 		 int min_gfn_in, int max_gfn_in,
 		 int interpolator_order);
 	void setup_hemisphere_patch_system
-		(int N_ghost_points, int N_overlap_points, fp delta_drho_dsigma,
+		(int N_ghost_points, int N_extend_points,
+		 int N_overlap_points, fp delta_drho_dsigma,
 		 int min_gfn_in, int max_gfn_in,
 		 int interpolator_order);
 	void setup_quadrant_patch_system
-		(int N_ghost_points, int N_overlap_points, fp delta_drho_dsigma,
+		(int N_ghost_points, int N_extend_points,
+		 int N_overlap_points, fp delta_drho_dsigma,
 		 int min_gfn_in, int max_gfn_in,
 		 int interpolator_order);
 	void setup_octant_patch_system
-		(int N_ghost_points, int N_overlap_points, fp delta_drho_dsigma,
+		(int N_ghost_points, int N_extend_points,
+		 int N_overlap_points, fp delta_drho_dsigma,
 		 int min_gfn_in, int max_gfn_in,
 		 int interpolator_order);