《并行程序设计导论》mpi中I/O处理

1.输出
几乎所有的MPI实现都允许MPI_COMM_WORLD里的所有进程都能访问标准输出stdout和标准错误输出stderr。但大部分的MPI实现并不提供对这些I/O设备访问的自动调度。也就是说如果多个进程试图写标准输出stdout，那么这些进程的输出是无法预测的，甚至会发生一个进程的输出被另一个输出打断的情况。
产生这一现象的原因是MPI进程在相互竞争，以取得对共享输出设备，标准输出stdout的访问。如果不希望进程的输出以随机顺序出现，那我们就应该按自己的想法去修改代码。例如，让除了0号进程以外的其它进程向0号进程发送他的输出，然后0号根据进程号的顺序打印输出结果。

if (my_rank!=0)
{
	sprintf(greeting,"Greetings from process %d of %d",my_rank,comm_sz);
	MPI_Send(greeting,strlen(greeting)+1,MPI_CHAR,0,0,MPI_COMM_WORLD);
}
else
{
	printf("greeting from process %d of %d",my_rank,comm_sz);
	for(int q=1;q<comm_sz;q++)
	{
		MPI_Recv(greeting,max_string,MPI_CHAR,q,0,MPI_COMM_WORLD,MPI_STATUS_IGNORE);
		printf("%s\n",greeting);
	}
}

2.输入
与输出不同，大部分的MPI实现只允许MPI_COMM_WORLD中的0号进程访问标准输入stdin。为了编写能够使用scanf的MPI程序，我们根据进程号来选取转移分支。0号进程负责读取数据，并将数据发送给其它进程。其余进程只需要recv一下。
例如，我想要从标准输入流中读取a,b,n,在main中调用一下

get_input(my_rank,comm_sz,*a,*b,*n);

他的函数内容为

void get_input(int my_rank,int comm_sz,double* a_p,double* b_p,int* n_p)
{
	int dest;
	if(my_rank==0)
	{
		printf("enter a,b, and n\n");
		scanf("%lf %lf %d",a_p,b_p,n_p);
		for(dest=1;dest<comm_sz;dest++)
		{
			MPI_Send(a_p,1,MPI_DOUBLE,dest,0,MPI_COMM_WORLD);
			MPI_Send(b_P,1,MPI_DOUBLE,dest,0,MPI_COMM_WORLD);
			MPI_Send(n_p,1,MPI_INT,dest,0,MPI_COMM_WORLD);
		}
	else
	{
		MPI_Recv(a_p,1,MPI_DOUBLE,0,0,MPI_COMM_WORLD,MPI+STATUS_IGNORE);
		MPI_Recv(b_p,1,MPI_DOUBLE,0,0,MPI_COMM_WORLD,MPI+STATUS_IGNORE);
		MPI_Recv(n_p,1,MPI_INT,0,0,MPI_COMM_WORLD,MPI+STATUS_IGNORE);
	}
}

	

这个输入和接受函数有点不好记啊！
MPI_Send函数是发送给几号进程
MPI_Recv函数是从几号进程接受