|
55 | 55 | *> where: |
56 | 56 | *> |
57 | 57 | *> P(K) is an N-by-N permutation matrix; |
58 | | -*> Q(K) is an M-by-M orthogonal matrix; |
| 58 | +*> Q(K) is an M-by-M unitary matrix; |
59 | 59 | *> R(K)_approx = ( R11(K), R12(K) ) is a rank K approximation of the |
60 | 60 | *> full rank factor R with K-by-K upper-triangular |
61 | 61 | *> R11(K) and K-by-N rectangular R12(K). The diagonal |
|
124 | 124 | *> d) RELMAXC2NRMK equals MAXC2NRMK divided by MAXC2NRM, the maximum |
125 | 125 | *> column 2-norm of the original matrix A, which is equal |
126 | 126 | *> to abs(R(1,1)), ( if K = min(M,N), RELMAXC2NRMK = 0.0 ); |
127 | | -*> e) Q(K)**H * B, the matrix B with the orthogonal |
| 127 | +*> e) Q(K)**H * B, the matrix B with the unitary |
128 | 128 | *> transformation Q(K)**H applied on the left. |
129 | 129 | *> |
130 | 130 | *> The N-by-N permutation matrix P(K) is stored in a compact form in |
131 | 131 | *> the integer array JPIV. For 1 <= j <= N, column j |
132 | 132 | *> of the matrix A was interchanged with column JPIV(j). |
133 | 133 | *> |
134 | | -*> The M-by-M orthogonal matrix Q is represented as a product |
| 134 | +*> The M-by-M unitary matrix Q is represented as a product |
135 | 135 | *> of elementary Householder reflectors |
136 | 136 | *> |
137 | 137 | *> Q(K) = H(1) * H(2) * . . . * H(K), |
|
300 | 300 | *> |
301 | 301 | *> 1. The elements below the diagonal of the subarray |
302 | 302 | *> A(1:M,1:K) together with TAU(1:K) represent the |
303 | | -*> orthogonal matrix Q(K) as a product of K Householder |
| 303 | +*> unitary matrix Q(K) as a product of K Householder |
304 | 304 | *> elementary reflectors. |
305 | 305 | *> |
306 | 306 | *> 2. The elements on and above the diagonal of |
|
579 | 579 | *> \verbatim |
580 | 580 | *> |
581 | 581 | *> November 2023, Igor Kozachenko, James Demmel, |
582 | | -*> Computer Science Division, |
583 | | -*> University of California, Berkeley |
| 582 | +*> EECS Department, |
| 583 | +*> University of California, Berkeley, USA. |
584 | 584 | *> |
585 | 585 | *> \endverbatim |
586 | 586 | * |
|
0 commit comments