MPI-jacobi迭代（雅克比迭代）C语言

最新推荐文章于 2024-01-25 11:28:21 发布
最新推荐文章于 2024-01-25 11:28:21 发布
阅读量1k 收藏 3
点赞数 1
分类专栏： 计算机体系结构 并行计算 C
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_44841518/article/details/127018504
版权

c语言 算法 矩阵

这两天完成了jacobi迭代的算法，在此记录一下心得
算法原理部分是参考的链接

主要有以下注意点：
1.这个算法是需要再来申请一个数组的，如果不申请数组的话求出来的数就是乱的，也就是代码中的B，这个B可以用N/size那么大，但是我还没改。
2.进程：每一个进程负责一部分，然后呢就要说到最后printf这个数组的时候，不能用某一个rank来printf，每个进程只改变了一块的数据，如果想输出的话，需要在每个进程中挨个输出，我之前就是最后想用0进程全部输出出来，但是发现只有0进程的数据进行了更新，其他的并没有更新，后来才改正确。
3.输出的时候最好是设置一个信号量，也就是代码中的x，让进程输出是按照顺序输出的，这样比较好检查。

这个代码中代码重复量比较大，接下来会将大量冗余代码封装成函数，然后传数据的时候采用指针来传送，还有就是，这个代码还没有测性能，这些都等着接下来再高吧！

代码如下，有大量的注释，数组用的是8*8的测试

#include<stdio.h>
#include<stdlib.h>
#include<mpi.h>
#include<time.h>

//int N = 8; // 这个代表的意思是本身的矩阵是400行的矩阵，然后在各个分块的上下各加一行，用来进程传递信息 
#define N 8
float B[N][N];


