?? gauss.cpp
字號:
#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);
}
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -