73 parameter ( nmax = 2 )
76 INTEGER i, info, j, nb
79 REAL a( nmax, nmax ), t( nmax, nmax ), w( nmax ),
80 $ c( nmax, nmax ), tau(nmax)
92 COMMON / infoc / infot, nout, ok, lerr
93 COMMON / srnamc / srnamt
101 WRITE( nout, fmt = * )
107 a( i, j ) = 1.d0 /
REAL( i+j )
108 c( i, j ) = 1.d0 /
REAL( i+j )
109 t( i, j ) = 1.d0 /
REAL( i+j )
121 CALL sgeqr( -1, 0, a, 1, tau, 1, w, 1, info )
122 CALL chkxer(
'SGEQR', infot, nout, lerr, ok )
124 CALL sgeqr( 0, -1, a, 1, tau, 1, w, 1, info )
125 CALL chkxer(
'SGEQR', infot, nout, lerr, ok )
127 CALL sgeqr( 1, 1, a, 0, tau, 1, w, 1, info )
128 CALL chkxer(
'SGEQR', infot, nout, lerr, ok )
130 CALL sgeqr( 3, 2, a, 3, tau, 1, w, 1, info )
131 CALL chkxer(
'SGEQR', infot, nout, lerr, ok )
133 CALL sgeqr( 3, 2, a, 3, tau, 7, w, 0, info )
134 CALL chkxer(
'SGEQR', infot, nout, lerr, ok )
143 CALL sgemqr(
'/',
'N', 0, 0, 0, a, 1, tau, 1, c, 1, w, 1,info)
144 CALL chkxer(
'SGEMQR', infot, nout, lerr, ok )
146 CALL sgemqr(
'L',
'/', 0, 0, 0, a, 1, tau, 1, c, 1, w, 1,info)
147 CALL chkxer(
'SGEMQR', infot, nout, lerr, ok )
149 CALL sgemqr(
'L',
'N', -1, 0, 0, a, 1, tau, 1, c, 1, w,1,info)
150 CALL chkxer(
'SGEMQR', infot, nout, lerr, ok )
152 CALL sgemqr(
'L',
'N', 0, -1, 0, a, 1, tau, 1, c, 1, w,1,info)
153 CALL chkxer(
'SGEMQR', infot, nout, lerr, ok )
155 CALL sgemqr(
'L',
'N', 0, 0, -1, a, 1, tau, 1, c, 1, w,1,info)
156 CALL chkxer(
'SGEMQR', infot, nout, lerr, ok )
158 CALL sgemqr(
'R',
'N', 0, 0, -1, a, 1, tau, 1, c, 1, w,1,info)
159 CALL chkxer(
'SGEMQR', infot, nout, lerr, ok )
161 CALL sgemqr(
'L',
'N', 2, 1, 0, a, 0, tau, 1, c, 1, w, 1,info)
162 CALL chkxer(
'SGEMQR', infot, nout, lerr, ok )
164 CALL sgemqr(
'R',
'N', 2, 2, 1, a, 2, tau, 0, c, 1, w, 1,info)
165 CALL chkxer(
'SGEMQR', infot, nout, lerr, ok )
167 CALL sgemqr(
'L',
'N', 2, 2, 1, a, 2, tau, 0, c, 1, w, 1,info)
168 CALL chkxer(
'SGEMQR', infot, nout, lerr, ok )
170 CALL sgemqr(
'L',
'N', 2, 1, 1, a, 2, tau, 6, c, 0, w, 1,info)
171 CALL chkxer(
'SGEMQR', infot, nout, lerr, ok )
173 CALL sgemqr(
'L',
'N', 2, 2, 1, a, 2, tau, 6, c, 2, w, 0,info)
174 CALL chkxer(
'SGEMQR', infot, nout, lerr, ok )
180 CALL sgelq( -1, 0, a, 1, tau, 1, w, 1, info )
181 CALL chkxer(
'SGELQ', infot, nout, lerr, ok )
183 CALL sgelq( 0, -1, a, 1, tau, 1, w, 1, info )
184 CALL chkxer(
'SGELQ', infot, nout, lerr, ok )
186 CALL sgelq( 1, 1, a, 0, tau, 1, w, 1, info )
187 CALL chkxer(
'SGELQ', infot, nout, lerr, ok )
189 CALL sgelq( 2, 3, a, 3, tau, 1, w, 1, info )
190 CALL chkxer(
'SGELQ', infot, nout, lerr, ok )
192 CALL sgelq( 2, 3, a, 3, tau, 7, w, 0, info )
193 CALL chkxer(
'SGELQ', infot, nout, lerr, ok )
202 CALL sgemlq(
'/',
'N', 0, 0, 0, a, 1, tau, 1, c, 1, w, 1,info)
203 CALL chkxer(
'SGEMLQ', infot, nout, lerr, ok )
205 CALL sgemlq(
'L',
'/', 0, 0, 0, a, 1, tau, 1, c, 1, w, 1,info)
206 CALL chkxer(
'SGEMLQ', infot, nout, lerr, ok )
208 CALL sgemlq(
'L',
'N', -1, 0, 0, a, 1, tau, 1, c, 1, w,1,info)
209 CALL chkxer(
'SGEMLQ', infot, nout, lerr, ok )
211 CALL sgemlq(
'L',
'N', 0, -1, 0, a, 1, tau, 1, c, 1, w,1,info)
212 CALL chkxer(
'SGEMLQ', infot, nout, lerr, ok )
214 CALL sgemlq(
'L',
'N', 0, 0, -1, a, 1, tau, 1, c, 1, w,1,info)
215 CALL chkxer(
'SGEMLQ', infot, nout, lerr, ok )
217 CALL sgemlq(
'R',
'N', 0, 0, -1, a, 1, tau, 1, c, 1, w,1,info)
218 CALL chkxer(
'SGEMLQ', infot, nout, lerr, ok )
220 CALL sgemlq(
'L',
'N', 1, 2, 0, a, 0, tau, 1, c, 1, w, 1,info)
221 CALL chkxer(
'SGEMLQ', infot, nout, lerr, ok )
223 CALL sgemlq(
'R',
'N', 2, 2, 1, a, 1, tau, 0, c, 1, w, 1,info)
224 CALL chkxer(
'SGEMLQ', infot, nout, lerr, ok )
226 CALL sgemlq(
'L',
'N', 2, 2, 1, a, 1, tau, 0, c, 1, w, 1,info)
227 CALL chkxer(
'SGEMLQ', infot, nout, lerr, ok )
229 CALL sgemlq(
'L',
'N', 1, 2, 1, a, 1, tau, 6, c, 0, w, 1,info)
230 CALL chkxer(
'SGEMLQ', infot, nout, lerr, ok )
232 CALL sgemlq(
'L',
'N', 2, 2, 1, a, 2, tau, 6, c, 2, w, 0,info)
233 CALL chkxer(
'SGEMLQ', infot, nout, lerr, ok )
237 CALL alaesm( path, ok, nout )
subroutine alaesm(PATH, OK, NOUT)
ALAESM
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
subroutine sgeqr(M, N, A, LDA, T, TSIZE, WORK, LWORK, INFO)
subroutine sgemqr(SIDE, TRANS, M, N, K, A, LDA, T, TSIZE, C, LDC, WORK, LWORK, INFO)
subroutine sgemlq(SIDE, TRANS, M, N, K, A, LDA, T, TSIZE, C, LDC, WORK, LWORK, INFO)
subroutine sgelq(M, N, A, LDA, T, TSIZE, WORK, LWORK, INFO)