ungl2#

cungl2

Functions

void cungl2(
    const INT           m,
    const INT           n,
    const INT           k,
          c64* restrict A,
    const INT           lda,
    const c64* restrict tau,
          c64* restrict work,
          INT*          info
);

void cungl2(const INT m, const INT n, const INT k, c64 *restrict A, const INT lda, const c64 *restrict tau, c64 *restrict work, INT *info)#

CUNGL2 generates an m-by-n complex matrix Q with orthonormal rows, which is defined as the first m rows of a product of k elementary reflectors of order n.

  Q  =  H(k)**H . . . H(2)**H H(1)**H

as returned by CGELQF.

Parameters

in

m

The number of rows of the matrix Q. m >= 0.

in

n

The number of columns of the matrix Q. n >= m.

in

k

The number of elementary reflectors whose product defines the matrix Q. m >= k >= 0.

inout

A

Single complex array, dimension (lda, n). On entry, the i-th row must contain the vector which defines the elementary reflector H(i), for i = 0,1,…,k-1, as returned by CGELQF in the first k rows of its array argument A. On exit, the m by n matrix Q.

in

lda

The first dimension of the array A. lda >= max(1, m).

in

tau

Single complex array, dimension (k). TAU(i) must contain the scalar factor of the elementary reflector H(i), as returned by CGELQF.

out

work

Single complex array, dimension (m).

out

info

= 0: successful exit
< 0: if info = -i, the i-th argument had an illegal value.

zungl2

Functions

void zungl2(
    const INT            m,
    const INT            n,
    const INT            k,
          c128* restrict A,
    const INT            lda,
    const c128* restrict tau,
          c128* restrict work,
          INT*           info
);

void zungl2(const INT m, const INT n, const INT k, c128 *restrict A, const INT lda, const c128 *restrict tau, c128 *restrict work, INT *info)#

ZUNGL2 generates an m-by-n complex matrix Q with orthonormal rows, which is defined as the first m rows of a product of k elementary reflectors of order n.

  Q  =  H(k)**H . . . H(2)**H H(1)**H

as returned by ZGELQF.