How can I use LU factorization in the most memory saving way?
조회 수: 3 (최근 30일)
이전 댓글 표시
Hello,
I would like to solve linear equations of the form Ax=b by means of LU decomposition. As the solution is seeked for several right hand side vectors b, I would like to store the decomposition of A such that it can be reused.
As the matrix A is very big, the otherwise very usefull package "factorize" from file exchange causes problems, as it returns a structure containing A, L and U as full matrices, which don't fit into the memory together.
The Matlab function lu can be called in a way that the L and U matrices are put together into one matrix Y = lu(A) with Y = U+L-eye(size(A)). This way, only A and one other matrix of the same size need to fit into the memory. However, in this case the permutation matrix P is not returned and I don't know how I can solve for x without it.
Do you have a suggestion how I can resolve this problem?
Thanks in advance, Amos
댓글 수: 0
답변 (0개)
참고 항목
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!