1 files changed, 337 insertions, 0 deletions
diff --git a/Carpet/Carpet/src/helpers.cc b/Carpet/Carpet/src/helpers.cc
new file mode 100644
index 000000000..70e892b4d
--- /dev/null
+++ b/Carpet/Carpet/src/helpers.cc
@@ -0,0 +1,337 @@
+#include <assert.h>
+#include <stdarg.h>
+#include <stdio.h>
+#include <stdlib.h>
+#include <string.h>
+
+#include <mpi.h>
+
+#include "cctk.h"
+#include "cctk_FortranString.h"
+#include "cctk_Parameters.h"
+
+#include "defs.hh"
+#include "dist.hh"
+#include "ggf.hh"
+
+#include "carpet.hh"
+
+extern "C" {
+  static const char* rcsid = "$Header: /home/eschnett/C/carpet/Carpet/Carpet/Carpet/src/helpers.cc,v 1.49 2004/07/08 12:43:34 tradke Exp $";
+  CCTK_FILEVERSION(Carpet_Carpet_helpers_cc);
+}
+
+
+
+namespace Carpet {
+
+  using namespace std;
+
+
+
+  // Get Carpet's GH extension
+  CarpetGH const * GetCarpetGH (const cGH * const cgh)
+  {
+    assert (cgh);
+    return &carpetGH;
+  }
+
+
+
+  // Enable or disable prolongating
+  CCTK_INT CarpetEnableProlongating (const CCTK_INT flag)
+  {
+    assert (flag==0 || flag==1);
+    do_prolongate = flag;
+    if (do_prolongate) {
+      Checkpoint ("Prolongating enabled");
+    } else {
+      Checkpoint ("Prolongating disabled");
+    }
+    return 0;
+  }
+
+
+
+  // Communication
+
+  int Barrier (const cGH* cgh)
+  {
+    void *dummy = &dummy;
+    dummy = &cgh;
+
+    MPI_Barrier (dist::comm);
+    return 0;
+  }
+
+
+
+  int Exit (cGH* cgh, int retval)
+  {
+    CCTK_Barrier (cgh);
+    dist::finalize();
+    exit (retval);
+    return -999;
+  }
+
+  int Abort (cGH* cgh, int retval)
+  {
+    void *dummy = &dummy;
+    dummy = &cgh;
+
+    MPI_Comm comm = dist::comm;
+    if (comm == MPI_COMM_NULL)
+    {
+      comm = MPI_COMM_WORLD;
+    }
+    MPI_Abort (comm, retval);
+    abort ();
+    return -999;
+  }
+
+
+
+  int MyProc (const cGH* cgh)
+  {
+    // if there is no cgh yet, assume nothing has been initialised
+    // yet, and don't use dist::comm
+    int rank;
+    MPI_Comm_rank (cgh ? dist::comm : MPI_COMM_WORLD, &rank);
+    return rank;
+  }
+
+  int nProcs (const cGH* cgh)
+  {
+    // if there is no cgh yet, assume nothing has been initialised
+    // yet, and don't use dist::comm
+    int size;
+    MPI_Comm_size (cgh ? dist::comm : MPI_COMM_WORLD, &size);
+    return size;
+  }
+
+
+
+  MPI_Comm CarpetMPIComm ()
+  {
+    return dist::comm;
+  }
+
+
+
+  // Datatypes
+
+  MPI_Datatype CarpetMPIDatatype (const int vartype)
+  {
+    switch (vartype) {
+#define TYPECASE(N,T)				\
+    case N: {					\
+      T dummy;					\
+      return dist::datatype(dummy);		\
+    }
+#include "typecase"
+#undef TYPECASE
+    default:
+      CCTK_VWarn (0, __LINE__, __FILE__, CCTK_THORNSTRING,
+		  "Carpet does not support the variable type %d.", vartype);
+    }
+    // notreached
+    return MPI_CHAR;
+  }
+
+  MPI_Datatype CarpetSimpleMPIDatatype (const int vartype)
+  {
+    switch (vartype) {
+#ifdef CARPET_COMPLEX
+    case CCTK_VARIABLE_COMPLEX:
+      return CarpetMPIDatatype (CCTK_VARIABLE_REAL);
+#endif
+#ifdef CARPET_COMPLEX8
+#  ifdef CCTK_REAL4
+    case CCTK_VARIABLE_COMPLEX8:
+      return CarpetMPIDatatype (CCTK_VARIABLE_REAL4);
+#  endif
+#endif
+#ifdef CARPET_COMPLEX16
+#  ifdef CCTK_REAL8
+    case CCTK_VARIABLE_COMPLEX16:
+      return CarpetMPIDatatype (CCTK_VARIABLE_REAL8);
+#  endif
+#endif
+#ifdef CARPET_COMPLEX32
+#  ifdef CCTK_REAL16
+    case CCTK_VARIABLE_COMPLEX32:
+      return CarpetMPIDatatype (CCTK_VARIABLE_REAL16);
+#  endif
+#endif
+    default:
+      return CarpetMPIDatatype (vartype);
+    }
+    // notreached
+    return MPI_CHAR;
+  }
+
+  int CarpetSimpleMPIDatatypeLength (const int vartype)
+  {
+    switch (vartype) {
+#ifdef CARPET_COMPLEX
+    case CCTK_VARIABLE_COMPLEX:
+#endif
+#ifdef CARPET_COMPLEX8
+#  ifdef CCTK_REAL4
+    case CCTK_VARIABLE_COMPLEX8:
+#  endif
+#endif
+#ifdef CARPET_COMPLEX16
+#  ifdef CCTK_REAL8
+    case CCTK_VARIABLE_COMPLEX16:
+#  endif
+#endif
+#ifdef CARPET_COMPLEX32
+#  ifdef CCTK_REAL16
+    case CCTK_VARIABLE_COMPLEX32:
+#  endif
+#endif
+      return 2;
+    default:
+      return 1;
+    }
+    // notreached
+    return 0;
+  }
+
+
+
+  // Timelevels
+
+  int mintl (const checktimes where, const int num_tl)
+  {
+    assert (num_tl>0);
+    switch (where) {
+    case currenttime:
+      return 0;
+    case currenttimebutnotifonly:
+      // don't include current time if there is only one time level
+      return num_tl>1 ? 0 : 1;
+    case previoustime:
+      return 1;
+    case allbutlasttime:
+      // do include current time if there is only one time level
+      return 0;
+    case allbutcurrenttime:
+      return 1;
+    case alltimes:
+      return 0;
+    default:
+      assert (0);
+    }
+    return -999;
+  }
+
+  int maxtl (const checktimes where, const int num_tl)
+  {
+    assert (num_tl>0);
+    switch (where) {
+    case currenttime:
+      return 0;
+    case currenttimebutnotifonly:
+      return 0;
+    case previoustime:
+      return num_tl>1 ? 1 : 0;
+    case allbutlasttime:
+      return num_tl-2;
+    case allbutcurrenttime:
+      return num_tl-1;
+    case alltimes:
+      return num_tl-1;
+    default:
+      assert (0);
+    }
+    return -999;
+  }
+
+
+
+  // Diagnostic output
+
+  static void prepend_id (char* const msg, size_t const msglen)
+  {
+    if (mglevel!=-1) {
+      snprintf (msg+strlen(msg), msglen-strlen(msg), "[%d]", mglevel);
+      if (reflevel!=-1) {
+        snprintf (msg+strlen(msg), msglen-strlen(msg), "[%d]", reflevel);
+        if (map!=-1) {
+          snprintf (msg+strlen(msg), msglen-strlen(msg), "[%d]", map);
+          if (component!=-1) {
+            snprintf (msg+strlen(msg), msglen-strlen(msg), "[%d]", component);
+          }
+        }
+      }
+      snprintf (msg+strlen(msg), msglen-strlen(msg), " ");
+    }
+  }
+
+  void Output (const char* fmt, ...)
+  {
+    DECLARE_CCTK_PARAMETERS;
+    va_list args;
+    char msg[1000];
+    snprintf (msg, sizeof msg, "%s", "");
+    prepend_id (msg + strlen(msg), sizeof msg - strlen(msg));
+    va_start (args, fmt);
+    vsnprintf (msg + strlen(msg), sizeof msg - strlen(msg), fmt, args);
+    va_end (args);
+    CCTK_INFO (msg);
+    if (barriers) {
+      MPI_Barrier (dist::comm);
+    }
+  }
+
+  void Waypoint (const char* fmt, ...)
+  {
+    DECLARE_CCTK_PARAMETERS;
+    if (verbose || veryverbose) {
+      va_list args;
+      char msg[1000];
+      snprintf (msg, sizeof msg, "%s", "");
+      prepend_id (msg + strlen(msg), sizeof msg - strlen(msg));
+      va_start (args, fmt);
+      vsnprintf (msg + strlen(msg), sizeof msg - strlen(msg), fmt, args);
+      va_end (args);
+      CCTK_INFO (msg);
+    }
+    if (barriers) {
+      MPI_Barrier (dist::comm);
+    }
+  }
+
+  void Checkpoint (const char* fmt, ...)
+  {
+    DECLARE_CCTK_PARAMETERS;
+    if (veryverbose) {
+      va_list args;
+      char msg[1000];
+      snprintf (msg, sizeof msg, "%s", "");
+      prepend_id (msg + strlen(msg), sizeof msg - strlen(msg));
+      va_start (args, fmt);
+      vsnprintf (msg + strlen(msg), sizeof msg - strlen(msg), fmt, args);
+      va_end (args);
+      CCTK_INFO (msg);
+    }
+    if (barriers) {
+      MPI_Barrier (dist::comm);
+    }
+  }
+
+
+
+  void UnsupportedVarType (const int vindex)
+  {
+    assert (vindex>=0 && vindex<CCTK_NumVars());
+    CCTK_VWarn
+      (0, __LINE__, __FILE__, CCTK_THORNSTRING,
+       "Carpet does not support the type of the variable \"%s\".\n"
+       "Either enable support for this type, "
+       "or change the type of this variable.", CCTK_FullName(vindex));
+  }
+
+} // namespace Carpet