For decades, traditional numerical methods have been used for load flow analysis of electric
power systems. From the historical point of view, first numerical algorithm employed was the Gauss-
Seidel method successfully solving the system of nonlinear algebraic equations with complex
coefficients due to fair convergence properties. Although the CPU time per iteration was relatively
small, strong dependency of network size on total number of iterations excluded the Gauss-Seidel
method from simulations of larger power systems