openmpi雅可比迭代法_雅可比迭代和高斯赛德尔迭代

算法伪代码如下，有错误的话请给我指出来：

//假设所有a[i][i]都不是0，为0的话还不知道怎么办

输入系数于数组a[][];

输入常数于数组b[];

//计算新系数于数组c[][]

for(int i=0; i

c[i][n-1] = b[i]/a[i][i];

}

//不是n-1行，经过商议不处理，只是列数减1

for(int i=0; i

for(int j=0; j

if(i

c[i][j] = -a[i][j+1]/a[i][i];

else if(i>j)

c[i][j] = -a[i][j]/a[i][i];

else

continue;

}

//迭代100次，自己加上精度判断吧

//从向量形式可以看出，初试解为全0的话第一次解为c[i][n-1]，如果为了省着一次时间的话可以先把

//这个复制给解数组 root[],再仔细看看会发现每次迭代的话都需要加上c[i][n-1]

while(T--) {

//k行

for(int j=0; j

//n-1列，最后的c[i][n-1]到最后加上

int temp = 0;

for(int k=0; k

if(j>k)

temp += c[j][k]*root[k];

else if(j

temp += c[j][k]*root[k+1];

else

continue;

}

//再仔细看看会发现每次迭代的话都需要加上c[i][n-1]

root[j] += temp;

}

openmpi雅可比迭代法_雅可比迭代和高斯赛德尔迭代相关推荐

openmpi雅可比迭代法_用雅可比迭代法求线性方程组的解的并行算法（MPI）
1 //=================================================================2 //Name : 使用雅可比迭代法求解线性方程组3 //A ...
雅可比迭代和高斯—赛德尔迭代法
一.雅可比迭代法对于线性方程组AX=b,我们首先将系数矩阵A分解为对角矩阵D.下三角矩阵L和上三角矩阵U: 1.1雅可比迭代法的matlab代码在这里,我们求解下面的带状方程(以下程序均是以求解该 ...
雅可比迭代与高斯-赛德尔迭代的C++实现，及判断收敛性
随机生成若干个 nnn 阶方阵与 nnn 阶向量构成 Ax=bAx=bAx=b 分别判断J法和GS法的收敛性是否能收敛报告中应生成部分不收敛的矩阵估计J法和GS法收敛速度哪个更快实现用J法和G ...
雅可比迭代，高斯-赛德尔迭代
#include <stdio.h> #include <math.h>#define max_error 0.000125double data[3][4] = {{8,-3 ...
运用雅可比和高斯赛德尔迭代公式求解方程组，并尝试将矩阵变为主对角占优矩阵
程序描述首先要求用户输入矩阵的大小n(默认不超过10),然后再提示用户输入大小为n的方阵.因为输入的方阵可能含有较多的0元素,因此用了数据结构上的矩阵的压缩方法来存储稀疏矩阵.矩阵的每一个非零元用一 ...
如何利用MATLAB求解线型方程组--雅可比迭代法、高斯赛德尔迭代法
文章目录前言 1 直接法 2 迭代法小结前言今天我们要说的就是数值微积分,赶紧看看他和高等数学中的微积分有什么区别吧.本文是科学计算与MATLAB语言专题六第2小节的学习笔记,如果大家有时间的 ...
matlab用雅可比迭代法,如何利用MATLAB求解线型方程组--雅可比迭代法、高斯赛德尔迭代法...
文章目录前言 1 直接法 2 迭代法小结前言今天我们要说的就是数值微积分,赶紧看看他和高等数学中的微积分有什么区别吧.本文是科学计算与MATLAB语言专题六第2小节的学习笔记,如果大家有时间的 ...
python解二元一次方程组迭代法_使用python实现高斯-赛德尔迭代法的简单运算
今天传热学老师说到高斯-赛德尔迭代法,我就想拿python写一个小程序来计算. 在网上找例子,发现居然没有人拿python写,都是C / C++ / Matlab的.没有参考答案,真是写得我焦头烂额啊 ...
计算方法：列主元消去法，LU分解法，雅可比迭代法，高斯塞德尔迭代法解线性方程（C++）
Matrix.h包括矩阵类Matrix的定义,Matrix.cpp包括该类成员函数的实现,LinearEqu.h包括线性方程类LinearEqu的定义,继承自Matrix类,其中solve()方法为列 ...

openmpi雅可比迭代法_雅可比迭代和高斯赛德尔迭代

openmpi雅可比迭代法_雅可比迭代和高斯赛德尔迭代相关推荐

最新文章

热门文章