#include "StdAfx.h"
#include "mpi.h"
#include<stdlib.h>
#include <stdio.h>
#include <math.h>
#include <time.h>
#include<iostream>
using namespace std;
#define N 16
#define numprocs 4
int myid, numproc;
MPI_Status status;
void display(int map[N][N])
{if (myid>0) return;
for(int i=0;i<N;i++)
{ for(int j=0;j<N;j++)
cout<<" "<<map[i][j];
cout<<endl;
}
}
void init(int map[N][N], int prob)
{ srand(time(NULL));
for(int i=0;i<N;i++)
for(int j=0;j<N;j++)
if (rand()%10<prob)
map[i][j]=1;
else map[i][j]=0;
}
void dostep(int map[N][N],int area[N+2][N/numprocs+2])
{
int send[N+2],recv[N+2];
int left_rank = (myid==0?MPI_PROC_NULL : myid-1);
int right_rank = (myid==numprocs-1?MPI_PROC_NULL :myid+1);
MPI_Recv(recv, N, MPI_INT,right_rank, right_rank,
MPI_COMM_WORLD,&status);//from right
for(int i=0;i<N+2;i++)
{
send[i]=area[i][1];
area[i][N/numprocs+2]=recv[i];
}
MPI_Send(send,N,MPI_INT,left_rank,myid,MPI_COMM_WORLD);//to left
for(int i=0;i<N+2;i++)
send[i]=area[i][N/numprocs+2];
MPI_Send(send,N,MPI_INT,right_rank,myid,MPI_COMM_WORLD);//to right
MPI_Recv(recv, N, MPI_INT,left_rank, left_rank,
MPI_COMM_WORLD,&status);//from left
for(int i=0;i<N+2;i++)
area[i][0]=recv[i];
MPI_Barrier(MPI_COMM_WORLD);
for(int i=1;i<N+1;i++)
for(int j=1;j<=N/numprocs;j++)
{int neighbors=0;
for (int k=i-1;k<=i+1; k++)
for (int l=j-1;l<=j+1; l++)
if (!(k==i&&l==j)&&area[k][l]>0)
neighbors++;
if (neighbors<2||neighbors>3)
资源评论
songchengye2015-02-05不够完善，没什么大用
hongweizou
- 粉丝: 0
- 资源: 3
最新资源
- uniapp vue3 自定义下拉刷新组件pullRefresh，带释放刷新状态、更新时间、加载动画
- WINDOWS 2003邮箱服务器搭建
- AI's prompts
- 该存储库将演示如何使用 OpenVINO 运行时 API 部署官方 YOLOv7 预训练模型.zip
- 该存储库包含使用 YOLOv9 对象检测模型和 DeepSORT 算法在视频中进行对象检测和跟踪的代码 .zip
- 论文《YOLO-ReT在边缘 GPU 上实现高精度实时物体检测》的实现.zip
- 让yolov6可以更方便的改变网络结构.zip
- springboot0桂林旅游景点导游平台(代码+数据库+LW)
- mmexportf3d00a398950f9982c0f132475da3f26_1732379945062.jpeg
- mmexport1732556836794.jpg
资源上传下载、课程学习等过程中有任何疑问或建议，欢迎提出宝贵意见哦~我们会及时处理！
点击此处反馈
