OpenMPI 实现CPU并行计算模拟GA算法

#include "mpi.h" /*MPI head file*/
#include<time.h>
#include <stdio.h>
#include<stdlib.h>
#define LEN 10000
#define MAXGEN 100


//qsort()使用的比较函数。
int comp(const void*a,const void*b)
{
	return (*(int*)b-*(int*)a);
}

int main( int argc, char** argv )
{
	int i, rank, size;
	//int in[LEN];
	int *in = (int *) malloc (sizeof(int)*LEN);
	int len2;
	int rand1,rand2;
	int generation = 0;
	//int sendave[LEN], recvave[LEN];
	int *sendave = (int *) malloc (sizeof(int)*LEN);
	int	*recvave = (int *) malloc (sizeof(int)*LEN);
	int flag = 1;	//标记数组是否已经初始化
	double sum=0.0;
	MPI_Datatype arrint;
	MPI_Status status;

	MPI_Init(&argc, &argv); /*initializing */
	MPI_Comm_rank(MPI_COMM_WORLD, &rank); /*Process#*/
	MPI_Comm_size(MPI_COMM_WORLD, &size); /*Total processes#*/
	//自定义进程间通信数据类型
	MPI_Type_contiguous( LEN, MPI_INT, &arrint );
    //提交数据类型
	MPI_Type_commit( &arrint );
	srand((unsigned 
  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
下面是使用openmpi实现矩阵相乘的简单示例代码: ```c #include <stdio.h> #include <stdlib.h> #include <mpi.h> #define N 1000 int main(int argc, char** argv) { int rank, size; MPI_Init(&argc, &argv); MPI_Comm_rank(MPI_COMM_WORLD, &rank); MPI_Comm_size(MPI_COMM_WORLD, &size); int n = N / size; // 每个进程处理的行数 int *A = NULL, *B = NULL, *C = NULL; if (rank == 0) { // 初始化矩阵 A 和 B A = (int*) malloc(N * N * sizeof(int)); B = (int*) malloc(N * N * sizeof(int)); C = (int*) malloc(N * N * sizeof(int)); for (int i = 0; i < N; i++) { for (int j = 0; j < N; j++) { A[i * N + j] = i * N + j; B[i * N + j] = i * N + j; } } } // 分发矩阵 A 和 B MPI_Bcast(B, N * N, MPI_INT, 0, MPI_COMM_WORLD); int *subA = (int*) malloc(n * N * sizeof(int)); MPI_Scatter(A, n * N, MPI_INT, subA, n * N, MPI_INT, 0, MPI_COMM_WORLD); // 计算各自的部分 int *subC = (int*) malloc(n * N * sizeof(int)); for (int i = 0; i < n; i++) { for (int j = 0; j < N; j++) { subC[i * N + j] = 0; for (int k = 0; k < N; k++) { subC[i * N + j] += subA[i * N + k] * B[k * N + j]; } } } // 收集各自的部分到进程 0 MPI_Gather(subC, n * N, MPI_INT, C, n * N, MPI_INT, 0, MPI_COMM_WORLD); if (rank == 0) { // 输出结果 for (int i = 0; i < N; i++) { for (int j = 0; j < N; j++) { printf("%d ", C[i * N + j]); } printf("\n"); } free(A); free(B); free(C); } free(subA); free(subC); MPI_Finalize(); return 0; } ``` 此代码中,每个进程先从矩阵 B 中接收数据,然后从矩阵 A 中接收对应的行,计算自己的部分并将结果收集到进程 0 中。需要注意的是,这里假设进程数能够整除矩阵的行数,否则需要对代码进行一些调整。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值