Report4 LU 분해법, Gauss-Jordan 소거법 with partial pivoting, 역행렬

#include <iostream> 
#include <fstream> 
#include <cmath> 
#include <iomanip>
using namespace std;
float A[3][3]={0,}; 
float B[3][1]={0,}; 
float b_[3][1]={0,}; 
float result[3][1]={0,};
float U[3][3]={0,}; 
float L[3][3]={1,0,0,0,1,0,0,0,1};
int main(void) 
{ 
	cout.setf(ios::fixed); 
	cout.precision(4); 
	ifstream fin; 
	fin.open("input4_1.txt");
	int i=0; 
	int j=0; 
	int p=0;
	int inputx=0; 
	int inputy=0;
	float x; 
	float y; 
	float z;
	fin>>inputx; 
	fin>>inputy; 
	for(int j=0;j<inputx;j++) 
	{ 
		for(int p=0;p<inputy;p++) 
		{ 
			if( p!=3) 
				fin>>A[j][p]; 
			else 
				fin>>B[j][0]; 
		} 
	}
	for(int i=0;i<3;i++) 
	{ 
		for(int j=0;j<3;j++) 
		{ 
			cout <<A[i][j]<<" "; 
		} 
		cout <<endl; 
	} 
	cout <<endl;
	float L21 = A[1][0]/A[0][0]; 
	L[1][0]=L21; 
	for(int i=0;i<3;i++) 
	{ 
		A[1][i] = A[1][i]-(L21*A[0][i]); 
	}
	for(int i=0;i<3;i++) 
	{ 
		for(int j=0;j<3;j++) 
		{ 
			cout <<A[i][j]<<" "; 
		} 
		cout <<endl; 
	} 
	cout <<endl;
	float L31 = A[2][0]/A[0][0]; 
	L[2][0]=L31; 
	for(int i=0;i<3;i++) 
	{ 
		A[2][i] = A[2][i]-(L31*A[0][i]); 
	}
	for(int i=0;i<3;i++) 
	{ 
		for(int j=0;j<3;j++) 
		{ 
			cout <<A[i][j]<<" "; 
		} 
		cout <<endl; 
	} 
	cout <<endl;
	float L32 = A[2][1]/A[1][1]; 
	L[2][1]=L32; 
	for(int i=0;i<3;i++) 
	{ 
		A[2][i] = A[2][i]-(L32*A[1][i]); 
	}
	cout<<"L"<<endl; 
	for(int i=0;i<3;i++) 
	{ 
		for(int j=0;j<3;j++) 
		{ 
			cout <<L[i][j]<<" "; 
		} 
		cout <<endl; 
	} 
	cout <<endl;
	cout<<"U"<<endl; 
	for(int i=0;i<3;i++) 
	{ 
		for(int j=0;j<3;j++) 
		{ 
			U[i][j]=A[i][j]; 
			cout <<U[i][j]<<" "; 
		} 
		cout <<endl; 
	} 
	cout <<endl;
	b_[0][0] = B[0][0]; 
	b_[1][0] = (B[1][0]-(L[1][0]*b_[0][0])); 
	b_[2][0] = B[2][0]-(L[2][0]*b_[0][0])-(L[2][1]*b_[2][0]);
	for(int i=0;i<3;i++) 
	{ 
		cout<<"b"<<i<<" : "<<b_[i][0]<<endl; 
	}
	result[2][0] = b_[2][0]/U[2][2]; 
	result[1][0] = (b_[1][0]-(U[1][2]*result[2][0]))/U[1][1]; 
	result[0][0] = (b_[0][0]-(U[0][1]*result[1][0])-(U[0][2]*result[2][0]))/U[0][0];
	for(int i=0;i<3;i++) 
	{ 
		cout<<"result"<<i<<" : "<<result[i][0]<<endl; 
	} 
	return 0; 
}
 </iomanip></cmath></fstream></iostream></textarea>
</div>

<img src="/assets/cfile29.uf.1825AA3B4D4E09682383D7.png" class="aligncenter" width="573" height="633" alt="" filename="제목 없음1.png" filemime="image/jpeg" /> 

