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 ctfttp(TRANSR, UPLO, N, ARF, AP, INFO)
CTFTTP copies a triangular matrix from the rectangular full packed format (TF) to the standard 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 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 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 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 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 chfrk(TRANSR, UPLO, TRANS, N, K, ALPHA, A, LDA, BETA, C)
CHFRK performs a Hermitian rank-k operation for matrix in RFP format.
subroutine cpftri(TRANSR, UPLO, N, A, INFO)
CPFTRI
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...