This paper presents design of generic Linear System Solver (LSS) for a class of large sparse symmetric matrices over real and complex numbers. These matrices correspond to either of the following 1) Symmetric Positive Definite (SPD) matrices, 2) Complex Hermitian matrices, 3) Complex matrices with SPD real and imaginary matrices. Such matrices arise in various power system analysis applications like load flow analysis and short circuit analysis. Template facility of C++ is used to write a generic program on float, double and complex data types. Design of algorithm guarantees numerical stability and efficient sparsity implementation. A reusable class SET is defined to cater to graph theoretic computations. LSS problems with matrices upto 20 000 nodes have been tested. Another feature of the proposed LSS is implementation of associative array, which allows subscripting an array with character strings, such as bus names. This helps in making the power system analysis software user friendly. The proposed LSS reflects an important development toward a truly object oriented power system analysis software.Index Terms-Graph theoretic applications, linear system solver, LU decomposition, object oriented programming, sparse matrix computations.