sor.h

来自「本实例通过ＶＣ＋＋实现了JACOBI 和SOR 算法」· C头文件 代码 · 共 34 行

#include<iostream.h>


void sor(double w,double **a,double *b,int n,int m)//假设初值和循环次数m给出，比较W不同时迭代的好坏。
{
	double *result=new double[n];


	int i;

	for(i=1;i<=n;i++)
		result[i-1]=1;

	for(int m1=1;m1<=m;m1++)
	{	
		for(i=1;i<=n;i++)
		{
			double temp1=0;
			double temp2=0;

			for(int j=1;j<=i-1;j++)
				temp1=temp1+a[i-1][j-1]*result[j-1];
	
			for(int k=i+1;k<=n;k++)
				temp2=temp2+a[i-1][k-1]*result[k-1];

			result[i-1]=(b[i-1]-temp1-temp2)/a[i-1][i-1];
		}
	}

	for(i=1;i<=n;i++)
		cout<<result[i-1]<<" ";
	cout<<endl;
}