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

I have a positive definite matrix of the form $Q+sI-\alpha J$ ($s>2, 0 < \alpha <1$ and $J$ is the all-ones matrix), where $Q$ is "nice", nonnegative and known. I'd like to know if there is a way to obtain an explicit expression for the Cholesky factorization of my matrix in this special case. Thanks!

share|improve this question

1 Answer 1

up vote 6 down vote accepted

The matrix $\alpha J$ is a rank one matrix, so there are simple update/downdate formulas for computing the Choleksy factorization of $Q+sI-\alpha J$ if you start with the factorization of $Q+sI$.

I'm not aware of any update formulas that get you from the Cholesky factorization of $Q$ to a Cholesky factorization of $Q+sI$.

share|improve this answer
    
Thanks! Can you give me a particular reference? My focus is on theory, not computation. –  Felix Goldberg Apr 25 '12 at 7:53
    
A classic reference (and the paper is available online as a free .pdf) is: P. Gill, G. Golub, W. Murray, and M. Saunders. Methods for modifying matrix factorizations. Mathematics of Computation, 126(28):505-535, 1974. stanford.edu/group/SOL/papers/ggms74.pdf –  Brian Borchers Apr 26 '12 at 1:02
    
Let me also mention that if $\alpha$ is fixed and you want to vary $s$, then you might find that a better way to go is to compute the eigenvalue decomposition of $Q-\alpha J$, and then adjust for $sI$ by adding $s$ to the eigenvalues. –  Brian Borchers Apr 26 '12 at 4:52
    
Well, actually it's $\alpha$ that's varying and $s$ is fixed. Thanks for the reference, I'll be sure to read. –  Felix Goldberg Apr 29 '12 at 14:08

Your Answer

 
discard

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

Not the answer you're looking for? Browse other questions tagged or ask your own question.