Drexel dragonThe Math ForumDonate to the Math Forum



Search All of the Math Forum:

Views expressed in these public forums are not endorsed by Drexel University or The Math Forum.


Math Forum » Discussions » sci.math.* » sci.math.num-analysis.independent

Topic: Woodbury formula in practice
Replies: 6   Last Post: Oct 2, 2012 3:20 PM

Advanced Search

Back to Topic List Back to Topic List Jump to Tree View Jump to Tree View   Messages: [ Previous | Next ]
srobinson3130@gmail.com

Posts: 4
Registered: 9/27/12
Woodbury formula in practice
Posted: Sep 27, 2012 11:08 AM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply

Hello,

I am trying to use the Woodbury formula [1] in order to apply corrections on the inverse of a matrix A I already calculated once as X (rather big: n=6000+). Once in a while I would like to compute the inverse of A with a "k" changes ranging from 2-3 to 20.

Applying the woodbury formula should be of big help there, and it works perfectly when I only apply 2 or 3 changes. But when I apply 20, one of the matrix I need to invert in the calculation (kind of known as the capacitance matrix [1]) becomes singular (determinant gets way too small). I am thus unable to apply the formula. Note that the corrections applied are valid, i.e. if I create the new matrix from scratch, it is non-singular.

My questions are the following:
- is it something I should expect?
- are there anyway to solve it?

I have been looking everywhere for practical information regarding woodbury's usage in a numerical analysis context but can not find much information. I would be very grateful for any guidance on that topic.


SR

[1] http://en.wikipedia.org/wiki/Woodbury_matrix_identity




Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© Drexel University 1994-2014. All Rights Reserved.
The Math Forum is a research and educational enterprise of the Drexel University School of Education.