49 assert(Ninstance == 1);
58 fprintf(stderr,
"Error at Communicator::setup(): inappropriate grid_size.\n");
68 const std::vector<int>& lattice_size,
69 std::vector<int>& grid_size,
72 int Ndim = lattice_size.size();
74 if (grid_size.size() != Ndim) {
75 grid_size.resize(Ndim, 1);
79 for (
int i = 0; i < Ndim; ++i) {
80 if (grid_size[i] != 1) {
86 printf(
"Error at Communicator: %s: unexpected grid_size.\n", __func__);
123 int Communicator::self_global()
156 for (
int i = 0; i < Ndim; ++i) {
168 for (
int i = 0; i < Ndim; ++i) {
191 int Communicator::sync_global()
243 memcpy(recv_buf, send_buf, size);
250 memcpy(recv_buf, send_buf,
sizeof(dcomplex) * count);
257 memcpy(recv_buf, send_buf,
sizeof(
double) * count);
264 memcpy(recv_buf, send_buf,
sizeof(
float) * count);
271 memcpy(recv_buf, send_buf,
sizeof(
int) * count);
279 memcpy(recv_buf, send_buf, size);
284 int Communicator::send_1to1(
int count, dcomplex *recv_buf, dcomplex *send_buf,
int send_to,
int recv_from,
int tag)
286 memcpy(recv_buf, send_buf,
sizeof(dcomplex) * count);
291 int Communicator::send_1to1(
int count,
double *recv_buf,
double *send_buf,
int send_to,
int recv_from,
int tag)
293 memcpy(recv_buf, send_buf,
sizeof(
double) * count);
300 memcpy(recv_buf, send_buf,
sizeof(
float) * count);
307 memcpy(recv_buf, send_buf,
sizeof(
int) * count);
315 memcpy(recv_buf, send_buf,
sizeof(dcomplex) * count);
322 memcpy(recv_buf, send_buf,
sizeof(
double) * count);
329 memcpy(recv_buf, send_buf,
sizeof(
float) * count);
336 memcpy(recv_buf, send_buf,
sizeof(
int) * count);
395 printf(
"Communicator Single\n");
406 if (gettimeofday(&now, (
struct timezone *)0) != 0) {
410 double sec = (double)now.tv_sec + ((
double)now.tv_usec) * 1.0e-6;
static int npe(const int dir)
logical grid extent
static int reduce_max(int count, double *recv_buf, double *send_buf, int pattern=0)
find a global maximum of an array of double over the communicator. pattern specifies the dimensions t...
static void abort()
terminate communicator
static int self()
rank within small world.
static bool is_primary_master()
check if the present node is primary in global communicator.
static int init(int *pargc, char ***pargv)
initialize communicator
static int grid_coord(int *grid_coord, const int rank)
find grid coordinate from rank number.
static int reduce_min(int count, double *recv_buf, double *send_buf, int pattern=0)
find a global minimum of an array of double over the communicator. pattern specifies the dimensions t...
static int broadcast(size_t size, void *data, int sender)
static int ipe(const int dir)
logical coordinate of current proc.
static int reduce_sum(int count, dcomplex *recv_buf, dcomplex *send_buf, int pattern=0)
make a global sum of an array of dcomplex over the communicator. pattern specifies the dimensions to ...
static int finalize()
finalize communicator
static int exchange(size_t size, void *recv_buf, void *send_buf, int idir, int ipm, int tag)
static int grid_dims(int *grid_dims)
find grid dimensions.
static int send_1to1(int count, dcomplex *recv_buf, dcomplex *send_buf, int p_to, int p_from, int tag)
send array of dcomplex from rank p_from to rank p_to. communication distinguished by tag...
static int grid_rank(int *rank, const int *grid_coord)
find rank number from grid coordinate.
static int setup(int ninstance=1)
setup communicator
static int sync_usleep()
synchronize within small world. (slow but no busy wait)
static int send_1to1(size_t size, void *recv_buf, void *send_buf, int send_to, int recv_from, int tag)
static int size()
size of small world.
static int broadcast(int count, dcomplex *data, int sender)
broadcast array of dcomplex from sender.
static int exchange(int count, dcomplex *recv_buf, dcomplex *send_buf, int idir, int ipm, int tag)
receive array of dcomplex from upstream specified by idir and ipm, and send array to downstream...
static double get_time()
obtain a wall-clock time.
static int sync()
synchronize within small world.
static bool is_primary()
check if the present node is primary in small communicator.