Sparse solver for non-square regressor matrix.

Question

Paulo Cerqueira 2020-3-5

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/509307-sparse-solver-for-non-square-regressor-matrix

编辑： Paulo Cerqueira 2020-3-13

When using the backslash operator to solve the least squares problem

for x, where A is highly sparse and non-square, MATLAB seems to apply a fast procedure that scales quite well with the size of A, as opposed to what one would expect from a standard QR-solver for full matrices.

However, the documentation is not quite clear with regard to what the specific procedure is when the regressor is sparse but not square, and I'd like to have a notion of its computational complexity.

I would like to know if someone could provide the article/book on which the procedure applied by MATLAB's backslash is based when the regressor A is sparse and not square.

Regards.

2 个评论
显示无隐藏无

David Goodmanson 2020-3-6

编辑：David Goodmanson 2020-3-6

Hi Paulo,

see 'doc mldivide' which is really a more obscure reference than it should be. The flowchart there says nonsquare sparse uses QR, although they do not say anything more. A lot of their core algorithms they play pretty close to the vest.

Paulo Cerqueira 2020-3-6

Thank you for your reponse =)

I had been through that page, and unfortunately I would have liked a more specific description of the algorithm or its at least of its computational complexity. They provide clear references for the procedures for square regressor matrices, but not for non-square ones.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Paulo Cerqueira 2020-3-13

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/509307-sparse-solver-for-non-square-regressor-matrix#answer_420061

编辑：Paulo Cerqueira 2020-3-13

If anyone's interested, "SuiteSparseQR: Multifrontal Multithreaded Rank-Revealing Sparse QR Factorization" by Timothy Davis details the sparse QR implementation implemented in MATLAB R2009a.

Unsure if this is the current one.