Check out PETSc. It is a library that can be called from C/C++ or Fortran, and provides several sparse solvers, as well as provides mechanisms to call other software libraries (lapack, superlu ...

Given a covariance matrix, we consider the problem of maximizing the variance explained by a particular linear combination of the input variables while constraining the number of nonzero coefficients ...