Untitled

 avatar
unknown
c_cpp
2 years ago
2.6 kB
4
Indexable
#include <mpi.h>
#include <stdio.h>
#include <stdlib.h>

#define ARRAY_SIZE 20

int main(int argc, char *argv[])
{
	/* Break into many part (3 part, init input, cal, 
							gather output like in the class, but actually they are the same in here)

	*/

	int rank, size;


	MPI_Init(&argc, &argv);
	MPI_Comm_rank(MPI_COMM_WORLD, &rank);
	MPI_Comm_size(MPI_COMM_WORLD, &size);

	int break_part_size = ARRAY_SIZE / size;


	if (rank == 0){
		// Init A, B, C 
		int *A,*B,*C, i;
		A = (int *)malloc(ARRAY_SIZE*sizeof(int));
		B = (int *)malloc(ARRAY_SIZE*sizeof(int));
		C = (int *)malloc(ARRAY_SIZE*sizeof(int));

		for (i=0;i < ARRAY_SIZE;i++){
			*(A+i) = i+1;
			*(B+i) = 2*(i+1);
		}

		for (i=0; i < break_part_size;i++ ){
			*(C+i)  = *(A+i) + *(B+i);
		}

		
		// Send all except last piece
		for (i=1; i<size-1; i++){
			MPI_Send(A + break_part_size * i,
					 break_part_size,
					 MPI_INT,
					 i,
					 i,
					 MPI_COMM_WORLD
					);

			MPI_Send(B + break_part_size * i,
					 break_part_size,
					 MPI_INT,
					 i,
					 i+1000,
					 MPI_COMM_WORLD
					);
		}
  
		// Send last piece
		MPI_Send(A + break_part_size * (size - 1),
				 ARRAY_SIZE - break_part_size * (size - 1), 
				 MPI_INT,
				 size-1,
				 size-1,
				 MPI_COMM_WORLD);


		MPI_Send(B + break_part_size * (size - 1),
				 ARRAY_SIZE - break_part_size * (size - 1), 
				 MPI_INT,
				 size-1,
				 size+999,
				 MPI_COMM_WORLD);



		for (i=1; i<size-1; i++){
			MPI_Recv(C + break_part_size * i,
					 break_part_size,
					 MPI_INT,
					 i,
					 0,
					 MPI_COMM_WORLD,
					 MPI_STATUS_IGNORE
					);

		}

		MPI_Recv(C + break_part_size * (size - 1),
				 ARRAY_SIZE - break_part_size * (size - 1), 
				 MPI_INT,
				 size-1,
				 0,
				 MPI_COMM_WORLD,
				 MPI_STATUS_IGNORE);

		
		// print C
		printf("From %d rank, C =\n", rank);
		for (i=0; i<ARRAY_SIZE;i++){
			printf("%d ", *(C+i));
		}
		printf("\n");
	}




	else {
		int *As, *Bs,*Cs, i, size_part;
		if (rank == size-1) {
			size_part =  ARRAY_SIZE - break_part_size * (size - 1);
		}

		else{
			size_part = break_part_size;
		}

		As = (int *)malloc(size_part*sizeof(int));
		Bs = (int *)malloc(size_part*sizeof(int));
		Cs = (int *)malloc(size_part*sizeof(int));

		

		MPI_Recv(As,
			size_part,
			MPI_INT,
			0,
			rank,
			MPI_COMM_WORLD,
			MPI_STATUS_IGNORE
			);

		MPI_Recv(Bs,
			size_part,
			MPI_INT,
			0,
			rank + 1000,
			MPI_COMM_WORLD,
			MPI_STATUS_IGNORE
			);

		for(i=0; i<size_part;i++){
			*(Cs+i)=  *(As+i) + *(Bs+i);
		}

		MPI_Send(Cs,
			size_part,
			MPI_INT,
			0,
			0,
			MPI_COMM_WORLD
			);

	}

	MPI_Finalize();
	return 0;
}
Editor is loading...