79 COMPLEX a( nmax, nmax ), t( nmax, nmax ), w( nmax ),
80 $
b( nmax, nmax ), c( nmax, 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)
111 w(
j ) = cmplx(0.0,0.0)
121 CALL
ctpqrt( -1, 1, 0, 1, a, 1,
b, 1, t, 1, w, info )
122 CALL
chkxer(
'CTPQRT', infot, nout, lerr, ok )
124 CALL
ctpqrt( 1, -1, 0, 1, a, 1,
b, 1, t, 1, w, info )
125 CALL
chkxer(
'CTPQRT', infot, nout, lerr, ok )
127 CALL
ctpqrt( 0, 1, -1, 1, a, 1,
b, 1, t, 1, w, info )
128 CALL
chkxer(
'CTPQRT', infot, nout, lerr, ok )
130 CALL
ctpqrt( 0, 1, 1, 1, a, 1,
b, 1, t, 1, w, info )
131 CALL
chkxer(
'CTPQRT', infot, nout, lerr, ok )
133 CALL
ctpqrt( 0, 1, 0, 0, a, 1,
b, 1, t, 1, w, info )
134 CALL
chkxer(
'CTPQRT', infot, nout, lerr, ok )
136 CALL
ctpqrt( 0, 1, 0, 2, a, 1,
b, 1, t, 1, w, info )
137 CALL
chkxer(
'CTPQRT', infot, nout, lerr, ok )
139 CALL
ctpqrt( 1, 2, 0, 2, a, 1,
b, 1, t, 1, w, info )
140 CALL
chkxer(
'CTPQRT', infot, nout, lerr, ok )
142 CALL
ctpqrt( 2, 1, 0, 1, a, 1,
b, 1, t, 1, w, info )
143 CALL
chkxer(
'CTPQRT', infot, nout, lerr, ok )
145 CALL
ctpqrt( 2, 2, 1, 2, a, 2,
b, 2, t, 1, w, info )
146 CALL
chkxer(
'CTPQRT', infot, nout, lerr, ok )
152 CALL
ctpqrt2( -1, 0, 0, a, 1,
b, 1, t, 1, info )
153 CALL
chkxer(
'CTPQRT2', infot, nout, lerr, ok )
155 CALL
ctpqrt2( 0, -1, 0, a, 1,
b, 1, t, 1, info )
156 CALL
chkxer(
'CTPQRT2', infot, nout, lerr, ok )
158 CALL
ctpqrt2( 0, 0, -1, a, 1,
b, 1, t, 1, info )
159 CALL
chkxer(
'CTPQRT2', infot, nout, lerr, ok )
161 CALL
ctpqrt2( 2, 2, 0, a, 1,
b, 2, t, 2, info )
162 CALL
chkxer(
'CTPQRT2', infot, nout, lerr, ok )
164 CALL
ctpqrt2( 2, 2, 0, a, 2,
b, 1, t, 2, info )
165 CALL
chkxer(
'CTPQRT2', infot, nout, lerr, ok )
167 CALL
ctpqrt2( 2, 2, 0, a, 2,
b, 2, t, 1, info )
168 CALL
chkxer(
'CTPQRT2', infot, nout, lerr, ok )
174 CALL
ctpmqrt(
'/',
'N', 0, 0, 0, 0, 1, a, 1, t, 1,
b, 1, c, 1,
176 CALL
chkxer(
'CTPMQRT', infot, nout, lerr, ok )
178 CALL
ctpmqrt(
'L',
'/', 0, 0, 0, 0, 1, a, 1, t, 1,
b, 1, c, 1,
180 CALL
chkxer(
'CTPMQRT', infot, nout, lerr, ok )
182 CALL
ctpmqrt(
'L',
'N', -1, 0, 0, 0, 1, a, 1, t, 1,
b, 1, c, 1,
184 CALL
chkxer(
'CTPMQRT', infot, nout, lerr, ok )
186 CALL
ctpmqrt(
'L',
'N', 0, -1, 0, 0, 1, a, 1, t, 1,
b, 1, c, 1,
188 CALL
chkxer(
'CTPMQRT', infot, nout, lerr, ok )
190 CALL
ctpmqrt(
'L',
'N', 0, 0, -1, 0, 1, a, 1, t, 1,
b, 1, c, 1,
193 CALL
ctpmqrt(
'L',
'N', 0, 0, 0, -1, 1, a, 1, t, 1,
b, 1, c, 1,
195 CALL
chkxer(
'CTPMQRT', infot, nout, lerr, ok )
197 CALL
ctpmqrt(
'L',
'N', 0, 0, 0, 0, 0, a, 1, t, 1,
b, 1, c, 1,
199 CALL
chkxer(
'CTPMQRT', infot, nout, lerr, ok )
201 CALL
ctpmqrt(
'R',
'N', 1, 2, 1, 1, 1, a, 1, t, 1,
b, 1, c, 1,
203 CALL
chkxer(
'CTPMQRT', infot, nout, lerr, ok )
205 CALL
ctpmqrt(
'L',
'N', 2, 1, 1, 1, 1, a, 1, t, 1,
b, 1, c, 1,
207 CALL
chkxer(
'CTPMQRT', infot, nout, lerr, ok )
209 CALL
ctpmqrt(
'R',
'N', 1, 1, 1, 1, 1, a, 1, t, 0,
b, 1, c, 1,
211 CALL
chkxer(
'CTPMQRT', infot, nout, lerr, ok )
213 CALL
ctpmqrt(
'L',
'N', 1, 1, 1, 1, 1, a, 1, t, 1,
b, 0, c, 1,
215 CALL
chkxer(
'CTPMQRT', infot, nout, lerr, ok )
217 CALL
ctpmqrt(
'L',
'N', 1, 1, 1, 1, 1, a, 1, t, 1,
b, 1, c, 0,
219 CALL
chkxer(
'CTPMQRT', infot, nout, lerr, ok )
223 CALL
alaesm( path, ok, nout )
subroutine cerrqrtp(PATH, NUNIT)
CERRQRTP
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real b(3) integer i
subroutine ctpqrt(M, N, L, NB, A, LDA, B, LDB, T, LDT, WORK, INFO)
CTPQRT
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
subroutine alaesm(PATH, OK, NOUT)
ALAESM
subroutine ctpqrt2(M, N, L, A, LDA, B, LDB, T, LDT, INFO)
CTPQRT2 computes a QR factorization of a real or complex "triangular-pentagonal" matrix, which is composed of a triangular block and a pentagonal block, using the compact WY representation for Q.
subroutine ctpmqrt(SIDE, TRANS, M, N, K, L, NB, V, LDV, T, LDT, A, LDA, B, LDB, WORK, INFO)
CTPMQRT