MPI编程----计算cosx的积分

原理：

将积分区间平均分成多份，每个进程计算相应的任务，计算完成后通过根进程收集并求和得出答案并计算出相应的误差。

代码如下：

#include <stdio.h>
#include "mpi.h"
#include <string.h>
#include<math.h>
#define PI 3.1415926
double my_abs(double x)
{if(x <= 0)x = x * -1;return x;
}
int main(int argc,char **argv)
{MPI_Comm comm = MPI_COMM_WORLD;MPI_Status status;int size,rank,dest;char str[100];MPI_Init(&argc,&argv);MPI_Comm_size(comm,&size);MPI_Comm_rank(comm,&rank);int n = 100000;        //将区间分成100000份double dx = PI / 2 / n;if(rank==0){double sum,sum1,sum2,sum3,sum4;double error_ ;MPI_Recv(&sum1,1,MPI_DOUBLE,1,1,comm,&status);     //接收子进程计算结果MPI_Recv(&sum2,1,MPI_DOUBLE,2,2,comm,&status);MPI_Recv(&sum3,1,MPI_DOUBLE,3,3,comm,&status);//MPI_Recv(&sum4,1,MPI_DOUBLE,4,4,comm,&status);sum = sum1 + sum2 + sum3;error_ = (1 - my_abs(1 - my_abs(sum))) * 100;       //计算误差printf("result is : %lf\n",sum);printf("right is : %lf%% \n",error_);}else if(rank==1){double aij = 0;double sum = 0;int i = 0;//while(aij + i * dx <= PI / 2)while(i <= n)                     //计算积分{sum += dx * cos(aij + dx / 2);i++;aij += dx;}MPI_Send(&sum,1,MPI_DOUBLE,0,1,comm);//将计算结果发给根进程}else if(rank==2){double aij = PI / 2;double sum = 0;int i = 0;while(i <= n){sum += dx * cos(aij + dx / 2);i++;aij += dx;}MPI_Send(&sum,1,MPI_DOUBLE,0,2,comm);}else if(rank==3){double aij = PI;double sum = 0;int i = 0;while(i <= n){sum += dx * cos(aij + dx / 2);i++;aij += dx;}MPI_Send(&sum,1,MPI_DOUBLE,0,3,comm);}/*else if(rank==4){double aij = 3*PI /2;double sum = 0;int i = 0;while(aij + i * dx <= PI*2){sum += dx * cos(aij + dx / 2);i++;}MPI_Send(&sum,1,MPI_DOUBLE,0,4,comm);}*/MPI_Finalize();return 0;
}

运行结果：

如图所示，区间分的份数越多，误差越小。

MPI编程----计算cosx的积分相关推荐

MPI求计算cosx的定积分
计算定积分设计思想计算cos(x)中从0-2π的定积分,基本思想,分成p块,每一块都是平均分成n个小区,然后n可以初始化定义也可以输入,n理论上越小,精确值越高,用长方形的面积累加表示区间面积,利 ...
MPI编程及性能优化
第1节 MPI简介 1.1 MPI及其历史与OpenMP相似,消息传递接口(Message Passing Interface,简称MPI)是一种编程接口标准,而不是一种具体的编程语言.该标准是由消 ...
高性能MPI编程实验
文章目录一.实验目的二.实验要求三.华为云平台的功能特性和实践关键步骤四.高性能MPI编程开发关键过程五.WordCount算法的MPI实现关键点六.MPI编程与MapReduce的比较分 ...
MPI编程实例（PI、AllGather、排序）
1.MPI 简介 MPI(Message Passing Interface)是目前最重要的一个基于消息传递的并行编程工具,它具有移植性好.功能强大.效率高等许多优点,而且有多种不同的免费.高效.实用 ...
编程计算1+2+3+4+...+n的值超过500的和的第一个值
<程序设计基础-c语言>杨莉刘鸿翔 ISBN-978-7-03-032903-5 p91 习题3 12.编程计算1+2+3+4+-+n的值超过500的第一个和值为多少:n的值为多少. # ...
【MPI编程】任意数节点的树形求和（高性能计算）
简述之前写过一个只能在2幂次数的节点环境下使用的树形结构求和 [MPI编程]2幂节点的树形求和(高性能计算) 稍微修改了一下对应代码,就可以使用了.其实就是在逻辑上扩充一下为2的幂数.也就是说接受节 ...
php单词出现频率,PHP编程计算文件或数组中单词出现频率的方法
本文实例讲述了PHP编程计算文件或数组中单词出现频率的方法.分享给大家供大家参考,具体如下: 如果是小文件,可以一次性读入到数组中,使用方便的数组计数函数进行词频统计(假设文件中内容都是空格隔开的单词 ...
MPI编程的常用接口速查
获取当前时间在插入MPI提供的头文件后,可以获得获取时间的函数. double MPI_Wtime(void) 取得当前时间, 计时的精度由 double MPI_Wtick(void) 取得作为 ...
[转]C#算法一对小兔子一年后长成大兔子；一对大兔子每半年生一对小兔子。大兔子的繁殖期为4年，兔子的寿命是6年。假定第一年年初投放了一对小兔子，试编程计算，第n年末总共会有多少对兔子...
/// <summary>/// 一对小兔子一年后长成大兔子:一对大兔子每半年生一对小兔子.大兔子的繁殖期为4年,兔子的寿命是6年./// 假定第一年年初投放了一对小兔子,试编程计算,第n ...

MPI编程----计算cosx的积分

原理：

代码如下：

运行结果：

MPI编程----计算cosx的积分相关推荐

最新文章

热门文章