Parallel computation of large power system networks using the multi-area Thévenin equivalents