* finish Jacobian support in src/util

* test driver for this in src/util/test_patch_system.cc
  --> IT WORKS!!!!!
* also fix comments on #include prerequisites to say ../jtutil/ not jt/


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

1 files changed, 4 insertions, 59 deletions

diff --git a/src/patch/patch.hh b/src/patch/patch.hh
index 6f812e5..b58882d 100644
--- a/src/patch/patch.hh
+++ b/src/patch/patch.hh
@@ -3,7 +3,6 @@
 //
 // ***** how patch boundaries are handled *****
 // patch - abstract base class to describe a coordinate/grid patch
-// patch_info - POD struct with summary info for patch construction
 //
 // z_patch - derived class for a +/- z patch
 // x_patch - derived class for a +/- x patch
@@ -15,10 +14,11 @@
 //	<stdio.h>
 //	<assert.h>
 //	<math.h>
-//	"jt/util.hh"
-//	"jt/array.hh"
-//	"jt/linear_map.hh"
+//	"stdc.h"
 //	"config.hh"
+//	"../jtutil/util.hh"
+//	"../jtutil/array.hh"
+//	"../jtutil/linear_map.hh"
 //	"coords.hh"
 //	"grid.hh"
 //	"fd_grid.hh"
@@ -425,61 +425,6 @@ private:
 //*****************************************************************************
 
 //
-// This (POD, and hence static-initializable) struct gives a minimal
-// set of information which varies from one patch to another.
-//
-// The member functions allow a patch to be constructed from this struct
-// together with the additional information provided by their arguments.
-// This doesn't allow the most general possible type of patch (eg it
-// constrains all ghost zones to have the same width), but does cover
-// the common casees.
-//
-// Arguments for member functions:
-// N_ghost_points = Width in grid points of all ghost zones.
-// N_extend_points = Number of grid points to extend each patch past
-//		     "just touching" so as to overlap neighboring patches.
-//		     Thus patches overlap by
-//			N_overlap_points = 2*N_extend_points + 1
-//		     grid points.  For example, with N_extend_points == 2,
-//		     here are the grid points of two neighboring patches:
-//			x   x   x   x   x   X   X
-//                                      |
-//			        O   O   o   o   o   o   o
-//		     Here | marks the "just touching" boundary,
-//		     x and o the grid points before this extension,
-//		     and X and O the extra grid points added by this
-//		     extension.
-// delta_drho_dsigma = Grid spacing (both rho and sigma) in degrees.
-//
-struct	patch_info
-	{
-	const char* name;
-	bool is_plus;
-	char ctype;
-	fp min_drho, max_drho;
-	fp min_dsigma, max_dsigma;
-
-	// compute and return reference to  struct grid_arrays::grid_array_pars
-	// ... 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_extend_points,
-			  fp delta_drho_dsigma)
-		const;
-
-	// compute and return reference to  struct grid::grid_pars
-	// ... result refers to internal static buffer;
-	//     the usual caveats about lifetimes/overwriting apply
-	const grid::grid_pars& grid_pars(int N_extend_points,
-					 fp delta_drho_dsigma)
-		const;
-	};
-
-//*****************************************************************************
-//*****************************************************************************
-//*****************************************************************************
-
-//
 // This class describes a +/- z patch.  It doesn't define any new
 // functions not already present in  class patch ; it "just" defines
 // non-virtual versions of all the pure virtual functions defined there.