blob: dac83076fce047439673c36b0eeec5596daeb494 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
|
/***************************************************************************
dist.hh - Helpers for distributed computing
-------------------
begin : Sun Jun 11 2000
copyright : (C) 2000 by Erik Schnetter
email : schnetter@astro.psu.edu
$Header: /home/eschnett/C/carpet/Carpet/Carpet/CarpetLib/src/dist.hh,v 1.3 2001/03/07 13:00:57 eschnett Exp $
***************************************************************************/
/***************************************************************************
* *
* This program is free software; you can redistribute it and/or modify *
* it under the terms of the GNU General Public License as published by *
* the Free Software Foundation; either version 2 of the License, or *
* (at your option) any later version. *
* *
***************************************************************************/
#ifndef DIST_HH
#define DIST_HH
#include <cassert>
#include <cstdio>
#include <cstdlib>
#include <complex>
#include <mpi.h>
#include "defs.hh"
namespace dist {
const int tag = 1;
extern MPI_Comm comm;
extern MPI_Datatype mpi_complex_float;
extern MPI_Datatype mpi_complex_double;
extern MPI_Datatype mpi_complex_long_double;
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
// This generic routine is only declared and not defined. Only
// specialised versions have definitions. Type errors will be
// caught by the compiler.
template<class T>
MPI_Datatype datatype (const T& dummy);
template<class T>
inline MPI_Datatype datatype (const T& dummy)
{ abort(); return -1; }
template<>
inline MPI_Datatype datatype (const char& dummy)
{ return MPI_CHAR; }
template<>
inline MPI_Datatype datatype (const signed char& dummy)
{ return MPI_UNSIGNED_CHAR; }
template<>
inline MPI_Datatype datatype (const unsigned char& dummy)
{ return MPI_BYTE; }
template<>
inline MPI_Datatype datatype (const short& dummy)
{ return MPI_SHORT; }
template<>
inline MPI_Datatype datatype (const unsigned short& dummy)
{ return MPI_UNSIGNED_SHORT; }
template<>
inline MPI_Datatype datatype (const int& dummy)
{ return MPI_INT; }
template<>
inline MPI_Datatype datatype (const unsigned int& dummy)
{ return MPI_UNSIGNED; }
template<>
inline MPI_Datatype datatype (const long& dummy)
{ return MPI_LONG; }
template<>
inline MPI_Datatype datatype (const unsigned long& dummy)
{ return MPI_UNSIGNED_LONG; }
template<>
inline MPI_Datatype datatype (const long long& dummy)
{ return MPI_LONG_LONG_INT; }
template<>
inline MPI_Datatype datatype (const float& dummy)
{ return MPI_FLOAT; }
template<>
inline MPI_Datatype datatype (const double& dummy)
{ return MPI_DOUBLE; }
template<>
inline MPI_Datatype datatype (const long double& dummy)
{ return MPI_LONG_DOUBLE; }
template<>
inline MPI_Datatype datatype (const complex<float>& dummy)
{ return mpi_complex_float; }
template<>
inline MPI_Datatype datatype (const complex<double>& dummy)
{ return mpi_complex_double; }
template<>
inline MPI_Datatype datatype (const complex<long double>& dummy)
{ return mpi_complex_long_double; }
} // namespace dist
#if defined(TMPL_IMPLICIT)
# include "dist.cc"
#endif
#endif // DIST_HH
|
