Regrid all grid functions at once. This is necessary to treate vector

Regrid all grid functions at once.  This is necessary to treate vector
grid functions correctly, where otherwise storage would be released
too early.  This requires the regridding interface between the dh and
ggf classes to change.  Another advantage is that this should now be
faster on multiple processors.

darcs-hash:20040419193833-07bb3-f7474744feaa57f8a1621aeca34de11b65eaeed8.gz

1 files changed, 11 insertions, 2 deletions

diff --git a/Carpet/CarpetLib/src/ggf.hh b/Carpet/CarpetLib/src/ggf.hh
index 645d047e2..3b7035bf1 100644
--- a/Carpet/CarpetLib/src/ggf.hh
+++ b/Carpet/CarpetLib/src/ggf.hh
@@ -1,4 +1,4 @@
-// $Header: /home/eschnett/C/carpet/Carpet/Carpet/CarpetLib/src/ggf.hh,v 1.22 2004/04/19 07:56:35 schnetter Exp $
+// $Header: /home/eschnett/C/carpet/Carpet/Carpet/CarpetLib/src/ggf.hh,v 1.23 2004/04/19 21:38:33 schnetter Exp $
 
 #ifndef GGF_HH
 #define GGF_HH
@@ -68,6 +68,9 @@ protected:
   int vectorindex;              // index of *this
   ggf* vectorleader;            // first vector element
   
+private:
+  fdata oldstorage;
+  
 public:
 
   // Constructors
@@ -87,7 +90,13 @@ public:
 
 
   // Modifiers
-  void recompose ();
+  // void recompose ();
+  void recompose_crop ();
+  void recompose_allocate (int rl);
+  void recompose_fill (comm_state<D>& state, int rl);
+  void recompose_free (int rl);
+  void recompose_bnd_prolongate (comm_state<D>& state, int rl);
+  void recompose_sync (comm_state<D>& state, int rl);
 
   // Cycle the time levels by rotating the data sets
   void cycle (int rl, int c, int ml);