le_sparseeuationtotalchoicegaussjordan.cpp

何光渝所著《计算机常用数值算法与程序(C++版)》的随书源代码。 是想学习C++科学计算的人士的必看代码。

//LE_SparseEuationTotalChoiceGaussJordan.cpp  
//求解大型稀疏方程组的全选主元高斯-约当消去法

#include <iostream>				//输入输出流头文件
#include "LinearEquation.h"		//线性方程(组)求解头文件

void main()
{
	int i;
    
	double a[8][8] = 
	{
		{0.0, 0.0, -1.0,   0.0,   0.0,   0.0,  2.0,  0.0},
		{0.0, 6.0,  0.0,   0.0,   0.0,  -6.0,  0.0,  0.0},
		{0.0, 0.0,  0.0,   2.0,   0.0,   0.0,  0.0, -4.0},
		{3.0, 0.0,  0.0,   0.0,  -2.0,   0.0,  1.0,  0.0},
		{0.0, 0.0,  6.0,   0.0,   0.0,   0.0,  5.0,  0.0},
		{1.0, 0.0,  0.0,   0.0,  -3.0,   0.0,  0.0,  2.0},
		{0.0, 4.0,  0.0,  -1.0,   0.0,   0.0,  0.0,  0.0},
		{0.0, 0.0,  1.0,   0.0,  -1.0,   0.0,  0.0, -2.0}
	};

	double b[8] = {4.0, 6.0, -8.0, -2.0, 27.0, -9.0, 2.0, -4.0};

	valarray<double> bb(b,8);
	matrix<double> aa(&a[0][0], 8, 8);
	    
	if (LE_SparseEuationTotalChoiceGaussJordan(aa, bb) != 0)
      for (i=0; i<8; i++)
	  {
		  cout << "x(" << i << ") = " << bb[i] << endl;
	  }
	  cout<<endl;       
}