Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Complex Cholesky Decomposition

  1. Jul 22, 2008 #1
    I am having trouble finding information about decomposing a complex symmetric positive definite matrix. I was wondering if the cholesky decomposition would change to accommodate complex numbers. I understand that multiplying, dividing, and taking the square root of complex numbers is not the same as with real numbers, but would any additional components need to be added to the algorithm?

    I only ask because I am writing a C program and need to implement the Cholesky function. I have written code that works properly for the real decomposition, but when I try to use complex numbers, I am not getting the correct answer.

    I did read this, however...
    "All the square roots appearing are real. So if one writes a computer programme
    implementing the Cholesky factorisation one can be sure that no
    complex numbers appear in the course of the computation."

    How would this be so? The answer has imaginary parts, not on the diagonals (is that what this means), but elsewhere.
  2. jcsd
  3. Jul 22, 2008 #2
  4. Jul 22, 2008 #3
    So how would that be incorporated into the equations in the attachments...

    Attached Files:

  5. Jul 22, 2008 #4
    You would use the equations in the attachment to find L. To find L*, you take the transpose of L, and then wherever there is something like 3-2i, change it to 3+2i (5+.9i -> 5-.9i, etc)
  6. Jul 22, 2008 #5
    Thank you for you help. I was truly over complicating it.

    Thanks once again.
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook