1 / 10

Linear Systems Pivoting in Gaussian Elim.

Linear Systems Pivoting in Gaussian Elim. CSE 541 Roger Crawfis. Limitations of Gaussian Elimation. The naïve implementation of Gaussian Elimination is not robust and can suffer from severe round-off errors due to: Dividing by zero Dividing by small numbers and adding.

lonna
Download Presentation

Linear Systems Pivoting in Gaussian Elim.

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Linear SystemsPivoting in Gaussian Elim. CSE 541 Roger Crawfis

  2. Limitations of Gaussian Elimation • The naïve implementation of Gaussian Elimination is not robust and can suffer from severe round-off errors due to: • Dividing by zero • Dividing by small numbers and adding. • Both can be solved with pivoting

  3. Factored Portion Row i Row j Partial Pivoting • What if at step i, Aii = 0? • Simple Fix: If Aii = 0 Find Aji 0 j > i Swap Row j with i

  4. Forward Elimination Example – Partial Pivoting

  5. Forward Elimination Rounded to 3 digits Example – Partial Pivoting

  6. Better Pivoting • Partial Pivoting to mitigate round-off error • Adds an O(n) search. Avoids Small Multipliers

  7. Partial Pivoting swap Forward Elimination Rounded to 3 digits

  8. k k k k Pivoting strategies • Partial Pivoting: • Only row interchange • Complete (Full) Pivoting • Row and Column interchange • Threshold Pivoting • Only if prospective pivot is found to be smaller than a certain threshold

  9. Pivoting With Permutations • Adding permutation matrices in the mix: • However, in Gaussian Elimination we will only swap rows or columns below the current pivot point. This implies a global reordering of the equations will work:

  10. Pivoting • Again, the pivoting is strictly a function of the matrix A, so once we determine P it is trivial to apply it to many problems bk. • For LU factorization we have: • LU = PA • Ly = Pb • Ux = y

More Related