single.cpp

// thread.cpp
// It runs one sender and one receiver but they are multithreaded, with MPI and openmp.
#include <mpi.h>
#include <omp.h>
#include <iostream>
#include <chrono>
#include <string>
#include <vector>
#include "timer.h"

using namespace std;

const double M=1000000;
int len=1073741824;
int *buf=NULL;
    
int main(int argc, char *argv[]){
    int rank, size;
    int num_thread=-1;
    char hostname[20];
    int name_len;

    MPI_Init(&argc, &argv);
    
    MPI_Comm_size(MPI_COMM_WORLD, &size);
    if(size!=2){
        cout << "# of processes is not 2\n";
        MPI_Abort(MPI_COMM_WORLD, -1);
    }
    MPI_Comm_rank(MPI_COMM_WORLD, &rank);

    // Every process spawns two threads.
    // Thread 0 sends.
    // Thread 1 receives.
    buf=new int[len];
    if(len%num_thread!=0){
        cout << "Buffer size cannot be divided\n";
        return -1;
    }

    MPI_Status recv_stat;
    MPI_Request recv_req;

    Timer t;
    MPI_Barrier(MPI_COMM_WORLD);
    t.start();

    if(rank==0){
        // Send.
        MPI_Send(buf, len, MPI_INT, 1, 0, MPI_COMM_WORLD);
    }
    else{
        // Receive.
        MPI_Recv(buf, len, MPI_INT, 0, 0, MPI_COMM_WORLD, &recv_stat);
    }

    MPI_Barrier(MPI_COMM_WORLD);
    if(rank==0){
	    cout << t.seconds() << endl;
    }
    

    delete[] buf;
    MPI_Finalize();
    return 0;
}