set gridfn min/max ipar in interpolator param table

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

1 files changed, 25 insertions, 4 deletions

diff --git a/src/patch/patch_interp.cc b/src/patch/patch_interp.cc
index 98f8284..be26ac6 100644
--- a/src/patch/patch_interp.cc
+++ b/src/patch/patch_interp.cc
@@ -90,11 +90,32 @@ status = Util_TableSetIntArray(interp_par_table_handle_,
 if (status < 0)
    then error_exit(ERROR_EXIT,
 "***** patch_interp::patch_interp():\n"
-"        can't set stride in interpolator parmameter table (handle=%d)!\n"
-"        error status=%d\n"
+"        can't set gridfn stride in interpolator parmameter table!\n"
+"        table handle=%d   error status=%d\n"
 ,
-		   interp_par_table_handle_,
-		   status);					/*NOTREACHED*/
+		   interp_par_table_handle_, status);		/*NOTREACHED*/
+
+// set up patch interpolation region in interpolator parameter table
+status = Util_TableSetIntArray(interp_par_table_handle_,
+			       N_dims, &min_ipar_,
+			       "input_array_min_subscripts");
+if (status < 0)
+   then error_exit(ERROR_EXIT,
+"***** patch_interp::patch_interp():\n"
+"        can't set min gridfn subscript in interpolator parmameter table!\n"
+"        table handle=%d   error status=%d\n"
+,
+		   interp_par_table_handle_, status);		/*NOTREACHED*/
+status = Util_TableSetIntArray(interp_par_table_handle_,
+			       N_dims, &max_ipar_,
+			       "input_array_max_subscripts");
+if (status < 0)
+   then error_exit(ERROR_EXIT,
+"***** patch_interp::patch_interp():\n"
+"        can't set max gridfn subscript in interpolator parmameter table!\n"
+"        table handle=%d   error status=%d\n"
+,
+		   interp_par_table_handle_, status);		/*NOTREACHED*/
 
 // set up gridfn type codes
 	for (int gfn = my_patch().ghosted_min_gfn() ;