Take the 2-minute tour ×
MathOverflow is a question and answer site for professional mathematicians. It's 100% free, no registration required.

I wish to solve for $x$ in

$$ (A+B)x=y $$

given square symmetric matrices $A$ and $B$. For certain reasons I have already computed the Cholesky decompositions for A and B:

$$ A = L^T L $$ $$ B = M^T M $$

Can I use these solve for $x$ more efficiently than by naively computing the Cholesky decomposition of $A+B$?

share|improve this question
    
For the Cholesky decomposition we should also assume that $A$ and $B$ are positive definite, I thought. –  Dietrich Burde May 22 '13 at 19:13
    
In general, no. One can make low-rank updates of factorization easily, but for anything different the usual answer is "you have to compute everything from scratch". Maybe you can use CG with $A$ (or $B$) as a preconditioner, but that will be numerically convenient only in some special cases. –  Federico Poloni May 22 '13 at 19:50

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Browse other questions tagged or ask your own question.