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?