第一題

<span style="color:#3333FF;"><span style="background-color: rgb(255, 255, 255);">#include <stdio.h>
#include "mpi.h"

#define row_A 10  // row of Matrix A
#define col_A 20  // column of Matrix A

void rowcolcomm( int myid , int col_P)
{
	int rows=row_A, cols;//每個處理器裏面矩陣的行數，列數
	int i,j;
	int A[ row_A ][ col_A ];

	if( myid < col_A % col_P)
		cols = col_A / col_P +1;
	else
		cols = col_A / col_P;

	sleep(myid);
	printf("Process%2d裏面的行數爲：rows=%d  列數爲:cols=%d  子矩陣爲：\n",myid,rows,cols); 
    	for( i = 0 ; i < rows; i++ )
	{
     		for( j = 0 ; j < cols ; j++ )
		{          
       			A[i][j] =  i + j*col_P + myid ;
			printf("%3d",A[i][j]); 
       		}
		printf("\n");
	}	    
    	printf("\n");
}

int main( int argc, char *argv[] )
{
    	int rank;
    	int size;
    
   	MPI_Init(&argc,&argv);
    	MPI_Comm_rank(MPI_COMM_WORLD, &rank);
    	MPI_Comm_size(MPI_COMM_WORLD, &size);

    	rowcolcomm( rank,size) ;

    	MPI_Finalize();
    	return 0;
}
</span></span>