Есть ли способ вычислить
\argmin_{x}\|Ax-b\|_2
на основе функции, которая вычисляет матричные векторные произведения Ax без явного сохранения, разреженного или неразреженного, A в памяти?
В Python я бы использовал для этого scipy.sparse.linalg.lsqr (несмотря на имя пакета, эта функция не требует разреженных матриц, но допускает LinearOperators.
