56 SUBROUTINE cerrqrt( PATH, NUNIT )
73 parameter ( nmax = 2 )
79 COMPLEX A( nmax, nmax ), T( nmax, nmax ), W( nmax ),
92 COMMON / infoc / infot, nout, ok, lerr
93 COMMON / srnamc / srnamt
96 INTRINSIC float, cmplx
101 WRITE( nout, fmt = * )
107 a( i, j ) = 1.0 / cmplx( float(i+j), 0.0 )
108 c( i, j ) = 1.0 / cmplx( float(i+j), 0.0 )
109 t( i, j ) = 1.0 / cmplx( float(i+j), 0.0 )
121 CALL cgeqrt( -1, 0, 1, a, 1, t, 1, w, info )
122 CALL chkxer(
'CGEQRT', infot, nout, lerr, ok )
124 CALL cgeqrt( 0, -1, 1, a, 1, t, 1, w, info )
125 CALL chkxer(
'CGEQRT', infot, nout, lerr, ok )
127 CALL cgeqrt( 0, 0, 0, a, 1, t, 1, w, info )
128 CALL chkxer(
'CGEQRT', infot, nout, lerr, ok )
130 CALL cgeqrt( 2, 1, 1, a, 1, t, 1, w, info )
131 CALL chkxer(
'CGEQRT', infot, nout, lerr, ok )
133 CALL cgeqrt( 2, 2, 2, a, 2, t, 1, w, info )
134 CALL chkxer(
'CGEQRT', infot, nout, lerr, ok )
140 CALL cgeqrt2( -1, 0, a, 1, t, 1, info )
141 CALL chkxer(
'CGEQRT2', infot, nout, lerr, ok )
143 CALL cgeqrt2( 0, -1, a, 1, t, 1, info )
144 CALL chkxer(
'CGEQRT2', infot, nout, lerr, ok )
146 CALL cgeqrt2( 2, 1, a, 1, t, 1, info )
147 CALL chkxer(
'CGEQRT2', infot, nout, lerr, ok )
149 CALL cgeqrt2( 2, 2, a, 2, t, 1, info )
150 CALL chkxer(
'CGEQRT2', infot, nout, lerr, ok )
156 CALL cgeqrt3( -1, 0, a, 1, t, 1, info )
157 CALL chkxer(
'CGEQRT3', infot, nout, lerr, ok )
159 CALL cgeqrt3( 0, -1, a, 1, t, 1, info )
160 CALL chkxer(
'CGEQRT3', infot, nout, lerr, ok )
162 CALL cgeqrt3( 2, 1, a, 1, t, 1, info )
163 CALL chkxer(
'CGEQRT3', infot, nout, lerr, ok )
165 CALL cgeqrt3( 2, 2, a, 2, t, 1, info )
166 CALL chkxer(
'CGEQRT3', infot, nout, lerr, ok )
172 CALL cgemqrt(
'/',
'N', 0, 0, 0, 1, a, 1, t, 1, c, 1, w, info )
173 CALL chkxer(
'CGEMQRT', infot, nout, lerr, ok )
175 CALL cgemqrt(
'L',
'/', 0, 0, 0, 1, a, 1, t, 1, c, 1, w, info )
176 CALL chkxer(
'CGEMQRT', infot, nout, lerr, ok )
178 CALL cgemqrt(
'L',
'N', -1, 0, 0, 1, a, 1, t, 1, c, 1, w, info )
179 CALL chkxer(
'CGEMQRT', infot, nout, lerr, ok )
181 CALL cgemqrt(
'L',
'N', 0, -1, 0, 1, a, 1, t, 1, c, 1, w, info )
182 CALL chkxer(
'CGEMQRT', infot, nout, lerr, ok )
184 CALL cgemqrt(
'L',
'N', 0, 0, -1, 1, a, 1, t, 1, c, 1, w, info )
185 CALL chkxer(
'CGEMQRT', infot, nout, lerr, ok )
187 CALL cgemqrt(
'R',
'N', 0, 0, -1, 1, a, 1, t, 1, c, 1, w, info )
188 CALL chkxer(
'CGEMQRT', infot, nout, lerr, ok )
190 CALL cgemqrt(
'L',
'N', 0, 0, 0, 0, a, 1, t, 1, c, 1, w, info )
191 CALL chkxer(
'CGEMQRT', infot, nout, lerr, ok )
193 CALL cgemqrt(
'R',
'N', 1, 2, 1, 1, a, 1, t, 1, c, 1, w, info )
194 CALL chkxer(
'CGEMQRT', infot, nout, lerr, ok )
196 CALL cgemqrt(
'L',
'N', 2, 1, 1, 1, a, 1, t, 1, c, 1, w, info )
197 CALL chkxer(
'CGEMQRT', infot, nout, lerr, ok )
199 CALL cgemqrt(
'R',
'N', 1, 1, 1, 1, a, 1, t, 0, c, 1, w, info )
200 CALL chkxer(
'CGEMQRT', infot, nout, lerr, ok )
202 CALL cgemqrt(
'L',
'N', 1, 1, 1, 1, a, 1, t, 1, c, 0, w, info )
203 CALL chkxer(
'CGEMQRT', infot, nout, lerr, ok )
207 CALL alaesm( path, ok, nout )
subroutine alaesm(PATH, OK, NOUT)
ALAESM
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
subroutine cgemqrt(SIDE, TRANS, M, N, K, NB, V, LDV, T, LDT, C, LDC, WORK, INFO)
CGEMQRT
subroutine cgeqrt2(M, N, A, LDA, T, LDT, INFO)
CGEQRT2 computes a QR factorization of a general real or complex matrix using the compact WY represen...
subroutine cerrqrt(PATH, NUNIT)
CERRQRT
recursive subroutine cgeqrt3(M, N, A, LDA, T, LDT, INFO)
CGEQRT3 recursively computes a QR factorization of a general real or complex matrix using the compac...
subroutine cgeqrt(M, N, NB, A, LDA, T, LDT, WORK, INFO)
CGEQRT