Definition
MPI_C_DOUBLE_COMPLEX is an MPI_Datatype that represents a double precision floating complex type in MPI, it corresponds to a double _Complex in C.
Datatypes
C
MPI_C_DOUBLE_COMPLEX is an MPI_Datatype that represents a double precision floating complex type in MPI, it corresponds to a double _Complex in C.
#include <stdio.h>
#include <stdlib.h>
#include <complex.h>
#include <mpi.h>
/**
* @brief Illustrate how to communicate a complex double precision float between
* 2 MPI processes.
* @details This application is meant to be run with 2 MPI processes: 1 sender
* and 1 receiver. The former sends a complex double precision float to the
* latter, which prints it.
**/
int main(int argc, char* argv[])
{
MPI_Init(&argc, &argv);
// Check that 2 MPI processes are used.
int size;
MPI_Comm_size(MPI_COMM_WORLD, &size);
if(size != 2)
{
printf("This application is meant to be run with 2 MPI processes.\n");
MPI_Abort(MPI_COMM_WORLD, EXIT_FAILURE);
}
// Get my rank and do the corresponding job.
enum role_ranks { SENDER, RECEIVER };
int my_rank;
MPI_Comm_rank(MPI_COMM_WORLD, &my_rank);
switch(my_rank)
{
case SENDER:
{
// Send the double complex
double _Complex complexToSend = I * I;
printf("[MPI process %d] I send double complex: %.1f+%.1fi.\n", my_rank, creal(complexToSend), cimag(complexToSend));
MPI_Ssend(&complexToSend, 1, MPI_C_DOUBLE_COMPLEX, RECEIVER, 0, MPI_COMM_WORLD);
break;
}
case RECEIVER:
{
// Receive the double complex
double _Complex complexReceived;
MPI_Recv(&complexReceived, 1, MPI_C_DOUBLE_COMPLEX, SENDER, 0, MPI_COMM_WORLD, MPI_STATUS_IGNORE);
printf("[MPI process %d] I received double complex: %.1f+%.1fi.\n", my_rank, creal(complexReceived), cimag(complexReceived));
break;
}
}
MPI_Finalize();
return EXIT_SUCCESS;
}
#include <stdio.h>
#include <stdlib.h>
#include <complex.h>
#include <mpi.h>
/**
* @brief Illustrate how to communicate an array of double precision floats
* between 2 MPI processes.
* @details This application is meant to be run with 2 MPI processes: 1 sender
* and 1 receiver. The former sends an array of double precision floats to the
* latter, which prints it.
**/
int main(int argc, char* argv[])
{
MPI_Init(&argc, &argv);
// Check that 2 MPI processes are used.
int size;
MPI_Comm_size(MPI_COMM_WORLD, &size);
if(size != 2)
{
printf("This application is meant to be run with 2 MPI processes.\n");
MPI_Abort(MPI_COMM_WORLD, EXIT_FAILURE);
}
// Get my rank and do the corresponding job.
enum role_ranks { SENDER, RECEIVER };
int my_rank;
MPI_Comm_rank(MPI_COMM_WORLD, &my_rank);
switch(my_rank)
{
case SENDER:
{
// Send the double complexes
double _Complex complexesToSend[2] = { I * I, 2 * I * I };
printf("[MPI process %d] I send double complexes: %.1f+%.1fi and %.1f+%.1fi.\n", my_rank, creal(complexesToSend[0]), cimag(complexesToSend[0]), creal(complexesToSend[1]), cimag(complexesToSend[1]));
MPI_Ssend(complexesToSend, 2, MPI_C_DOUBLE_COMPLEX, RECEIVER, 0, MPI_COMM_WORLD);
break;
}
case RECEIVER:
{
// Receive the double complexes
double _Complex complexesReceived[2];
MPI_Recv(complexesReceived, 2, MPI_C_DOUBLE_COMPLEX, SENDER, 0, MPI_COMM_WORLD, MPI_STATUS_IGNORE);
printf("[MPI process %d] I received double complexes: %.1f+%.1fi and %.1f+%.1fi.\n", my_rank, creal(complexesReceived[0]), cimag(complexesReceived[0]), creal(complexesReceived[1]), cimag(complexesReceived[1]));
break;
}
}
MPI_Finalize();
return EXIT_SUCCESS;
}