Conclusion
Iterative methods are used when an exact answer is not computable or needed.
Gauss-Seidel converges faster than Jacobi but parallelism is trickier.
Finite element codes are simply systems of equations
- Solve with either Jacobi or Gauss-Seidel