PtLU Factorization
- 최초 등록일
- 2005.01.01
- 최종 저작일
- 2005.01
- 2페이지/ 한컴오피스
- 가격 1,000원
소개글
PtLU Factorization (Decomposition) 수치해석 기법입니다. 과제하실때 활용하세요.
컴파일 실행환경
Visual C++ 6.0
1.5G 노트북
본문내용
/*
PtLU factorization.
*/
#include<stdio.h>
main(){
int i,j,k;
float M[3][3]={{0,2,3},{1,1,-1},{0,-1,1}};
float PT[3][3]={{1,0,0},{0,1,0},{0,0,1}};
float L[3][3]={{0,0,0},{0,0,0},{0,0,0}},t;
for(i=0;i<2;i++){
if(M[i][i]==0){
for(j=i+1;j<3;j++){
if(M[i][i]==0 && M[j][i]!=0){
for(k=0;k<3;k++){
t = M[i][k];
M[i][k] = M[j][k];
M[j][k] = t;
}
for(k=0;k<3;k++){
t = PT[k][i];
PT[k][i] = PT[k][j];
PT[k][j] = t;
}
}
}
}
for(j=i+1;j<3;j++){
t=M[j][i]/M[i][i];
for(k=i;k<3;k++){
M[j][k]=M[j][k]-M[i][k]*t;
}
L[j][i]=t;
}
}
for(i=0;i<3;i++)L[i][i]=1;
printf("Pt Matrix\n");
for(j=0;j<3;j++){
for(k=0;k<3;k++){
printf("%3.3f ",PT[j][k]);
}
printf("\n");
}
printf("\nL Matrix\n");
for(j=0;j<3;j++){
for(k=0;k<3;k++){
printf("%3.3f ",L[j][k]);
}
printf("\n");
}
printf("\nU Matrix\n");
for(j=0;j<3;j++){
참고 자료
없음