72 COMPLEX A( 1, 1), B( 1, 1)
88 COMMON / infoc / infot, nout, ok, lerr
89 COMMON / srnamc / srnamt
95 a( 1, 1 ) = cmplx( 1.d0 , 1.d0 )
96 b( 1, 1 ) = cmplx( 1.d0 , 1.d0 )
97 alpha = cmplx( 1.d0 , 1.d0 )
98 beta = cmplx( 1.d0 , 1.d0 )
102 CALL
cpftrf(
'/',
'U', 0, a, info )
103 CALL
chkxer(
'CPFTRF', infot, nout, lerr, ok )
105 CALL
cpftrf(
'N',
'/', 0, a, info )
106 CALL
chkxer(
'CPFTRF', infot, nout, lerr, ok )
108 CALL
cpftrf(
'N',
'U', -1, a, info )
109 CALL
chkxer(
'CPFTRF', infot, nout, lerr, ok )
113 CALL
cpftrs(
'/',
'U', 0, 0, a, b, 1, info )
114 CALL
chkxer(
'CPFTRS', infot, nout, lerr, ok )
116 CALL
cpftrs(
'N',
'/', 0, 0, a, b, 1, info )
117 CALL
chkxer(
'CPFTRS', infot, nout, lerr, ok )
119 CALL
cpftrs(
'N',
'U', -1, 0, a, b, 1, info )
120 CALL
chkxer(
'CPFTRS', infot, nout, lerr, ok )
122 CALL
cpftrs(
'N',
'U', 0, -1, a, b, 1, info )
123 CALL
chkxer(
'CPFTRS', infot, nout, lerr, ok )
125 CALL
cpftrs(
'N',
'U', 0, 0, a, b, 0, info )
126 CALL
chkxer(
'CPFTRS', infot, nout, lerr, ok )
130 CALL
cpftri(
'/',
'U', 0, a, info )
131 CALL
chkxer(
'CPFTRI', infot, nout, lerr, ok )
133 CALL
cpftri(
'N',
'/', 0, a, info )
134 CALL
chkxer(
'CPFTRI', infot, nout, lerr, ok )
136 CALL
cpftri(
'N',
'U', -1, a, info )
137 CALL
chkxer(
'CPFTRI', infot, nout, lerr, ok )
141 CALL
ctfsm(
'/',
'L',
'U',
'C',
'U', 0, 0, alpha, a, b, 1 )
142 CALL
chkxer(
'CTFSM ', infot, nout, lerr, ok )
144 CALL
ctfsm(
'N',
'/',
'U',
'C',
'U', 0, 0, alpha, a, b, 1 )
145 CALL
chkxer(
'CTFSM ', infot, nout, lerr, ok )
147 CALL
ctfsm(
'N',
'L',
'/',
'C',
'U', 0, 0, alpha, a, b, 1 )
148 CALL
chkxer(
'CTFSM ', infot, nout, lerr, ok )
150 CALL
ctfsm(
'N',
'L',
'U',
'/',
'U', 0, 0, alpha, a, b, 1 )
151 CALL
chkxer(
'CTFSM ', infot, nout, lerr, ok )
153 CALL
ctfsm(
'N',
'L',
'U',
'C',
'/', 0, 0, alpha, a, b, 1 )
154 CALL
chkxer(
'CTFSM ', infot, nout, lerr, ok )
156 CALL
ctfsm(
'N',
'L',
'U',
'C',
'U', -1, 0, alpha, a, b, 1 )
157 CALL
chkxer(
'CTFSM ', infot, nout, lerr, ok )
159 CALL
ctfsm(
'N',
'L',
'U',
'C',
'U', 0, -1, alpha, a, b, 1 )
160 CALL
chkxer(
'CTFSM ', infot, nout, lerr, ok )
162 CALL
ctfsm(
'N',
'L',
'U',
'C',
'U', 0, 0, alpha, a, b, 0 )
163 CALL
chkxer(
'CTFSM ', infot, nout, lerr, ok )
167 CALL
ctftri(
'/',
'L',
'N', 0, a, info )
168 CALL
chkxer(
'CTFTRI', infot, nout, lerr, ok )
170 CALL
ctftri(
'N',
'/',
'N', 0, a, info )
171 CALL
chkxer(
'CTFTRI', infot, nout, lerr, ok )
173 CALL
ctftri(
'N',
'L',
'/', 0, a, info )
174 CALL
chkxer(
'CTFTRI', infot, nout, lerr, ok )
176 CALL
ctftri(
'N',
'L',
'N', -1, a, info )
177 CALL
chkxer(
'CTFTRI', infot, nout, lerr, ok )
181 CALL
ctfttr(
'/',
'U', 0, a, b, 1, info )
182 CALL
chkxer(
'CTFTTR', infot, nout, lerr, ok )
184 CALL
ctfttr(
'N',
'/', 0, a, b, 1, info )
185 CALL
chkxer(
'CTFTTR', infot, nout, lerr, ok )
187 CALL
ctfttr(
'N',
'U', -1, a, b, 1, info )
188 CALL
chkxer(
'CTFTTR', infot, nout, lerr, ok )
190 CALL
ctfttr(
'N',
'U', 0, a, b, 0, info )
191 CALL
chkxer(
'CTFTTR', infot, nout, lerr, ok )
195 CALL
ctrttf(
'/',
'U', 0, a, 1, b, info )
196 CALL
chkxer(
'CTRTTF', infot, nout, lerr, ok )
198 CALL
ctrttf(
'N',
'/', 0, a, 1, b, info )
199 CALL
chkxer(
'CTRTTF', infot, nout, lerr, ok )
201 CALL
ctrttf(
'N',
'U', -1, a, 1, b, info )
202 CALL
chkxer(
'CTRTTF', infot, nout, lerr, ok )
204 CALL
ctrttf(
'N',
'U', 0, a, 0, b, info )
205 CALL
chkxer(
'CTRTTF', infot, nout, lerr, ok )
209 CALL
ctfttp(
'/',
'U', 0, a, b, info )
210 CALL
chkxer(
'CTFTTP', infot, nout, lerr, ok )
212 CALL
ctfttp(
'N',
'/', 0, a, b, info )
213 CALL
chkxer(
'CTFTTP', infot, nout, lerr, ok )
215 CALL
ctfttp(
'N',
'U', -1, a, b, info )
216 CALL
chkxer(
'CTFTTP', infot, nout, lerr, ok )
220 CALL
ctpttf(
'/',
'U', 0, a, b, info )
221 CALL
chkxer(
'CTPTTF', infot, nout, lerr, ok )
223 CALL
ctpttf(
'N',
'/', 0, a, b, info )
224 CALL
chkxer(
'CTPTTF', infot, nout, lerr, ok )
226 CALL
ctpttf(
'N',
'U', -1, a, b, info )
227 CALL
chkxer(
'CTPTTF', infot, nout, lerr, ok )
231 CALL
ctrttp(
'/', 0, a, 1, b, info )
232 CALL
chkxer(
'CTRTTP', infot, nout, lerr, ok )
234 CALL
ctrttp(
'U', -1, a, 1, b, info )
235 CALL
chkxer(
'CTRTTP', infot, nout, lerr, ok )
237 CALL
ctrttp(
'U', 0, a, 0, b, info )
238 CALL
chkxer(
'CTRTTP', infot, nout, lerr, ok )
242 CALL
ctpttr(
'/', 0, a, b, 1, info )
243 CALL
chkxer(
'CTPTTR', infot, nout, lerr, ok )
245 CALL
ctpttr(
'U', -1, a, b, 1, info )
246 CALL
chkxer(
'CTPTTR', infot, nout, lerr, ok )
248 CALL
ctpttr(
'U', 0, a, b, 0, info )
249 CALL
chkxer(
'CTPTTR', infot, nout, lerr, ok )
253 CALL
chfrk(
'/',
'U',
'N', 0, 0, alpha, a, 1, beta, b )
254 CALL
chkxer(
'CHFRK ', infot, nout, lerr, ok )
256 CALL
chfrk(
'N',
'/',
'N', 0, 0, alpha, a, 1, beta, b )
257 CALL
chkxer(
'CHFRK ', infot, nout, lerr, ok )
259 CALL
chfrk(
'N',
'U',
'/', 0, 0, alpha, a, 1, beta, b )
260 CALL
chkxer(
'CHFRK ', infot, nout, lerr, ok )
262 CALL
chfrk(
'N',
'U',
'N', -1, 0, alpha, a, 1, beta, b )
263 CALL
chkxer(
'CHFRK ', infot, nout, lerr, ok )
265 CALL
chfrk(
'N',
'U',
'N', 0, -1, alpha, a, 1, beta, b )
266 CALL
chkxer(
'CHFRK ', infot, nout, lerr, ok )
268 CALL
chfrk(
'N',
'U',
'N', 0, 0, alpha, a, 0, beta, b )
269 CALL
chkxer(
'CHFRK ', infot, nout, lerr, ok )
274 WRITE( nout, fmt = 9999 )
276 WRITE( nout, fmt = 9998 )
279 9999
FORMAT( 1x,
'COMPLEX RFP routines passed the tests of the ',
281 9998
FORMAT(
' *** RFP routines failed the tests of the error ',
subroutine cpftrf(TRANSR, UPLO, N, A, INFO)
CPFTRF
subroutine ctrttf(TRANSR, UPLO, N, A, LDA, ARF, INFO)
CTRTTF copies a triangular matrix from the standard full format (TR) to the rectangular full packed f...
subroutine cpftrs(TRANSR, UPLO, N, NRHS, A, B, LDB, INFO)
CPFTRS
subroutine ctpttr(UPLO, N, AP, A, LDA, INFO)
CTPTTR copies a triangular matrix from the standard packed format (TP) to the standard full format (T...
subroutine cerrrfp(NUNIT)
CERRRFP
subroutine ctfsm(TRANSR, SIDE, UPLO, TRANS, DIAG, M, N, ALPHA, A, B, LDB)
CTFSM solves a matrix equation (one operand is a triangular matrix in RFP format).
subroutine ctftri(TRANSR, UPLO, DIAG, N, A, INFO)
CTFTRI
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
subroutine ctpttf(TRANSR, UPLO, N, AP, ARF, INFO)
CTPTTF copies a triangular matrix from the standard packed format (TP) to the rectangular full packed...
subroutine ctrttp(UPLO, N, A, LDA, AP, INFO)
CTRTTP copies a triangular matrix from the standard full format (TR) to the standard packed format (T...
subroutine ctfttr(TRANSR, UPLO, N, ARF, A, LDA, INFO)
CTFTTR copies a triangular matrix from the rectangular full packed format (TF) to the standard full f...
subroutine cpftri(TRANSR, UPLO, N, A, INFO)
CPFTRI
subroutine chfrk(TRANSR, UPLO, TRANS, N, K, ALPHA, A, LDA, BETA, C)
CHFRK performs a Hermitian rank-k operation for matrix in RFP format.
subroutine ctfttp(TRANSR, UPLO, N, ARF, AP, INFO)
CTFTTP copies a triangular matrix from the rectangular full packed format (TF) to the standard packed...