How can I use "symbfact" to speed up several consecutive Cholesky factorizations?

1 view (last 30 days)
Hi all,
I just came across the "symbfact" command. It seems powerful, but I'm wondering exactly how take advantage of it. I need to do several hundred consecutive Cholesky factorizations of matrices that have the same sparsity structure (different entries). Symbfact tells me the sparsity structure of the resulting Cholesky factors, but how can I take advantage of that?
>> % For nxn Q, n~20000
>> L = chol(Q);
>> S = symbfact(L);
>> % How can I use S to efficiently calculate subsequent L's?
Thanks!

Accepted Answer

Richard Brown
Richard Brown on 21 Jun 2013
You can't (unless something has changed recently). When I had to do this, I ended up working with the C code in CSparse directly.
  5 Comments
Alec Jacobson
Alec Jacobson on 5 Dec 2016
Has anything changed on this front? Is it still impossible to use the symbolic factorization from symbfact for repeated solves?

Sign in to comment.

More Answers (0)

Categories

Find more on Mathematics in Help Center and File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!