Many large-scale numerical simulations can be broken down into common mathematical routines. While the applications may differ, the need to perform functions such as matrix solves, Fourier transforms, or eigenvalue analysis routinely arise. Consequently, targeting fast, efficient implementations of these methods will benefit a large number of applications. Graphics Processing Units (GPUs) are emerging as an attractive platform to perform these types of simulations. There FLOPS/Watt and FLOPS/dollar figures are far below competing alternatives. In previous work, EM Photonics has implemented dense matrix solvers using a hybrid GPU/multicore microprocessor approach. This has shown the ability to significantly outperform either platform when used independently. In this project, we will develop a complimentary library focused on performing routines on sparse matrices. This will be extremely valuable to a wide set of users including those doing finite-element analysis and computational fluid dynamics. Using GPUs, users are able to build single workstations with an excess of four teraFLOPS of computational power as well as create large, high-performance computing systems that are efficient in terms of both cost and power. By leveraging libraries such as the ones we will develop for this project, the user is shielded from the intricacies of GPU programming while still able to access their computational performance.