69 COMPLEX A( 1, 1), B( 1, 1)
85 COMMON / infoc / infot, nout, ok, lerr
86 COMMON / srnamc / srnamt
92 a( 1, 1 ) = cmplx( 1.d0 , 1.d0 )
93 b( 1, 1 ) = cmplx( 1.d0 , 1.d0 )
94 alpha = cmplx( 1.d0 , 1.d0 )
95 beta = cmplx( 1.d0 , 1.d0 )
99 CALL cpftrf(
'/',
'U', 0, a, info )
100 CALL chkxer(
'CPFTRF', infot, nout, lerr, ok )
102 CALL cpftrf(
'N',
'/', 0, a, info )
103 CALL chkxer(
'CPFTRF', infot, nout, lerr, ok )
105 CALL cpftrf(
'N',
'U', -1, a, info )
106 CALL chkxer(
'CPFTRF', infot, nout, lerr, ok )
110 CALL cpftrs(
'/',
'U', 0, 0, a, b, 1, info )
111 CALL chkxer(
'CPFTRS', infot, nout, lerr, ok )
113 CALL cpftrs(
'N',
'/', 0, 0, a, b, 1, info )
114 CALL chkxer(
'CPFTRS', infot, nout, lerr, ok )
116 CALL cpftrs(
'N',
'U', -1, 0, a, b, 1, info )
117 CALL chkxer(
'CPFTRS', infot, nout, lerr, ok )
119 CALL cpftrs(
'N',
'U', 0, -1, a, b, 1, info )
120 CALL chkxer(
'CPFTRS', infot, nout, lerr, ok )
122 CALL cpftrs(
'N',
'U', 0, 0, a, b, 0, info )
123 CALL chkxer(
'CPFTRS', infot, nout, lerr, ok )
127 CALL cpftri(
'/',
'U', 0, a, info )
128 CALL chkxer(
'CPFTRI', infot, nout, lerr, ok )
130 CALL cpftri(
'N',
'/', 0, a, info )
131 CALL chkxer(
'CPFTRI', infot, nout, lerr, ok )
133 CALL cpftri(
'N',
'U', -1, a, info )
134 CALL chkxer(
'CPFTRI', infot, nout, lerr, ok )
138 CALL ctfsm(
'/',
'L',
'U',
'C',
'U', 0, 0, alpha, a, b, 1 )
139 CALL chkxer(
'CTFSM ', infot, nout, lerr, ok )
141 CALL ctfsm(
'N',
'/',
'U',
'C',
'U', 0, 0, alpha, a, b, 1 )
142 CALL chkxer(
'CTFSM ', infot, nout, lerr, ok )
144 CALL ctfsm(
'N',
'L',
'/',
'C',
'U', 0, 0, alpha, a, b, 1 )
145 CALL chkxer(
'CTFSM ', infot, nout, lerr, ok )
147 CALL ctfsm(
'N',
'L',
'U',
'/',
'U', 0, 0, alpha, a, b, 1 )
148 CALL chkxer(
'CTFSM ', infot, nout, lerr, ok )
150 CALL ctfsm(
'N',
'L',
'U',
'C',
'/', 0, 0, alpha, a, b, 1 )
151 CALL chkxer(
'CTFSM ', infot, nout, lerr, ok )
153 CALL ctfsm(
'N',
'L',
'U',
'C',
'U', -1, 0, alpha, a, b, 1 )
154 CALL chkxer(
'CTFSM ', infot, nout, lerr, ok )
156 CALL ctfsm(
'N',
'L',
'U',
'C',
'U', 0, -1, alpha, a, b, 1 )
157 CALL chkxer(
'CTFSM ', infot, nout, lerr, ok )
159 CALL ctfsm(
'N',
'L',
'U',
'C',
'U', 0, 0, alpha, a, b, 0 )
160 CALL chkxer(
'CTFSM ', infot, nout, lerr, ok )
164 CALL ctftri(
'/',
'L',
'N', 0, a, info )
165 CALL chkxer(
'CTFTRI', infot, nout, lerr, ok )
167 CALL ctftri(
'N',
'/',
'N', 0, a, info )
168 CALL chkxer(
'CTFTRI', infot, nout, lerr, ok )
170 CALL ctftri(
'N',
'L',
'/', 0, a, info )
171 CALL chkxer(
'CTFTRI', infot, nout, lerr, ok )
173 CALL ctftri(
'N',
'L',
'N', -1, a, info )
174 CALL chkxer(
'CTFTRI', infot, nout, lerr, ok )
178 CALL ctfttr(
'/',
'U', 0, a, b, 1, info )
179 CALL chkxer(
'CTFTTR', infot, nout, lerr, ok )
181 CALL ctfttr(
'N',
'/', 0, a, b, 1, info )
182 CALL chkxer(
'CTFTTR', infot, nout, lerr, ok )
184 CALL ctfttr(
'N',
'U', -1, a, b, 1, info )
185 CALL chkxer(
'CTFTTR', infot, nout, lerr, ok )
187 CALL ctfttr(
'N',
'U', 0, a, b, 0, info )
188 CALL chkxer(
'CTFTTR', infot, nout, lerr, ok )
192 CALL ctrttf(
'/',
'U', 0, a, 1, b, info )
193 CALL chkxer(
'CTRTTF', infot, nout, lerr, ok )
195 CALL ctrttf(
'N',
'/', 0, a, 1, b, info )
196 CALL chkxer(
'CTRTTF', infot, nout, lerr, ok )
198 CALL ctrttf(
'N',
'U', -1, a, 1, b, info )
199 CALL chkxer(
'CTRTTF', infot, nout, lerr, ok )
201 CALL ctrttf(
'N',
'U', 0, a, 0, b, info )
202 CALL chkxer(
'CTRTTF', infot, nout, lerr, ok )
206 CALL ctfttp(
'/',
'U', 0, a, b, info )
207 CALL chkxer(
'CTFTTP', infot, nout, lerr, ok )
209 CALL ctfttp(
'N',
'/', 0, a, b, info )
210 CALL chkxer(
'CTFTTP', infot, nout, lerr, ok )
212 CALL ctfttp(
'N',
'U', -1, a, b, info )
213 CALL chkxer(
'CTFTTP', infot, nout, lerr, ok )
217 CALL ctpttf(
'/',
'U', 0, a, b, info )
218 CALL chkxer(
'CTPTTF', infot, nout, lerr, ok )
220 CALL ctpttf(
'N',
'/', 0, a, b, info )
221 CALL chkxer(
'CTPTTF', infot, nout, lerr, ok )
223 CALL ctpttf(
'N',
'U', -1, a, b, info )
224 CALL chkxer(
'CTPTTF', infot, nout, lerr, ok )
228 CALL ctrttp(
'/', 0, a, 1, b, info )
229 CALL chkxer(
'CTRTTP', infot, nout, lerr, ok )
231 CALL ctrttp(
'U', -1, a, 1, b, info )
232 CALL chkxer(
'CTRTTP', infot, nout, lerr, ok )
234 CALL ctrttp(
'U', 0, a, 0, b, info )
235 CALL chkxer(
'CTRTTP', infot, nout, lerr, ok )
239 CALL ctpttr(
'/', 0, a, b, 1, info )
240 CALL chkxer(
'CTPTTR', infot, nout, lerr, ok )
242 CALL ctpttr(
'U', -1, a, b, 1, info )
243 CALL chkxer(
'CTPTTR', infot, nout, lerr, ok )
245 CALL ctpttr(
'U', 0, a, b, 0, info )
246 CALL chkxer(
'CTPTTR', infot, nout, lerr, ok )
250 CALL chfrk(
'/',
'U',
'N', 0, 0, alpha, a, 1, beta, b )
251 CALL chkxer(
'CHFRK ', infot, nout, lerr, ok )
253 CALL chfrk(
'N',
'/',
'N', 0, 0, alpha, a, 1, beta, b )
254 CALL chkxer(
'CHFRK ', infot, nout, lerr, ok )
256 CALL chfrk(
'N',
'U',
'/', 0, 0, alpha, a, 1, beta, b )
257 CALL chkxer(
'CHFRK ', infot, nout, lerr, ok )
259 CALL chfrk(
'N',
'U',
'N', -1, 0, alpha, a, 1, beta, b )
260 CALL chkxer(
'CHFRK ', infot, nout, lerr, ok )
262 CALL chfrk(
'N',
'U',
'N', 0, -1, alpha, a, 1, beta, b )
263 CALL chkxer(
'CHFRK ', infot, nout, lerr, ok )
265 CALL chfrk(
'N',
'U',
'N', 0, 0, alpha, a, 0, beta, b )
266 CALL chkxer(
'CHFRK ', infot, nout, lerr, ok )
271 WRITE( nout, fmt = 9999 )
273 WRITE( nout, fmt = 9998 )
276 9999
FORMAT( 1x,
'COMPLEX RFP routines passed the tests of the ',
278 9998
FORMAT(
' *** RFP routines failed the tests of the error ',
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
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 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 ctftri(TRANSR, UPLO, DIAG, N, A, INFO)
CTFTRI
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 cpftrs(TRANSR, UPLO, N, NRHS, A, B, LDB, INFO)
CPFTRS
subroutine cpftri(TRANSR, UPLO, N, A, INFO)
CPFTRI
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 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 cpftrf(TRANSR, UPLO, N, A, INFO)
CPFTRF
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 ctfttp(TRANSR, UPLO, N, ARF, AP, INFO)
CTFTTP copies a triangular matrix from the rectangular full packed format (TF) to the standard packed...
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...