#LnormInf corresponds to the absolute value of the greatest element of the vector. In this article, we are going to develop algorithm for Gauss Seidel method. For guaranteed convergence, system must be in Diagonally Dominant Form. In this method, first given system of linear equations are arranged in diagonally dominant form. Print ("The solution vector in iteration", iter1, "is:", x) The Gauss-Seidel Method is a specific iterative method, that is always using the latest estimated value for each elements in x. Gauss Seidel method is iterative approach for solving system of linear equations. def gauss_seidel(A, b, tolerance, max_iterations, x): Instead I created my own little function that with the help of a permutation matrix as seen in another answer of mine permutation matrix will produce the solution (x vector) for any square matrix, including those with zeros on the diagonal. I know this is old but, I haven't found any pre existing library in python for gauss - seidel.
0 Comments
Leave a Reply. |
Details
AuthorWrite something about yourself. No need to be fancy, just an overview. ArchivesCategories |