|Apr13-12, 11:50 AM||#1|
Gauss-Seidel Fortran 90/95
Hi, I have to program using Fortran 90/95 simultaneous linear equations. My question is, is it possible to program this? First I did the inverse matrix to find the values of the solution of a 4X4 and 3X3 matrix. I set up [A][X] = [X]' and did the inverse matrix to find the values for [X]. So I got the solutions for x1 x2 x3 x4, ( numbers are subscripts ).
Now that I know that solutions exist, can I program this with the Gauss Seidel iteration method using F 90/95? May you please explain how to set up a tolerance (say 0.000001) between two iterations because I don't understand how to since I have so many variables, or what I'm looking for? I am not asking for the specific code, just a few words to clear up how to approach this program. Thanks for your help.
|Jul30-12, 09:04 AM||#2|
About the code, sure you can.
About the tolerance, the most common stopping criteria is |rk|/|r0|≤ε where ε is the desired precision of the approximation and |rk|,|r0| the norms of the residual r=b-Ax on each iteration.
|Similar Threads for: Gauss-Seidel Fortran 90/95|
|Gauss-Seidel Method (Matlab)||Math & Science Software||2|
|Jacobi, Gauss-Seidel, SOR question||Calculus & Beyond Homework||0|
|Jacobi Vs gauss-seidel||Calculus||2|
|gauss seidel method of convergence||General Math||2|
|Gauss-Seidel Method||General Math||0|