Алгоритм неотрицательного множественного линейного наименьших квадратов?

В настоящее время я использую в своем коде алгоритм множественного линейного наименьшего квадрата, но я достиг точки, в которой я действительно хотел бы наложить ограничение неотрицательности на фитинг.

я знать что это другая математическая задача, где были проведены различные исследования, и не может быть простого «лучшего» ответа. Тем не менее, я надеюсь получить помощь здесь, в SO, — несколько указателей на соответствующий код [c ++], который я мог бы использовать. Я в первую очередь не ищу решающие библиотеки, но хочу изучить и понять некоторый фрагмент кода, но любые предложения приветствуются.