orgtr#

Functions

void sorgtr(
    const char*          uplo,
    const INT            n,
          f32*  restrict A,
    const INT            lda,
    const f32*  restrict tau,
          f32*  restrict work,
    const INT            lwork,
          INT*           info
);
void sorgtr(const char *uplo, const INT n, f32 *restrict A, const INT lda, const f32 *restrict tau, f32 *restrict work, const INT lwork, INT *info)#

SORGTR generates a real orthogonal matrix Q which is defined as the product of n-1 elementary reflectors of order N, as returned by SSYTRD:

if UPLO = ‘U’, Q = H(n-1) … H(2) H(1),

if UPLO = ‘L’, Q = H(1) H(2) … H(n-1).

Parameters

in
uplo

= ‘U’: Upper triangle of A contains elementary reflectors from SSYTRD; = ‘L’: Lower triangle of A contains elementary reflectors from SSYTRD.

in
n

The order of the matrix Q. n >= 0.

inout
A

Double precision array, dimension (lda, n). On entry, the vectors which define the elementary reflectors, as returned by SSYTRD. On exit, the N-by-N orthogonal matrix Q.

in
lda

The leading dimension of the array A. lda >= max(1, n).

in
tau

Double precision array, dimension (n-1). TAU(i) must contain the scalar factor of the elementary reflector H(i), as returned by SSYTRD.

out
work

Double precision array, dimension (max(1, lwork)). On exit, if info = 0, work[0] returns the optimal lwork.

in
lwork

The dimension of the array work. lwork >= max(1, n-1). For optimum performance lwork >= (n-1)*nb, where nb is the optimal blocksize. If lwork = -1, then a workspace query is assumed.

out
info

  • = 0: successful exit

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

Functions

void dorgtr(
    const char*          uplo,
    const INT            n,
          f64*  restrict A,
    const INT            lda,
    const f64*  restrict tau,
          f64*  restrict work,
    const INT            lwork,
          INT*           info
);
void dorgtr(const char *uplo, const INT n, f64 *restrict A, const INT lda, const f64 *restrict tau, f64 *restrict work, const INT lwork, INT *info)#

DORGTR generates a real orthogonal matrix Q which is defined as the product of n-1 elementary reflectors of order N, as returned by DSYTRD:

if UPLO = ‘U’, Q = H(n-1) … H(2) H(1),

if UPLO = ‘L’, Q = H(1) H(2) … H(n-1).

Parameters

in
uplo

= ‘U’: Upper triangle of A contains elementary reflectors from DSYTRD; = ‘L’: Lower triangle of A contains elementary reflectors from DSYTRD.

in
n

The order of the matrix Q. n >= 0.

inout
A

Double precision array, dimension (lda, n). On entry, the vectors which define the elementary reflectors, as returned by DSYTRD. On exit, the N-by-N orthogonal matrix Q.

in
lda

The leading dimension of the array A. lda >= max(1, n).

in
tau

Double precision array, dimension (n-1). TAU(i) must contain the scalar factor of the elementary reflector H(i), as returned by DSYTRD.

out
work

Double precision array, dimension (max(1, lwork)). On exit, if info = 0, work[0] returns the optimal lwork.

in
lwork

The dimension of the array work. lwork >= max(1, n-1). For optimum performance lwork >= (n-1)*nb, where nb is the optimal blocksize. If lwork = -1, then a workspace query is assumed.

out
info

  • = 0: successful exit

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