gauss.cpp

利用gauss分解求解线性方程组

#include<stdio.h>
#include<math.h>
#define n 150
#define eps 1.0e-40
#include"down_triangle.h"
#include"up_triangle.h"
#include"Gauss_C.h"
#include"Gauss.h"
#include"Cholesky.h"
#include"improved_Cholesky.h"
#include"Guass_resolution.h"
#include"Guass_C_resolution.h"
#include"Cholesky_resolution.h"
#include"improved_Cholesky_resolution.h"
void main()
{
int i,j,u[n];
double A[n][n],B[n][n],bb[n];
//FILE *fp1,*fp2;
//fp1=fopen("F:/myfiles/Guass_equation/solution1.txt","w+");
//fp2=fopen("F:/myfiles/Guass_equation/solution2.txt","w+");

//input A, b, u   u is a vector that record the trasformation
for(i=0;i<n;i++)
{	
	for(j=0;j<n;j++)
	{
		if(j==i)
		{
		  A[i][j]=6;
		  B[i][j]=10;
		}
		else if(j==i+1)
		{
			A[i][j]=1;
			B[i][j]=1;
		}
		else if(j==i-1)
		{
			A[i][j]=8;
			B[i][j]=1;
		}
		else 
		{
			A[i][j]=0;
			B[i][j]=0;
		}
	}
}
for(i=0;i<n;i++)
 u[i]=i;
for(i=1;i<n-1;i++)
bb[i]=12;
bb[0]=11;
bb[n-1]=11;

//for(i=0;i<n;i++)
 // printf("bb[%d]=%f\n",i,bb[i]);

//Cholesky_resolution(B,bb);
//improved_Cholesky(B);
improved_Cholesky_resolution(B,bb);
//Guass_C_resolution(A,bb);
//Guass_resolution(A,bb);
//Gauss_C(A,u);

/*for(i=0;i<n;i++)
{
for(j=0;j<n;j++)
 printf("%f\t",B[i][j]);
 printf("\n");
}*/

/*for(i=0;i<n;i++)
{
  printf("u[%d]=%d\t",i,u[i]);
}*/

printf("\n");
	
for(i=0;i<n;i++)
printf("bb[%d]=%f\t",i,bb[i]);

printf("\n");

//for(i=0;i<n;i++)
//  fprintf(fp1,"bb[%d]=%f\n",i,bb[i]);


//for(i=0;i<n;i++)
//  fprintf(fp2,"bb[%d]=%f\n",i,bb[i]);
 //fclose(fp1);
// fclose(fp2);
}