sopgtr - generate a real orthogonal matrix Q which is defined as the product of n-1 elementary reflectors H(i)
of order n, as returned by SSPTRD using packed storage
SUBROUTINE SOPGTR( UPLO, N, AP, TAU, Q, LDQ, WORK, INFO) CHARACTER * 1 UPLO INTEGER N, LDQ, INFO REAL AP(*), TAU(*), Q(LDQ,*), WORK(*)
SUBROUTINE SOPGTR_64( UPLO, N, AP, TAU, Q, LDQ, WORK, INFO) CHARACTER * 1 UPLO INTEGER*8 N, LDQ, INFO REAL AP(*), TAU(*), Q(LDQ,*), WORK(*)
SUBROUTINE OPGTR( UPLO, [N], AP, TAU, Q, [LDQ], [WORK], [INFO]) CHARACTER(LEN=1) :: UPLO INTEGER :: N, LDQ, INFO REAL, DIMENSION(:) :: AP, TAU, WORK REAL, DIMENSION(:,:) :: Q
SUBROUTINE OPGTR_64( UPLO, [N], AP, TAU, Q, [LDQ], [WORK], [INFO]) CHARACTER(LEN=1) :: UPLO INTEGER(8) :: N, LDQ, INFO REAL, DIMENSION(:) :: AP, TAU, WORK REAL, DIMENSION(:,:) :: Q
#include <sunperf.h>
void sopgtr(char uplo, int n, float *ap, float *tau, float *q, int ldq, int *info);
void sopgtr_64(char uplo, long n, float *ap, float *tau, float *q, long ldq, long *info);
sopgtr generates a real orthogonal matrix Q which is defined as the
product of n-1 elementary reflectors H(i)
of order n, as returned by
SSPTRD using packed storage:
if UPLO = 'U', Q = H(n-1)
. . . H(2)
H(1),
if UPLO = 'L', Q = H(1)
H(2)
. . . H(n-1).
= 'U': Upper triangular packed storage used in previous call to SSPTRD; = 'L': Lower triangular packed storage used in previous call to SSPTRD.
TAU(i)
must contain the scalar factor of the elementary
reflector H(i), as returned by SSPTRD.
dimension(N-1)
= 0: successful exit
< 0: if INFO = -i, the i-th argument had an illegal value