next up previous contents
Next: Problem with Allocatable Arrays Up: Procedures Previous: Solution

 

Pivotal LU Decomposition

In order to make your LU Decomposition algorithm secure partial pivoting must be added.

Pivoting is acheived by ensuring that the diagonal element of the row under consideration is the largest available. This means that before each row operation, the column below is searched and the row containing the largest element is swapped with the current row. A record must be kept of the new permutation so that the RHS elements can be related to the correct rows.


next up previous contents
Next: Problem with Allocatable Arrays Up: Procedures Previous: Solution

Adam Marshall ©University of Liverpool, 1996
Fri Dec 6 14:10:26 GMT 1996