An iterative algorithm for large size least-squares constrained regularization problems