Last active
January 17, 2024 00:55
-
-
Save ehamberg/1263868 to your computer and use it in GitHub Desktop.
MPI_Scatterv example
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
#include <mpi.h> | |
#include <stdio.h> | |
#include <stdlib.h> | |
#include <math.h> | |
#define SIZE 4 | |
int main(int argc, char *argv[]) | |
{ | |
int rank, size; // for storing this process' rank, and the number of processes | |
MPI_Init(&argc, &argv); | |
MPI_Comm_rank(MPI_COMM_WORLD, &rank); | |
MPI_Comm_size(MPI_COMM_WORLD, &size); | |
int *sendcounts; // array describing how many elements to send to each process | |
int *displs; // array describing the displacements where each segment begins | |
int rem = (SIZE*SIZE)%size; // elements remaining after division among processes | |
int sum = 0; // Sum of counts. Used to calculate displacements | |
char rec_buf[100]; // buffer where the received data should be stored | |
// the data to be distributed | |
char data[SIZE][SIZE] = { | |
{'a', 'b', 'c', 'd'}, | |
{'e', 'f', 'g', 'h'}, | |
{'i', 'j', 'k', 'l'}, | |
{'m', 'n', 'o', 'p'} | |
}; | |
sendcounts = malloc(sizeof(int)*size); | |
displs = malloc(sizeof(int)*size); | |
// calculate send counts and displacements | |
for (int i = 0; i < size; i++) { | |
sendcounts[i] = (SIZE*SIZE)/size; | |
if (rem > 0) { | |
sendcounts[i]++; | |
rem--; | |
} | |
displs[i] = sum; | |
sum += sendcounts[i]; | |
} | |
// print calculated send counts and displacements for each process | |
if (0 == rank) { | |
for (int i = 0; i < size; i++) { | |
printf("sendcounts[%d] = %d\tdispls[%d] = %d\n", i, sendcounts[i], i, displs[i]); | |
} | |
} | |
// divide the data among processes as described by sendcounts and displs | |
MPI_Scatterv(&data, sendcounts, displs, MPI_CHAR, &rec_buf, 100, MPI_CHAR, 0, MPI_COMM_WORLD); | |
// print what each process received | |
printf("%d: ", rank); | |
for (int i = 0; i < sendcounts[rank]; i++) { | |
printf("%c\t", rec_buf[i]); | |
} | |
printf("\n"); | |
MPI_Finalize(); | |
free(sendcounts); | |
free(displs); | |
return 0; | |
} |
Thank you for the code, very helpfull
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Thanks a lot, very helpful ❤️