complete interpatch-interpolation Jacobian code

add 1st draft of code in src/util/test_patch_system.cc to test this


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

1 files changed, 7 insertions, 7 deletions

diff --git a/src/patch/ghost_zone.cc b/src/patch/ghost_zone.cc
index 8516241..3197e2e 100644
--- a/src/patch/ghost_zone.cc
+++ b/src/patch/ghost_zone.cc
@@ -175,7 +175,7 @@ interpatch_ghost_zone::interpatch_ghost_zone(const patch_edge& my_edge_in,
           min_ipar_used_(NULL),	max_ipar_used_(NULL),
 	  other_par_(NULL),
 	  interp_result_buffer_(NULL),
-	  Jacobian_oipar_posn_(NULL), Jacobian_buffer_(NULL) // no comma
+	  Jacobian_y_ipar_posn_(NULL), Jacobian_buffer_(NULL) // no comma
 {
 //
 // verify that we have the expected relationships between
@@ -330,7 +330,7 @@ other_iperp_ = new jtutil::cpm_map<fp>(min_iperp(), max_iperp(),
 interpatch_ghost_zone::~interpatch_ghost_zone()
 {
 delete Jacobian_buffer_;
-delete Jacobian_oipar_posn_;
+delete Jacobian_y_ipar_posn_;
 delete interp_result_buffer_;
 delete other_par_;
 delete max_ipar_used_;
@@ -589,18 +589,18 @@ assert(other_patch_interp_ != NULL);
 
 other_patch_interp_->verify_Jacobian_sparsity_pattern_ok();
 
-other_patch_interp_->molecule_minmax_ipar_m(min_oipar_m_, max_oipar_m_);
+other_patch_interp_->molecule_minmax_ipar_m(min_y_ipar_m_, max_y_ipar_m_);
 
-if (Jacobian_oipar_posn_ == NULL)
-   then Jacobian_oipar_posn_ = new jtutil::array2d<CCTK_INT>
+if (Jacobian_y_ipar_posn_ == NULL)
+   then Jacobian_y_ipar_posn_ = new jtutil::array2d<CCTK_INT>
 				     (min_other_iperp_, max_other_iperp_,
 				      extreme_min_ipar_, extreme_max_ipar_);
-other_patch_interp_->molecule_posn(*Jacobian_oipar_posn_);
+other_patch_interp_->molecule_posn(*Jacobian_y_ipar_posn_);
 
 if (Jacobian_buffer_ == NULL)
    then Jacobian_buffer_ = new jtutil::array3d<fp>
 				 (min_other_iperp_, max_other_iperp_,
 				  extreme_min_ipar_, extreme_max_ipar_,
-				  min_oipar_m_, max_oipar_m_);
+				  min_y_ipar_m_, max_y_ipar_m_);
 other_patch_interp_->Jacobian(*Jacobian_buffer_);
 }