69 DOUBLE PRECISION alpha, beta
73 COMPLEX*16 a( 1, 1),
b( 1, 1)
89 COMMON / infoc / infot, nout, ok, lerr
90 COMMON / srnamc / srnamt
96 a( 1, 1 ) = dcmplx( 1.0d0 , 1.0d0 )
97 b( 1, 1 ) = dcmplx( 1.0d0 , 1.0d0 )
99 calpha = dcmplx( 1.0d0 , 1.0d0 )
104 CALL
zpftrf(
'/',
'U', 0, a, info )
105 CALL
chkxer(
'ZPFTRF', infot, nout, lerr, ok )
107 CALL
zpftrf(
'N',
'/', 0, a, info )
108 CALL
chkxer(
'ZPFTRF', infot, nout, lerr, ok )
110 CALL
zpftrf(
'N',
'U', -1, a, info )
111 CALL
chkxer(
'ZPFTRF', infot, nout, lerr, ok )
115 CALL
zpftrs(
'/',
'U', 0, 0, a,
b, 1, info )
116 CALL
chkxer(
'ZPFTRS', infot, nout, lerr, ok )
118 CALL
zpftrs(
'N',
'/', 0, 0, a,
b, 1, info )
119 CALL
chkxer(
'ZPFTRS', infot, nout, lerr, ok )
121 CALL
zpftrs(
'N',
'U', -1, 0, a,
b, 1, info )
122 CALL
chkxer(
'ZPFTRS', infot, nout, lerr, ok )
124 CALL
zpftrs(
'N',
'U', 0, -1, a,
b, 1, info )
125 CALL
chkxer(
'ZPFTRS', infot, nout, lerr, ok )
127 CALL
zpftrs(
'N',
'U', 0, 0, a,
b, 0, info )
128 CALL
chkxer(
'ZPFTRS', infot, nout, lerr, ok )
132 CALL
zpftri(
'/',
'U', 0, a, info )
133 CALL
chkxer(
'ZPFTRI', infot, nout, lerr, ok )
135 CALL
zpftri(
'N',
'/', 0, a, info )
136 CALL
chkxer(
'ZPFTRI', infot, nout, lerr, ok )
138 CALL
zpftri(
'N',
'U', -1, a, info )
139 CALL
chkxer(
'ZPFTRI', infot, nout, lerr, ok )
143 CALL
ztfsm(
'/',
'L',
'U',
'C',
'U', 0, 0, calpha, a,
b, 1 )
144 CALL
chkxer(
'ZTFSM ', infot, nout, lerr, ok )
146 CALL
ztfsm(
'N',
'/',
'U',
'C',
'U', 0, 0, calpha, a,
b, 1 )
147 CALL
chkxer(
'ZTFSM ', infot, nout, lerr, ok )
149 CALL
ztfsm(
'N',
'L',
'/',
'C',
'U', 0, 0, calpha, a,
b, 1 )
150 CALL
chkxer(
'ZTFSM ', infot, nout, lerr, ok )
152 CALL
ztfsm(
'N',
'L',
'U',
'/',
'U', 0, 0, calpha, a,
b, 1 )
153 CALL
chkxer(
'ZTFSM ', infot, nout, lerr, ok )
155 CALL
ztfsm(
'N',
'L',
'U',
'C',
'/', 0, 0, calpha, a,
b, 1 )
156 CALL
chkxer(
'ZTFSM ', infot, nout, lerr, ok )
158 CALL
ztfsm(
'N',
'L',
'U',
'C',
'U', -1, 0, calpha, a,
b, 1 )
159 CALL
chkxer(
'ZTFSM ', infot, nout, lerr, ok )
161 CALL
ztfsm(
'N',
'L',
'U',
'C',
'U', 0, -1, calpha, a,
b, 1 )
162 CALL
chkxer(
'ZTFSM ', infot, nout, lerr, ok )
164 CALL
ztfsm(
'N',
'L',
'U',
'C',
'U', 0, 0, calpha, a,
b, 0 )
165 CALL
chkxer(
'ZTFSM ', infot, nout, lerr, ok )
169 CALL
ztftri(
'/',
'L',
'N', 0, a, info )
170 CALL
chkxer(
'ZTFTRI', infot, nout, lerr, ok )
172 CALL
ztftri(
'N',
'/',
'N', 0, a, info )
173 CALL
chkxer(
'ZTFTRI', infot, nout, lerr, ok )
175 CALL
ztftri(
'N',
'L',
'/', 0, a, info )
176 CALL
chkxer(
'ZTFTRI', infot, nout, lerr, ok )
178 CALL
ztftri(
'N',
'L',
'N', -1, a, info )
179 CALL
chkxer(
'ZTFTRI', infot, nout, lerr, ok )
183 CALL
ztfttr(
'/',
'U', 0, a,
b, 1, info )
184 CALL
chkxer(
'ZTFTTR', infot, nout, lerr, ok )
186 CALL
ztfttr(
'N',
'/', 0, a,
b, 1, info )
187 CALL
chkxer(
'ZTFTTR', infot, nout, lerr, ok )
189 CALL
ztfttr(
'N',
'U', -1, a,
b, 1, info )
190 CALL
chkxer(
'ZTFTTR', infot, nout, lerr, ok )
192 CALL
ztfttr(
'N',
'U', 0, a,
b, 0, info )
193 CALL
chkxer(
'ZTFTTR', infot, nout, lerr, ok )
197 CALL
ztrttf(
'/',
'U', 0, a, 1,
b, info )
198 CALL
chkxer(
'ZTRTTF', infot, nout, lerr, ok )
200 CALL
ztrttf(
'N',
'/', 0, a, 1,
b, info )
201 CALL
chkxer(
'ZTRTTF', infot, nout, lerr, ok )
203 CALL
ztrttf(
'N',
'U', -1, a, 1,
b, info )
204 CALL
chkxer(
'ZTRTTF', infot, nout, lerr, ok )
206 CALL
ztrttf(
'N',
'U', 0, a, 0,
b, info )
207 CALL
chkxer(
'ZTRTTF', infot, nout, lerr, ok )
211 CALL
ztfttp(
'/',
'U', 0, a,
b, info )
212 CALL
chkxer(
'ZTFTTP', infot, nout, lerr, ok )
214 CALL
ztfttp(
'N',
'/', 0, a,
b, info )
215 CALL
chkxer(
'ZTFTTP', infot, nout, lerr, ok )
217 CALL
ztfttp(
'N',
'U', -1, a,
b, info )
218 CALL
chkxer(
'ZTFTTP', infot, nout, lerr, ok )
222 CALL
ztpttf(
'/',
'U', 0, a,
b, info )
223 CALL
chkxer(
'ZTPTTF', infot, nout, lerr, ok )
225 CALL
ztpttf(
'N',
'/', 0, a,
b, info )
226 CALL
chkxer(
'ZTPTTF', infot, nout, lerr, ok )
228 CALL
ztpttf(
'N',
'U', -1, a,
b, info )
229 CALL
chkxer(
'ZTPTTF', infot, nout, lerr, ok )
233 CALL
ztrttp(
'/', 0, a, 1,
b, info )
234 CALL
chkxer(
'ZTRTTP', infot, nout, lerr, ok )
236 CALL
ztrttp(
'U', -1, a, 1,
b, info )
237 CALL
chkxer(
'ZTRTTP', infot, nout, lerr, ok )
239 CALL
ztrttp(
'U', 0, a, 0,
b, info )
240 CALL
chkxer(
'ZTRTTP', infot, nout, lerr, ok )
244 CALL
ztpttr(
'/', 0, a,
b, 1, info )
245 CALL
chkxer(
'ZTPTTR', infot, nout, lerr, ok )
247 CALL
ztpttr(
'U', -1, a,
b, 1, info )
248 CALL
chkxer(
'ZTPTTR', infot, nout, lerr, ok )
250 CALL
ztpttr(
'U', 0, a,
b, 0, info )
251 CALL
chkxer(
'ZTPTTR', infot, nout, lerr, ok )
255 CALL
zhfrk(
'/',
'U',
'N', 0, 0, alpha, a, 1, beta,
b )
256 CALL
chkxer(
'ZHFRK ', infot, nout, lerr, ok )
258 CALL
zhfrk(
'N',
'/',
'N', 0, 0, alpha, a, 1, beta,
b )
259 CALL
chkxer(
'ZHFRK ', infot, nout, lerr, ok )
261 CALL
zhfrk(
'N',
'U',
'/', 0, 0, alpha, a, 1, beta,
b )
262 CALL
chkxer(
'ZHFRK ', infot, nout, lerr, ok )
264 CALL
zhfrk(
'N',
'U',
'N', -1, 0, alpha, a, 1, beta,
b )
265 CALL
chkxer(
'ZHFRK ', infot, nout, lerr, ok )
267 CALL
zhfrk(
'N',
'U',
'N', 0, -1, alpha, a, 1, beta,
b )
268 CALL
chkxer(
'ZHFRK ', infot, nout, lerr, ok )
270 CALL
zhfrk(
'N',
'U',
'N', 0, 0, alpha, a, 0, beta,
b )
271 CALL
chkxer(
'ZHFRK ', infot, nout, lerr, ok )
276 WRITE( nout, fmt = 9999 )
278 WRITE( nout, fmt = 9998 )
281 9999
FORMAT( 1
x,
'COMPLEX*16 RFP routines passed the tests of the ',
283 9998
FORMAT(
' *** RFP routines failed the tests of the error ',
subroutine ztfttp(TRANSR, UPLO, N, ARF, AP, INFO)
ZTFTTP copies a triangular matrix from the rectangular full packed format (TF) to the standard packed...
subroutine ztftri(TRANSR, UPLO, DIAG, N, A, INFO)
ZTFTRI
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
subroutine ztfsm(TRANSR, SIDE, UPLO, TRANS, DIAG, M, N, ALPHA, A, B, LDB)
ZTFSM solves a matrix equation (one operand is a triangular matrix in RFP format).
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real b(3) integer i
subroutine zerrrfp(NUNIT)
ZERRRFP
subroutine ztpttf(TRANSR, UPLO, N, AP, ARF, INFO)
ZTPTTF copies a triangular matrix from the standard packed format (TP) to the rectangular full packed...
subroutine zpftrs(TRANSR, UPLO, N, NRHS, A, B, LDB, INFO)
ZPFTRS
subroutine zhfrk(TRANSR, UPLO, TRANS, N, K, ALPHA, A, LDA, BETA, C)
ZHFRK performs a Hermitian rank-k operation for matrix in RFP format.
subroutine ztrttf(TRANSR, UPLO, N, A, LDA, ARF, INFO)
ZTRTTF copies a triangular matrix from the standard full format (TR) to the rectangular full packed f...
subroutine zpftrf(TRANSR, UPLO, N, A, INFO)
ZPFTRF
subroutine ztpttr(UPLO, N, AP, A, LDA, INFO)
ZTPTTR copies a triangular matrix from the standard packed format (TP) to the standard full format (T...
subroutine ztrttp(UPLO, N, A, LDA, AP, INFO)
ZTRTTP copies a triangular matrix from the standard full format (TR) to the standard packed format (T...
subroutine zpftri(TRANSR, UPLO, N, A, INFO)
ZPFTRI
subroutine ztfttr(TRANSR, UPLO, N, ARF, A, LDA, INFO)
ZTFTTR copies a triangular matrix from the rectangular full packed format (TF) to the standard full f...