Fast Newton load flow

The Newton-Raphson method is widely used to solve load flow problems. Traditionally a direct solver is used to solve the linear systems within this method. In this paper we explore the use of an iterative method to solve the linear systems, leading to an inexact Newton-Krylov method. The main parameters of this method are the preconditioner and the forcing terms. Several candidate choices for these parameters are discussed and tested. With the proper preconditioner, and forcing terms, the inexact Newton-Krylov method is shown to greatly improve on using a direct solver.