Introduction to MPI

Aveuh

326.6K views

GitHub

Edit

Open Source Your Knowledge, Become a Contributor

Technology knowledge has to be shared and made accessible for free. Join the movement.

Create Content

Broadcasting - exercise 1

Statement

In this exercise, we are going to setup a simple program with 3 processes. Each process will have a different role. Process #0 will read data (a buffer of doubles) on stdin, and send this data to the other processes. This will be done in two steps : first the process will send the number of elements it read from stdin, then the actual buffer. We do this, so processes #1 and #2 are able to create the reception buffers with the right size. Then every process will sum some elements and print the resulting value on stdout :

Process 0 will sum and display all the elements of the buffer
Process 1 will sum and display the positive elements
Process 2 will sum and display the negative elements

Broadcasting

The basic broadcasting function in MPI is MPI_Bcast, which is defined as follows :

int MPI_Bcast(void* buffer, int count, MPI_Datatype datatype, int root, MPI_Comm comm);

There is also a non-blocking version that works as the non-blocking p2p communications (and that should be followed by MPI_Test and MPI_Wait)

int MPI_Ibcast(void *buffer, int count, MPI_Datatype datatype, int root, MPI_Comm comm, MPI_Request *request)

All parameters to these functions should be pretty obvious now. root is the id of the process sending the data, all other processes will be receivers and, as such, the values in buffer will be overwritten.

Broadcasting, exercise 1

Open Source Your Knowledge: become a Contributor and help others learn. Create New Content

#include <iostream>
#include <iomanip>
#include <mpi.h>
int main(int argc, char **argv) {
  MPI_Init(&argc, &argv);
  int rank, size;
  MPI_Comm_rank(MPI_COMM_WORLD, &rank);
  MPI_Comm_size(MPI_COMM_WORLD, &size);
  // Input reading for process 0
  int n_elements;
  double *buffer;
  if (rank == 0) {
    std::cin >> n_elements;
    buffer = new double[n_elements];
    
    for (int i=0; i < n_elements; ++i)
      std::cin >> buffer[i];
  }
  // 1- Broadcast the value n_elements to all processes
  // [...]
  // Here we create the buffers for the non-root processes now that
  // we have n_elements
  if (rank != 0)
    buffer = new double[n_elements];

Open Source Your Knowledge, Become a Contributor

16/30 Broadcasting - exercise

Broadcasting - exercise 1

Statement

Broadcasting

PYTHON: BEGINNER QUIZ (10 Questions)

Simple Python Test

PYTHON STRING METHODS: center(), rjust() and ljust().

Stored Procedures con T-SQL