Let , which is a nonzero symmetric orthogonal projection: and . Then is symmetric and
Since , this equals the identity exactly when . Thus
If , then and the defining formula immediately gives . For arbitrary ,
Thus the Householder reflection fixes the hyperplane perpendicular to pointwise and reverses the component along its normal.
Let . Choose if and if , and set . This opposite-sign choice ensures . Because ,
Therefore
The denominator is positive, including when already lies on the first coordinate axis.
Apply successive Householder reflections supported on smaller trailing coordinate sets. At step , let be the entries in rows of column of the current symmetric matrix. If , part (c) gives an orthogonal reflection taking it to ; if , use the identity. Embed this as and replace the current matrix by .
The first columns are unchanged by right multiplication. Left multiplication turns the tail of column into a vector supported only in row . Previously treated columns have zero entries throughout this trailing coordinate set, so their zeros are preserved. Symmetry makes the matching row zeros appear as well. After the final step all entries below the first subdiagonal, and by symmetry above the first superdiagonal, vanish.
The product is orthogonal, and
Restricting each transformation to its trailing coordinates is what preserves earlier zero patterns.

Articles by others on the same topic (0)

There are currently no matching articles.