backup/C-Plus-Plus
1
0
Fork 0
mirror of https://github.com/TheAlgorithms/C-Plus-Plus.git synced 2026-04-03 02:29:53 +08:00
Code Issues Projects Releases Wiki Activity

attempt to document gaussian elimination

Browse Source
This commit is contained in:
Krishna Vedala
2020-05-28 23:08:24 -04:00
parent 7af2ad1eaa
commit bf3ec59ec8
1 changed files with 13 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

15
computer_oriented_statistical_methods/gaussian_elimination.cpp
View File

@@ -1,17 +1,26 @@
/**
* \file
* \brief [Gaussian elimination
* method](https://en.wikipedia.org/wiki/Gaussian_elimination)
*/
#include <iostream>
/** Main function */
int main() {
int mat_size, i, j, step;
std::cout << "Matrix size: ";
std::cin >> mat_size;
// create a 2D matrix by dynamic memory allocation
double **mat = new double *[mat_size + 1], **x = new double *[mat_size];
for (i = 0; i <= mat_size; i++) {
mat[i] = new double[mat_size + 1];
if (i < mat_size) x[i] = new double[mat_size + 1];
if (i < mat_size)
x[i] = new double[mat_size + 1];
}
// get the matrix elements from user
std::cout << std::endl << "Enter value of the matrix: " << std::endl;
for (i = 0; i < mat_size; i++) {
for (j = 0; j <= mat_size; j++) {
@@ -20,6 +29,7 @@ int main() {
}
}
// perform Gaussian elimination
for (step = 0; step < mat_size - 1; step++) {
for (i = step; i < mat_size - 1; i++) {
double a = (mat[i + 1][step] / mat[step][step]);
@@ -56,7 +66,8 @@ int main() {
for (i = 0; i <= mat_size; i++) {
delete[] mat[i];
if (i < mat_size) delete[] x[i];
if (i < mat_size)
delete[] x[i];
}
delete[] mat;
delete[] x;