The method has a good efficiency and convergence behavior by using power iteration.
Load flow is calculated with real-number matrix operations.
Accuracy level equals to Newton–Raphson method and loop theory based method.
Link branch, slack bus and PV node can be processed together with a combined matrix.