MPI 实现雅可比迭代法

使用逼近的方法求方程组的解 只使用第i个方程计算第i个x的值 必须是对角占优的矩阵(否则结果会是正负无穷)

每个线程负责某个x的计算 线程对整个矩阵都能完全直接访问 计算结果分别输出

#include "mpi.h"
#include 
#include 

//broadcast x
void bcastx(float *xs,int size) {
	for(int i=0;ithd?1:0;
		//update all xs
		bcastx(xs,size);
	}
	printf("%d,%f \n",self,nx);
	free(xs);
	MPI_Finalize();
	return 0;
}


你可能感兴趣的:(做并行,float)