并行计算之π的计算

本是古典何须时尚 2022-08-09 18:52 115阅读 0赞

#include <stdio.h>
    #include <mpi.h>
    
    double f(double x) {
    	return 4/(1+x*x);
    }
    
    double Trap (double left_endpt, double right_endpt, int trap_count, double base_len) {
    	double estimate, x;
    	int i;
    
    	estimate = (f(left_endpt) + f(right_endpt))/2.0;
    	for (i = 1; i<=trap_count-1; i++) {
    		x = left_endpt + i*base_len;
    		estimate+=f(x);
    	}
    	estimate = estimate*base_len;
    
    	return estimate;
    }
    
    int main(void) {
    	int my_rank, comm_sz, n=100000000, local_n;
    	double a=0.0, b=1.0, h, local_a, local_b;
    	double local_int, total_int;
    	int source;
    
    	MPI_Init(NULL, NULL);
    	MPI_Comm_rank(MPI_COMM_WORLD, &my_rank);
    	MPI_Comm_size(MPI_COMM_WORLD, &comm_sz);
    
    	h=(b-a)/n;
    	local_n=n/comm_sz;
    
    	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);
    
    	if (my_rank!=0) {
    		MPI_Send(&local_int, 1, MPI_DOUBLE, 0, 0, MPI_COMM_WORLD);
    	} else {
    		total_int = local_int;
    		for (source=1; source<comm_sz; source++) {
    			MPI_Recv(&local_int, 1, MPI_DOUBLE, source, 0, MPI_COMM_WORLD, MPI_STATUS_IGNORE);
    			total_int += local_int;
    		}
    	}
    
    	if(my_rank==0)  {
    		printf("with n = %d trapezoids, or estimate\n", n);
    		printf("of the integral form %f to %f = %.15e\n", a, b, total_int);
    	}
    	MPI_Finalize();
    	return 0;
    }

MPI运行效果图

![20151104141244796][]

VC6运行效果图

![20151104141308695][]

[20151104141244796]: /images/20220731/0f933798eade45118050450ea30248ac.png
[20151104141308695]: /images/20220731/d0be3d3d434643e69fa4c1440766cd20.png

发表评论取消回复

表情：

评论列表（有 0 条评论，115人围观）

还没有评论，来说两句吧...

相关阅读

相关 C++与并行计算：利用并行计算加速程序运行

目录 C++与并行计算：利用并行计算加速程序运行什么是并行计算 C++中的并行计算工具并行计算的注意事项结论 -------------------- C

系统管理员/ 2024年02月05日 12:04/ 0 赞/ 20 阅读

相关并行计算与分布式计算原理

![图片][0b9c84e94b45e12555e051f784dd2ae3.png] ![图片][7febf7e122189fc0f44196feaf4dfdfd.gif]

骑猪看日落/ 2022年10月08日 05:39/ 0 赞/ 217 阅读

相关并行计算之梯形积分

方法1 include<stdio.h> include"mpi.h" double Trap(double left_endpt,double ri

女爷i/ 2022年08月10日 08:34/ 0 赞/ 190 阅读

相关并行计算之向量的计算

include <stdio.h> include <stdlib.h> include <mpi.h> void Read_n(in

谁践踏了优雅/ 2022年08月10日 00:58/ 0 赞/ 104 阅读

相关并行计算之π的计算

include <stdio.h> include <mpi.h> double f(double x) { return 4/(1

本是古典何须时尚/ 2022年08月09日 18:52/ 0 赞/ 116 阅读

相关并行计算之HelloWorld

include <stdio.h> include "mpi.h" int main(int argc, char argv[]) {

「爱情、让人受尽委屈。」/ 2022年08月09日 18:52/ 0 赞/ 164 阅读

相关课程资源并行计算

华东师范大学数学系潘建瑜老师讲的[并行计算][Link 1] [Link 1]: http://math.ecnu.edu.cn/~jypan/Teaching/Par

男娘i/ 2022年07月21日 11:13/ 0 赞/ 185 阅读

相关 CUDA并行计算框架编程+矩阵相乘并行计算

当下的GPGPU(General Purpose GPU(Graphic Process Unit))—（CUDA: Compute Unified Device Archit

清疚/ 2022年05月25日 21:06/ 0 赞/ 214 阅读

相关分布式并行计算MapReduce

1、用自己的话阐明Hadoop平台上HDFS和MapReduce的功能、工作原理和工作过程。 MapReduce是我们再进行离线大数据处理的时候经常要使用的计算模型，Ma

落日映苍穹つ/ 2021年12月10日 13:53/ 0 赞/ 292 阅读

相关分布式并行计算MapReduce

作业要求来自于：https://edu.cnblogs.com/campus/gzcc/GZCC-16SE2/homework/3319 1. 用自己的话阐明Hadoop

谁践踏了优雅/ 2021年12月09日 03:01/ 0 赞/ 338 阅读