并行计算求解三对角线性方程组

本程序是C++实现三对角线性方程组的求解，用的是OpenMP的并行算法！

#include 
#include  
//函数声明
void Speed_up（int ndouble **Adouble *Bdouble *x）;
template 
T** Allocation2D（int Mint N）;
template 
void Delete2D（int MT **p）;

void main（）//Hermite配点法解吸模型
{
	int ijn;
	n=3;//矩阵阶数
	double a[3][3]={{110}{111}{011}};
	double b[3]={011};
	double *x;
	x=new double [n];
	
	//Speed_up（）函数要求输入动态数组因此
	double **A*B;
	A=Allocation2D（nn）;
	B=new double [n];	
	for（i=0;i	{
		for（j=0;j		B[i]=b[i];
	}

	//Speed_up（）函数调用
	Speed_up（nABx）;

	//输出结果
	for（i=0;i	{
		cout<<“x“<	}
	Delete2D（nA）;
	delete [] B;
}

void Speed_up（int ndouble **Adouble *Bdouble *x）
{
	int i;
	double **a;
	a=Allocation2D（n3）;
	for（i=0;i	{
		x[i]=B[i];
		if（i>0）
		{
			a[i][0]=A[i][0];
		    a[i][1]=A[i][1];
			if（i		}	
	}
	a[0][0]=0; a[0][1]=A[0][0]; a[0][2]=A[0][1]; a[n-1][2]=0;
    for（i=1;i	{
        a[i][0]/=a[i-1][1];
        a[i][1]-=a[i][0]*a[i-1][2];
        if（fabs（a[i][1]）==0） 
		{
			Delete2D（na）;
		    return;
		}
        x[i]-=a[i][0]*x[i-1];
    }
    x[n-1]/=a[n-1][1];
    for（i=n-2;i>=0;i--） x[i]=（x[i]-a[i][2]*x[i+1]）/a[i][1];	
	Delete2D（na）;
}

template 
T** Allocation2D（int Mint N） 
{
    T **a;
    a=new T*[M];
    for （int i=0;i	{
		a[i]=new T[N];
	}
    return a;
}

template 
void Delete2D（int MT **p）
{
	int i;
	for（i=0;i	{
		delete [] p[i];
	}
	delete [] p;
}

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2011-11-15 17:42  各种三对角线性方程组的求解\
     文件        2307  2004-09-25 13:29  各种三对角线性方程组的求解\16.txt
     目录           0  2011-11-15 17:41  各种三对角线性方程组的求解\CH1\
     文件        1660  1996-02-09 17:48  各种三对角线性方程组的求解\CH1\ABAND.C
     文件         660  1996-02-09 17:55  各种三对角线性方程组的求解\CH1\ABAND0.C
     文件         416  1996-02-09 17:51  各种三对角线性方程组的求解\CH1\ABAND0.DAT
     文件        1198  1996-02-11 00:58  各种三对角线性方程组的求解\CH1\ABINT.C
     文件         603  1996-02-11 01:01  各种三对角线性方程组的求解\CH1\ABINT0.C
     文件          84  1996-02-11 01:02  各种三对角线性方程组的求解\CH1\ABINT0.DAT
     文件        2529  1996-02-09 17:14  各种三对角线性方程组的求解\CH1\ACGAS.C
     文件         692  1996-02-09 17:19  各种三对角线性方程组的求解\CH1\ACGAS0.C
     文件         152  1996-02-09 17:21  各种三对角线性方程组的求解\CH1\ACGAS0.DAT
     文件        1394  1996-02-09 17:56  各种三对角线性方程组的求解\CH1\ACHOL.C
     文件         418  1996-02-09 17:58  各种三对角线性方程组的求解\CH1\ACHOL0.C
     文件         148  1996-02-09 17:58  各种三对角线性方程组的求解\CH1\ACHOL0.DAT
     文件        2405  1996-02-09 17:25  各种三对角线性方程组的求解\CH1\ACJDN.C
     文件         873  1996-02-09 17:28  各种三对角线性方程组的求解\CH1\ACJDN0.C
     文件         284  1996-02-09 17:28  各种三对角线性方程组的求解\CH1\ACJDN0.DAT
     文件        1712  1996-02-09 16:48  各种三对角线性方程组的求解\CH1\AGAUS.C
     文件         431  1996-02-09 16:49  各种三对角线性方程组的求解\CH1\AGAUS0.C
     文件          72  1996-02-09 16:52  各种三对角线性方程组的求解\CH1\AGAUS0.DAT
     文件        1622  1996-02-09 18:00  各种三对角线性方程组的求解\CH1\AGGJE.C
     文件         736  1996-02-09 18:01  各种三对角线性方程组的求解\CH1\AGGJE0.C
     文件         160  1996-02-09 18:01  各种三对角线性方程组的求解\CH1\AGGJE0.DAT
     文件        1769  1996-02-09 17:08  各种三对角线性方程组的求解\CH1\AGJDN.C
     文件         496  1996-02-09 17:20  各种三对角线性方程组的求解\CH1\AGJDN0.C
     文件         144  1996-02-09 17:20  各种三对角线性方程组的求解\CH1\AGJDN0.DAT
     文件         346  1996-02-10 19:05  各种三对角线性方程组的求解\CH1\AGMIV.C
     文件         744  1996-02-10 19:06  各种三对角线性方程组的求解\CH1\AGMIV0.C
     文件         245  1996-02-10 19:07  各种三对角线性方程组的求解\CH1\AGMIV0.DAT
     文件         579  1996-02-10 19:00  各种三对角线性方程组的求解\CH1\AGMQR.C
............此处省略377个文件信息