Lets say we have a block matrix $ M =\left( \begin{array}{ccc} A & B\\\\ B^{*} & C \end{array} \right)$ where M is positive definite. (A, and C are also pos def) There is a formula for carrying out block Cholesky decomposition. See http://en.wikipedia.org/wiki/Block_LU_decomposition. Summarising we have the following result. The matrix $<math>\begin{matrix}M = LU\end{matrix}</math>$ can be decomposed in an algebraic manner into $<math>L = \begin{pmatrix} A^{\frac{1}{2}} & 0 \\\\ B^{\*} A^{-\frac{*}{2}} & Q^{\frac{1}{2}} \end{pmatrix}$ where $\begin{matrix} Q = C - B^{*} A^{-1} B \end{matrix}$ $\*$ indicates transpose in this case Now lets say we have already carried out the cholesky decomposition for A, and C. So we have already calculated $A^{1/2}$, and $C^{1/2}$ (It is therefore straightforward to calculate the inverses $A^{-1/2}$, and $C^{-1/2}$ using forward substitution). Rewriting the Q in terms of these quantities we now have. $Q = Q^{1/2}Q^{\*/2} = C^{1/2} C^{\*/2} - (B^{*} A^{-\*/2})(A^{-1/2} B)$ My question is this: Given this set up is it possible to algebraicly calculate $Q^{1/2}$ without having to apply cholesky decomposition to $Q$. Or in other words can I use $C^{1/2}$ to help me in the calculation of $Q^{1/2}$. Thanks in advance for any replies. Matt.