@mengxr I don't think so.  Regularisation is equivalent to tossing a new block of constraints that just say "lambda x = 0" in this context.  So it's already rolled into ata and atb before it makes it to the NNLS solver.