The Gauss-Seidel method, introduced in Sec. 5.2 , provides convergent solutions for many problems in CFD. It is most eﬀective when a modest reduction in residual is required, e.g. as part of a steady-state solution described in Sec. 5.12 .
When the Gauss-Seidel method requires a lot of sweeps (e.g. over 10) to converge to a suitable tolerance, alternative methods may be more eﬃcient. Descent methods provide alternative matrix solvers that are often used in CFD.
Equating the gradient to zero, , corresponds to a minimum in the quadratic function. At the same time, it is the solution to . The method is therefore concerned with ﬁnding the minimum of the quadratic form eﬃciently.
For this method to work, the quadratic form must have a minimum, which requires that is symmetric and positive-deﬁnite. A positive-deﬁnite matrix is hard to visualise, but for a 2-value function it ensures the quadratic function is a paraboloid.
Diagonal dominance is the convergence condition for the Gauss-Seidel method, discussed in Sec. 5.3 . Importantly, a symmetric matrix that is diagonally dominant, and has positive diagonal coeﬃcients, is also positive-deﬁnite.