56 SUBROUTINE derrlqt( 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 dgelqt( -1, 0, 1, a, 1, t, 1, w, info )
122 CALL chkxer(
'DGELQT', infot, nout, lerr, ok )
124 CALL dgelqt( 0, -1, 1, a, 1, t, 1, w, info )
125 CALL chkxer(
'DGELQT', infot, nout, lerr, ok )
127 CALL dgelqt( 0, 0, 0, a, 1, t, 1, w, info )
128 CALL chkxer(
'DGELQT', infot, nout, lerr, ok )
130 CALL dgelqt( 2, 1, 1, a, 1, t, 1, w, info )
131 CALL chkxer(
'DGELQT', infot, nout, lerr, ok )
133 CALL dgelqt( 2, 2, 2, a, 2, t, 1, w, info )
134 CALL chkxer(
'DGELQT', infot, nout, lerr, ok )
140 CALL dgelqt3( -1, 0, a, 1, t, 1, info )
141 CALL chkxer(
'DGELQT3', infot, nout, lerr, ok )
143 CALL dgelqt3( 0, -1, a, 1, t, 1, info )
144 CALL chkxer(
'DGELQT3', infot, nout, lerr, ok )
146 CALL dgelqt3( 2, 2, a, 1, t, 1, info )
147 CALL chkxer(
'DGELQT3', infot, nout, lerr, ok )
149 CALL dgelqt3( 2, 2, a, 2, t, 1, info )
150 CALL chkxer(
'DGELQT3', infot, nout, lerr, ok )
156 CALL dgemlqt(
'/',
'N', 0, 0, 0, 1, a, 1, t, 1, c, 1, w, info )
157 CALL chkxer(
'DGEMLQT', infot, nout, lerr, ok )
159 CALL dgemlqt(
'L',
'/', 0, 0, 0, 1, a, 1, t, 1, c, 1, w, info )
160 CALL chkxer(
'DGEMLQT', infot, nout, lerr, ok )
162 CALL dgemlqt(
'L',
'N', -1, 0, 0, 1, a, 1, t, 1, c, 1, w, info )
163 CALL chkxer(
'DGEMLQT', infot, nout, lerr, ok )
165 CALL dgemlqt(
'L',
'N', 0, -1, 0, 1, a, 1, t, 1, c, 1, w, info )
166 CALL chkxer(
'DGEMLQT', infot, nout, lerr, ok )
168 CALL dgemlqt(
'L',
'N', 0, 0, -1, 1, a, 1, t, 1, c, 1, w, info )
169 CALL chkxer(
'DGEMLQT', infot, nout, lerr, ok )
171 CALL dgemlqt(
'R',
'N', 0, 0, -1, 1, a, 1, t, 1, c, 1, w, info )
172 CALL chkxer(
'DGEMLQT', infot, nout, lerr, ok )
174 CALL dgemlqt(
'L',
'N', 0, 0, 0, 0, a, 1, t, 1, c, 1, w, info )
175 CALL chkxer(
'DGEMLQT', infot, nout, lerr, ok )
177 CALL dgemlqt(
'R',
'N', 2, 2, 2, 1, a, 1, t, 1, c, 1, w, info )
178 CALL chkxer(
'DGEMLQT', infot, nout, lerr, ok )
180 CALL dgemlqt(
'L',
'N', 2, 2, 2, 1, a, 1, t, 1, c, 1, w, info )
181 CALL chkxer(
'DGEMLQT', infot, nout, lerr, ok )
183 CALL dgemlqt(
'R',
'N', 1, 1, 1, 1, a, 1, t, 0, c, 1, w, info )
184 CALL chkxer(
'DGEMLQT', infot, nout, lerr, ok )
186 CALL dgemlqt(
'L',
'N', 1, 1, 1, 1, a, 1, t, 1, c, 0, w, info )
187 CALL chkxer(
'DGEMLQT', infot, nout, lerr, ok )
191 CALL alaesm( path, ok, nout )
subroutine derrlqt(PATH, NUNIT)
DERLQT
subroutine alaesm(PATH, OK, NOUT)
ALAESM
subroutine dgemlqt(SIDE, TRANS, M, N, K, MB, V, LDV, T, LDT, C, LDC, WORK, INFO)
DGEMLQT
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
recursive subroutine dgelqt3(M, N, A, LDA, T, LDT, INFO)
DGELQT3 recursively computes a LQ factorization of a general real or complex matrix using the compact...
subroutine dgelqt(M, N, MB, A, LDA, T, LDT, WORK, INFO)
DGELQT