开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 并行计算 > 查看源码
including.c
资源名称：mpi.zip　[点击查看]
上传用户：hhyinxing
上传日期：2013-09-10
资源大小：833k
文件大小：3k
源码类别：

并行计算

开发平台：

Unix_Linux

including.c：源码内容
  1. #include <mpi.h>
  2. #include <stdio.h>
  4. int n;
  5. double xtemp[20],ytemp[20];
  6. double x,y;
  7. int s,mys;
  8. int group_size,my_rank;
  10. /*判断射线与线段是否有交点*/
  11. int cal_inter(int number,int i,double x,double y)
  12. {
  13. double x1,y1,x2,y2,temp;
  14.    int result;
  15.    result=0;
  17.    if(number+i*group_size>=n)
  18. return result;
  19. x1=xtemp[number+i*group_size];
  20.    y1=ytemp[number+i*group_size];
  21.    x2=xtemp[(number+1+i*group_size)%n];
  22.    y2=ytemp[(number+1+i*group_size)%n];
  23.  
  24.   if(y1>y2) 
  25.    {
  26.     temp=x1;
  27.     x1=x2;
  28.     x2=temp;
  29.     temp=y1;
  30.     y1=y2;
  31.     y2=temp;
  32.    }
  34.     /*判断竖直边的情况*/
  35.    if(x1==x2)
  36.    {
  37.     if((x>x1)&&(y<=y2)&&(y>y1))
  38.       result=1;
  39.     else
  40. result=0;
  41.         /*点在竖直边上，应该对result赋一个比较的大的值，这里是100*/
  42. if((x==x1)&&((y-y1)*(y2-y)>=0))
  43. result=100;
  44. }
  45. else
  46. {
  47. /*非竖直边，非水平边*/
  48.         if (y1!=y2)
  49. {
  50. temp=x2+(y-y2)*(x2-x1)/(y2-y1);
  51. /*交点刚好在边上，且不为下顶点*/
  52.       if((temp<x)&&(y<=y2)&&(y>y1))
  53.          result=1;
  54.       else
  55.          result=0;
  57.       /*点在边上，应该对result赋一个比较的大的值，这里是100*/
  58. if((temp==x)&&((y-y2)*(y1-y)>=0))
  59. result=100;
  60. }
  61. else
  62. {
  63. /*点在水平边上，应该对result赋一个比较的大的值，这里是100*/
  64.       if((y==y1)&&((x1-x)*(x-x2)>=0)) result=100;
  65. }
  66. }
  67. return result;
  68. }
  70. main(int argc,char* argv[])
  71. {
  72. int i;
  74.     MPI_Init(&argc,&argv);  
  75.     MPI_Comm_rank(MPI_COMM_WORLD,&my_rank);
  76. MPI_Comm_size(MPI_COMM_WORLD,&group_size);
  78.    /*各处理器计数器初始化*/
  79.    mys=0;
  81.     /*主处理器读入多边形顶点和要判断的点的坐标*/
  82.    if(my_rank==0)
  83.    {
  84.    printf("请输入点的个数:");
  85.   scanf("%d",&n);
  86.    printf("请输入各点的坐标n");
  87.    for(i=0;i<n;i++)
  88.     {
  89.     printf("%d:",i);
  90.     scanf("%lf",&xtemp[i]);
  91.     scanf("%lf",&ytemp[i]);
  92.     }
  93.    printf("请输入要判断点的坐标n");
  94.    scanf("%lf %lf",&x,&y);
  95.    }
  96. MPI_Barrier(MPI_COMM_WORLD);
  97. /*把多边形的顶点数、顶点坐标与要判别的点的坐标播送给所有进程*/
  98.    MPI_Bcast(&n,1,MPI_INT,0,MPI_COMM_WORLD);
  99.    MPI_Bcast(&x,1,MPI_DOUBLE,0,MPI_COMM_WORLD);
  100.    MPI_Bcast(&y,1,MPI_DOUBLE,0,MPI_COMM_WORLD);
  101.    MPI_Bcast(xtemp,n,MPI_DOUBLE,0,MPI_COMM_WORLD);
  102.    MPI_Bcast(ytemp,n,MPI_DOUBLE,0,MPI_COMM_WORLD);
  104.    MPI_Barrier(MPI_COMM_WORLD);
  105.     /*每一个处理器处理n/group_size条边上的情况并求和*/
  106.    for(i=0;i<n/group_size+1;i++)
  107.     {
  108.      mys+=cal_inter(my_rank,i,x,y); 
  109. }
  111.    MPI_Barrier(MPI_COMM_WORLD);
  113.     /*把mys的值规约到s*/
  114.     MPI_Reduce(&mys,&s,1,MPI_INT,MPI_SUM,0,MPI_COMM_WORLD);
  116.     /*根据s值确定输出结果*/
  117. if(my_rank==0)
  118.     {
  119. if(s>=100)  
  120. printf("vertex p is in polygonn");
  121.      else
  122.         if(s%2==1)
  123.          printf("vertex p is in polygonn");
  124.         else
  125.           printf("vertex p is out of polygonn");
  126. }
  127.     MPI_Finalize();
  128. }