OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
zerrrfp.f
Go to the documentation of this file.
1*> \brief \b ZERRRFP
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 ZERRRFP( NUNIT )
12*
13* .. Scalar Arguments ..
14* INTEGER NUNIT
15* ..
16*
17*
18*> \par Purpose:
19* =============
20*>
21*> \verbatim
22*>
23*> ZERRRFP tests the error exits for the COMPLEX*16 driver routines
24*> for solving linear systems of equations.
25*>
26*> ZDRVRFP tests the COMPLEX*16 LAPACK RFP routines:
27*> ZTFSM, ZTFTRI, ZHFRK, ZTFTTP, ZTFTTR, ZPFTRF, ZPFTRS, ZTPTTF,
28*> ZTPTTR, ZTRTTF, and ZTRTTP
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 complex16_lin
49*
50* =====================================================================
51 SUBROUTINE zerrrfp( 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 DOUBLE PRECISION ALPHA, BETA
67 COMPLEX*16 CALPHA
68* ..
69* .. Local Arrays ..
70 COMPLEX*16 A( 1, 1), B( 1, 1)
71* ..
72* .. External Subroutines ..
73 EXTERNAL chkxer, ztfsm, ztftri, zhfrk, ztfttp, ztfttr,
74 + zpftri, zpftrf, zpftrs, ztpttf, ztpttr, ztrttf,
75 + ztrttp
76* ..
77* .. Scalars in Common ..
78 LOGICAL LERR, OK
79 CHARACTER*32 SRNAMT
80 INTEGER INFOT, NOUT
81* ..
82* .. Intrinsic Functions ..
83 INTRINSIC dcmplx
84* ..
85* .. Common blocks ..
86 COMMON / infoc / infot, nout, ok, lerr
87 COMMON / srnamc / srnamt
88* ..
89* .. Executable Statements ..
90*
91 nout = nunit
92 ok = .true.
93 a( 1, 1 ) = dcmplx( 1.0d0 , 1.0d0 )
94 b( 1, 1 ) = dcmplx( 1.0d0 , 1.0d0 )
95 alpha = 1.0d0
96 calpha = dcmplx( 1.0d0 , 1.0d0 )
97 beta = 1.0d0
98*
99 srnamt = 'ZPFTRF'
100 infot = 1
101 CALL zpftrf( '/', 'U', 0, a, info )
102 CALL chkxer( 'ZPFTRF', infot, nout, lerr, ok )
103 infot = 2
104 CALL zpftrf( 'N', '/', 0, a, info )
105 CALL chkxer( 'ZPFTRF', infot, nout, lerr, ok )
106 infot = 3
107 CALL zpftrf( 'N', 'U', -1, a, info )
108 CALL chkxer( 'ZPFTRF', infot, nout, lerr, ok )
109*
110 srnamt = 'ZPFTRS'
111 infot = 1
112 CALL zpftrs( '/', 'U', 0, 0, a, b, 1, info )
113 CALL chkxer( 'ZPFTRS', infot, nout, lerr, ok )
114 infot = 2
115 CALL zpftrs( 'N', '/', 0, 0, a, b, 1, info )
116 CALL chkxer( 'ZPFTRS', infot, nout, lerr, ok )
117 infot = 3
118 CALL zpftrs( 'N', 'U', -1, 0, a, b, 1, info )
119 CALL chkxer( 'ZPFTRS', infot, nout, lerr, ok )
120 infot = 4
121 CALL zpftrs( 'N', 'U', 0, -1, a, b, 1, info )
122 CALL chkxer( 'ZPFTRS', infot, nout, lerr, ok )
123 infot = 7
124 CALL zpftrs( 'N', 'U', 0, 0, a, b, 0, info )
125 CALL chkxer( 'ZPFTRS', infot, nout, lerr, ok )
126*
127 srnamt = 'ZPFTRI'
128 infot = 1
129 CALL zpftri( '/', 'U', 0, a, info )
130 CALL chkxer( 'ZPFTRI', infot, nout, lerr, ok )
131 infot = 2
132 CALL zpftri( 'N', '/', 0, a, info )
133 CALL chkxer( 'ZPFTRI', infot, nout, lerr, ok )
134 infot = 3
135 CALL zpftri( 'N', 'U', -1, a, info )
136 CALL chkxer( 'ZPFTRI', infot, nout, lerr, ok )
137*
138 srnamt = 'ZTFSM '
139 infot = 1
140 CALL ztfsm( '/', 'L', 'U', 'C', 'U', 0, 0, calpha, a, b, 1 )
141 CALL chkxer( 'ZTFSM ', infot, nout, lerr, ok )
142 infot = 2
143 CALL ztfsm( 'N', '/', 'U', 'C', 'U', 0, 0, calpha, a, b, 1 )
144 CALL chkxer( 'ZTFSM ', infot, nout, lerr, ok )
145 infot = 3
146 CALL ztfsm( 'n', 'l', '/', 'c', 'u', 0, 0, CALPHA, A, B, 1 )
147 CALL CHKXER( 'ztfsm ', INFOT, NOUT, LERR, OK )
148 INFOT = 4
149 CALL ZTFSM( 'n', 'l', 'u', '/', 'u', 0, 0, CALPHA, A, B, 1 )
150 CALL CHKXER( 'ztfsm ', INFOT, NOUT, LERR, OK )
151 INFOT = 5
152 CALL ZTFSM( 'n', 'l', 'u', 'c', '/', 0, 0, CALPHA, A, B, 1 )
153 CALL CHKXER( 'ztfsm ', INFOT, NOUT, LERR, OK )
154 INFOT = 6
155 CALL ZTFSM( 'n', 'l', 'u', 'c', 'u', -1, 0, CALPHA, A, B, 1 )
156 CALL CHKXER( 'ztfsm ', INFOT, NOUT, LERR, OK )
157 INFOT = 7
158 CALL ZTFSM( 'n', 'l', 'u', 'c', 'u', 0, -1, CALPHA, A, B, 1 )
159 CALL CHKXER( 'ztfsm ', INFOT, NOUT, LERR, OK )
160 INFOT = 11
161 CALL ZTFSM( 'n', 'l', 'u', 'c', 'u', 0, 0, CALPHA, A, B, 0 )
162 CALL CHKXER( 'ztfsm ', INFOT, NOUT, LERR, OK )
163*
164 SRNAMT = 'ztftri'
165 INFOT = 1
166 CALL ZTFTRI( '/', 'l', 'n', 0, A, INFO )
167 CALL CHKXER( 'ztftri', INFOT, NOUT, LERR, OK )
168 INFOT = 2
169 CALL ZTFTRI( 'n', '/', 'n', 0, A, INFO )
170 CALL CHKXER( 'ztftri', INFOT, NOUT, LERR, OK )
171 INFOT = 3
172 CALL ZTFTRI( 'n', 'l', '/', 0, A, INFO )
173 CALL CHKXER( 'ztftri', INFOT, NOUT, LERR, OK )
174 INFOT = 4
175 CALL ZTFTRI( 'n', 'l', 'n', -1, A, INFO )
176 CALL CHKXER( 'ztftri', INFOT, NOUT, LERR, OK )
177*
178 SRNAMT = 'ztfttr'
179 INFOT = 1
180 CALL ZTFTTR( '/', 'u', 0, A, B, 1, INFO )
181 CALL CHKXER( 'ztfttr', INFOT, NOUT, LERR, OK )
182 INFOT = 2
183 CALL ZTFTTR( 'n', '/', 0, A, B, 1, INFO )
184 CALL CHKXER( 'ztfttr', INFOT, NOUT, LERR, OK )
185 INFOT = 3
186 CALL ZTFTTR( 'n', 'u', -1, A, B, 1, INFO )
187 CALL CHKXER( 'ztfttr', INFOT, NOUT, LERR, OK )
188 INFOT = 6
189 CALL ZTFTTR( 'n', 'u', 0, A, B, 0, INFO )
190 CALL CHKXER( 'ztfttr', INFOT, NOUT, LERR, OK )
191*
192 SRNAMT = 'ztrttf'
193 INFOT = 1
194 CALL ZTRTTF( '/', 'u', 0, A, 1, B, INFO )
195 CALL CHKXER( 'ztrttf', INFOT, NOUT, LERR, OK )
196 INFOT = 2
197 CALL ZTRTTF( 'n', '/', 0, A, 1, B, INFO )
198 CALL CHKXER( 'ztrttf', INFOT, NOUT, LERR, OK )
199 INFOT = 3
200 CALL ZTRTTF( 'n', 'u', -1, A, 1, B, INFO )
201 CALL CHKXER( 'ztrttf', INFOT, NOUT, LERR, OK )
202 INFOT = 5
203 CALL ZTRTTF( 'n', 'u', 0, A, 0, B, INFO )
204 CALL CHKXER( 'ztrttf', INFOT, NOUT, LERR, OK )
205*
206 SRNAMT = 'ztfttp'
207 INFOT = 1
208 CALL ZTFTTP( '/', 'u', 0, A, B, INFO )
209 CALL CHKXER( 'ztfttp', INFOT, NOUT, LERR, OK )
210 INFOT = 2
211 CALL ZTFTTP( 'n', '/', 0, A, B, INFO )
212 CALL CHKXER( 'ztfttp', infot, nout, lerr, ok )
213 infot = 3
214 CALL ztfttp( 'N', 'U', -1, a, b, info )
215 CALL chkxer( 'ZTFTTP', infot, nout, lerr, ok )
216*
217 srnamt = 'ZTPTTF'
218 infot = 1
219 CALL ztpttf( '/', 'U', 0, a, b, info )
220 CALL chkxer( 'ZTPTTF', infot, nout, lerr, ok )
221 infot = 2
222 CALL ztpttf( 'N', '/', 0, a, b, info )
223 CALL chkxer( 'ZTPTTF', infot, nout, lerr, ok )
224 infot = 3
225 CALL ztpttf( 'N', 'U', -1, a, b, info )
226 CALL chkxer( 'ZTPTTF', infot, nout, lerr, ok )
227*
228 srnamt = 'ZTRTTP'
229 infot = 1
230 CALL ztrttp( '/', 0, a, 1, b, info )
231 CALL chkxer( 'ZTRTTP', infot, nout, lerr, ok )
232 infot = 2
233 CALL ztrttp( 'U', -1, a, 1, b, info )
234 CALL chkxer( 'ZTRTTP', infot, nout, lerr, ok )
235 infot = 4
236 CALL ztrttp( 'U', 0, a, 0, b, info )
237 CALL chkxer( 'ZTRTTP', infot, nout, lerr, ok )
238*
239 srnamt = 'ZTPTTR'
240 infot = 1
241 CALL ztpttr( '/', 0, a, b, 1, info )
242 CALL chkxer( 'ZTPTTR', infot, nout, lerr, ok )
243 infot = 2
244 CALL ztpttr( 'U', -1, a, b, 1, info )
245 CALL chkxer( 'ZTPTTR', infot, nout, lerr, ok )
246 infot = 5
247 CALL ztpttr( 'U', 0, a, b, 0, info )
248 CALL chkxer( 'ZTPTTR', infot, nout, lerr, ok )
249*
250 srnamt = 'ZHFRK '
251 infot = 1
252 CALL zhfrk( '/', 'U', 'N', 0, 0, alpha, a, 1, beta, b )
253 CALL chkxer( 'ZHFRK ', infot, nout, lerr, ok )
254 infot = 2
255 CALL zhfrk( 'N', '/', 'N', 0, 0, alpha, a, 1, beta, b )
256 CALL chkxer( 'ZHFRK ', infot, nout, lerr, ok )
257 infot = 3
258 CALL zhfrk( 'N', 'U', '/', 0, 0, alpha, a, 1, beta, b )
259 CALL chkxer( 'ZHFRK ', infot, nout, lerr, ok )
260 infot = 4
261 CALL zhfrk( 'N', 'U', 'N', -1, 0, alpha, a, 1, beta, b )
262 CALL chkxer( 'ZHFRK ', infot, nout, lerr, ok )
263 infot = 5
264 CALL zhfrk( 'N', 'U', 'N', 0, -1, alpha, a, 1, beta, b )
265 CALL chkxer( 'ZHFRK ', infot, nout, lerr, ok )
266 infot = 8
267 CALL zhfrk( 'N', 'U', 'N', 0, 0, alpha, a, 0, beta, b )
268 CALL chkxer( 'ZHFRK ', infot, nout, lerr, ok )
269*
270* Print a summary line.
271*
272 IF( ok ) THEN
273 WRITE( nout, fmt = 9999 )
274 ELSE
275 WRITE( nout, fmt = 9998 )
276 END IF
277*
278 9999 FORMAT( 1x, 'COMPLEX*16 RFP routines passed the tests of the ',
279 $ 'error exits' )
280 9998 FORMAT( ' *** RFP routines failed the tests of the error ',
281 $ 'exits ***' )
282 RETURN
283*
284* End of ZERRRFP
285*
286 END
chkxer
subroutine chkxer(srnamt, infot, nout, lerr, ok)
Definition cblat2.f:3196
ztpttr
subroutine ztpttr(uplo, n, ap, a, lda, info)
ZTPTTR copies a triangular matrix from the standard packed format (TP) to the standard full format (T...
Definition ztpttr.f:104
ztfsm
subroutine ztfsm(transr, side, uplo, trans, diag, m, n, alpha, a, b, ldb)
ZTFSM solves a matrix equation (one operand is a triangular matrix in RFP format).
Definition ztfsm.f:298
ztftri
subroutine ztftri(transr, uplo, diag, n, a, info)
ZTFTRI
Definition ztftri.f:221
ztfttp
subroutine ztfttp(transr, uplo, n, arf, ap, info)
ZTFTTP copies a triangular matrix from the rectangular full packed format (TF) to the standard packed...
Definition ztfttp.f:208
ztrttf
subroutine ztrttf(transr, uplo, n, a, lda, arf, info)
ZTRTTF copies a triangular matrix from the standard full format (TR) to the rectangular full packed f...
Definition ztrttf.f:216
zpftri
subroutine zpftri(transr, uplo, n, a, info)
ZPFTRI
Definition zpftri.f:212
ztrttp
subroutine ztrttp(uplo, n, a, lda, ap, info)
ZTRTTP copies a triangular matrix from the standard full format (TR) to the standard packed format (T...
Definition ztrttp.f:104
zpftrs
subroutine zpftrs(transr, uplo, n, nrhs, a, b, ldb, info)
ZPFTRS
Definition zpftrs.f:220
ztpttf
subroutine ztpttf(transr, uplo, n, ap, arf, info)
ZTPTTF copies a triangular matrix from the standard packed format (TP) to the rectangular full packed...
Definition ztpttf.f:207
zpftrf
subroutine zpftrf(transr, uplo, n, a, info)
ZPFTRF
Definition zpftrf.f:211
ztfttr
subroutine ztfttr(transr, uplo, n, arf, a, lda, info)
ZTFTTR copies a triangular matrix from the rectangular full packed format (TF) to the standard full f...
Definition ztfttr.f:216
zhfrk
subroutine zhfrk(transr, uplo, trans, n, k, alpha, a, lda, beta, c)
ZHFRK performs a Hermitian rank-k operation for matrix in RFP format.
Definition zhfrk.f:168
zerrrfp
subroutine zerrrfp(nunit)
ZERRRFP
Definition zerrrfp.f:52