diff options
Diffstat (limited to 'src/patch/patch_system.hh')
-rw-r--r-- | src/patch/patch_system.hh | 27 |
1 files changed, 22 insertions, 5 deletions
diff --git a/src/patch/patch_system.hh b/src/patch/patch_system.hh index a4777f0..a1b4227 100644 --- a/src/patch/patch_system.hh +++ b/src/patch/patch_system.hh @@ -27,7 +27,8 @@ //****************************************************************************** // -// a patch_system object describes a system of interlinked patches. +// A patch_system object describes a system of interlinked patches. +// Its const qualifiers refer to the gridfn data. // class patch_system @@ -245,13 +246,26 @@ public: void synchronize(int ghosted_min_gfn_to_sync, int ghosted_max_gfn_to_sync); + // ... do this for all ghosted gridfns + void synchronize() + { synchronize(ghosted_min_gfn(), ghosted_max_gfn()); } + // // compute Jacobian of synchronize() into internal buffers, // taking into account synchronize()'s full 3-phase algorithm // void compute_synchronize_Jacobian(int ghosted_min_gfn_to_sync, - int ghosted_max_gfn_to_sync); + int ghosted_max_gfn_to_sync) + const; + + // ... do this for all ghosted gridfns + void compute_synchronize_Jacobian() + const + { + compute_synchronize_Jacobian(ghosted_min_gfn(), + ghosted_max_gfn()); + } // // The following functions access the Jacobian computed by @@ -266,7 +280,8 @@ public: const { return xgz.Jacobian_y_patch(); } const patch_edge& synchronize_Jacobian_y_edge (const ghost_zone& xgz) - const { return xgz.Jacobian_y_edge(); } + const + { return xgz.Jacobian_y_edge(); } // what is the [min,max] range of m for a Jacobian row? int synchronize_Jacobian_min_y_ipar_m(const ghost_zone& xgz) @@ -283,7 +298,8 @@ public: // what is the posn value of the y points in this Jacobian row? int synchronize_Jacobian_y_ipar_posn(const ghost_zone& xgz, - int x_iperp, int x_ipar) const; + int x_iperp, int x_ipar) + const; // what is the Jacobian // partial synchronize() gridfn(ghosted_gfn, px, x_iperp, x_ipar) @@ -294,7 +310,8 @@ public: // y_posn = Jacobian_y_ipar_posn(x_iperp, x_ipar) // taking into account synchronize()'s full 3-phase algorithm fp synchronize_Jacobian(const ghost_zone& xgz, int x_iperp, int x_ipar, - int y_ipar_m) const; + int y_ipar_m) + const; // |