Many changes that accumulated while Cactus and Carpet diverged.

Many changes that accumulated while Cactus and Carpet diverged.
Add processor splitting mechanism "along-dir" that splits along a
specified direction.

Rename group PostRestrict to bin POSTRESTRICT.

Prolongate initial data only when desired.  This saves much time.
Sorry, Ian.

Fix bug in time level cycling of grid arrays.  (Note: grid arrays
should not have time levels.)

Fix time_t bug on IRIX.

Make sure that there is no integer overflow when there are many
refinement levels.  Always put parentheses around
(maxreflevelfact/reflevelfact).

Fix typo in Carpet verbose output.

Add debug output in processor splitting.

Communicate in three stages: Irecv, (work), Isend, Wait.  This might
be more efficient.  Much more, potentially.

Fix bug in processor layout of grid arrays.  Sorry, Ian.

Make the interpolator interpolate between time levels.  Untested.

Fix bug in processor communication in interpolator.  Sorry, Ian.

Rewrite prolongation operators to make them twice as fast.  There you
are, Ian.

Move prolongation operator kind handling from data to gdata.

Add official hyperslabbing interfaces to CarpetSlab.

Adapt to new cGH * handling.

darcs-hash:20031105151837-07bb3-758a87ff0355dba053269df4b7d7d79bea018669.gz

1 files changed, 6 insertions, 2 deletions

diff --git a/CarpetAttic/CarpetIOFlexIO/src/ioflexio.cc b/CarpetAttic/CarpetIOFlexIO/src/ioflexio.cc
index 7456091de..2a00ebe5c 100644
--- a/CarpetAttic/CarpetIOFlexIO/src/ioflexio.cc
+++ b/CarpetAttic/CarpetIOFlexIO/src/ioflexio.cc
@@ -305,7 +305,9 @@ namespace CarpetIOFlexIO {
       
       gdata<dim>* const tmp = data->make_typed (n);
       tmp->allocate (ext, 0);
-      tmp->copy_from (data, ext);
+      for (comm_state<dim> state; !state.done(); state.step()) {
+        tmp->copy_from (state, data, ext);
+      }
       
       // Write data
       if (CCTK_MyProc(cgh)==0) {
@@ -600,7 +602,9 @@ namespace CarpetIOFlexIO {
         }
         
         // Copy into grid function
-        data->copy_from (tmp, ext);
+        for (comm_state<dim> state; !state.done(); state.step()) {
+          data->copy_from (state, tmp, ext);
+        }
         
         // Delete temporary copy
         delete tmp;