76 DOUBLE PRECISION A( NMAX, NMAX ), T( NMAX, NMAX ), W( NMAX ),
89 COMMON / infoc / infot, nout, ok, lerr
90 COMMON / srnamc / srnamt
98 WRITE( nout, fmt = * )
104 a( i, j ) = 1.d0 / dble( i+j )
105 c( i, j ) = 1.d0 / dble( i+j )
106 t( i, j ) = 1.d0 / dble( i+j )
118 CALL dgeqrt( -1, 0, 1, a, 1, t, 1, w, info )
119 CALL chkxer(
'DGEQRT', infot, nout, lerr, ok )
121 CALL dgeqrt( 0, -1, 1, a, 1, t, 1, w, info )
122 CALL chkxer(
'DGEQRT', infot, nout, lerr, ok )
124 CALL dgeqrt( 0, 0, 0, a, 1, t, 1, w, info )
125 CALL chkxer(
'DGEQRT', infot, nout, lerr, ok )
127 CALL dgeqrt( 2, 1, 1, a, 1, t, 1, w, info )
128 CALL chkxer(
'DGEQRT', infot, nout, lerr, ok )
130 CALL dgeqrt( 2, 2, 2, a, 2, t, 1, w, info )
131 CALL chkxer(
'DGEQRT', infot, nout, lerr, ok )
137 CALL dgeqrt2( -1, 0, a, 1, t, 1, info )
138 CALL chkxer(
'DGEQRT2', infot, nout, lerr, ok )
140 CALL dgeqrt2( 0, -1, a, 1, t, 1, info )
141 CALL chkxer(
'DGEQRT2', infot, nout, lerr, ok )
143 CALL dgeqrt2( 2, 1, a, 1, t, 1, info )
144 CALL chkxer(
'DGEQRT2', infot, nout, lerr, ok )
146 CALL dgeqrt2( 2, 2, a, 2, t, 1, info )
147 CALL chkxer(
'DGEQRT2', infot, nout, lerr, ok )
153 CALL dgeqrt3( -1, 0, a, 1, t, 1, info )
154 CALL chkxer(
'DGEQRT3', infot, nout, lerr, ok )
156 CALL dgeqrt3( 0, -1, a, 1, t, 1, info )
157 CALL chkxer(
'DGEQRT3', infot, nout, lerr, ok )
159 CALL dgeqrt3( 2, 1, a, 1, t, 1, info )
160 CALL chkxer(
'DGEQRT3', infot, nout, lerr, ok )
162 CALL dgeqrt3( 2, 2, a, 2, t, 1, info )
163 CALL chkxer(
'DGEQRT3', infot, nout, lerr, ok )
169 CALL dgemqrt(
'/',
'N', 0, 0, 0, 1, a, 1, t, 1, c, 1, w, info )
170 CALL chkxer(
'DGEMQRT', infot, nout, lerr, ok )
172 CALL dgemqrt(
'L',
'/', 0, 0, 0, 1, a, 1, t, 1, c, 1, w, info )
173 CALL chkxer(
'DGEMQRT', infot, nout, lerr, ok )
175 CALL dgemqrt(
'L',
'N', -1, 0, 0, 1, a, 1, t, 1, c, 1, w, info )
176 CALL chkxer(
'DGEMQRT', infot, nout, lerr, ok )
178 CALL dgemqrt(
'L',
'N', 0, -1, 0, 1, a, 1, t, 1, c, 1, w, info )
179 CALL chkxer(
'DGEMQRT', infot, nout, lerr, ok )
181 CALL dgemqrt(
'L',
'N', 0, 0, -1, 1, a, 1, t, 1, c, 1, w, info )
182 CALL chkxer(
'DGEMQRT', infot, nout, lerr, ok )
184 CALL dgemqrt(
'R',
'N', 0, 0, -1, 1, a, 1, t, 1, c, 1, w, info )
185 CALL chkxer(
'DGEMQRT', infot, nout, lerr, ok )
187 CALL dgemqrt(
'L',
'N', 0, 0, 0, 0, a, 1, t, 1, c, 1, w, info )
188 CALL chkxer(
'DGEMQRT', infot, nout, lerr, ok )
190 CALL dgemqrt(
'R',
'N', 1, 2, 1, 1, a, 1, t, 1, c, 1, w, info )
191 CALL chkxer(
'DGEMQRT', infot, nout, lerr, ok )
193 CALL dgemqrt(
'L',
'N', 2, 1, 1, 1, a, 1, t, 1, c, 1, w, info )
194 CALL chkxer(
'DGEMQRT', infot, nout, lerr, ok )
196 CALL dgemqrt(
'R',
'N', 1, 1, 1, 1, a, 1, t, 0, c, 1, w, info )
197 CALL chkxer(
'DGEMQRT', infot, nout, lerr, ok )
199 CALL dgemqrt(
'L',
'N', 1, 1, 1, 1, a, 1, t, 1, c, 0, w, info )
200 CALL chkxer(
'DGEMQRT', infot, nout, lerr, ok )
204 CALL alaesm( path, ok, nout )
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
subroutine alaesm(PATH, OK, NOUT)
ALAESM
subroutine derrqrt(PATH, NUNIT)
DERRQRT
recursive subroutine dgeqrt3(M, N, A, LDA, T, LDT, INFO)
DGEQRT3 recursively computes a QR factorization of a general real or complex matrix using the compact...
subroutine dgeqrt(M, N, NB, A, LDA, T, LDT, WORK, INFO)
DGEQRT
subroutine dgeqrt2(M, N, A, LDA, T, LDT, INFO)
DGEQRT2 computes a QR factorization of a general real or complex matrix using the compact WY represen...
subroutine dgemqrt(SIDE, TRANS, M, N, K, NB, V, LDV, T, LDT, C, LDC, WORK, INFO)
DGEMQRT