Nowadays, unprecedented computing power of desktop personal computers and efficient computational methodologies, such as the global gradient algorithm (GGA), make large water distribution system modeling feasible. However, many network analysis applications, such as optimization models, require running numerous hydraulic simulations with modified input parameters. Therefore, a methodology that could reduce the computational burden of network analysis, and still provide the required model accuracy, is needed. This paper presents a matrix transformation approach to convert the classic GGA, which is implemented within the widely