For a nonzero real vector , choose with and put . When , the Householder reflectionsatisfies : every component except the first becomes zero. Here is the first Cartesian basis vector. The calculation is short:If , the vector is already at the selected target, so the identity matrix suffices.
In floating-point arithmetic, a customary choice isThen , avoiding subtraction of nearly equal positive quantities in the first component. Either sign can be valid in exact arithmetic; the numerical choice protects the construction from catastrophic cancellation. [ref]
https://math.mit.edu/classes/18.085/summer2016/handouts/qr.pdf. Olver, Orthogonal Bases and the QR Algorithm: Householder factorization and sign selection
For , this choice gives and . Dividing by its Euclidean norm gives , exactly the normal in the geometric example.
Articles by others on the same topic
There are currently no matching articles.