icpi.c

MPICH2+VC6.0 并行计算程序

//组通信
//各进程分别计算一部分分计算区间
//然后高用MPI_REDUCE(聚合通信-归约)将各进和的结果加起来


#include "mpi.h"
#include <stdio.h>
#include <math.h>

double f(double);

double f(double a)
{
    return (4.0 / (1.0 + a*a));
}

int main(int argc,char *argv[])     //将命令行参数,传入主函数,argc,argv分别是命令行参数的个数和参数数组的指针
{
    int done = 0, n, myid, numprocs, i;
    double PI25DT = 3.141592653589793238462643;
    double mypi, pi, h, sum, x;
    double startwtime = 0.0, endwtime;
    int  namelen;                        //CPU NAME length
    char processor_name[MPI_MAX_PROCESSOR_NAME];//MPI_MAX_PROCESSOR_NAME返回处理器名称最大长度
    //*******进入并行环境
    MPI_Init(&argc,&argv);//时入并行环境,启动所有进程,构建初始的MPI通信环境和通信器MPI_COMM_WORLD,MPI_COMM_SELF
    MPI_Comm_size(MPI_COMM_WORLD,&numprocs);//获得通信器中的进程数,并传入numprocs
    MPI_Comm_rank(MPI_COMM_WORLD,&myid);//获得本进程的进程号,并传入myid
    MPI_Get_processor_name(processor_name,&namelen);//获得处理器名称,存入processor_name,长度存入namelen


    printf("%d %d %d\n",namelen,numprocs,myid);

    //*******主体并行任务
    while (!done) {
        if (myid == 0) {//如果是0号进程(MASTER进程),输入区间个数 interval
            fprintf(stdout, "Enter the number of intervals: (0 /CTRL+Z quits) ");//把格式化输出写入指定流
	        fflush(stdout);//刷新指定文件
            if (scanf("%d",&n) != 1) {
		    fprintf( stdout, "No number entered; quitting\n" );
		    n = 0;
			}
	        startwtime = MPI_Wtime();//返回当前时钟,从某时刻起开始秒数为单位
		}
        MPI_Bcast(&n, 1, MPI_INT, 0, MPI_COMM_WORLD); //通信器将O号进程中的&n内容发送给通信器中的所有其他进程
		                                              //O号为主进程,1 为发送数据个数,MPI_INI为数据类型
        if (n == 0)
            done = 1;
        else {
            h   = 1.0 / (double) n;
            sum = 0.0;
            for (i = myid + 1; i <= n; i += numprocs) //将n个区间依次分配到numprocs个进程上计算
			{
                x = h * ((double)i - 0.5);
                sum += f(x);
            }//end of for
            mypi = h * sum;
			//MPI_Reduce,聚合通信_归约运算,在分布在不同进程的数据间进行指定的运算,次处为求和
			//在指写进程(根进程,此处为0号进程)中返回归约运算结果
			//CONT=1,MPI_TYPE为MPI_DOUBLE,MPI_SUM为运算符
            MPI_Reduce(&mypi, &pi, 1, MPI_DOUBLE, MPI_SUM, 0, MPI_COMM_WORLD);
            if (myid == 0) {
                printf("pi is approximately %.16f, Error is %.16f\n",pi,fabs(pi - PI25DT));
		        endwtime = MPI_Wtime();//返回当前时钟,从某时刻起开始秒数为单位
		        printf("wall clock time = %f\n", endwtime-startwtime);	       
		        fflush( stdout );
			}//end of if
		}//end of else 
    }
	//********退出并行环境
    MPI_Finalize();//MPI程序中最后一个被调用的MPI函数,正常退出必须调用该函数
    return 0;
}