mpi my_rank全为0而comm_sz全为1的解决办法

今天在编写mpi并行程序时,出现了一个莫名其妙的问题。

MPI_Comm_rank(MPI_COMM_WORLD,&my_rank);
MPI_COmm_size(MPI_COMM_WORLD,&comm_sz);

这两行语句分别获取通信域MPI_COMM_WORLD中的进程号和进程数量,然而我的程序在执行这两条语句后,my_rank全部变为0,comm_sz全部变为1。

在网上查阅相关问题后,找到了https://segmentfault.com/a/1190000000485292。原来是因为一时手贱,在安装了mpich后,又安装了openmpi,导致两个冲突了。于是乎,卸载掉任意一个就好了。

我卸载的是openmpi,卸载方法如下:

dpkg -l | grep openmpi      //列出有关openmpi的包
sudo apt-get remove 包名      //将列出的包全部卸载

至此,重新编译程序,程序可以正常运行了。

以下是将 MPI 梯形积分程序修改为 Pthread 梯形积分程序的代码: ```c #include <stdio.h> #include <stdlib.h> #include <pthread.h> #include <assert.h> #define intT long #define MAX_THREADS 1024 double f(double x) { return x * x; } double trap(double left_endpt, double right_endpt, intT trap_count, double base_len) { double estimate, x; estimate = (f(left_endpt) + f(right_endpt)) / 2.0; for (intT i = 1; i < trap_count; i++) { x = left_endpt + i * base_len; estimate += f(x); } estimate *= base_len; return estimate; } void* calculate_trap(void* arg) { intT my_rank = *((intT*) arg); double a = 0.0, b = 2.0, h, local_a, local_b; double local_int = 0.0, total_int = 0.0; intT n = 300000000, local_n; assert(n % MAX_THREADS == 0); h = (b - a) / n; local_n = n / MAX_THREADS; local_a = a + my_rank * local_n * h; local_b = local_a + local_n * h; local_int = trap(local_a, local_b, local_n, h); total_int += local_int; return (void*) total_int; } int main(void) { pthread_t threads[MAX_THREADS]; intT thread_args[MAX_THREADS]; void* thread_return[MAX_THREADS]; double total_int = 0.0; for (intT i = 0; i < MAX_THREADS; i++) { thread_args[i] = i; pthread_create(&threads[i], NULL, calculate_trap, (void*) &thread_args[i]); } for (intT i = 0; i < MAX_THREADS; i++) { pthread_join(threads[i], &thread_return[i]); total_int += *((double*) thread_return[i]); } printf("With a = %lf, b = %lf, n = %ld, total_int = %.10lf\n", 0.0, 2.0, 300000000, total_int); return 0; } ``` 在这个代码中,我们创建了 `MAX_THREADS` 个线程,并且每个线程都会计算一部分积分值,最后将所有线程的积分值求和得到最终的结果。由于 Pthread 不像 MPI 那样需要初始化和结束通信,因此代码更加简洁。
评论 5
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值