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.
Fortran solver for the Sparse LSE problem
sparse;fortran;constrained optimization;least squares