computes A*Q, where Q is the product of Householder transformations
Given an M by N matrix A, this function computes AQ where Q is the product of 2(N - 1) transformations
GV(N-1)*...*GV(1)*GW(1)*...*GW(N-1)
and GV(I), GW(I) are Givens rotations in the (I,N) plane which eliminate elements in the I-th and N-th planes, respectively. Q itself is not given, rather the information to recover the GV, GW rotations is supplied.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| integer | :: | m |
In - number of rows of A |
|||
| integer | :: | n |
In - number of columns of A |
|||
| real(kind=dp) | :: | a(lda,n) |
In/Out - the M by N array. On input, the matrix A to be postmultiplied by the orthogonal matrix Q. On output, the value of A*Q. |
|||
| integer | :: | lda |
In - leading dimension of A, which must not be less than M. |
|||
| real(kind=dp) | :: | v(n) |
In - contain the information necessary to recover the Givens rotations GV |
|||
| real(kind=dp) | :: | w(n) |
In - contain the information necessary to recover the Givens rotations GW |