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