# Fortran solver for the Sparse LSE problem

I was wondering if there is a Fortran library that contains a solver for the Sparse LSE(linear equality-constrained least squares) problem $$min_{x}\|Cx-d\|^2 \text{ subject to } Ax=b$$ where $A$ and $C$ are large and sparse.

For dense linear systems, LAPACK has the subroutine DGGLSE, but I am looking for a solver for sparse systems.