1 files changed, 120 insertions, 0 deletions
diff --git a/Carpet/CarpetLib/src/dist.hh b/Carpet/CarpetLib/src/dist.hh
new file mode 100644
index 000000000..ab89ae580
--- /dev/null
+++ b/Carpet/CarpetLib/src/dist.hh
@@ -0,0 +1,120 @@
+// $Header: /home/eschnett/C/carpet/Carpet/Carpet/CarpetLib/src/dist.hh,v 1.10 2004/03/01 19:43:39 schnetter Exp $
+
+#ifndef DIST_HH
+#define DIST_HH
+
+#include <assert.h>
+#include <stdio.h>
+#include <stdlib.h>
+
+#include <complex>
+
+#include <mpi.h>
+
+#include "cctk.h"
+
+#include "defs.hh"
+
+using namespace std;
+
+
+
+namespace dist {
+  
+  extern MPI_Comm comm;
+  
+#if 0
+  extern MPI_Datatype mpi_complex_float;
+  extern MPI_Datatype mpi_complex_double;
+  extern MPI_Datatype mpi_complex_long_double;
+#else
+  extern MPI_Datatype mpi_complex8;
+  extern MPI_Datatype mpi_complex16;
+  extern MPI_Datatype mpi_complex32;
+#endif
+  
+  void init (int& argc, char**& argv);
+  void pseudoinit ();
+  void finalize ();
+  
+  // Debugging output
+#define CHECKPOINT dist::checkpoint(__FILE__, __LINE__)
+  void checkpoint (const char* file, int line);
+  
+  
+  
+  // Datatype helpers
+  inline MPI_Datatype datatype (const char& dummy)
+  { return MPI_CHAR; }
+  
+  inline MPI_Datatype datatype (const signed char& dummy)
+  { return MPI_UNSIGNED_CHAR; }
+  
+  inline MPI_Datatype datatype (const unsigned char& dummy)
+  { return MPI_BYTE; }
+  
+  inline MPI_Datatype datatype (const short& dummy)
+  { return MPI_SHORT; }
+  
+  inline MPI_Datatype datatype (const unsigned short& dummy)
+  { return MPI_UNSIGNED_SHORT; }
+  
+  inline MPI_Datatype datatype (const int& dummy)
+  { return MPI_INT; }
+  
+  inline MPI_Datatype datatype (const unsigned int& dummy)
+  { return MPI_UNSIGNED; }
+  
+  inline MPI_Datatype datatype (const long& dummy)
+  { return MPI_LONG; }
+  
+  inline MPI_Datatype datatype (const unsigned long& dummy)
+  { return MPI_UNSIGNED_LONG; }
+  
+  inline MPI_Datatype datatype (const long long& dummy)
+  { return MPI_LONG_LONG_INT; }
+  
+  inline MPI_Datatype datatype (const float& dummy)
+  { return MPI_FLOAT; }
+  
+  inline MPI_Datatype datatype (const double& dummy)
+  { return MPI_DOUBLE; }
+  
+  inline MPI_Datatype datatype (const long double& dummy)
+  { return MPI_LONG_DOUBLE; }
+  
+#if 0
+  
+  inline MPI_Datatype datatype (const complex<float>& dummy)
+  { return mpi_complex_float; }
+  
+  inline MPI_Datatype datatype (const complex<double>& dummy)
+  { return mpi_complex_double; }
+  
+  inline MPI_Datatype datatype (const complex<long double>& dummy)
+  { return mpi_complex_long_double; }
+  
+#else
+  
+#  ifdef CCTK_REAL4
+  inline MPI_Datatype datatype (const CCTK_COMPLEX8& dummy)
+  { return mpi_complex8; }
+#  endif
+  
+#  ifdef CCTK_REAL8
+  inline MPI_Datatype datatype (const CCTK_COMPLEX16& dummy)
+  { return mpi_complex16; }
+#  endif
+  
+#  ifdef CCTK_REAL16
+  inline MPI_Datatype datatype (const CCTK_COMPLEX32& dummy)
+  { return mpi_complex32; }
+#  endif
+  
+#endif
+  
+} // namespace dist
+
+
+
+#endif // DIST_HH