int main(int arg, char* argv[]){
	
	MPI_Init(&arg, &argv);
	int size, rank;
	MPI_Comm_size(MPI_COMM_WORLD, &size);
	MPI_Comm_rank(MPI_COMM_WORLD, &rank);
	MPI_Status status;
	
	//首先创建一下二维数组
	int i, j, k, x, y;
	float A[N + 8][N];
	int size2 = (N+8)/4;
//	float (*p)[N + 8] = A; 等完事了用指针试一下 ,直接传一个指针 
	srand(time(0));
	for(i = 0; i < N + 8; i++){
		if(i != 0 && i != size2 - 1 && i != size2 && i != 2*size2 - 1 && i != 2*size2 && i != 3*size2 - 1 && i != 3*size2 && i != 4*size2 - 1) 
			for(j = 0; j < N; j++){
				A[i][j] = rand() % 10;
			}
		else{
			for(k = 0; k < N; k++){
				A[i][k] = 0.0;
			}
		}
	}
	
	if(rank == 0){
		for(i = 0; i < N + 8;i++){
			if(i != 0 && i != size2 - 1 && i != size2 && i != 2*size2 - 1 && i != 2*size2 && i != 3*size2 - 1 && i != 3*size2 && i != 4*size2 - 1){
				printf("\n");
				for(j = 0;j < N;j++){
					printf("%f ", A[i][j]);
				}
		}
		}
	} 
	
	printf("\n");
	
	int m = 0;
	//int pr = N / 4; //因为本机是4核的， 所以就用4个进程
	if(rank == 0){ // 0号进程处理的是矩阵中第0行到第N/4 - 1行 
		for(k = 0; k < N; k++){
			MPI_Send(&A[size2 - 2][k], 1, MPI_FLOAT, rank + 1, k, MPI_COMM_WORLD); 
			
		} 
		
		for(k = 0; k < N; k++){// 接收1号进程发来的消息 
			MPI_Recv(&A[size2 - 1][k], 1, MPI_FLOAT, rank + 1, k, MPI_COMM_WORLD, &status);
		} 
		
		// 赋值+计算
		for(i = 1; i < size2 - 1;i++){
			if(i == 1){
				for(j = 0;j < N;j++){
					B[i-1][j] = A[i][j];
				}
			}
			else{
				for(j = 0;j < N;j++){
					if(j == 0 || j == N-1){
						
						B[m][j] = A[i][j];
					} 
					else{
//						B[i-1][j] = 0.25 * (A[i-1][j] + A[i][j-1] + A[i+1][j] + A[i][j+1]);
						B[m][j] = 0.25 * (A[i-1][j] + A[i][j-1] + A[i+1][j] + A[i][j+1]);
						
					}
				}
			}
			m += 1;
		}
		
		// 给1发消息阻塞一下 
		MPI_Send(&y, 1, MPI_INT, 1, 200, MPI_COMM_WORLD);
		for(i = 0;i < N/size;i++){
			printf("\n");
			for(j = 0; j < N;j++){
				printf("%f ", B[i][j]);
			}
		}
		
	}
	else if(rank == 1){
		
		for(k = 0; k < N; k++){//不整一些花里胡哨，直接把 N/4 - 1行全传过去	
//			if(k == 0){
//				printf("\n  %d \n", A[size2][k]);
//			}
			// 给2号进程发消息
			MPI_Send(&A[2 * size2 - 2][k], 1, MPI_FLOAT, rank + 1, k, MPI_COMM_WORLD); 
			// 给0号进程发消息
			MPI_Send(&A[size2 + 1][k], 1, MPI_FLOAT, rank - 1, k, MPI_COMM_WORLD); 	
			// 收0号进程给发过来的消息
			MPI_Recv(&A[size2][k], 1, MPI_FLOAT, rank - 1, k, MPI_COMM_WORLD, &status); 
			
		}
		
//		printf("\n=================================================\n");
//		for(i = 0; i < N;i++){
//			printf("%f ",A[size2][i] );
//		}	
		// 接收2号进程发来的消息 
		for(k = 0; k < N; k++){
			MPI_Recv(&A[2*size2 - 1][k], 1, MPI_FLOAT, rank + 1, k, MPI_COMM_WORLD, &status);
			
		}
		
//		printf("\n=================================================\n");
//		
//		for(i = 0; i < N;i++){
//			printf("%f ",A[2*size2 - 1][i] );
//		}	
		
		// 赋值+计算
		for(i = size2 + 1; i < 2*size2 - 1;i++){
			for(j = 0;j < N;j++){
				if(j == 0 || j == N-1){
					B[m][j] = A[i][j];
				} 
				else{
//					B[i-2 * rank - 1][j] = 0.25 * (A[i-1][j] + A[i][j-1] + A[i+1][j] + A[i][j+1]);
					B[m][j] = 0.25 * (A[i-1][j] + A[i][j-1] + A[i+1][j] + A[i][j+1]);
					
				}
			}
			m += 1;
		}
		
		MPI_Recv(&x, 1, MPI_INT, 0, 200, MPI_COMM_WORLD, &status);
		
		
//		for(i = N/size;i < 2 * N/size - 1;i++){
//			printf("\n");
//			for(j = 0; j < N;j++){
//				printf("%f ", B[i][j]);
//			}
//		}
		for(i = 0;i < N/size;i++){
			printf("\n");
			for(j = 0; j < N;j++){
				printf("%f ", B[i][j]);
			}
		}
		MPI_Send(&x, 1, MPI_INT, 2, 200, MPI_COMM_WORLD);
		
		
	
		
		
	}
	else if(rank == 2){
		// 给3号进程和2号进程发消息 
		for(k = 0; k < N; k++){//不整一些花里胡哨，直接把 N/4 - 1行全传过去
			// 接收1号进程发来的消息 
			
			MPI_Send(&A[3 * size2 - 2][k], 1, MPI_FLOAT, rank + 1, k, MPI_COMM_WORLD); 
			MPI_Send(&A[2*size2 + 1][k], 1, MPI_FLOAT, rank - 1, k, MPI_COMM_WORLD); 
			MPI_Recv(&A[2*size2][k], 1, MPI_FLOAT, rank - 1, k, MPI_COMM_WORLD, &status);
			
		} 
//		printf("\n=================================================\n");
//		for(i = 0; i < N;i++){
//			printf("%f ",A[2*size2][i] );
//		}	
		
		// 接收3号进程发来的消息 
		for(k = 0; k < N; k++){
			MPI_Recv(&A[3*size2 - 1][k], 1, MPI_FLOAT, rank + 1, k, MPI_COMM_WORLD, &status); 
		}
		
//		printf("\n=================================================\n");
//		
//		for(i = 0; i < N;i++){
//			printf("%f ",A[3*size2 - 1][i] );
//		}	
		
		// 赋值+计算
		for(i = 2*size2 + 1; i < 3*size2 - 1;i++){
			for(j = 0;j < N;j++){
				if(j == 0 || j == N-1){
					B[m][j] = A[i][j];
				} 
				else{
//					B[i-3 * rank - 1][j] = 0.25 * (A[i-1][j] + A[i][j-1] + A[i+1][j] + A[i][j+1]);
					B[m][j] = 0.25 * (A[i-1][j] + A[i][j-1] + A[i+1][j] + A[i][j+1]);
					
				}
			}
			m += 1;
		}
		
		
		MPI_Recv(&x, 1, MPI_INT, 1, 200, MPI_COMM_WORLD, &status);
		
		
//		for(i = 2*N/size + 1;i < 3*N/size - 1; i++){
//			printf("\n");
//			for(j = 0; j < N;j++){
//				printf("%f ", B[i][j]);
//			}
//		}
		for(i = 0;i < N/size;i++){
			printf("\n");
			for(j = 0; j < N;j++){
				printf("%f ", B[i][j]);
			}
		}
		MPI_Send(&x, 1, MPI_INT, 3, 200, MPI_COMM_WORLD);
		
	
		
	}
	else{
		// 给2号进程发消息 
		for(k = 0; k < N; k++){
			MPI_Send(&A[3*size2 + 1][k], 1, MPI_FLOAT, rank - 1, k, MPI_COMM_WORLD);
		}
		
		// 接收2号进程发来的消息 
		for(k = 0; k < N; k++){
			MPI_Recv(&A[3*size2][k], 1, MPI_FLOAT, rank - 1, k, MPI_COMM_WORLD, &status);
		}
		
		// 赋值+计算
		for(i = 3*size2 + 1; i < 4*size2 - 1;i++){
			if( i == 4*size2 - 2){
				for(j = 0;j < N;j++){
					B[i-4 * rank - 1][j] = A[i][j];
				}
			}
			else{
				for(j = 0;j < N;j++){
					if(j == 0 || j == N-1){
						B[m][j] = A[i][j];
					} 
					else{
//						B[i-4 * rank - 1][j] = 0.25 * (A[i-1][j] + A[i][j-1] + A[i+1][j] + A[i][j+1]);
						B[m][j] = 0.25 * (A[i-1][j] + A[i][j-1] + A[i+1][j] + A[i][j+1]);
						
					}
				}
			}
			m += 1;
		}
		
		MPI_Recv(&x, 1, MPI_INT, 2, 200, MPI_COMM_WORLD, &status);
		
//		for(i = 3* N/size + 1;i < 4*N/size - 1; i++){
//			printf("\n");
//			for(j = 0; j < N;j++){
//				printf("%f ", B[i][j]);
//			}
//		}
		for(i = 0;i < N/size;i++){
			printf("\n");
			for(j = 0; j < N;j++){
				printf("%f ", B[i][j]);
			}
		}
	}
	printf("\n");
	MPI_Barrier(MPI_COMM_WORLD);
	MPI_Finalize();
	return 0;
	
}
zHgZ
关注
  • 1
    点赞
  • 踩
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
用C语言实现雅可比迭代法.txt
02-15
用C语言实现雅可比迭代法
C语言实现雅克比迭代法求根
qq_41982200的博客
10-22 9197
C语言实现雅克比迭代法求根 雅克比迭代法求根C语言实现雅克比迭代法求根问题描述算法思想C语言程序实验结果 问题描述 设方程组Ax=bAx = bAx=b的系数矩阵AAA非奇异 ，且aii≠0{a_{ii}} \ne 0aii​​=0将AAA分裂为：A=D+L+UA = D + L + UA=D+L+U D=diag(a11,a22,⋯ann)D = diag({a_{11}},{a_{22}}, \cdots {a_{nn}})D=diag(a11​,a22​,⋯ann​), LLL和UUU分别是将AAA
使用C语言编写一个Jacobi迭代程序
weixin_42610671的博客
02-14 843
Jacobi迭代是一种解决线性方程组的数值方法，用C语言编写Jacobi迭代程序的基本步骤如下： 定义矩阵和向量：在C语言中，可以使用二维数组存储矩阵，一维数组存储向量。 输入矩阵和向量：让用户输入矩阵和向量的元素值。 初始化解向量：选择一个初始解向量，它的大小应该和方程的未知数数量相同。 进行迭代：对于每一个未知数，利用Jacobi迭代公式计算新的近似解。 检验迭代是否收敛：检查新的解...
史上最完整的MPI实现Jacobi(雅可比)迭代算法[C++]
yxg520s的博客
12-03 7649
史上最详细的MPI 雅可比迭代算法实现(c++)
雅可比迭代法v.s.高斯-赛德尔迭代的C语言实现(数值分析经典算法)
lois107的博客
12-23 8214
雅可比迭代v.s.高斯-赛德尔迭代的C语言实现算法理论雅可比迭代高斯-赛德尔迭代C语言程序接口定义及说明解读（思路）代码 算法理论 求解给定线性方程组：Ax⃗=b⃗A\vec{x}=\vec{b}Ax=b 核心思想：Ax⃗=b⃗→x⃗(k+1)=Bx⃗(k)+fA\vec{x}=\vec{b} \to \vec{x}^{(k+1)}=B\vec{x}^{(k)}+fAx=b→x(k+1)=Bx(k...
MPI并行程序雅克比迭代
02-28
这是大学课程 mpi并行程序设计实例教程上课中 讲解的雅克比迭代 的代码。
基于MPI的并行计算实现Jacobi迭代
12-13
### 基于MPI的并行计算实现Jacobi迭代 #### 概述 本文主要讨论了如何使用消息传递接口（Message Passing Interface, MPI）来实现并行计算中的Jacobi迭代法。Jacobi迭代法是一种广泛应用于求解线性方程组的数值方法...
JOCOBI_Jacobi迭代法_线性方程求解_
10-03
** Jacobi 迭代法详解 ** 在数值分析领域，Jacobi 迭代法是一种用于求解大型稀疏线性方程组的有效方法。这种方法基于矩阵的分解，尤其适用于那些对角占优或者分块对角占优的矩阵。线性方程组通常表示为 Ax = b 的...
线性方程组Jacobi迭代的并行求解
06-24
在附件中，`jacobi.c`是一个可能的C语言实现，它实现了Jacobi迭代算法。程序可能包含以下关键部分： 1. 读取输入文件`dataIn.txt`，获取矩阵A和向量b的信息。 2. 初始化迭代变量x^(0)，通常设置为零向量。 3. 循环...
mesh:MPI-IS网格处理库
05-14
感知系统网格封装该软件包包含用于操纵网格和使其可视化的核心功能。 它需要Python 3.5+并且在Linux和macOS操作系统上受支持。 Mesh处理库支持我们的一些项目，例如要求您首先需要安装Boost <http> _库。...
MPI C++ jacobi迭代
zx1772874955的博客
09-14 932
MPI C++ 实现jacobi迭代，
雅可比迭代法程序c语言,求雅可比迭代法解方程组的C\C++程序
weixin_34487308的博客
05-18 1544
满意答案singleycf2013.07.05采纳率：54%等级：13已帮助：7908人#include #include Jacobidiedai(int n, double *a, double *b,double *x){int i,j;double *x0,m=0,eps;x0= (double *) malloc(n*sizeof(double));for(i=0;ix0[i]=...
基于Jacobi迭代的OpenMP 和MPI优化并行
forrestguang的博客
06-09 1451
问题描述 求解微分方程： φxx+φyy−αφ=F.\varphi_{xx} + \varphi_{yy} - \alpha \varphi = F.φxx​+φyy​−αφ=F. 其中测试问题的时候选取φ=(1−x2)(1−y2)\varphi = (1 - x^2)(1 - y^2)φ=(1−x2)(1−y2)，计算区域是[−1,1]2[-1,1]^2[−1,1]2使用五点差分格式来求解问题: φi−1,j−2φi,j+φi+1,jh2+φi,j−1−2φi,j+φi,j+1h2−αφi,j=Fi,j.
雅可比迭代c语言实验报告,雅克比迭代法求解线性方程组的C语言程序?
weixin_35686386的博客
05-18 1090
导航：网站首页 >雅克比迭代法求解线性方程组的C语言程序?时间：2018-10-29雅克比迭代法求解线性方程组的C语言程序?相关问题:匿名网友:void Solve ( double dCoef [] ,double dY [] ,unsigned int iOrder ,double dErr){//用Jacobi迭代法解方程组,dCoef[]系数阵,Y[]向量,iOrder给出方程阶数,...
C语言线性方程组求解：LU分解、雅可比迭代、高斯-赛德尔迭代
最新发布
jk_lovers的博客
01-25 1451
（注：该程序C和C++通用）我们知道在Matlab中，矩阵的求解是非常方便的，对于Ax = b，只需要做x = A \ b的运算就可以得到x的结果，然而对于C语言，我们就需要寻找一些算法来实现线性矩阵的求解。由于本学期刚学完数值分析，并且目前做的项目中刚好有需要求解线性方程组的需求，于是将求解线性方程组的三种方法（其实还有更多，我这里仅给出三种）：LU分解、雅可比迭代以及高斯-赛德尔迭代。这三个算法的基本原理就不过多赘述了，大家想进一步学习的可以网上搜，一大堆教程，而且原理也较为简单。
用c++代码实现使用雅可比迭代法法解4阶线性代数方程的解
wgq2020的博客
11-09 302
注：该代码仅供参考，具体实现方式可能存在不同，如矩阵存储方式等。
雅可比高斯迭代法原理及其实际用法介绍举例（附C++代码）
qq_60811855的博客
07-19 7055
雅可比迭代法与高斯迭代法主要用于求解线性方程组，经过n次迭代收敛于精确值达到求解近似解的目的。本文忽视了复杂的证明过程，使读者可以更加有效的了解雅可比与高斯迭代算法的实际应用。...
C语言实现雅可比迭代（计算方法）
qq_43452446的博客
06-01 9907
#include <stdio.h> #include <stdlib.h> #define MAX 10 double A[MAX][MAX]; double b[MAX]; double X[MAX]; double Y[MAX]; double X0[MAX]; int NUM; double epsilon=1e-5; void Input_Matrix()//...
jacobi的c语言程序,用C语言实现雅可比迭代法.doc
weixin_36029647的博客
05-17 1626
用C语言实现雅可比迭代法， 源程序代码如下：#include "stdlib.h"#include "stdio.h"#include "conio.h"#include "string.h"#include "math.h"#define N 100float Table(int n,float a[N][N],float b[N]){int i,j;float c[N][N];printf("...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值