SUBROUTINE DOPGTR( UPLO, N, AP, TAU, Q, LDQ, WORK, INFO) CHARACTER * 1 UPLO INTEGER N, LDQ, INFO DOUBLE PRECISION AP(*), TAU(*), Q(LDQ,*), WORK(*) SUBROUTINE DOPGTR_64( UPLO, N, AP, TAU, Q, LDQ, WORK, INFO) CHARACTER * 1 UPLO INTEGER*8 N, LDQ, INFO DOUBLE PRECISION AP(*), TAU(*), Q(LDQ,*), WORK(*)
SUBROUTINE OPGTR( UPLO, [N], AP, TAU, Q, [LDQ], [WORK], [INFO]) CHARACTER(LEN=1) :: UPLO INTEGER :: N, LDQ, INFO REAL(8), DIMENSION(:) :: AP, TAU, WORK REAL(8), DIMENSION(:,:) :: Q SUBROUTINE OPGTR_64( UPLO, [N], AP, TAU, Q, [LDQ], [WORK], [INFO]) CHARACTER(LEN=1) :: UPLO INTEGER(8) :: N, LDQ, INFO REAL(8), DIMENSION(:) :: AP, TAU, WORK REAL(8), DIMENSION(:,:) :: Q
void dopgtr(char uplo, int n, double *ap, double *tau, double *q, int ldq, int *info);
void dopgtr_64(char uplo, long n, double *ap, double *tau, double *q, long ldq, long *info);
if UPLO = 'U', Q = H(n-1) . . . H(2) H(1),
if UPLO = 'L', Q = H(1) H(2) . . . H(n-1).