Basic HTML version of Foils prepared 6 December 96

Foil 40 Sequential LU Algorithm

From Full Matrices - December 4, 1995 CPS615 Basic Simulation Track for Computational Science -- Fall Semester 95. by Geoffrey C. Fox, Nancy J. McCracken


The algorithm proceeds in M steps.
At the start of step k, we identify the row, r, containing the largest value of |Ai,k| for k £ i < M. If r ¹ k then rows r and k are exchanged. This is called partial pivoting, and is done to improve the numerical stability. After the exchange the element that is now Ak,k is called the pivot.



© on Tue Oct 7 1997