72 REAL A( 1, 1), B( 1, 1)
85 COMMON / infoc / infot, nout, ok, lerr
86 COMMON / srnamc / srnamt
99 CALL spftrf(
'/',
'U', 0, a, info )
100 CALL chkxer(
'SPFTRF', infot, nout, lerr, ok )
102 CALL spftrf(
'N',
'/', 0, a, info )
103 CALL chkxer(
'SPFTRF', infot, nout, lerr, ok )
105 CALL spftrf(
'N',
'U', -1, a, info )
106 CALL chkxer(
'SPFTRF', infot, nout, lerr, ok )
110 CALL spftrs(
'/',
'U', 0, 0, a, b, 1, info )
111 CALL chkxer(
'SPFTRS', infot, nout, lerr, ok )
113 CALL spftrs(
'N',
'/', 0, 0, a, b, 1, info )
114 CALL chkxer(
'SPFTRS', infot, nout, lerr, ok )
116 CALL spftrs(
'N',
'U', -1, 0, a, b, 1, info )
117 CALL chkxer(
'SPFTRS', infot, nout, lerr, ok )
119 CALL spftrs(
'N',
'U', 0, -1, a, b, 1, info )
120 CALL chkxer(
'SPFTRS', infot, nout, lerr, ok )
122 CALL spftrs(
'N',
'U', 0, 0, a, b, 0, info )
123 CALL chkxer(
'SPFTRS', infot, nout, lerr, ok )
127 CALL spftri(
'/',
'U', 0, a, info )
128 CALL chkxer(
'SPFTRI', infot, nout, lerr, ok )
130 CALL spftri(
'N',
'/', 0, a, info )
131 CALL chkxer(
'SPFTRI', infot, nout, lerr, ok )
133 CALL spftri(
'N',
'U', -1, a, info )
134 CALL chkxer(
'SPFTRI', infot, nout, lerr, ok )
138 CALL stfsm(
'/',
'L',
'U',
'T',
'U', 0, 0, alpha, a, b, 1 )
139 CALL chkxer(
'STFSM ', infot, nout, lerr, ok )
141 CALL stfsm(
'N',
'/',
'U',
'T',
'U', 0, 0, alpha, a, b, 1 )
142 CALL chkxer(
'STFSM ', infot, nout, lerr, ok )
144 CALL stfsm(
'N',
'L',
'/',
'T',
'U', 0, 0, alpha, a, b, 1 )
145 CALL chkxer(
'STFSM ', infot, nout, lerr, ok )
147 CALL stfsm(
'N',
'L',
'U',
'/',
'U', 0, 0, alpha, a, b, 1 )
148 CALL chkxer(
'STFSM ', infot, nout, lerr, ok )
150 CALL stfsm(
'N',
'L',
'U',
'T',
'/', 0, 0, alpha, a, b, 1 )
151 CALL chkxer(
'STFSM ', infot, nout, lerr, ok )
153 CALL stfsm(
'N',
'L',
'U',
'T',
'U', -1, 0, alpha, a, b, 1 )
154 CALL chkxer(
'STFSM ', infot, nout, lerr, ok )
156 CALL stfsm(
'N',
'L',
'U',
'T',
'U', 0, -1, alpha, a, b, 1 )
157 CALL chkxer(
'STFSM ', infot, nout, lerr, ok )
159 CALL stfsm(
'N',
'L',
'U',
'T',
'U', 0, 0, alpha, a, b, 0 )
160 CALL chkxer(
'STFSM ', infot, nout, lerr, ok )
164 CALL stftri(
'/',
'L',
'N', 0, a, info )
165 CALL chkxer(
'STFTRI', infot, nout, lerr, ok )
167 CALL stftri(
'N',
'/',
'N', 0, a, info )
168 CALL chkxer(
'STFTRI', infot, nout, lerr, ok )
170 CALL stftri(
'N',
'L',
'/', 0, a, info )
171 CALL chkxer(
'STFTRI', infot, nout, lerr, ok )
173 CALL stftri(
'N',
'L',
'N', -1, a, info )
174 CALL chkxer(
'STFTRI', infot, nout, lerr, ok )
178 CALL stfttr(
'/',
'U', 0, a, b, 1, info )
179 CALL chkxer(
'STFTTR', infot, nout, lerr, ok )
181 CALL stfttr(
'N',
'/', 0, a, b, 1, info )
182 CALL chkxer(
'STFTTR', infot, nout, lerr, ok )
184 CALL stfttr(
'N',
'U', -1, a, b, 1, info )
185 CALL chkxer(
'STFTTR', infot, nout, lerr, ok )
187 CALL stfttr(
'N',
'U', 0, a, b, 0, info )
188 CALL chkxer(
'STFTTR', infot, nout, lerr, ok )
192 CALL strttf(
'/',
'U', 0, a, 1, b, info )
193 CALL chkxer(
'STRTTF', infot, nout, lerr, ok )
195 CALL strttf(
'N',
'/', 0, a, 1, b, info )
196 CALL chkxer(
'STRTTF', infot, nout, lerr, ok )
198 CALL strttf(
'N',
'U', -1, a, 1, b, info )
199 CALL chkxer(
'STRTTF', infot, nout, lerr, ok )
201 CALL strttf(
'N',
'U', 0, a, 0, b, info )
202 CALL chkxer(
'STRTTF', infot, nout, lerr, ok )
206 CALL stfttp(
'/',
'U', 0, a, b, info )
207 CALL chkxer(
'STFTTP', infot, nout, lerr, ok )
209 CALL stfttp(
'N',
'/', 0, a, b, info )
210 CALL chkxer(
'STFTTP', infot, nout, lerr, ok )
212 CALL stfttp(
'N',
'U', -1, a, b, info )
213 CALL chkxer(
'STFTTP', infot, nout, lerr, ok )
217 CALL stpttf(
'/',
'U', 0, a, b, info )
218 CALL chkxer(
'STPTTF', infot, nout, lerr, ok )
220 CALL stpttf(
'N',
'/', 0, a, b, info )
221 CALL chkxer(
'STPTTF', infot, nout, lerr, ok )
223 CALL stpttf(
'N',
'U', -1, a, b, info )
224 CALL chkxer(
'STPTTF', infot, nout, lerr, ok )
228 CALL strttp(
'/', 0, a, 1, b, info )
229 CALL chkxer(
'STRTTP', infot, nout, lerr, ok )
231 CALL strttp(
'U', -1, a, 1, b, info )
232 CALL chkxer(
'STRTTP', infot, nout, lerr, ok )
234 CALL strttp(
'U', 0, a, 0, b, info )
235 CALL chkxer(
'STRTTP', infot, nout, lerr, ok )
239 CALL stpttr(
'/', 0, a, b, 1, info )
240 CALL chkxer(
'STPTTR', infot, nout, lerr, ok )
242 CALL stpttr(
'U', -1, a, b, 1, info )
243 CALL chkxer(
'STPTTR', infot, nout, lerr, ok )
245 CALL stpttr(
'U', 0, a, b, 0, info )
246 CALL chkxer(
'STPTTR', infot, nout, lerr, ok )
250 CALL ssfrk(
'/',
'U',
'N', 0, 0, alpha, a, 1, beta, b )
251 CALL chkxer(
'SSFRK ', infot, nout, lerr, ok )
253 CALL ssfrk(
'N',
'/',
'N', 0, 0, alpha, a, 1, beta, b )
254 CALL chkxer(
'SSFRK ', infot, nout, lerr, ok )
256 CALL ssfrk(
'N',
'U',
'/', 0, 0, alpha, a, 1, beta, b )
257 CALL chkxer(
'SSFRK ', infot, nout, lerr, ok )
259 CALL ssfrk(
'N',
'U',
'N', -1, 0, alpha, a, 1, beta, b )
260 CALL chkxer(
'SSFRK ', infot, nout, lerr, ok )
262 CALL ssfrk(
'N',
'U',
'N', 0, -1, alpha, a, 1, beta, b )
263 CALL chkxer(
'SSFRK ', infot, nout, lerr, ok )
265 CALL ssfrk(
'N',
'U',
'N', 0, 0, alpha, a, 0, beta, b )
266 CALL chkxer(
'SSFRK ', infot, nout, lerr, ok )
271 WRITE( nout, fmt = 9999 )
273 WRITE( nout, fmt = 9998 )
276 9999
FORMAT( 1x,
'REAL RFP routines passed the tests of ',
277 $
'the error exits' )
278 9998
FORMAT(
' *** RFP routines failed the tests of the error ',
subroutine ssfrk(TRANSR, UPLO, TRANS, N, K, ALPHA, A, LDA, BETA, C)
SSFRK performs a symmetric rank-k operation for matrix in RFP format.
subroutine strttf(TRANSR, UPLO, N, A, LDA, ARF, INFO)
STRTTF copies a triangular matrix from the standard full format (TR) to the rectangular full packed f...
subroutine stfttp(TRANSR, UPLO, N, ARF, AP, INFO)
STFTTP copies a triangular matrix from the rectangular full packed format (TF) to the standard packed...
subroutine stpttr(UPLO, N, AP, A, LDA, INFO)
STPTTR copies a triangular matrix from the standard packed format (TP) to the standard full format (T...
subroutine chkxer(SRNAMT, INFOT, NOUT, LERR, OK)
subroutine stftri(TRANSR, UPLO, DIAG, N, A, INFO)
STFTRI
subroutine spftrf(TRANSR, UPLO, N, A, INFO)
SPFTRF
subroutine stfttr(TRANSR, UPLO, N, ARF, A, LDA, INFO)
STFTTR copies a triangular matrix from the rectangular full packed format (TF) to the standard full f...
subroutine serrrfp(NUNIT)
SERRRFP
subroutine stpttf(TRANSR, UPLO, N, AP, ARF, INFO)
STPTTF copies a triangular matrix from the standard packed format (TP) to the rectangular full packed...
subroutine spftri(TRANSR, UPLO, N, A, INFO)
SPFTRI
subroutine spftrs(TRANSR, UPLO, N, NRHS, A, B, LDB, INFO)
SPFTRS
subroutine stfsm(TRANSR, SIDE, UPLO, TRANS, DIAG, M, N, ALPHA, A, B, LDB)
STFSM solves a matrix equation (one operand is a triangular matrix in RFP format).
subroutine strttp(UPLO, N, A, LDA, AP, INFO)
STRTTP copies a triangular matrix from the standard full format (TR) to the standard packed format (T...