add setup_gridfn_storage() calls in ctor,

prototype in class declaration (.hh)
--> still need to write the actual function (.cc)
VS: ----------------------------------------------------------------------


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

2 files changed, 18 insertions, 4 deletions

diff --git a/src/patch/patch_system.cc b/src/patch/patch_system.cc
index 969adaf..66d9ef6 100644
--- a/src/patch/patch_system.cc
+++ b/src/patch/patch_system.cc
@@ -71,8 +71,8 @@ patch_system::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_extend_points,
 			   fp delta_drho_dsigma,
-			   int min_gfn_in, int max_gfn_in,
-			   int ghosted_min_gfn_in, int ghosted_max_gfn_in,
+			   int min_gfn, int max_gfn,
+			   int ghosted_min_gfn, int ghosted_max_gfn,
 			   int interp_handle_in, int interp_par_table_handle_in)
 
 	: global_coords_(origin_x_in, origin_y_in, origin_z_in),
@@ -89,6 +89,8 @@ case full_sphere_patch_system:
 	construct_patches(patch_system_info::full_sphere::patch_info_array,
 			  N_ghost_points, N_extend_points,
 			  delta_drho_dsigma);
+	setup_gridfn_storage(min_gfn, max_gfn,
+			     ghosted_min_gfn, ghosted_max_gfn);
 	interlink_full_sphere_patch_system(N_ghost_points, N_extend_points,
 				       N_overlap_points, delta_drho_dsigma,
 				       min_gfn_in, max_gfn_in,
@@ -99,6 +101,8 @@ case plus_z_hemisphere_patch_system:
 					   ::patch_info_array,
 			  N_ghost_points, N_extend_points,
 			  delta_drho_dsigma);
+	setup_gridfn_storage(min_gfn, max_gfn,
+			     ghosted_min_gfn, ghosted_max_gfn);
 	setup_plus_z_hemisphere_patch_system
 		(N_ghost_points, N_extend_points,
 		 N_overlap_points, delta_drho_dsigma,
@@ -109,6 +113,8 @@ case plus_xy_quadrant_patch_system:
 	construct_patches(patch_system_info::plus_xy_quadrant::patch_info_array,
 			  N_ghost_points, N_extend_points,
 			  delta_drho_dsigma);
+	setup_gridfn_storage(min_gfn, max_gfn,
+			     ghosted_min_gfn, ghosted_max_gfn);
 	setup_plus_xy_quadrant_patch_system(N_ghost_points, N_extend_points,
 					    N_overlap_points, delta_drho_dsigma,
 					    min_gfn_in, max_gfn_in,
@@ -118,6 +124,8 @@ case plus_xz_quadrant_patch_system:
 	construct_patches(patch_system_info::plus_xz_quadrant::patch_info_array,
 			  N_ghost_points, N_extend_points,
 			  delta_drho_dsigma);
+	setup_gridfn_storage(min_gfn, max_gfn,
+			     ghosted_min_gfn, ghosted_max_gfn);
 	setup_plus_xz_quadrant_patch_system(N_ghost_points, N_extend_points,
 					    N_overlap_points, delta_drho_dsigma,
 					    min_gfn_in, max_gfn_in,
@@ -127,6 +135,8 @@ case plus_xyz_octant_patch_system:
 	construct_patches(patch_system_info::plus_xyz_octant::patch_info_array,
 			  N_ghost_points, N_extend_points,
 			  delta_drho_dsigma);
+	setup_gridfn_storage(min_gfn, max_gfn,
+			     ghosted_min_gfn, ghosted_max_gfn);
 	setup_plus_xyz_octant_patch_system(N_ghost_points, N_extend_points,
 				  N_overlap_points, delta_drho_dsigma,
 				  min_gfn_in, max_gfn_in,
diff --git a/src/patch/patch_system.hh b/src/patch/patch_system.hh
index fa3c4d1..046e30f 100644
--- a/src/patch/patch_system.hh
+++ b/src/patch/patch_system.hh
@@ -144,8 +144,8 @@ public:
 		     enum patch_system_type type_in,
 		     int N_ghost_points, int N_extend_points,
 		     fp delta_drho_dsigma,
-		     int         min_gfn_in, int         max_gfn_in,
-		     int ghosted_min_gfn_in, int ghosted_max_gfn_in,
+		     int min_gfn, max_gfn,
+		     int ghosted_min_gfn, int ghosted_max_gfn);
 		     int interp_handle_in, int interp_par_table_handle_in);
 	~patch_system();
 
@@ -213,6 +213,10 @@ private:
 	// of all patches are fully setup
 	void assert_all_ghost_zones_fully_setup() const;
 
+	// setup all gridfns with contiguous-across-patches storage
+	void setup_gridfn_storage(int min_gfn, int max_gfn,
+				  int ghosted_min_gfn, int ghosted_max_gfn);
+
 private:
         // we forbid copying and passing by value
         // by declaring the copy constructor and assignment operator