73 parameter ( nmax = 2 )
79 REAL 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.0 / float( i+j )
108 c( i, j ) = 1.0 / float( i+j )
109 t( i, j ) = 1.0 / float( i+j )
121 CALL sgeqrt( -1, 0, 1, a, 1, t, 1, w, info )
122 CALL chkxer(
'SGEQRT', infot, nout, lerr, ok )
124 CALL sgeqrt( 0, -1, 1, a, 1, t, 1, w, info )
125 CALL chkxer(
'SGEQRT', infot, nout, lerr, ok )
127 CALL sgeqrt( 0, 0, 0, a, 1, t, 1, w, info )
128 CALL chkxer(
'SGEQRT', infot, nout, lerr, ok )
130 CALL sgeqrt( 2, 1, 1, a, 1, t, 1, w, info )
131 CALL chkxer(
'SGEQRT', infot, nout, lerr, ok )
133 CALL sgeqrt( 2, 2, 2, a, 2, t, 1, w, info )
134 CALL chkxer(
'SGEQRT', infot, nout, lerr, ok )
140 CALL sgeqrt2( -1, 0, a, 1, t, 1, info )
141 CALL chkxer(
'SGEQRT2', infot, nout, lerr, ok )
143 CALL sgeqrt2( 0, -1, a, 1, t, 1, info )
144 CALL chkxer(
'SGEQRT2', infot, nout, lerr, ok )
146 CALL sgeqrt2( 2, 1, a, 1, t, 1, info )
147 CALL chkxer(
'SGEQRT2', infot, nout, lerr, ok )
149 CALL sgeqrt2( 2, 2, a, 2, t, 1, info )
150 CALL chkxer(
'SGEQRT2', infot, nout, lerr, ok )
156 CALL sgeqrt3( -1, 0, a, 1, t, 1, info )
157 CALL chkxer(
'SGEQRT3', infot, nout, lerr, ok )
159 CALL sgeqrt3( 0, -1, a, 1, t, 1, info )
160 CALL chkxer(
'SGEQRT3', infot, nout, lerr, ok )
162 CALL sgeqrt3( 2, 1, a, 1, t, 1, info )
163 CALL chkxer(
'SGEQRT3', infot, nout, lerr, ok )
165 CALL sgeqrt3( 2, 2, a, 2, t, 1, info )
166 CALL chkxer(
'SGEQRT3', infot, nout, lerr, ok )
172 CALL sgemqrt(
'/',
'N', 0, 0, 0, 1, a, 1, t, 1, c, 1, w, info )
173 CALL chkxer(
'SGEMQRT', infot, nout, lerr, ok )
175 CALL sgemqrt(
'L',
'/', 0, 0, 0, 1, a, 1, t, 1, c, 1, w, info )
176 CALL chkxer(
'SGEMQRT', infot, nout, lerr, ok )
178 CALL sgemqrt(
'L',
'N', -1, 0, 0, 1, a, 1, t, 1, c, 1, w, info )
179 CALL chkxer(
'SGEMQRT', infot, nout, lerr, ok )
181 CALL sgemqrt(
'L',
'N', 0, -1, 0, 1, a, 1, t, 1, c, 1, w, info )
182 CALL chkxer(
'SGEMQRT', infot, nout, lerr, ok )
184 CALL sgemqrt(
'L',
'N', 0, 0, -1, 1, a, 1, t, 1, c, 1, w, info )
185 CALL chkxer(
'SGEMQRT', infot, nout, lerr, ok )
187 CALL sgemqrt(
'R',
'N', 0, 0, -1, 1, a, 1, t, 1, c, 1, w, info )
188 CALL chkxer(
'SGEMQRT', infot, nout, lerr, ok )
190 CALL sgemqrt(
'L',
'N', 0, 0, 0, 0, a, 1, t, 1, c, 1, w, info )
191 CALL chkxer(
'SGEMQRT', infot, nout, lerr, ok )
193 CALL sgemqrt(
'R',
'N', 1, 2, 1, 1, a, 1, t, 1, c, 1, w, info )
194 CALL chkxer(
'SGEMQRT', infot, nout, lerr, ok )
196 CALL sgemqrt(
'L',
'N', 2, 1, 1, 1, a, 1, t, 1, c, 1, w, info )
197 CALL chkxer(
'SGEMQRT', infot, nout, lerr, ok )
199 CALL sgemqrt(
'R',
'N', 1, 1, 1, 1, a, 1, t, 0, c, 1, w, info )
200 CALL chkxer(
'SGEMQRT', infot, nout, lerr, ok )
202 CALL sgemqrt(
'L',
'N', 1, 1, 1, 1, a, 1, t, 1, c, 0, w, info )
203 CALL chkxer(
'SGEMQRT', infot, nout, lerr, ok )
207 CALL alaesm( path, ok, nout )
subroutine sgeqrt2(M, N, A, LDA, T, LDT, INFO)
SGEQRT2 computes a QR factorization of a general real or complex matrix using the compact WY represen...
subroutine sgeqrt(M, N, NB, A, LDA, T, LDT, WORK, INFO)
SGEQRT
subroutine alaesm(PATH, OK, NOUT)
ALAESM
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
subroutine sgemqrt(SIDE, TRANS, M, N, K, NB, V, LDV, T, LDT, C, LDC, WORK, INFO)
SGEMQRT
recursive subroutine sgeqrt3(M, N, A, LDA, T, LDT, INFO)
SGEQRT3 recursively computes a QR factorization of a general real or complex matrix using the compact...