Page 36 - Matrix Analysis & Applied Linear Algebra
P. 36

28               Chapter 1                                            Linear Equations

                                        Experience has shown that the following strategy for combining row scaling
                                    with column scaling usually works reasonably well.

                                                       Practical Scaling Strategy

                                         1.   Choose units that are natural to the problem and do not dis-
                                              tort the relationships between the sizes of things. These natural
                                              units are usually self-evident, and further column scaling past
                                              this point is not ordinarily attempted.
                                         2.   Row scale the system [A|b] so that the coefficient of maximum
                                              magnitude in each row of A is equal to 1. That is, divide each
                                              equation by the coefficient of maximum magnitude.
                                       Partial pivoting together with the scaling strategy described above
                                       makes Gaussian elimination with back substitution an extremely effec-
                                       tive tool. Over the course of time, this technique has proven to be reliable
                                       for solving a majority of linear systems encountered in practical work.


                                        Although it is not extensively used, there is an extension of partial pivoting
                                    known as complete pivoting which, in some special cases, can be more effective
                                    than partial pivoting in helping to control the effects of roundoff error.


                                                           Complete Pivoting
                                       If [A|b] is the augmented matrix at the k th  step of Gaussian elimina-
                                       tion, then search the pivotal position together with every position in A
                                       that is below or to the right of the pivotal position for the coefficient
                                       of maximum magnitude. If necessary, perform the appropriate row and
                                       column interchanges to bring the coefficient of maximum magnitude into
                                       the pivotal position. Shown below is the third step in a typical situation:
                                                                               
                                                          ∗∗     ∗  ∗   ∗     ∗
                                                         0  ∗   ∗  ∗   ∗     ∗ 
                                                                               
                                                         00        S   S     ∗ 
                                                                 S
                                                          00    S   S   S     ∗
                                                                               
                                                          00    S   S   S     ∗
                                       Search the positions marked “ S ” for the coefficient of maximal magni-
                                       tude. If necessary, interchange rows and columns to bring this maximal
                                       coefficient into the circled pivotal position. Recall from Exercise 1.2.13
                                       that the effect of a column interchange in A is equivalent to permuting
                                       (or renaming) the associated unknowns.
   31   32   33   34   35   36   37   38   39   40   41