<div>
<textarea name="code" class="language" style="color: rgb(51, 51, 51); font-family: 굴림; font-size: 9pt; line-height: 1.5; "> 
#include <iostream> 
#include <fstream> 
#include <cmath> 
#include <iomanip>
using namespace std;
float A[3][4]={0,}; 
float result[3][1]={0,};
int main(void) 
{ 
	cout.setf(ios::fixed); 
	cout.precision(4); 
	ifstream fin; 
	fin.open("input4_1.txt");
	int i=0; 
	int j=0; 
	int p=0;
	int inputx=0; 
	int inputy=0;
	fin>>inputx; 
	fin>>inputy; 
	for(int j=0;j<inputx;j++) 
	{ 
		for(int p=0;p<inputy;p++) 
		{ 
				fin>>A[j][p]; 
				cout<<A[j][p]<<" "; 
		} 
		cout<<endl; 
	} 
	cout<<endl;
	for(int i=inputx-1;i>=0;i--) 
	{ 
		if(A[i][i]!=1) 
		{ 
			cout<<i<<endl; 
			float temp = A[i][i]; 
			for(int j=0;j<inputy;j++) 
			{ 
				A[i][j] = A[i][j]/temp; 
			}
			for(int p=0;p<inputx;p++) 
			{ 
				for(int o=0;o<inputy;o++) 
				{ 
					cout <<A[p][o]<<" "; 
				} 
				cout<<endl; 
			} 
			cout<<endl; 
		}
		for(int j=0;j<inputx;j++) 
		{ 
			if(i!=j) 
			{ 
				float temp = A[j][i]; 
				for(int p=0;p<inputy;p++) 
				{ 
					A[j][p]=A[j][p]-A[i][p]*temp; 
				} 
			} 
		}
		for(int j=0;j<inputx;j++) 
		{ 
			for(int p=0;p<inputy;p++) 
			{ 
				cout<<A[j][p]<<" "; 
			} 
			cout<<endl; 
		} 
		cout<<endl; 
	}
	for(int j=0;j<inputx;j++) 
	{ 
		for(int p=3;p<inputy;p++) 
		{ 
			cout<<A[j][p]<<" "; 
		} 
		cout<<endl; 
	} 
	cout<<endl;
	return 0; 
}
 </iomanip></cmath></fstream></iostream></textarea>
</div>
<img src="/assets/cfile22.uf.1925AA3B4D4E096824D276.png" class="aligncenter" width="576" height="610" alt="" filename="제목 없음2.png" filemime="image/jpeg" />
</div>
<div>
<textarea name="code" class="language" style="color: rgb(51, 51, 51); font-family: 굴림; font-size: 9pt; line-height: 1.5; "> 
#include <iostream> 
#include <fstream> 
#include <cmath> 
#include <iomanip>
using namespace std;
float A[4][8]={0,};
int main(void) 
{ 
	cout.setf(ios::fixed); 
	cout.precision(4); 
	ifstream fin; 
	fin.open("input4_2.txt");
	int i=0; 
	int j=0; 
	int p=0;
	int inputx=0; 
	int inputy=0;
	fin>>inputx; 
	fin>>inputy; 
	for(int j=0;j<inputx;j++) 
	{ 
		for(int p=0;p<inputy;p++) 
		{ 
				fin>>A[j][p]; 
		} 
	}
	for(int j=0;j<inputx;j++) 
	{ 
		for(int p=0;p<inputy;p++) 
		{ 
			cout<<A[j][p]<<" "; 
		} 
		cout<<endl; 
	} 
	cout<<endl;
	for(int i=0;i<inputx;i++) 
	{ 
		if(A[i][i]!=1) 
		{ 
			float temp = A[i][i]; 
			for(int j=0;j<inputy;j++) 
			{ 
				A[i][j] = A[i][j]/temp; 
			}
			for(int p=0;p<inputx;p++) 
			{ 
				for(int o=0;o<inputy;o++) 
				{ 
					cout <<A[p][o]<<" "; 
				} 
				cout<<endl; 
			} 
			cout<<endl; 
		}
		for(int j=0;j<inputx;j++) 
		{ 
			if(i!=j) 
			{ 
				float temp = A[j][i]; 
				for(int p=0;p<inputy;p++) 
				{ 
					A[j][p]=A[j][p]-A[i][p]*temp; 
				} 
			} 
		}
		for(int j=0;j<inputx;j++) 
		{ 
			for(int p=0;p<inputy;p++) 
			{ 
				cout<<A[j][p]<<" "; 
			} 
			cout<<endl; 
		} 
		cout<<endl; 
	}
	for(int j=0;j<inputx;j++) 
	{ 
		for(int p=4;p<inputy;p++) 
		{ 
			cout<<A[j][p]<<" "; 
		} 
		cout<<endl; 
	} 
	cout<<endl;
	return 0; 
}
 </iomanip></cmath></fstream></iostream></textarea>
</div>
</iomanip></cmath></fstream></iostream>