This paper presents the implementation details of a power flow method for unbalanced three-phase distribution networks. The proposed technique utilizes the current injection (CI) method for power flow calculation and considers the modeling of unbalanced loads, load tap-changing (LTC) transformers, and shunts capacitors. Compensation currents are employed for modeling changes in the settings of LTC transformers and shunt capacitors. The compensation currents eliminate the need to re-factorize the admittance matrix during the iterative solution of the power flow equations and provide a computationally efficient and fast method. The proposed method has been tested on the modified IEEE 13 node test feeder, and the presented results show that the proposed method leads to an accurate solution with a reduced computational complexity when compared to the conventional formulation.