#LnormInf corresponds to the absolute value of the greatest element of the vector. You can rate examples to help us improve the quality of examples. Print ('The solution vector in iteration', iter1, 'is:', x) These are the top rated real world Python examples of pyamgrelaxation.gaussseidel extracted from open source projects.
The Poisson Equation in Any Space Dimensions def sor_method(A, b, omega, initial_guess, tolerance, max_iterations): See section 3 on the paper The Optimal Relaxation Parameter for the SOR Method Applied to Obviously, with higher omega values the number of iterations should decrease.Īs for a working algorithm on SOR this is what I have computed, where best convergence is reached when the optimal omega is used. Phi = sor_solver(A, b, omega, initial_guess, residual_convergence)įor an extended answer on omega and its uses please refer to my other answer SOR method as what is quoted below is not accurate. #An example case that mirrors the one in the Wikipedia article Here I have some python script, which solves the system of linear equations using Gauss-Seidel method: import numpy as np
