A solver for large-scale sparse reconstruction