diff options
author | jthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5> | 2002-07-08 13:52:34 +0000 |
---|---|---|
committer | jthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5> | 2002-07-08 13:52:34 +0000 |
commit | 75a52836516dc22f60133fa132fe5f63b07a7309 (patch) | |
tree | 8ccb2d2c6bd4d54c262a7bb138babca0083aa9c1 /src/patch/ghost_zone.cc | |
parent | d7f9c48e97fef4ee571f0e5ae5b493cce36bc981 (diff) |
* make want_*_corner arguments optional (default = true) for interpatch ghost
zone synchronize and Jacobian fns
* rename Jacobian_y_min_ipar_m() --> Jacobian_min_y_ipar_m()
etc etc
* change interpatch_ghost_zone::Jacobian()
to only accept in-range y_ipar_m coords (assert failure if out of range)
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@612 f88db872-0e4f-0410-b76b-b9085cfa78c5
Diffstat (limited to 'src/patch/ghost_zone.cc')
-rw-r--r-- | src/patch/ghost_zone.cc | 30 |
1 files changed, 16 insertions, 14 deletions
diff --git a/src/patch/ghost_zone.cc b/src/patch/ghost_zone.cc index 0b7825e..72d9738 100644 --- a/src/patch/ghost_zone.cc +++ b/src/patch/ghost_zone.cc @@ -122,9 +122,9 @@ delete iperp_map_; // we want. // void symmetry_ghost_zone::synchronize(int ghosted_min_gfn, int ghosted_max_gfn, - bool want_min_par_corner, - bool want_non_corner, - bool want_max_par_corner) + bool want_min_par_corner = true, + bool want_non_corner = true, + bool want_max_par_corner = true) { for (int gfn = ghosted_min_gfn ; gfn <= ghosted_max_gfn ; ++gfn) { @@ -537,9 +537,9 @@ if (! (want_min_par_corner && want_non_corner && want_max_par_corner)) // void interpatch_ghost_zone::synchronize (int ghosted_min_gfn, int ghosted_max_gfn, - bool want_min_par_corner, - bool want_non_corner, - bool want_max_par_corner) + bool want_min_par_corner = true, + bool want_non_corner = true, + bool want_max_par_corner = true) { verify_caller_wants_all_of_ghost_zone(want_min_par_corner, want_non_corner, @@ -582,9 +582,9 @@ other_patch_interp_->interpolate(ghosted_min_gfn, ghosted_max_gfn, // void interpatch_ghost_zone::compute_Jacobian (int ghosted_min_gfn, int ghosted_max_gfn, - bool want_min_par_corner, - bool want_non_corner, - bool want_max_par_corner) + bool want_min_par_corner = true, + bool want_non_corner = true, + bool want_max_par_corner = true) { verify_caller_wants_all_of_ghost_zone(want_min_par_corner, want_non_corner, @@ -593,7 +593,8 @@ assert(other_patch_interp_ != NULL); other_patch_interp_->verify_Jacobian_sparsity_pattern_ok(); -other_patch_interp_->molecule_minmax_ipar_m(min_y_ipar_m_, max_y_ipar_m_); +other_patch_interp_->molecule_minmax_ipar_m(Jacobian_min_y_ipar_m_, + Jacobian_max_y_ipar_m_); if (Jacobian_y_ipar_posn_ == NULL) then Jacobian_y_ipar_posn_ = new jtutil::array2d<CCTK_INT> @@ -602,9 +603,10 @@ if (Jacobian_y_ipar_posn_ == NULL) 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_y_ipar_m_, max_y_ipar_m_); + then Jacobian_buffer_ + = new jtutil::array3d<fp> + (min_other_iperp_, max_other_iperp_, + extreme_min_ipar(), extreme_max_ipar(), + Jacobian_min_y_ipar_m_, Jacobian_max_y_ipar_m_); other_patch_interp_->Jacobian(*Jacobian_buffer_); } |