OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
zerrpox.f
Go to the documentation of this file.
1*> \brief \b ZERRPOX
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 ZERRPO( PATH, NUNIT )
12*
13* .. Scalar Arguments ..
14* CHARACTER*3 PATH
15* INTEGER NUNIT
16* ..
17*
18*
19*> \par Purpose:
20* =============
21*>
22*> \verbatim
23*>
24*> ZERRPO tests the error exits for the COMPLEX*16 routines
25*> for Hermitian positive definite matrices.
26*>
27*> Note that this file is used only when the XBLAS are available,
28*> otherwise zerrpo.f defines this subroutine.
29*> \endverbatim
30*
31* Arguments:
32* ==========
33*
34*> \param[in] PATH
35*> \verbatim
36*> PATH is CHARACTER*3
37*> The LAPACK path name for the routines to be tested.
38*> \endverbatim
39*>
40*> \param[in] NUNIT
41*> \verbatim
42*> NUNIT is INTEGER
43*> The unit number for output.
44*> \endverbatim
45*
46* Authors:
47* ========
48*
49*> \author Univ. of Tennessee
50*> \author Univ. of California Berkeley
51*> \author Univ. of Colorado Denver
52*> \author NAG Ltd.
53*
54*> \ingroup complex16_lin
55*
56* =====================================================================
57 SUBROUTINE zerrpo( PATH, NUNIT )
58*
59* -- LAPACK test routine --
60* -- LAPACK is a software package provided by Univ. of Tennessee, --
61* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
62*
63* .. Scalar Arguments ..
64 CHARACTER*3 PATH
65 INTEGER NUNIT
66* ..
67*
68* =====================================================================
69*
70* .. Parameters ..
71 INTEGER NMAX
72 parameter( nmax = 4 )
73* ..
74* .. Local Scalars ..
75 CHARACTER EQ
76 CHARACTER*2 C2
77 INTEGER I, INFO, J, N_ERR_BNDS, NPARAMS
78 DOUBLE PRECISION ANRM, RCOND, BERR
79* ..
80* .. Local Arrays ..
81 DOUBLE PRECISION S( NMAX ), R( NMAX ), R1( NMAX ), R2( NMAX ),
82 $ ERR_BNDS_N( NMAX, 3 ), ERR_BNDS_C( NMAX, 3 ),
83 $ PARAMS( 1 )
84 COMPLEX*16 A( NMAX, NMAX ), AF( NMAX, NMAX ), B( NMAX ),
85 $ W( 2*NMAX ), X( NMAX )
86* ..
87* .. External Functions ..
88 LOGICAL LSAMEN
89 EXTERNAL lsamen
90* ..
91* .. External Subroutines ..
92 EXTERNAL alaesm, chkxer, zpbcon, zpbequ, zpbrfs, zpbtf2,
93 $ zpbtrf, zpbtrs, zpocon, zpoequ, zporfs, zpotf2,
94 $ zpotrf, zpotri, zpotrs, zppcon, zppequ, zpprfs,
95 $ zpptrf, zpptri, zpptrs, zpoequb, zporfsx
96* ..
97* .. Scalars in Common ..
98 LOGICAL LERR, OK
99 CHARACTER*32 SRNAMT
100 INTEGER INFOT, NOUT
101* ..
102* .. Common blocks ..
103 COMMON / infoc / infot, nout, ok, lerr
104 COMMON / srnamc / srnamt
105* ..
106* .. Intrinsic Functions ..
107 INTRINSIC dble, dcmplx
108* ..
109* .. Executable Statements ..
110*
111 nout = nunit
112 WRITE( nout, fmt = * )
113 c2 = path( 2: 3 )
114*
115* Set the variables to innocuous values.
116*
117 DO 20 j = 1, nmax
118 DO 10 i = 1, nmax
119 a( i, j ) = dcmplx( 1.d0 / dble( i+j ),
120 $ -1.d0 / dble( i+j ) )
121 af( i, j ) = dcmplx( 1.d0 / dble( i+j ),
122 $ -1.d0 / dble( i+j ) )
123 10 CONTINUE
124 b( j ) = 0.d0
125 r1( j ) = 0.d0
126 r2( j ) = 0.d0
127 w( j ) = 0.d0
128 x( j ) = 0.d0
129 s( j ) = 0.d0
130 20 CONTINUE
131 anrm = 1.d0
132 ok = .true.
133*
134* Test error exits of the routines that use the Cholesky
135* decomposition of a Hermitian positive definite matrix.
136*
137 IF( lsamen( 2, c2, 'po' ) ) THEN
138*
139* ZPOTRF
140*
141 SRNAMT = 'zpotrf'
142 INFOT = 1
143 CALL ZPOTRF( '/', 0, A, 1, INFO )
144 CALL CHKXER( 'zpotrf', INFOT, NOUT, LERR, OK )
145 INFOT = 2
146 CALL ZPOTRF( 'u', -1, A, 1, INFO )
147 CALL CHKXER( 'zpotrf', INFOT, NOUT, LERR, OK )
148 INFOT = 4
149 CALL ZPOTRF( 'u', 2, A, 1, INFO )
150 CALL CHKXER( 'zpotrf', INFOT, NOUT, LERR, OK )
151*
152* ZPOTF2
153*
154 SRNAMT = 'zpotf2'
155 INFOT = 1
156 CALL ZPOTF2( '/', 0, A, 1, INFO )
157 CALL CHKXER( 'zpotf2', INFOT, NOUT, LERR, OK )
158 INFOT = 2
159 CALL ZPOTF2( 'u', -1, A, 1, INFO )
160 CALL CHKXER( 'zpotf2', INFOT, NOUT, LERR, OK )
161 INFOT = 4
162 CALL ZPOTF2( 'u', 2, A, 1, INFO )
163 CALL CHKXER( 'zpotf2', INFOT, NOUT, LERR, OK )
164*
165* ZPOTRI
166*
167 SRNAMT = 'zpotri'
168 INFOT = 1
169 CALL ZPOTRI( '/', 0, A, 1, INFO )
170 CALL CHKXER( 'zpotri', INFOT, NOUT, LERR, OK )
171 INFOT = 2
172 CALL ZPOTRI( 'u', -1, A, 1, INFO )
173 CALL CHKXER( 'zpotri', INFOT, NOUT, LERR, OK )
174 INFOT = 4
175 CALL ZPOTRI( 'u', 2, A, 1, INFO )
176 CALL CHKXER( 'zpotri', INFOT, NOUT, LERR, OK )
177*
178* ZPOTRS
179*
180 SRNAMT = 'zpotrs'
181 INFOT = 1
182 CALL ZPOTRS( '/', 0, 0, A, 1, B, 1, INFO )
183 CALL CHKXER( 'zpotrs', INFOT, NOUT, LERR, OK )
184 INFOT = 2
185 CALL ZPOTRS( 'u', -1, 0, A, 1, B, 1, INFO )
186 CALL CHKXER( 'zpotrs', INFOT, NOUT, LERR, OK )
187 INFOT = 3
188 CALL ZPOTRS( 'u', 0, -1, A, 1, B, 1, INFO )
189 CALL CHKXER( 'zpotrs', INFOT, NOUT, LERR, OK )
190 INFOT = 5
191 CALL ZPOTRS( 'u', 2, 1, A, 1, B, 2, INFO )
192 CALL CHKXER( 'zpotrs', INFOT, NOUT, LERR, OK )
193 INFOT = 7
194 CALL ZPOTRS( 'u', 2, 1, A, 2, B, 1, INFO )
195 CALL CHKXER( 'zpotrs', INFOT, NOUT, LERR, OK )
196*
197* ZPORFS
198*
199 SRNAMT = 'zporfs'
200 INFOT = 1
201 CALL ZPORFS( '/', 0, 0, A, 1, AF, 1, B, 1, X, 1, R1, R2, W, R,
202 $ INFO )
203 CALL CHKXER( 'zporfs', INFOT, NOUT, LERR, OK )
204 INFOT = 2
205 CALL ZPORFS( 'u', -1, 0, A, 1, AF, 1, B, 1, X, 1, R1, R2, W, R,
206 $ INFO )
207 CALL CHKXER( 'zporfs', INFOT, NOUT, LERR, OK )
208 INFOT = 3
209 CALL ZPORFS( 'u', 0, -1, A, 1, AF, 1, B, 1, X, 1, R1, R2, W, R,
210 $ INFO )
211 CALL CHKXER( 'zporfs', INFOT, NOUT, LERR, OK )
212 INFOT = 5
213 CALL ZPORFS( 'u', 2, 1, A, 1, AF, 2, B, 2, X, 2, R1, R2, W, R,
214 $ INFO )
215 CALL CHKXER( 'zporfs', INFOT, NOUT, LERR, OK )
216 INFOT = 7
217 CALL ZPORFS( 'u', 2, 1, A, 2, AF, 1, B, 2, X, 2, R1, R2, W, R,
218 $ INFO )
219 CALL CHKXER( 'zporfs', INFOT, NOUT, LERR, OK )
220 INFOT = 9
221 CALL ZPORFS( 'u', 2, 1, A, 2, AF, 2, B, 1, X, 2, R1, R2, W, R,
222 $ INFO )
223 CALL CHKXER( 'zporfs', INFOT, NOUT, LERR, OK )
224 INFOT = 11
225 CALL ZPORFS( 'u', 2, 1, A, 2, AF, 2, B, 2, X, 1, R1, R2, W, R,
226 $ INFO )
227 CALL CHKXER( 'zporfs', INFOT, NOUT, LERR, OK )
228*
229* ZPORFSX
230*
231 N_ERR_BNDS = 3
232 NPARAMS = 0
233 SRNAMT = 'zporfsx'
234 INFOT = 1
235 CALL ZPORFSX( '/', EQ, 0, 0, A, 1, AF, 1, S, B, 1, X, 1,
236 $ RCOND, BERR, N_ERR_BNDS, ERR_BNDS_N, ERR_BNDS_C, NPARAMS,
237 $ PARAMS, W, R, INFO )
238 CALL CHKXER( 'zporfsx', INFOT, NOUT, LERR, OK )
239 INFOT = 2
240 CALL ZPORFSX( 'u', "/", -1, 0, A, 1, AF, 1, S, B, 1, X, 1,
241 $ RCOND, BERR, N_ERR_BNDS, ERR_BNDS_N, ERR_BNDS_C, NPARAMS,
242 $ PARAMS, W, R, INFO )
243 CALL CHKXER( 'zporfsx', INFOT, NOUT, LERR, OK )
244 EQ = 'n'
245 INFOT = 3
246 CALL ZPORFSX( 'u', EQ, -1, 0, A, 1, AF, 1, S, B, 1, X, 1,
247 $ RCOND, BERR, N_ERR_BNDS, ERR_BNDS_N, ERR_BNDS_C, NPARAMS,
248 $ PARAMS, W, R, INFO )
249 CALL CHKXER( 'zporfsx', INFOT, NOUT, LERR, OK )
250 INFOT = 4
251 CALL ZPORFSX( 'u', EQ, 0, -1, A, 1, AF, 1, S, B, 1, X, 1,
252 $ RCOND, BERR, N_ERR_BNDS, ERR_BNDS_N, ERR_BNDS_C, NPARAMS,
253 $ PARAMS, W, R, INFO )
254 CALL CHKXER( 'zporfsx', INFOT, NOUT, LERR, OK )
255 INFOT = 6
256 CALL ZPORFSX( 'u', EQ, 2, 1, A, 1, AF, 2, S, B, 2, X, 2,
257 $ RCOND, BERR, N_ERR_BNDS, ERR_BNDS_N, ERR_BNDS_C, NPARAMS,
258 $ PARAMS, W, R, INFO )
259 CALL CHKXER( 'zporfsx', INFOT, NOUT, LERR, OK )
260 INFOT = 8
261 CALL ZPORFSX( 'u', EQ, 2, 1, A, 2, AF, 1, S, B, 2, X, 2,
262 $ RCOND, BERR, N_ERR_BNDS, ERR_BNDS_N, ERR_BNDS_C, NPARAMS,
263 $ PARAMS, W, R, INFO )
264 CALL CHKXER( 'zporfsx', INFOT, NOUT, LERR, OK )
265 INFOT = 11
266 CALL ZPORFSX( 'u', EQ, 2, 1, A, 2, AF, 2, S, B, 1, X, 2,
267 $ RCOND, BERR, N_ERR_BNDS, ERR_BNDS_N, ERR_BNDS_C, NPARAMS,
268 $ PARAMS, W, R, INFO )
269 CALL CHKXER( 'zporfsx', INFOT, NOUT, LERR, OK )
270 INFOT = 13
271 CALL ZPORFSX( 'u', EQ, 2, 1, A, 2, AF, 2, S, B, 2, X, 1,
272 $ RCOND, BERR, N_ERR_BNDS, ERR_BNDS_N, ERR_BNDS_C, NPARAMS,
273 $ PARAMS, W, R, INFO )
274 CALL CHKXER( 'zporfsx', INFOT, NOUT, LERR, OK )
275*
276* ZPOCON
277*
278 SRNAMT = 'zpocon'
279 INFOT = 1
280 CALL ZPOCON( '/', 0, A, 1, ANRM, RCOND, W, R, INFO )
281 CALL CHKXER( 'zpocon', INFOT, NOUT, LERR, OK )
282 INFOT = 2
283 CALL ZPOCON( 'u', -1, A, 1, ANRM, RCOND, W, R, INFO )
284 CALL CHKXER( 'zpocon', INFOT, NOUT, LERR, OK )
285 INFOT = 4
286 CALL ZPOCON( 'u', 2, A, 1, ANRM, RCOND, W, R, INFO )
287 CALL CHKXER( 'zpocon', INFOT, NOUT, LERR, OK )
288 INFOT = 5
289 CALL ZPOCON( 'u', 1, A, 1, -ANRM, RCOND, W, R, INFO )
290 CALL CHKXER( 'zpocon', INFOT, NOUT, LERR, OK )
291*
292* ZPOEQU
293*
294 SRNAMT = 'zpoequ'
295 INFOT = 1
296 CALL ZPOEQU( -1, A, 1, R1, RCOND, ANRM, INFO )
297 CALL CHKXER( 'zpoequ', INFOT, NOUT, LERR, OK )
298 INFOT = 3
299 CALL ZPOEQU( 2, A, 1, R1, RCOND, ANRM, INFO )
300 CALL CHKXER( 'zpoequ', INFOT, NOUT, LERR, OK )
301*
302* ZPOEQUB
303*
304 SRNAMT = 'zpoequb'
305 INFOT = 1
306 CALL ZPOEQUB( -1, A, 1, R1, RCOND, ANRM, INFO )
307 CALL CHKXER( 'zpoequb', INFOT, NOUT, LERR, OK )
308 INFOT = 3
309 CALL ZPOEQUB( 2, A, 1, R1, RCOND, ANRM, INFO )
310 CALL CHKXER( 'zpoequb', INFOT, NOUT, LERR, OK )
311*
312* Test error exits of the routines that use the Cholesky
313* decomposition of a Hermitian positive definite packed matrix.
314*
315 ELSE IF( LSAMEN( 2, C2, 'pp' ) ) THEN
316*
317* ZPPTRF
318*
319 SRNAMT = 'zpptrf'
320 INFOT = 1
321 CALL ZPPTRF( '/', 0, A, INFO )
322 CALL CHKXER( 'zpptrf', INFOT, NOUT, LERR, OK )
323 INFOT = 2
324 CALL ZPPTRF( 'u', -1, A, INFO )
325 CALL CHKXER( 'zpptrf', INFOT, NOUT, LERR, OK )
326*
327* ZPPTRI
328*
329 SRNAMT = 'zpptri'
330 INFOT = 1
331 CALL ZPPTRI( '/', 0, A, INFO )
332 CALL CHKXER( 'zpptri', INFOT, NOUT, LERR, OK )
333 INFOT = 2
334 CALL ZPPTRI( 'u', -1, A, INFO )
335 CALL CHKXER( 'zpptri', INFOT, NOUT, LERR, OK )
336*
337* ZPPTRS
338*
339 SRNAMT = 'zpptrs'
340 INFOT = 1
341 CALL ZPPTRS( '/', 0, 0, A, B, 1, INFO )
342 CALL CHKXER( 'zpptrs', INFOT, NOUT, LERR, OK )
343 INFOT = 2
344 CALL ZPPTRS( 'u', -1, 0, A, B, 1, INFO )
345 CALL CHKXER( 'zpptrs', INFOT, NOUT, LERR, OK )
346 INFOT = 3
347 CALL ZPPTRS( 'u', 0, -1, A, B, 1, INFO )
348 CALL CHKXER( 'zpptrs', INFOT, NOUT, LERR, OK )
349 INFOT = 6
350 CALL ZPPTRS( 'u', 2, 1, A, B, 1, INFO )
351 CALL CHKXER( 'zpptrs', INFOT, NOUT, LERR, OK )
352*
353* ZPPRFS
354*
355 SRNAMT = 'zpprfs'
356 INFOT = 1
357 CALL ZPPRFS( '/', 0, 0, A, AF, B, 1, X, 1, R1, R2, W, R, INFO )
358 CALL CHKXER( 'zpprfs', INFOT, NOUT, LERR, OK )
359 INFOT = 2
360 CALL ZPPRFS( 'u', -1, 0, A, AF, B, 1, X, 1, R1, R2, W, R,
361 $ INFO )
362 CALL CHKXER( 'zpprfs', INFOT, NOUT, LERR, OK )
363 INFOT = 3
364 CALL ZPPRFS( 'u', 0, -1, A, AF, B, 1, X, 1, R1, R2, W, R,
365 $ INFO )
366 CALL CHKXER( 'zpprfs', INFOT, NOUT, LERR, OK )
367 INFOT = 7
368 CALL ZPPRFS( 'u', 2, 1, A, AF, B, 1, X, 2, R1, R2, W, R, INFO )
369 CALL CHKXER( 'zpprfs', INFOT, NOUT, LERR, OK )
370 INFOT = 9
371 CALL ZPPRFS( 'u', 2, 1, A, AF, B, 2, X, 1, R1, R2, W, R, INFO )
372 CALL CHKXER( 'zpprfs', INFOT, NOUT, LERR, OK )
373*
374* ZPPCON
375*
376 SRNAMT = 'zppcon'
377 INFOT = 1
378 CALL ZPPCON( '/', 0, A, ANRM, RCOND, W, R, INFO )
379 CALL CHKXER( 'zppcon', INFOT, NOUT, LERR, OK )
380 INFOT = 2
381 CALL ZPPCON( 'u', -1, A, ANRM, RCOND, W, R, INFO )
382 CALL CHKXER( 'zppcon', INFOT, NOUT, LERR, OK )
383 INFOT = 4
384 CALL ZPPCON( 'u', 1, A, -ANRM, RCOND, W, R, INFO )
385 CALL CHKXER( 'zppcon', INFOT, NOUT, LERR, OK )
386*
387* ZPPEQU
388*
389 SRNAMT = 'zppequ'
390 INFOT = 1
391 CALL ZPPEQU( '/', 0, A, R1, RCOND, ANRM, INFO )
392 CALL CHKXER( 'zppequ', INFOT, NOUT, LERR, OK )
393 INFOT = 2
394 CALL ZPPEQU( 'u', -1, A, R1, RCOND, ANRM, INFO )
395 CALL CHKXER( 'zppequ', INFOT, NOUT, LERR, OK )
396*
397* Test error exits of the routines that use the Cholesky
398* decomposition of a Hermitian positive definite band matrix.
399*
400 ELSE IF( LSAMEN( 2, C2, 'pb' ) ) THEN
401*
402* ZPBTRF
403*
404 SRNAMT = 'zpbtrf'
405 INFOT = 1
406 CALL ZPBTRF( '/', 0, 0, A, 1, INFO )
407 CALL CHKXER( 'zpbtrf', INFOT, NOUT, LERR, OK )
408 INFOT = 2
409 CALL ZPBTRF( 'u', -1, 0, A, 1, INFO )
410 CALL CHKXER( 'zpbtrf', INFOT, NOUT, LERR, OK )
411 INFOT = 3
412 CALL ZPBTRF( 'u', 1, -1, A, 1, INFO )
413 CALL CHKXER( 'zpbtrf', INFOT, NOUT, LERR, OK )
414 INFOT = 5
415 CALL ZPBTRF( 'u', 2, 1, A, 1, INFO )
416 CALL CHKXER( 'zpbtrf', INFOT, NOUT, LERR, OK )
417*
418* ZPBTF2
419*
420 SRNAMT = 'zpbtf2'
421 INFOT = 1
422 CALL ZPBTF2( '/', 0, 0, A, 1, INFO )
423 CALL CHKXER( 'zpbtf2', INFOT, NOUT, LERR, OK )
424 INFOT = 2
425 CALL ZPBTF2( 'u', -1, 0, A, 1, INFO )
426 CALL CHKXER( 'zpbtf2', INFOT, NOUT, LERR, OK )
427 INFOT = 3
428 CALL ZPBTF2( 'u', 1, -1, A, 1, INFO )
429 CALL CHKXER( 'zpbtf2', INFOT, NOUT, LERR, OK )
430 INFOT = 5
431 CALL ZPBTF2( 'u', 2, 1, A, 1, INFO )
432 CALL CHKXER( 'zpbtf2', INFOT, NOUT, LERR, OK )
433*
434* ZPBTRS
435*
436 SRNAMT = 'zpbtrs'
437 INFOT = 1
438 CALL ZPBTRS( '/', 0, 0, 0, A, 1, B, 1, INFO )
439 CALL CHKXER( 'zpbtrs', INFOT, NOUT, LERR, OK )
440 INFOT = 2
441 CALL ZPBTRS( 'u', -1, 0, 0, A, 1, B, 1, INFO )
442 CALL CHKXER( 'zpbtrs', INFOT, NOUT, LERR, OK )
443 INFOT = 3
444 CALL ZPBTRS( 'u', 1, -1, 0, A, 1, B, 1, INFO )
445 CALL CHKXER( 'zpbtrs', INFOT, NOUT, LERR, OK )
446 INFOT = 4
447 CALL ZPBTRS( 'u', 0, 0, -1, A, 1, B, 1, INFO )
448 CALL CHKXER( 'zpbtrs', INFOT, NOUT, LERR, OK )
449 INFOT = 6
450 CALL ZPBTRS( 'u', 2, 1, 1, A, 1, B, 1, INFO )
451 CALL CHKXER( 'zpbtrs', INFOT, NOUT, LERR, OK )
452 INFOT = 8
453 CALL ZPBTRS( 'u', 2, 0, 1, A, 1, B, 1, INFO )
454 CALL CHKXER( 'zpbtrs', INFOT, NOUT, LERR, OK )
455*
456* ZPBRFS
457*
458 SRNAMT = 'zpbrfs'
459 INFOT = 1
460 CALL ZPBRFS( '/', 0, 0, 0, A, 1, AF, 1, B, 1, X, 1, R1, R2, W,
461 $ R, INFO )
462 CALL CHKXER( 'zpbrfs', INFOT, NOUT, LERR, OK )
463 INFOT = 2
464 CALL ZPBRFS( 'u', -1, 0, 0, A, 1, AF, 1, B, 1, X, 1, R1, R2, W,
465 $ R, INFO )
466 CALL CHKXER( 'zpbrfs', INFOT, NOUT, LERR, OK )
467 INFOT = 3
468 CALL ZPBRFS( 'u', 1, -1, 0, A, 1, AF, 1, B, 1, X, 1, R1, R2, W,
469 $ R, INFO )
470 CALL CHKXER( 'zpbrfs', INFOT, NOUT, LERR, OK )
471 INFOT = 4
472 CALL ZPBRFS( 'u', 0, 0, -1, A, 1, AF, 1, B, 1, X, 1, R1, R2, W,
473 $ R, INFO )
474 CALL CHKXER( 'zpbrfs', INFOT, NOUT, LERR, OK )
475 INFOT = 6
476 CALL ZPBRFS( 'u', 2, 1, 1, A, 1, AF, 2, B, 2, X, 2, R1, R2, W,
477 $ R, INFO )
478 CALL CHKXER( 'zpbrfs', INFOT, NOUT, LERR, OK )
479 INFOT = 8
480 CALL ZPBRFS( 'u', 2, 1, 1, A, 2, AF, 1, B, 2, X, 2, R1, R2, W,
481 $ R, INFO )
482 CALL CHKXER( 'zpbrfs', INFOT, NOUT, LERR, OK )
483 INFOT = 10
484 CALL ZPBRFS( 'u', 2, 0, 1, A, 1, AF, 1, B, 1, X, 2, R1, R2, W,
485 $ R, INFO )
486 CALL CHKXER( 'zpbrfs', INFOT, NOUT, LERR, OK )
487 INFOT = 12
488 CALL ZPBRFS( 'u', 2, 0, 1, A, 1, AF, 1, B, 2, X, 1, R1, R2, W,
489 $ R, INFO )
490 CALL CHKXER( 'zpbrfs', INFOT, NOUT, LERR, OK )
491*
492* ZPBCON
493*
494 SRNAMT = 'zpbcon'
495 INFOT = 1
496 CALL ZPBCON( '/', 0, 0, A, 1, ANRM, RCOND, W, R, INFO )
497 CALL CHKXER( 'zpbcon', INFOT, NOUT, LERR, OK )
498 INFOT = 2
499 CALL ZPBCON( 'u', -1, 0, A, 1, ANRM, RCOND, W, R, INFO )
500 CALL CHKXER( 'zpbcon', INFOT, NOUT, LERR, OK )
501 INFOT = 3
502 CALL ZPBCON( 'u', 1, -1, A, 1, ANRM, RCOND, W, R, INFO )
503 CALL CHKXER( 'zpbcon', INFOT, NOUT, LERR, OK )
504 INFOT = 5
505 CALL ZPBCON( 'u', 2, 1, A, 1, ANRM, RCOND, W, R, INFO )
506 CALL CHKXER( 'zpbcon', INFOT, NOUT, LERR, OK )
507 INFOT = 6
508 CALL ZPBCON( 'u', 1, 0, A, 1, -ANRM, RCOND, W, R, INFO )
509 CALL CHKXER( 'zpbcon', INFOT, NOUT, LERR, OK )
510*
511* ZPBEQU
512*
513 SRNAMT = 'zpbequ'
514 INFOT = 1
515 CALL ZPBEQU( '/', 0, 0, A, 1, R1, RCOND, ANRM, INFO )
516 CALL CHKXER( 'zpbequ', INFOT, NOUT, LERR, OK )
517 INFOT = 2
518 CALL ZPBEQU( 'u', -1, 0, A, 1, R1, RCOND, ANRM, INFO )
519 CALL CHKXER( 'zpbequ', INFOT, NOUT, LERR, OK )
520 INFOT = 3
521 CALL ZPBEQU( 'u', 1, -1, A, 1, R1, RCOND, ANRM, INFO )
522 CALL CHKXER( 'zpbequ', INFOT, NOUT, LERR, OK )
523 INFOT = 5
524 CALL ZPBEQU( 'u', 2, 1, A, 1, R1, RCOND, ANRM, INFO )
525 CALL CHKXER( 'zpbequ', INFOT, NOUT, LERR, OK )
526 END IF
527*
528* Print a summary line.
529*
530 CALL ALAESM( PATH, OK, NOUT )
531*
532 RETURN
533*
534* End of ZERRPOX
535*
536 END
chkxer
subroutine chkxer(srnamt, infot, nout, lerr, ok)
Definition cblat2.f:3196
alaesm
subroutine alaesm(path, ok, nout)
ALAESM
Definition alaesm.f:63
zpbtrs
subroutine zpbtrs(uplo, n, kd, nrhs, ab, ldab, b, ldb, info)
ZPBTRS
Definition zpbtrs.f:121
zpbtrf
subroutine zpbtrf(uplo, n, kd, ab, ldab, info)
ZPBTRF
Definition zpbtrf.f:142
zpptrs
subroutine zpptrs(uplo, n, nrhs, ap, b, ldb, info)
ZPPTRS
Definition zpptrs.f:108
zpbcon
subroutine zpbcon(uplo, n, kd, ab, ldab, anorm, rcond, work, rwork, info)
ZPBCON
Definition zpbcon.f:133
zppequ
subroutine zppequ(uplo, n, ap, s, scond, amax, info)
ZPPEQU
Definition zppequ.f:117
zpbequ
subroutine zpbequ(uplo, n, kd, ab, ldab, s, scond, amax, info)
ZPBEQU
Definition zpbequ.f:130
zpprfs
subroutine zpprfs(uplo, n, nrhs, ap, afp, b, ldb, x, ldx, ferr, berr, work, rwork, info)
ZPPRFS
Definition zpprfs.f:171
zpbtf2
subroutine zpbtf2(uplo, n, kd, ab, ldab, info)
ZPBTF2 computes the Cholesky factorization of a symmetric/Hermitian positive definite band matrix (un...
Definition zpbtf2.f:142
zppcon
subroutine zppcon(uplo, n, ap, anorm, rcond, work, rwork, info)
ZPPCON
Definition zppcon.f:118
zpptri
subroutine zpptri(uplo, n, ap, info)
ZPPTRI
Definition zpptri.f:93
zpbrfs
subroutine zpbrfs(uplo, n, kd, nrhs, ab, ldab, afb, ldafb, b, ldb, x, ldx, ferr, berr, work, rwork, info)
ZPBRFS
Definition zpbrfs.f:189
zpptrf
subroutine zpptrf(uplo, n, ap, info)
ZPPTRF
Definition zpptrf.f:119
zpotf2
subroutine zpotf2(uplo, n, a, lda, info)
ZPOTF2 computes the Cholesky factorization of a symmetric/Hermitian positive definite matrix (unblock...
Definition zpotf2.f:109
zporfsx
subroutine zporfsx(uplo, equed, n, nrhs, a, lda, af, ldaf, s, b, ldb, x, ldx, rcond, berr, n_err_bnds, err_bnds_norm, err_bnds_comp, nparams, params, work, rwork, info)
ZPORFSX
Definition zporfsx.f:393
zpoequb
subroutine zpoequb(n, a, lda, s, scond, amax, info)
ZPOEQUB
Definition zpoequb.f:119
zpocon
subroutine zpocon(uplo, n, a, lda, anorm, rcond, work, rwork, info)
ZPOCON
Definition zpocon.f:121
zporfs
subroutine zporfs(uplo, n, nrhs, a, lda, af, ldaf, b, ldb, x, ldx, ferr, berr, work, rwork, info)
ZPORFS
Definition zporfs.f:183
zpoequ
subroutine zpoequ(n, a, lda, s, scond, amax, info)
ZPOEQU
Definition zpoequ.f:113
zpotri
subroutine zpotri(uplo, n, a, lda, info)
ZPOTRI
Definition zpotri.f:95
zpotrs
subroutine zpotrs(uplo, n, nrhs, a, lda, b, ldb, info)
ZPOTRS
Definition zpotrs.f:110
zerrpo
subroutine zerrpo(path, nunit)
ZERRPO
Definition zerrpo.f:55
zpotrf
subroutine zpotrf(uplo, n, a, lda, info)
ZPOTRF VARIANT: right looking block version of the algorithm, calling Level 3 BLAS.
Definition zpotrf.f:102