jacobi.cpp

本程序采用并行方法实现jacobi迭代算法

// jacobi.cpp : 定义控制台应用程序的入口点。
//

#include "stdafx.h"


int _tmain(int argc, _TCHAR* argv[])
{
	int thread_n = 2,k = 1;
	float x1 = 1.0,x2 = 1.0,x3 = 1.0;
	float x1_n,x2_n,x3_n;
	float a[3][3] = {{4,3,0},{3,4,-1},{0,-1,4}};
	float b[3]={24,30,-24};

    while(k)
	{
		x1_n = (float)(24 - 3*x2)/4;
		x2_n = (float)(30 - 3*x1 + x3)/4;
		x3_n = (float)( -24 + x2)/4;

        k--;

		x1 = x1_n;
		x2 = x2_n;
		x3 = x3_n;
	}

	printf("x1 = %d,x2 = %d,x3 = %d\n",x1,x2,x3);

	return 0;
}