Nonnegative constraints on fitted coefficients?

[tomwenseleers]

Would it be possible by any chance to also allow constraints on the fitted coefficients to be specified (e.g. nonnegativity), as one can do in glmnet using options lower.limits and upper.limits ?