56 SUBROUTINE derrqrt( PATH, NUNIT )
73 parameter ( nmax = 2 )
79 DOUBLE PRECISION A( nmax, nmax ), T( nmax, nmax ), W( nmax ),
92 COMMON / infoc / infot, nout, ok, lerr
93 COMMON / srnamc / srnamt
101 WRITE( nout, fmt = * )
107 a( i, j ) = 1.d0 / dble( i+j )
108 c( i, j ) = 1.d0 / dble( i+j )
109 t( i, j ) = 1.d0 / dble( i+j )
121 CALL dgeqrt( -1, 0, 1, a, 1, t, 1, w, info )
122 CALL chkxer(
'DGEQRT', infot, nout, lerr, ok )
124 CALL dgeqrt( 0, -1, 1, a, 1, t, 1, w, info )
125 CALL chkxer(
'DGEQRT', infot, nout, lerr, ok )
127 CALL dgeqrt( 0, 0, 0, a, 1, t, 1, w, info )
128 CALL chkxer(
'DGEQRT', infot, nout, lerr, ok )
130 CALL dgeqrt( 2, 1, 1, a, 1, t, 1, w, info )
131 CALL chkxer(
'DGEQRT', infot, nout, lerr, ok )
133 CALL dgeqrt( 2, 2, 2, a, 2, t, 1, w, info )
134 CALL chkxer(
'DGEQRT', infot, nout, lerr, ok )
140 CALL dgeqrt2( -1, 0, a, 1, t, 1, info )
141 CALL chkxer(
'DGEQRT2', infot, nout, lerr, ok )
143 CALL dgeqrt2( 0, -1, a, 1, t, 1, info )
144 CALL chkxer(
'DGEQRT2', infot, nout, lerr, ok )
146 CALL dgeqrt2( 2, 1, a, 1, t, 1, info )
147 CALL chkxer(
'DGEQRT2', infot, nout, lerr, ok )
149 CALL dgeqrt2( 2, 2, a, 2, t, 1, info )
150 CALL chkxer(
'DGEQRT2', infot, nout, lerr, ok )
156 CALL dgeqrt3( -1, 0, a, 1, t, 1, info )
157 CALL chkxer(
'DGEQRT3', infot, nout, lerr, ok )
159 CALL dgeqrt3( 0, -1, a, 1, t, 1, info )
160 CALL chkxer(
'DGEQRT3', infot, nout, lerr, ok )
162 CALL dgeqrt3( 2, 1, a, 1, t, 1, info )
163 CALL chkxer(
'DGEQRT3', infot, nout, lerr, ok )
165 CALL dgeqrt3( 2, 2, a, 2, t, 1, info )
166 CALL chkxer(
'DGEQRT3', infot, nout, lerr, ok )
172 CALL dgemqrt(
'/',
'N', 0, 0, 0, 1, a, 1, t, 1, c, 1, w, info )
173 CALL chkxer(
'DGEMQRT', infot, nout, lerr, ok )
175 CALL dgemqrt(
'L',
'/', 0, 0, 0, 1, a, 1, t, 1, c, 1, w, info )
176 CALL chkxer(
'DGEMQRT', infot, nout, lerr, ok )
178 CALL dgemqrt(
'L',
'N', -1, 0, 0, 1, a, 1, t, 1, c, 1, w, info )
179 CALL chkxer(
'DGEMQRT', infot, nout, lerr, ok )
181 CALL dgemqrt(
'L',
'N', 0, -1, 0, 1, a, 1, t, 1, c, 1, w, info )
182 CALL chkxer(
'DGEMQRT', infot, nout, lerr, ok )
184 CALL dgemqrt(
'L',
'N', 0, 0, -1, 1, a, 1, t, 1, c, 1, w, info )
185 CALL chkxer(
'DGEMQRT', infot, nout, lerr, ok )
187 CALL dgemqrt(
'R',
'N', 0, 0, -1, 1, a, 1, t, 1, c, 1, w, info )
188 CALL chkxer(
'DGEMQRT', infot, nout, lerr, ok )
190 CALL dgemqrt(
'L',
'N', 0, 0, 0, 0, a, 1, t, 1, c, 1, w, info )
191 CALL chkxer(
'DGEMQRT', infot, nout, lerr, ok )
193 CALL dgemqrt(
'R',
'N', 1, 2, 1, 1, a, 1, t, 1, c, 1, w, info )
194 CALL chkxer(
'DGEMQRT', infot, nout, lerr, ok )
196 CALL dgemqrt(
'L',
'N', 2, 1, 1, 1, a, 1, t, 1, c, 1, w, info )
197 CALL chkxer(
'DGEMQRT', infot, nout, lerr, ok )
199 CALL dgemqrt(
'R',
'N', 1, 1, 1, 1, a, 1, t, 0, c, 1, w, info )
200 CALL chkxer(
'DGEMQRT', infot, nout, lerr, ok )
202 CALL dgemqrt(
'L',
'N', 1, 1, 1, 1, a, 1, t, 1, c, 0, w, info )
203 CALL chkxer(
'DGEMQRT', infot, nout, lerr, ok )
207 CALL alaesm( path, ok, nout )
subroutine dgemqrt(SIDE, TRANS, M, N, K, NB, V, LDV, T, LDT, C, LDC, WORK, INFO)
DGEMQRT
subroutine derrqrt(PATH, NUNIT)
DERRQRT
subroutine alaesm(PATH, OK, NOUT)
ALAESM
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 chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
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...