Computational Methods for Large Sparse Power Systems Analysis