OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
cerrrfp.f
Go to the documentation of this file.
1*> \brief \b CERRRFP
2*
3* =========== DOCUMENTATION ===========
4*
5* Online html documentation available at
6* http://www.netlib.org/lapack/explore-html/
7*
8* Definition:
9* ===========
10*
11* SUBROUTINE CERRRFP( NUNIT )
12*
13* .. Scalar Arguments ..
14* INTEGER NUNIT
15* ..
16*
17*
18*> \par Purpose:
19* =============
20*>
21*> \verbatim
22*>
23*> CERRRFP tests the error exits for the COMPLEX driver routines
24*> for solving linear systems of equations.
25*>
26*> CDRVRFP tests the COMPLEX LAPACK RFP routines:
27*> CTFSM, CTFTRI, CHFRK, CTFTTP, CTFTTR, CPFTRF, CPFTRS, CTPTTF,
28*> CTPTTR, CTRTTF, and CTRTTP
29*> \endverbatim
30*
31* Arguments:
32* ==========
33*
34*> \param[in] NUNIT
35*> \verbatim
36*> NUNIT is INTEGER
37*> The unit number for output.
38*> \endverbatim
39*
40* Authors:
41* ========
42*
43*> \author Univ. of Tennessee
44*> \author Univ. of California Berkeley
45*> \author Univ. of Colorado Denver
46*> \author NAG Ltd.
47*
48*> \ingroup complex_lin
49*
50* =====================================================================
51 SUBROUTINE cerrrfp( NUNIT )
52*
53* -- LAPACK test routine --
54* -- LAPACK is a software package provided by Univ. of Tennessee, --
55* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
56*
57* .. Scalar Arguments ..
58 INTEGER NUNIT
59* ..
60*
61* =====================================================================
62*
63* ..
64* .. Local Scalars ..
65 INTEGER INFO
66 COMPLEX ALPHA, BETA
67* ..
68* .. Local Arrays ..
69 COMPLEX A( 1, 1), B( 1, 1)
70* ..
71* .. External Subroutines ..
72 EXTERNAL chkxer, ctfsm, ctftri, chfrk, ctfttp, ctfttr,
73 + cpftri, cpftrf, cpftrs, ctpttf, ctpttr, ctrttf,
74 + ctrttp
75* ..
76* .. Scalars in Common ..
77 LOGICAL LERR, OK
78 CHARACTER*32 SRNAMT
79 INTEGER INFOT, NOUT
80* ..
81* .. Intrinsic Functions ..
82 INTRINSIC cmplx
83* ..
84* .. Common blocks ..
85 COMMON / infoc / infot, nout, ok, lerr
86 COMMON / srnamc / srnamt
87* ..
88* .. Executable Statements ..
89*
90 nout = nunit
91 ok = .true.
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 )
96*
97 srnamt = 'CPFTRF'
98 infot = 1
99 CALL cpftrf( '/', 'U', 0, a, info )
100 CALL chkxer( 'CPFTRF', infot, nout, lerr, ok )
101 infot = 2
102 CALL cpftrf( 'N', '/', 0, a, info )
103 CALL chkxer( 'CPFTRF', infot, nout, lerr, ok )
104 infot = 3
105 CALL cpftrf( 'N', 'U', -1, a, info )
106 CALL chkxer( 'CPFTRF', infot, nout, lerr, ok )
107*
108 srnamt = 'CPFTRS'
109 infot = 1
110 CALL cpftrs( '/', 'U', 0, 0, a, b, 1, info )
111 CALL chkxer( 'CPFTRS', infot, nout, lerr, ok )
112 infot = 2
113 CALL cpftrs( 'N', '/', 0, 0, a, b, 1, info )
114 CALL chkxer( 'CPFTRS', infot, nout, lerr, ok )
115 infot = 3
116 CALL cpftrs( 'N', 'U', -1, 0, a, b, 1, info )
117 CALL chkxer( 'CPFTRS', infot, nout, lerr, ok )
118 infot = 4
119 CALL cpftrs( 'N', 'U', 0, -1, a, b, 1, info )
120 CALL chkxer( 'CPFTRS', infot, nout, lerr, ok )
121 infot = 7
122 CALL cpftrs( 'N', 'U', 0, 0, a, b, 0, info )
123 CALL chkxer( 'CPFTRS', infot, nout, lerr, ok )
124*
125 srnamt = 'CPFTRI'
126 infot = 1
127 CALL cpftri( '/', 'u', 0, A, INFO )
128 CALL CHKXER( 'cpftri', INFOT, NOUT, LERR, OK )
129 INFOT = 2
130 CALL CPFTRI( 'n', '/', 0, A, INFO )
131 CALL CHKXER( 'cpftri', INFOT, NOUT, LERR, OK )
132 INFOT = 3
133 CALL CPFTRI( 'n', 'u', -1, A, INFO )
134 CALL CHKXER( 'cpftri', INFOT, NOUT, LERR, OK )
135*
136 SRNAMT = 'ctfsm '
137 INFOT = 1
138 CALL CTFSM( '/', 'l', 'u', 'c', 'u', 0, 0, ALPHA, A, B, 1 )
139 CALL CHKXER( 'ctfsm ', INFOT, NOUT, LERR, OK )
140 INFOT = 2
141 CALL CTFSM( 'n', '/', 'u', 'c', 'u', 0, 0, ALPHA, A, B, 1 )
142 CALL CHKXER( 'ctfsm ', INFOT, NOUT, LERR, OK )
143 INFOT = 3
144 CALL CTFSM( 'n', 'l', '/', 'c', 'u', 0, 0, ALPHA, A, B, 1 )
145 CALL CHKXER( 'ctfsm ', INFOT, NOUT, LERR, OK )
146 INFOT = 4
147 CALL CTFSM( 'n', 'l', 'u', '/', 'u', 0, 0, ALPHA, A, B, 1 )
148 CALL CHKXER( 'ctfsm ', INFOT, NOUT, LERR, OK )
149 INFOT = 5
150 CALL CTFSM( 'n', 'l', 'u', 'c', '/', 0, 0, ALPHA, A, B, 1 )
151 CALL CHKXER( 'ctfsm ', INFOT, NOUT, LERR, OK )
152 INFOT = 6
153 CALL CTFSM( 'n', 'l', 'u', 'c', 'u', -1, 0, ALPHA, A, B, 1 )
154 CALL CHKXER( 'ctfsm ', INFOT, NOUT, LERR, OK )
155 INFOT = 7
156 CALL CTFSM( 'n', 'l', 'u', 'c', 'u', 0, -1, ALPHA, A, B, 1 )
157 CALL CHKXER( 'ctfsm ', INFOT, NOUT, LERR, OK )
158 INFOT = 11
159 CALL CTFSM( 'n', 'l', 'u', 'c', 'u', 0, 0, ALPHA, A, B, 0 )
160 CALL CHKXER( 'ctfsm ', INFOT, NOUT, LERR, OK )
161*
162 SRNAMT = 'ctftri'
163 INFOT = 1
164 CALL CTFTRI( '/', 'l', 'n', 0, A, INFO )
165 CALL CHKXER( 'ctftri', INFOT, NOUT, LERR, OK )
166 INFOT = 2
167 CALL CTFTRI( 'n', '/', 'n', 0, A, INFO )
168 CALL CHKXER( 'ctftri', INFOT, NOUT, LERR, OK )
169 INFOT = 3
170 CALL CTFTRI( 'n', 'l', '/', 0, A, INFO )
171 CALL CHKXER( 'ctftri', INFOT, NOUT, LERR, OK )
172 INFOT = 4
173 CALL CTFTRI( 'n', 'l', 'n', -1, A, INFO )
174 CALL CHKXER( 'ctftri', INFOT, NOUT, LERR, OK )
175*
176 SRNAMT = 'ctfttr'
177 INFOT = 1
178 CALL CTFTTR( '/', 'u', 0, A, B, 1, INFO )
179 CALL CHKXER( 'ctfttr', INFOT, NOUT, LERR, OK )
180 INFOT = 2
181 CALL CTFTTR( 'n', '/', 0, A, B, 1, INFO )
182 CALL CHKXER( 'ctfttr', INFOT, NOUT, LERR, OK )
183 INFOT = 3
184 CALL CTFTTR( 'n', 'u', -1, A, B, 1, INFO )
185 CALL CHKXER( 'ctfttr', INFOT, NOUT, LERR, OK )
186 INFOT = 6
187 CALL CTFTTR( 'n', 'u', 0, A, B, 0, INFO )
188 CALL CHKXER( 'ctfttr', INFOT, NOUT, LERR, OK )
189*
190 SRNAMT = 'ctrttf'
191 INFOT = 1
192 CALL CTRTTF( '/', 'u', 0, A, 1, B, INFO )
193 CALL CHKXER( 'ctrttf', INFOT, NOUT, LERR, OK )
194 INFOT = 2
195 CALL CTRTTF( 'n', '/', 0, A, 1, B, INFO )
196 CALL CHKXER( 'ctrttf', INFOT, NOUT, LERR, OK )
197 INFOT = 3
198 CALL CTRTTF( 'n', 'u', -1, A, 1, B, INFO )
199 CALL CHKXER( 'ctrttf', INFOT, NOUT, LERR, OK )
200 INFOT = 5
201 CALL CTRTTF( 'n', 'u', 0, A, 0, B, INFO )
202 CALL CHKXER( 'ctrttf', INFOT, NOUT, LERR, OK )
203*
204 SRNAMT = 'ctfttp'
205 INFOT = 1
206 CALL CTFTTP( '/', 'u', 0, A, B, INFO )
207 CALL CHKXER( 'ctfttp', INFOT, NOUT, LERR, OK )
208 INFOT = 2
209 CALL CTFTTP( 'n', '/', 0, A, B, INFO )
210 CALL CHKXER( 'ctfttp', INFOT, NOUT, LERR, OK )
211 INFOT = 3
212 CALL CTFTTP( 'n', 'u', -1, A, B, INFO )
213 CALL CHKXER( 'ctfttp', INFOT, NOUT, LERR, OK )
214*
215 SRNAMT = 'ctpttf'
216 INFOT = 1
217 CALL CTPTTF( '/', 'u', 0, A, B, INFO )
218 CALL CHKXER( 'ctpttf', INFOT, NOUT, LERR, OK )
219 INFOT = 2
220 CALL CTPTTF( 'n', '/', 0, A, B, INFO )
221 CALL CHKXER( 'ctpttf', INFOT, NOUT, LERR, OK )
222 INFOT = 3
223 CALL CTPTTF( 'n', 'u', -1, A, B, INFO )
224 CALL CHKXER( 'ctpttf', INFOT, NOUT, LERR, OK )
225*
226 SRNAMT = 'ctrttp'
227 INFOT = 1
228 CALL CTRTTP( '/', 0, A, 1, B, INFO )
229 CALL CHKXER( 'ctrttp', INFOT, NOUT, LERR, OK )
230 INFOT = 2
231 CALL CTRTTP( 'u', -1, A, 1, B, INFO )
232 CALL CHKXER( 'ctrttp', INFOT, NOUT, LERR, OK )
233 INFOT = 4
234 CALL CTRTTP( 'u', 0, A, 0, B, INFO )
235 CALL CHKXER( 'ctrttp', INFOT, NOUT, LERR, OK )
236*
237 SRNAMT = 'ctpttr'
238 INFOT = 1
239 CALL CTPTTR( '/', 0, A, B, 1, INFO )
240 CALL CHKXER( 'ctpttr', INFOT, NOUT, LERR, OK )
241 INFOT = 2
242 CALL CTPTTR( 'u', -1, A, B, 1, INFO )
243 CALL CHKXER( 'ctpttr', INFOT, NOUT, LERR, OK )
244 INFOT = 5
245 CALL CTPTTR( 'u', 0, A, B, 0, INFO )
246 CALL CHKXER( 'ctpttr', INFOT, NOUT, LERR, OK )
247*
248 SRNAMT = 'chfrk '
249 INFOT = 1
250 CALL CHFRK( '/', 'u', 'n', 0, 0, ALPHA, A, 1, BETA, B )
251 CALL CHKXER( 'chfrk ', INFOT, NOUT, LERR, OK )
252 INFOT = 2
253 CALL CHFRK( 'n', '/', 'n', 0, 0, ALPHA, A, 1, BETA, B )
254 CALL CHKXER( 'chfrk ', INFOT, NOUT, LERR, OK )
255 INFOT = 3
256 CALL CHFRK( 'n', 'u', '/', 0, 0, ALPHA, A, 1, BETA, B )
257 CALL CHKXER( 'chfrk ', INFOT, NOUT, LERR, OK )
258 INFOT = 4
259 CALL CHFRK( 'n', 'u', 'n', -1, 0, ALPHA, A, 1, BETA, B )
260 CALL CHKXER( 'chfrk ', INFOT, NOUT, LERR, OK )
261 INFOT = 5
262 CALL CHFRK( 'n', 'u', 'n', 0, -1, ALPHA, A, 1, BETA, B )
263 CALL CHKXER( 'chfrk ', INFOT, NOUT, LERR, OK )
264 INFOT = 8
265 CALL CHFRK( 'n', 'u', 'n', 0, 0, ALPHA, A, 0, BETA, B )
266 CALL CHKXER( 'chfrk ', INFOT, NOUT, LERR, OK )
267*
268* Print a summary line.
269*
270 IF( OK ) THEN
271 WRITE( NOUT, FMT = 9999 )
272 ELSE
273 WRITE( NOUT, FMT = 9998 )
274 END IF
275*
276 9999 FORMAT( 1X, 'COMPLEX RFP routines passed the tests of the ',
277 $ 'error exits' )
278 9998 FORMAT( ' *** RFP routines failed the tests of the error ',
279 $ 'exits ***' )
280 RETURN
281*
282* End of CERRRFP
283*
284 END
cmplx
float cmplx[2]
Definition pblas.h:136
chkxer
subroutine chkxer(srnamt, infot, nout, lerr, ok)
Definition cblat2.f:3196
ctfttr
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...
Definition ctfttr.f:216
ctpttr
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...
Definition ctpttr.f:104
cpftrs
subroutine cpftrs(transr, uplo, n, nrhs, a, b, ldb, info)
CPFTRS
Definition cpftrs.f:220
cpftrf
subroutine cpftrf(transr, uplo, n, a, info)
CPFTRF
Definition cpftrf.f:211
chfrk
subroutine chfrk(transr, uplo, trans, n, k, alpha, a, lda, beta, c)
CHFRK performs a Hermitian rank-k operation for matrix in RFP format.
Definition chfrk.f:168
ctpttf
subroutine ctpttf(transr, uplo, n, ap, arf, info)
CTPTTF copies a triangular matrix from the standard packed format (TP) to the rectangular full packed...
Definition ctpttf.f:207
ctrttf
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...
Definition ctrttf.f:216
cpftri
subroutine cpftri(transr, uplo, n, a, info)
CPFTRI
Definition cpftri.f:212
ctftri
subroutine ctftri(transr, uplo, diag, n, a, info)
CTFTRI
Definition ctftri.f:221
ctfsm
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).
Definition ctfsm.f:298
ctfttp
subroutine ctfttp(transr, uplo, n, arf, ap, info)
CTFTTP copies a triangular matrix from the rectangular full packed format (TF) to the standard packed...
Definition ctfttp.f:208
ctrttp
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...
Definition ctrttp.f:104
cerrrfp
subroutine cerrrfp(nunit)
CERRRFP
Definition cerrrfp.f:52