An efficient method to include equality constraints in branch current distribution system state estimation