An interior point method for large constrained discrete ill-posed problems