Diagonal matrix size reduction

Question

Octavian 2015-3-30

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/195158-diagonal-matrix-size-reduction

评论： Octavian 2015-3-30

Dear All,

I have a loop written for GPU (gpuarrays) with multiple matrix inputs, one of which is a diagonal matrix A with size m (A=eye(m)), with m from 10000-100000. m proves limiting as I get GPU memory errors with higher m values. A is involved in several left and right multiplication steps, and is updated (fills) with each iteration. Is there any matrix decomposition/factorization implementation in matlab that will allow rendering of A as A= a*b (or a*b*c and so on) where a size is (m,n), b size (n,m) etc (n << m) so that I can introduce intermediate code steps using a, b etc in these operations and reduce memory cost? Thank you, as always,

Octavian

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Edric Ellis 2015-3-30

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/195158-diagonal-matrix-size-reduction#answer_173168

In R2015a, sparse support was added to gpuArray - perhaps this might help you? mtimes is one of the methods that is implemented for sparse gpuArray.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Octavian 2015-3-30

Thank you Edric, I will try what you suggested. My only concern is that, even if A enters the loop as speye(m), as A fills (the number of nonzero elements increases with every iteration), its memory allocation with drastically increase with the last iterations. If there is a quick way to decompose a nonsparse nonsymmetric square matrix to reduce its size, please let me know. Thank you again,

Octavian

请先登录，再进行评论。