A Three-Phase Unbalanced Load Flow Solver for Large-Scale Distribution Power Systems