Разреженный линейный решатель наименьших квадратов
Это отличный ответ ТАК указывает на хороший разреженный решатель дляAx=b
, но у меня есть ограничения наx
так, что каждый элемент вx
является>=0
<=N
.
Также,A
являетсяогромный (около 2e6x2e6) но очень редко<=4
элементов в строке.
Есть идеи / рекомендации? Я ищу что-то вроде MATLABlsqlin
но с огромными разреженными матрицами.
Я по сути пытаюсь решить крупномасштабнуюограниченная переменная задача наименьших квадратов на разреженных матрицах:
РЕДАКТИРОВАТЬ: ВCVX:
cvx_begin
variable x(n)
minimize( norm(A*x-b) );
subject to
x <= N;
x >= 0;
cvx_end