Skip to main content
11 events
when toggle format what by license comment
Apr 5, 2022 at 16:48 answer added José C Ferreira timeline score: 1
Mar 26, 2022 at 16:28 comment added wlad Uniqueness proof: $QR = Q'R'$ implies $Q^TQ' = R(R')^{-1}$. The LHS is orthogonal, and the RHS is upper triangular. By the orthogonality property, the inverse of the RHS is its transpose. But the inverse of an upper triangular matrix is upper triangular. Therefore the transpose of the RHS is upper triangular. Therefore the RHS is diagonal. A diagonal matrix which is orthogonal has $\pm 1$ along its diagonal. Therefore $R = R' \operatorname{diag}(\pm 1, \pm 1, \dotsc, \pm 1)$. Recall the claim about what $R$ is unique up to. From this, uniqueness of $Q$ is easy to show under the same conditions
Mar 26, 2022 at 16:16 comment added wlad @GottfriedHelms Make sure that the matrix R you get from your different algorithms has positive entries along its diagonal.
Mar 26, 2022 at 16:00 comment added wlad @GottfriedHelms The QR decomposition is only unique subject to the requirement that the diagonal entries of R be positive. This should be easy to fix in Matlab or Numpy, or whatever software you're using
Mar 21, 2022 at 4:51 history became hot network question
Mar 21, 2022 at 0:17 comment added Todd Trimble Does anyone know how to describe in formal terms how the "Francis function" is well-behaved over PSD matrices?
Mar 20, 2022 at 22:31 answer added Federico Poloni timeline score: 3
Mar 20, 2022 at 21:15 comment added Paul Tupper Is your Francis function uniquely defined if $M$ is singular?
Mar 20, 2022 at 21:06 history edited wlad
edited tags
Mar 20, 2022 at 20:57 history edited wlad CC BY-SA 4.0
added 42 characters in body
Mar 20, 2022 at 20:51 history asked wlad CC BY-SA 4.0