Modifying a Gaussian Elimination Algorithm to Perform Gauss-Jordan E.

AI Thread Summary
The discussion focuses on modifying a Gaussian elimination algorithm to achieve Gauss-Jordan elimination, specifically to reduce a matrix to reduced row-echelon form. The current algorithm successfully transforms the matrix into an upper triangular form but struggles with eliminating additional variables and scaling the diagonal entries to one. The user seeks guidance on adjusting the indices within the loops to facilitate this process. After achieving a matrix with a diagonal of ones, the next step involves eliminating the entries above the diagonal. The conversation emphasizes the need for clarity in modifying the existing algorithm to complete the Gauss-Jordan elimination process.
Illania
Messages
26
Reaction score
0

Homework Statement



I have an algorithm that implements Gaussian elimination. According to the text, with some modification of the indices and their in the loops, I should be able to have this algorithm perform Gauss-Jordan elimination. I also have to reduce the matrix to reduced row-echelon form, but for now I cannot figure out how I would go about modifying the indices to perform Gauss-Jordan elimination.

Homework Equations


The input is a matrix A[1...n, 1...n] and column vector b[1...n]
The output is an upper triangular matrix.
Code:
for i←1 to n do A[i, n+1]←b[i]
for i←1 to n−1 do
     pivot ← i
     for j←i+1 to n do
           if |A[j, i]| > |A[pivot, i]|, pivot←j
     for k←i to n+1 do
           swap(A[i, k],A[pivot, k])
     for j←i+1 to n do
           scale←A[j, i]/A[i, i]
     for k←i to n+1 do
           A[j, k]←A[j, k]−A[i, k]∗scale

The Attempt at a Solution



I have performed Gaussian Elimination on a maxtrix with 3 rows and 4 columns. This leaves a matrix with the form:
Code:
x[SUB]1[/SUB] y[SUB]1[/SUB] z[SUB]1[/SUB] b[SUB]1[/SUB]
0  y[SUB]2[/SUB] z[SUB]2[/SUB] b[SUB]2[/SUB]
0  0  z[SUB]3[/SUB] b[SUB]3[/SUB]

I understand that y1, z1, and z2 also need to be eliminated, but I can't see how to do this with the current algorithm. Could someone kindly give me a push in the right direction here?
 
Physics news on Phys.org
Once the matrix is in upper triangular form, the entries on the main diagonal need to be modified so that they all equal 1.

After the main diagonal is all 1's, then one can repeat the elimination procedure to eliminate the entries above the main diagonal, leaving the main diagonal as the only non-zero entries (with the exception of the augmentation columns).
 
Back
Top