Parallel power flow solutions using a biconjugate gradient algorithm and a Newton method: A GPU-based approach