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
|
/* Copyright (c) 2009-2018. The SimGrid Team.
* All rights reserved. */
/* This program is free software; you can redistribute it and/or modify it
* under the terms of the license (GNU LGPL) which comes with this package. */
#include <stdlib.h>
#include <stdio.h>
#include <string.h>
#include <errno.h>
#include "mpi.h"
int main(int argc, char *argv[])
{
int i;
int rank;
int size;
int status;
MPI_Init(&argc, &argv);
MPI_Comm_rank(MPI_COMM_WORLD, &rank);
MPI_Comm_size(MPI_COMM_WORLD, &size);
int* recv_counts = (int *) xbt_malloc(size * sizeof(int));
int* recv_disps = (int *) xbt_malloc(size * sizeof(int));
int recv_sb_size = 0;
for (i = 0; i < size; i++) {
recv_counts[i] = i + 1;
recv_disps[i] = recv_sb_size;
recv_sb_size += i + 1;
}
int* sb = (int *) xbt_malloc(recv_counts[rank] * sizeof(int));
int* rb = (int *) xbt_malloc(recv_sb_size * sizeof(int));
printf("[%d] sndbuf=[", rank);
for (i = 0; i < recv_counts[rank]; i++){
sb[i] = recv_disps[rank] + i;
printf("%d ", sb[i]);
}
printf("]\n");
for (i = 0; i < recv_sb_size; i++)
rb[i] = -1;
status = MPI_Allgatherv(sb, recv_counts[rank], MPI_INT, rb, recv_counts, recv_disps, MPI_INT, MPI_COMM_WORLD);
printf("[%d] rcvbuf=[", rank);
for (i = 0; i < recv_sb_size; i++)
printf("%d ", rb[i]);
printf("]\n");
if (rank == 0 && status != MPI_SUCCESS) {
printf("allgatherv returned %d\n", status);
fflush(stdout);
}
xbt_free(sb);
xbt_free(rb);
xbt_free(recv_counts);
xbt_free(recv_disps);
MPI_Finalize();
return (EXIT_SUCCESS);
}
|
