66 DOUBLE PRECISION ALPHA, BETA
70 COMPLEX*16 A( 1, 1), B( 1, 1)
86 COMMON / infoc / infot, nout, ok, lerr
87 COMMON / srnamc / srnamt
93 a( 1, 1 ) = dcmplx( 1.0d0 , 1.0d0 )
94 b( 1, 1 ) = dcmplx( 1.0d0 , 1.0d0 )
96 calpha = dcmplx( 1.0d0 , 1.0d0 )
101 CALL zpftrf(
'/',
'U', 0, a, info )
102 CALL chkxer(
'ZPFTRF', infot, nout, lerr, ok )
104 CALL zpftrf(
'N',
'/', 0, a, info )
105 CALL chkxer(
'ZPFTRF', infot, nout, lerr, ok )
107 CALL zpftrf(
'N',
'U', -1, a, info )
108 CALL chkxer(
'ZPFTRF', infot, nout, lerr, ok )
112 CALL zpftrs(
'/',
'U', 0, 0, a, b, 1, info )
113 CALL chkxer(
'ZPFTRS', infot, nout, lerr, ok )
115 CALL zpftrs(
'N',
'/', 0, 0, a, b, 1, info )
116 CALL chkxer(
'ZPFTRS', infot, nout, lerr, ok )
118 CALL zpftrs(
'N',
'U', -1, 0, a, b, 1, info )
119 CALL chkxer(
'ZPFTRS', infot, nout, lerr, ok )
121 CALL zpftrs(
'N',
'U', 0, -1, a, b, 1, info )
122 CALL chkxer(
'ZPFTRS', infot, nout, lerr, ok )
124 CALL zpftrs(
'N',
'U', 0, 0, a, b, 0, info )
125 CALL chkxer(
'ZPFTRS', infot, nout, lerr, ok )
129 CALL zpftri(
'/',
'U', 0, a, info )
130 CALL chkxer(
'ZPFTRI', infot, nout, lerr, ok )
132 CALL zpftri(
'N',
'/', 0, a, info )
133 CALL chkxer(
'ZPFTRI', infot, nout, lerr, ok )
135 CALL zpftri(
'N',
'U', -1, a, info )
136 CALL chkxer(
'ZPFTRI', infot, nout, lerr, ok )
140 CALL ztfsm(
'/',
'L',
'U',
'C',
'U', 0, 0, calpha, a, b, 1 )
141 CALL chkxer(
'ZTFSM ', infot, nout, lerr, ok )
143 CALL ztfsm(
'N',
'/',
'U',
'C',
'U', 0, 0, calpha, a, b, 1 )
144 CALL chkxer(
'ZTFSM ', infot, nout, lerr, ok )
146 CALL ztfsm(
'N',
'L',
'/',
'C',
'U', 0, 0, calpha, a, b, 1 )
147 CALL chkxer(
'ZTFSM ', infot, nout, lerr, ok )
149 CALL ztfsm(
'N',
'L',
'U',
'/',
'U', 0, 0, calpha, a, b, 1 )
150 CALL chkxer(
'ZTFSM ', infot, nout, lerr, ok )
152 CALL ztfsm(
'N',
'L',
'U',
'C',
'/', 0, 0, calpha, a, b, 1 )
153 CALL chkxer(
'ZTFSM ', infot, nout, lerr, ok )
155 CALL ztfsm(
'N',
'L',
'U',
'C',
'U', -1, 0, calpha, a, b, 1 )
156 CALL chkxer(
'ZTFSM ', infot, nout, lerr, ok )
158 CALL ztfsm(
'N',
'L',
'U',
'C',
'U', 0, -1, calpha, a, b, 1 )
159 CALL chkxer(
'ZTFSM ', infot, nout, lerr, ok )
161 CALL ztfsm(
'N',
'L',
'U',
'C',
'U', 0, 0, calpha, a, b, 0 )
162 CALL chkxer(
'ZTFSM ', infot, nout, lerr, ok )
166 CALL ztftri(
'/',
'L',
'N', 0, a, info )
167 CALL chkxer(
'ZTFTRI', infot, nout, lerr, ok )
169 CALL ztftri(
'N',
'/',
'N', 0, a, info )
170 CALL chkxer(
'ZTFTRI', infot, nout, lerr, ok )
172 CALL ztftri(
'N',
'L',
'/', 0, a, info )
173 CALL chkxer(
'ZTFTRI', infot, nout, lerr, ok )
175 CALL ztftri(
'N',
'L',
'N', -1, a, info )
176 CALL chkxer(
'ZTFTRI', infot, nout, lerr, ok )
180 CALL ztfttr(
'/',
'U', 0, a, b, 1, info )
181 CALL chkxer(
'ZTFTTR', infot, nout, lerr, ok )
183 CALL ztfttr(
'N',
'/', 0, a, b, 1, info )
184 CALL chkxer(
'ZTFTTR', infot, nout, lerr, ok )
186 CALL ztfttr(
'N',
'U', -1, a, b, 1, info )
187 CALL chkxer(
'ZTFTTR', infot, nout, lerr, ok )
189 CALL ztfttr(
'N',
'U', 0, a, b, 0, info )
190 CALL chkxer(
'ZTFTTR', infot, nout, lerr, ok )
194 CALL ztrttf(
'/',
'U', 0, a, 1, b, info )
195 CALL chkxer(
'ZTRTTF', infot, nout, lerr, ok )
197 CALL ztrttf(
'N',
'/', 0, a, 1, b, info )
198 CALL chkxer(
'ZTRTTF', infot, nout, lerr, ok )
200 CALL ztrttf(
'N',
'U', -1, a, 1, b, info )
201 CALL chkxer(
'ZTRTTF', infot, nout, lerr, ok )
203 CALL ztrttf(
'N',
'U', 0, a, 0, b, info )
204 CALL chkxer(
'ZTRTTF', infot, nout, lerr, ok )
208 CALL ztfttp(
'/',
'U', 0, a, b, info )
209 CALL chkxer(
'ZTFTTP', infot, nout, lerr, ok )
211 CALL ztfttp(
'N',
'/', 0, a, b, info )
212 CALL chkxer(
'ZTFTTP', infot, nout, lerr, ok )
214 CALL ztfttp(
'N',
'U', -1, a, b, info )
215 CALL chkxer(
'ZTFTTP', infot, nout, lerr, ok )
219 CALL ztpttf(
'/',
'U', 0, a, b, info )
220 CALL chkxer(
'ZTPTTF', infot, nout, lerr, ok )
222 CALL ztpttf(
'N',
'/', 0, a, b, info )
223 CALL chkxer(
'ZTPTTF', infot, nout, lerr, ok )
225 CALL ztpttf(
'N',
'U', -1, a, b, info )
226 CALL chkxer(
'ZTPTTF', infot, nout, lerr, ok )
230 CALL ztrttp(
'/', 0, a, 1, b, info )
231 CALL chkxer(
'ZTRTTP', infot, nout, lerr, ok )
233 CALL ztrttp(
'U', -1, a, 1, b, info )
234 CALL chkxer(
'ZTRTTP', infot, nout, lerr, ok )
236 CALL ztrttp(
'U', 0, a, 0, b, info )
237 CALL chkxer(
'ZTRTTP', infot, nout, lerr, ok )
241 CALL ztpttr(
'/', 0, a, b, 1, info )
242 CALL chkxer(
'ZTPTTR', infot, nout, lerr, ok )
244 CALL ztpttr(
'U', -1, a, b, 1, info )
245 CALL chkxer(
'ZTPTTR', infot, nout, lerr, ok )
247 CALL ztpttr(
'U', 0, a, b, 0, info )
248 CALL chkxer(
'ZTPTTR', infot, nout, lerr, ok )
252 CALL zhfrk(
'/',
'U',
'N', 0, 0, alpha, a, 1, beta, b )
253 CALL chkxer(
'ZHFRK ', infot, nout, lerr, ok )
255 CALL zhfrk(
'N',
'/',
'N', 0, 0, alpha, a, 1, beta, b )
256 CALL chkxer(
'ZHFRK ', infot, nout, lerr, ok )
258 CALL zhfrk(
'N',
'U',
'/', 0, 0, alpha, a, 1, beta, b )
259 CALL chkxer(
'ZHFRK ', infot, nout, lerr, ok )
261 CALL zhfrk(
'N',
'U',
'N', -1, 0, alpha, a, 1, beta, b )
262 CALL chkxer(
'ZHFRK ', infot, nout, lerr, ok )
264 CALL zhfrk(
'N',
'U',
'N', 0, -1, alpha, a, 1, beta, b )
265 CALL chkxer(
'ZHFRK ', infot, nout, lerr, ok )
267 CALL zhfrk(
'N',
'U',
'N', 0, 0, alpha, a, 0, beta, b )
268 CALL chkxer(
'ZHFRK ', infot, nout, lerr, ok )
273 WRITE( nout, fmt = 9999 )
275 WRITE( nout, fmt = 9998 )
278 9999
FORMAT( 1x,
'COMPLEX*16 RFP routines passed the tests of the ',
280 9998
FORMAT(
' *** RFP routines failed the tests of the error ',
subroutine chkxer(srnamt, infot, nout, lerr, ok)
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 zpftrf(transr, uplo, n, a, info)
ZPFTRF
subroutine zpftri(transr, uplo, n, a, info)
ZPFTRI
subroutine zpftrs(transr, uplo, n, nrhs, a, b, ldb, info)
ZPFTRS
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).
subroutine ztftri(transr, uplo, diag, n, a, info)
ZTFTRI
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 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...
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 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 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 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 zerrrfp(nunit)
ZERRRFP