OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
derrpox.f
Go to the documentation of this file.
1*> \brief \b DERRPOX
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 DERRPO( 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*> DERRPO tests the error exits for the DOUBLE PRECISION routines
25*> for symmetric positive definite matrices.
26*>
27*> Note that this file is used only when the XBLAS are available,
28*> otherwise derrpo.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 double_lin
55*
56* =====================================================================
57 SUBROUTINE derrpo( 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 INTEGER IW( NMAX )
82 DOUBLE PRECISION A( NMAX, NMAX ), AF( NMAX, NMAX ), B( NMAX ),
83 $ R1( NMAX ), R2( NMAX ), W( 3*NMAX ), X( NMAX ),
84 $ S( NMAX ), ERR_BNDS_N( NMAX, 3 ),
85 $ ERR_BNDS_C( NMAX, 3), PARAMS( 1 )
86* ..
87* .. External Functions ..
88 LOGICAL LSAMEN
89 EXTERNAL lsamen
90* ..
91* .. External Subroutines ..
92 EXTERNAL alaesm, chkxer, dpbcon, dpbequ, dpbrfs, dpbtf2,
93 $ dpbtrf, dpbtrs, dpocon, dpoequ, dporfs, dpotf2,
94 $ dpotrf, dpotri, dpotrs, dppcon, dppequ, dpprfs,
95 $ dpptrf, dpptri, dpptrs, dpoequb, dporfsx
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
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 ) = 1.d0 / dble( i+j )
120 af( i, j ) = 1.d0 / dble( i+j )
121 10 CONTINUE
122 b( j ) = 0.d0
123 r1( j ) = 0.d0
124 r2( j ) = 0.d0
125 w( j ) = 0.d0
126 x( j ) = 0.d0
127 s( j ) = 0.d0
128 iw( j ) = j
129 20 CONTINUE
130 ok = .true.
131*
132 IF( lsamen( 2, c2, 'PO' ) ) THEN
133*
134* Test error exits of the routines that use the Cholesky
135* decomposition of a symmetric positive definite matrix.
136*
137* DPOTRF
138*
139 srnamt = 'DPOTRF'
140 infot = 1
141 CALL dpotrf( '/', 0, a, 1, info )
142 CALL chkxer( 'DPOTRF', infot, nout, lerr, ok )
143 infot = 2
144 CALL dpotrf( 'U', -1, a, 1, info )
145 CALL chkxer( 'DPOTRF', infot, nout, lerr, ok )
146 infot = 4
147 CALL dpotrf( 'U', 2, a, 1, info )
148 CALL chkxer( 'DPOTRF', infot, nout, lerr, ok )
149*
150* DPOTF2
151*
152 srnamt = 'DPOTF2'
153 infot = 1
154 CALL dpotf2( '/', 0, a, 1, info )
155 CALL chkxer( 'DPOTF2', infot, nout, lerr, ok )
156 infot = 2
157 CALL dpotf2( 'U', -1, a, 1, info )
158 CALL chkxer( 'DPOTF2', infot, nout, lerr, ok )
159 infot = 4
160 CALL dpotf2( 'U', 2, a, 1, info )
161 CALL chkxer( 'DPOTF2', infot, nout, lerr, ok )
162*
163* DPOTRI
164*
165 srnamt = 'DPOTRI'
166 infot = 1
167 CALL dpotri( '/', 0, a, 1, info )
168 CALL chkxer( 'DPOTRI', infot, nout, lerr, ok )
169 infot = 2
170 CALL dpotri( 'U', -1, a, 1, info )
171 CALL chkxer( 'DPOTRI', infot, nout, lerr, ok )
172 infot = 4
173 CALL dpotri( 'U', 2, a, 1, info )
174 CALL chkxer( 'DPOTRI', infot, nout, lerr, ok )
175*
176* DPOTRS
177*
178 srnamt = 'DPOTRS'
179 infot = 1
180 CALL dpotrs( '/', 0, 0, a, 1, b, 1, info )
181 CALL chkxer( 'DPOTRS', infot, nout, lerr, ok )
182 infot = 2
183 CALL dpotrs( 'U', -1, 0, a, 1, b, 1, info )
184 CALL chkxer( 'DPOTRS', infot, nout, lerr, ok )
185 infot = 3
186 CALL dpotrs( 'U', 0, -1, a, 1, b, 1, info )
187 CALL chkxer( 'DPOTRS', infot, nout, lerr, ok )
188 infot = 5
189 CALL dpotrs( 'U', 2, 1, a, 1, b, 2, info )
190 CALL chkxer( 'DPOTRS', infot, nout, lerr, ok )
191 infot = 7
192 CALL dpotrs( 'U', 2, 1, a, 2, b, 1, info )
193 CALL chkxer( 'DPOTRS', infot, nout, lerr, ok )
194*
195* DPORFS
196*
197 srnamt = 'DPORFS'
198 infot = 1
199 CALL dporfs( '/', 0, 0, a, 1, af, 1, b, 1, x, 1, r1, r2, w, iw,
200 $ info )
201 CALL chkxer( 'DPORFS', infot, nout, lerr, ok )
202 infot = 2
203 CALL dporfs( 'U', -1, 0, a, 1, af, 1, b, 1, x, 1, r1, r2, w,
204 $ iw, info )
205 CALL chkxer( 'DPORFS', infot, nout, lerr, ok )
206 infot = 3
207 CALL dporfs( 'U', 0, -1, a, 1, af, 1, b, 1, x, 1, r1, r2, w,
208 $ iw, info )
209 CALL chkxer( 'DPORFS', infot, nout, lerr, ok )
210 infot = 5
211 CALL dporfs( 'U', 2, 1, a, 1, af, 2, b, 2, x, 2, r1, r2, w, iw,
212 $ info )
213 CALL chkxer( 'DPORFS', infot, nout, lerr, ok )
214 infot = 7
215 CALL dporfs( 'U', 2, 1, a, 2, af, 1, b, 2, x, 2, r1, r2, w, iw,
216 $ info )
217 CALL chkxer( 'DPORFS', infot, nout, lerr, ok )
218 infot = 9
219 CALL dporfs( 'U', 2, 1, a, 2, af, 2, b, 1, x, 2, r1, r2, w, iw,
220 $ info )
221 CALL chkxer( 'DPORFS', infot, nout, lerr, ok )
222 infot = 11
223 CALL dporfs( 'U', 2, 1, a, 2, af, 2, b, 2, x, 1, r1, r2, w, iw,
224 $ info )
225 CALL chkxer( 'DPORFS', infot, nout, lerr, ok )
226*
227* DPORFSX
228*
229 n_err_bnds = 3
230 nparams = 0
231 srnamt = 'DPORFSX'
232 infot = 1
233 CALL dporfsx( '/', eq, 0, 0, a, 1, af, 1, s, b, 1, x, 1,
234 $ rcond, berr, n_err_bnds, err_bnds_n, err_bnds_c, nparams,
235 $ params, w, iw, info )
236 CALL chkxer( 'DPORFSX', infot, nout, lerr, ok )
237 infot = 2
238 CALL dporfsx( 'U', "/", -1, 0, a, 1, af, 1, s, b, 1, x, 1,
239 $ rcond, berr, n_err_bnds, err_bnds_n, err_bnds_c, nparams,
240 $ params, w, iw, info )
241 CALL chkxer( 'DPORFSX', infot, nout, lerr, ok )
242 eq = 'N'
243 infot = 3
244 CALL dporfsx( 'U', eq, -1, 0, a, 1, af, 1, s, b, 1, x, 1,
245 $ rcond, berr, n_err_bnds, err_bnds_n, err_bnds_c, nparams,
246 $ params, w, iw, info )
247 CALL chkxer( 'DPORFSX', infot, nout, lerr, ok )
248 infot = 4
249 CALL dporfsx( 'U', eq, 0, -1, a, 1, af, 1, s, b, 1, x, 1,
250 $ rcond, berr, n_err_bnds, err_bnds_n, err_bnds_c, nparams,
251 $ params, w, iw, info )
252 CALL chkxer( 'DPORFSX', infot, nout, lerr, ok )
253 infot = 6
254 CALL dporfsx( 'U', eq, 2, 1, a, 1, af, 2, s, b, 2, x, 2,
255 $ rcond, berr, n_err_bnds, err_bnds_n, err_bnds_c, nparams,
256 $ params, w, iw, info )
257 CALL chkxer( 'DPORFSX', infot, nout, lerr, ok )
258 infot = 8
259 CALL dporfsx( 'U', eq, 2, 1, a, 2, af, 1, s, b, 2, x, 2,
260 $ rcond, berr, n_err_bnds, err_bnds_n, err_bnds_c, nparams,
261 $ params, w, iw, info )
262 CALL chkxer( 'DPORFSX', infot, nout, lerr, ok )
263 infot = 11
264 CALL dporfsx( 'U', eq, 2, 1, a, 2, af, 2, s, b, 1, x, 2,
265 $ rcond, berr, n_err_bnds, err_bnds_n, err_bnds_c, nparams,
266 $ params, w, iw, info )
267 CALL chkxer( 'DPORFSX', infot, nout, lerr, ok )
268 infot = 13
269 CALL dporfsx( 'U', eq, 2, 1, a, 2, af, 2, s, b, 2, x, 1,
270 $ rcond, berr, n_err_bnds, err_bnds_n, err_bnds_c, nparams,
271 $ params, w, iw, info )
272 CALL chkxer( 'DPORFSX', infot, nout, lerr, ok )
273*
274* DPOCON
275*
276 srnamt = 'DPOCON'
277 infot = 1
278 CALL dpocon( '/', 0, a, 1, anrm, rcond, w, iw, info )
279 CALL chkxer( 'DPOCON', infot, nout, lerr, ok )
280 infot = 2
281 CALL dpocon( 'U', -1, a, 1, anrm, rcond, w, iw, info )
282 CALL chkxer( 'DPOCON', infot, nout, lerr, ok )
283 infot = 4
284 CALL dpocon( 'U', 2, a, 1, anrm, rcond, w, iw, info )
285 CALL chkxer( 'DPOCON', infot, nout, lerr, ok )
286*
287* DPOEQU
288*
289 srnamt = 'DPOEQU'
290 infot = 1
291 CALL dpoequ( -1, a, 1, r1, rcond, anrm, info )
292 CALL chkxer( 'DPOEQU', infot, nout, lerr, ok )
293 infot = 3
294 CALL dpoequ( 2, a, 1, r1, rcond, anrm, info )
295 CALL chkxer( 'DPOEQU', infot, nout, lerr, ok )
296*
297* DPOEQUB
298*
299 srnamt = 'dpoequb'
300 INFOT = 1
301 CALL DPOEQUB( -1, A, 1, R1, RCOND, ANRM, INFO )
302 CALL CHKXER( 'dpoequb', INFOT, NOUT, LERR, OK )
303 INFOT = 3
304 CALL DPOEQUB( 2, A, 1, R1, RCOND, ANRM, INFO )
305 CALL CHKXER( 'dpoequb', INFOT, NOUT, LERR, OK )
306*
307 ELSE IF( LSAMEN( 2, C2, 'pp' ) ) THEN
308*
309* Test error exits of the routines that use the Cholesky
310* decomposition of a symmetric positive definite packed matrix.
311*
312* DPPTRF
313*
314 SRNAMT = 'dpptrf'
315 INFOT = 1
316 CALL DPPTRF( '/', 0, A, INFO )
317 CALL CHKXER( 'dpptrf', INFOT, NOUT, LERR, OK )
318 INFOT = 2
319 CALL DPPTRF( 'u', -1, A, INFO )
320 CALL CHKXER( 'dpptrf', INFOT, NOUT, LERR, OK )
321*
322* DPPTRI
323*
324 SRNAMT = 'dpptri'
325 INFOT = 1
326 CALL DPPTRI( '/', 0, A, INFO )
327 CALL CHKXER( 'dpptri', INFOT, NOUT, LERR, OK )
328 INFOT = 2
329 CALL DPPTRI( 'u', -1, A, INFO )
330 CALL CHKXER( 'dpptri', INFOT, NOUT, LERR, OK )
331*
332* DPPTRS
333*
334 SRNAMT = 'dpptrs'
335 INFOT = 1
336 CALL DPPTRS( '/', 0, 0, A, B, 1, INFO )
337 CALL CHKXER( 'dpptrs', INFOT, NOUT, LERR, OK )
338 INFOT = 2
339 CALL DPPTRS( 'u', -1, 0, A, B, 1, INFO )
340 CALL CHKXER( 'dpptrs', INFOT, NOUT, LERR, OK )
341 INFOT = 3
342 CALL DPPTRS( 'u', 0, -1, A, B, 1, INFO )
343 CALL CHKXER( 'dpptrs', INFOT, NOUT, LERR, OK )
344 INFOT = 6
345 CALL DPPTRS( 'u', 2, 1, A, B, 1, INFO )
346 CALL CHKXER( 'dpptrs', INFOT, NOUT, LERR, OK )
347*
348* DPPRFS
349*
350 SRNAMT = 'dpprfs'
351 INFOT = 1
352 CALL DPPRFS( '/', 0, 0, A, AF, B, 1, X, 1, R1, R2, W, IW,
353 $ INFO )
354 CALL CHKXER( 'dpprfs', INFOT, NOUT, LERR, OK )
355 INFOT = 2
356 CALL DPPRFS( 'u', -1, 0, A, AF, B, 1, X, 1, R1, R2, W, IW,
357 $ INFO )
358 CALL CHKXER( 'dpprfs', INFOT, NOUT, LERR, OK )
359 INFOT = 3
360 CALL DPPRFS( 'u', 0, -1, A, AF, B, 1, X, 1, R1, R2, W, IW,
361 $ INFO )
362 CALL CHKXER( 'dpprfs', INFOT, NOUT, LERR, OK )
363 INFOT = 7
364 CALL DPPRFS( 'u', 2, 1, A, AF, B, 1, X, 2, R1, R2, W, IW,
365 $ INFO )
366 CALL CHKXER( 'dpprfs', INFOT, NOUT, LERR, OK )
367 INFOT = 9
368 CALL DPPRFS( 'u', 2, 1, A, AF, B, 2, X, 1, R1, R2, W, IW,
369 $ INFO )
370 CALL CHKXER( 'dpprfs', INFOT, NOUT, LERR, OK )
371*
372* DPPCON
373*
374 SRNAMT = 'dppcon'
375 INFOT = 1
376 CALL DPPCON( '/', 0, A, ANRM, RCOND, W, IW, INFO )
377 CALL CHKXER( 'dppcon', INFOT, NOUT, LERR, OK )
378 INFOT = 2
379 CALL DPPCON( 'u', -1, A, ANRM, RCOND, W, IW, INFO )
380 CALL CHKXER( 'dppcon', INFOT, NOUT, LERR, OK )
381*
382* DPPEQU
383*
384 SRNAMT = 'dppequ'
385 INFOT = 1
386 CALL DPPEQU( '/', 0, A, R1, RCOND, ANRM, INFO )
387 CALL CHKXER( 'dppequ', INFOT, NOUT, LERR, OK )
388 INFOT = 2
389 CALL DPPEQU( 'u', -1, A, R1, RCOND, ANRM, INFO )
390 CALL CHKXER( 'dppequ', INFOT, NOUT, LERR, OK )
391*
392 ELSE IF( LSAMEN( 2, C2, 'pb' ) ) THEN
393*
394* Test error exits of the routines that use the Cholesky
395* decomposition of a symmetric positive definite band matrix.
396*
397* DPBTRF
398*
399 SRNAMT = 'dpbtrf'
400 INFOT = 1
401 CALL DPBTRF( '/', 0, 0, A, 1, INFO )
402 CALL CHKXER( 'dpbtrf', INFOT, NOUT, LERR, OK )
403 INFOT = 2
404 CALL DPBTRF( 'u', -1, 0, A, 1, INFO )
405 CALL CHKXER( 'dpbtrf', INFOT, NOUT, LERR, OK )
406 INFOT = 3
407 CALL DPBTRF( 'u', 1, -1, A, 1, INFO )
408 CALL CHKXER( 'dpbtrf', INFOT, NOUT, LERR, OK )
409 INFOT = 5
410 CALL DPBTRF( 'u', 2, 1, A, 1, INFO )
411 CALL CHKXER( 'dpbtrf', INFOT, NOUT, LERR, OK )
412*
413* DPBTF2
414*
415 SRNAMT = 'dpbtf2'
416 INFOT = 1
417 CALL DPBTF2( '/', 0, 0, A, 1, INFO )
418 CALL CHKXER( 'dpbtf2', INFOT, NOUT, LERR, OK )
419 INFOT = 2
420 CALL DPBTF2( 'u', -1, 0, A, 1, INFO )
421 CALL CHKXER( 'dpbtf2', INFOT, NOUT, LERR, OK )
422 INFOT = 3
423 CALL DPBTF2( 'u', 1, -1, A, 1, INFO )
424 CALL CHKXER( 'dpbtf2', INFOT, NOUT, LERR, OK )
425 INFOT = 5
426 CALL DPBTF2( 'u', 2, 1, A, 1, INFO )
427 CALL CHKXER( 'dpbtf2', INFOT, NOUT, LERR, OK )
428*
429* DPBTRS
430*
431 SRNAMT = 'dpbtrs'
432 INFOT = 1
433 CALL DPBTRS( '/', 0, 0, 0, A, 1, B, 1, INFO )
434 CALL CHKXER( 'dpbtrs', INFOT, NOUT, LERR, OK )
435 INFOT = 2
436 CALL DPBTRS( 'u', -1, 0, 0, A, 1, B, 1, INFO )
437 CALL CHKXER( 'dpbtrs', INFOT, NOUT, LERR, OK )
438 INFOT = 3
439 CALL DPBTRS( 'u', 1, -1, 0, A, 1, B, 1, INFO )
440 CALL CHKXER( 'dpbtrs', INFOT, NOUT, LERR, OK )
441 INFOT = 4
442 CALL DPBTRS( 'u', 0, 0, -1, A, 1, B, 1, INFO )
443 CALL CHKXER( 'dpbtrs', INFOT, NOUT, LERR, OK )
444 INFOT = 6
445 CALL DPBTRS( 'u', 2, 1, 1, A, 1, B, 1, INFO )
446 CALL CHKXER( 'dpbtrs', INFOT, NOUT, LERR, OK )
447 INFOT = 8
448 CALL DPBTRS( 'u', 2, 0, 1, A, 1, B, 1, INFO )
449 CALL CHKXER( 'dpbtrs', INFOT, NOUT, LERR, OK )
450*
451* DPBRFS
452*
453 SRNAMT = 'dpbrfs'
454 INFOT = 1
455 CALL DPBRFS( '/', 0, 0, 0, A, 1, AF, 1, B, 1, X, 1, R1, R2, W,
456 $ IW, INFO )
457 CALL CHKXER( 'dpbrfs', INFOT, NOUT, LERR, OK )
458 INFOT = 2
459 CALL DPBRFS( 'u', -1, 0, 0, A, 1, AF, 1, B, 1, X, 1, R1, R2, W,
460 $ IW, INFO )
461 CALL CHKXER( 'dpbrfs', INFOT, NOUT, LERR, OK )
462 INFOT = 3
463 CALL DPBRFS( 'u', 1, -1, 0, A, 1, AF, 1, B, 1, X, 1, R1, R2, W,
464 $ IW, INFO )
465 CALL CHKXER( 'dpbrfs', INFOT, NOUT, LERR, OK )
466 INFOT = 4
467 CALL DPBRFS( 'u', 0, 0, -1, A, 1, AF, 1, B, 1, X, 1, R1, R2, W,
468 $ IW, INFO )
469 CALL CHKXER( 'dpbrfs', INFOT, NOUT, LERR, OK )
470 INFOT = 6
471 CALL DPBRFS( 'u', 2, 1, 1, A, 1, AF, 2, B, 2, X, 2, R1, R2, W,
472 $ IW, INFO )
473 CALL CHKXER( 'dpbrfs', INFOT, NOUT, LERR, OK )
474 INFOT = 8
475 CALL DPBRFS( 'u', 2, 1, 1, A, 2, AF, 1, B, 2, X, 2, R1, R2, W,
476 $ IW, INFO )
477 CALL CHKXER( 'dpbrfs', INFOT, NOUT, LERR, OK )
478 INFOT = 10
479 CALL DPBRFS( 'u', 2, 0, 1, A, 1, AF, 1, B, 1, X, 2, R1, R2, W,
480 $ IW, INFO )
481 CALL CHKXER( 'dpbrfs', INFOT, NOUT, LERR, OK )
482 INFOT = 12
483 CALL DPBRFS( 'u', 2, 0, 1, A, 1, AF, 1, B, 2, X, 1, R1, R2, W,
484 $ IW, INFO )
485 CALL CHKXER( 'dpbrfs', INFOT, NOUT, LERR, OK )
486*
487* DPBCON
488*
489 SRNAMT = 'dpbcon'
490 INFOT = 1
491 CALL DPBCON( '/', 0, 0, A, 1, ANRM, RCOND, W, IW, INFO )
492 CALL CHKXER( 'dpbcon', INFOT, NOUT, LERR, OK )
493 INFOT = 2
494 CALL DPBCON( 'u', -1, 0, A, 1, ANRM, RCOND, W, IW, INFO )
495 CALL CHKXER( 'dpbcon', INFOT, NOUT, LERR, OK )
496 INFOT = 3
497 CALL DPBCON( 'u', 1, -1, A, 1, ANRM, RCOND, W, IW, INFO )
498 CALL CHKXER( 'dpbcon', INFOT, NOUT, LERR, OK )
499 INFOT = 5
500 CALL DPBCON( 'u', 2, 1, A, 1, ANRM, RCOND, W, IW, INFO )
501 CALL CHKXER( 'dpbcon', INFOT, NOUT, LERR, OK )
502*
503* DPBEQU
504*
505 SRNAMT = 'dpbequ'
506 INFOT = 1
507 CALL DPBEQU( '/', 0, 0, A, 1, R1, RCOND, ANRM, INFO )
508 CALL CHKXER( 'dpbequ', INFOT, NOUT, LERR, OK )
509 INFOT = 2
510 CALL DPBEQU( 'u', -1, 0, A, 1, R1, RCOND, ANRM, INFO )
511 CALL CHKXER( 'dpbequ', INFOT, NOUT, LERR, OK )
512 INFOT = 3
513 CALL DPBEQU( 'u', 1, -1, A, 1, R1, RCOND, ANRM, INFO )
514 CALL CHKXER( 'dpbequ', INFOT, NOUT, LERR, OK )
515 INFOT = 5
516 CALL DPBEQU( 'u', 2, 1, A, 1, R1, RCOND, ANRM, INFO )
517 CALL CHKXER( 'dpbequ', INFOT, NOUT, LERR, OK )
518 END IF
519*
520* Print a summary line.
521*
522 CALL ALAESM( PATH, OK, NOUT )
523*
524 RETURN
525*
526* End of DERRPOX
527*
528 END
chkxer
subroutine chkxer(srnamt, infot, nout, lerr, ok)
Definition cblat2.f:3196
alaesm
subroutine alaesm(path, ok, nout)
ALAESM
Definition alaesm.f:63
dpbtf2
subroutine dpbtf2(uplo, n, kd, ab, ldab, info)
DPBTF2 computes the Cholesky factorization of a symmetric/Hermitian positive definite band matrix (un...
Definition dpbtf2.f:142
dpbtrs
subroutine dpbtrs(uplo, n, kd, nrhs, ab, ldab, b, ldb, info)
DPBTRS
Definition dpbtrs.f:121
dppcon
subroutine dppcon(uplo, n, ap, anorm, rcond, work, iwork, info)
DPPCON
Definition dppcon.f:118
dpbequ
subroutine dpbequ(uplo, n, kd, ab, ldab, s, scond, amax, info)
DPBEQU
Definition dpbequ.f:129
dpptrs
subroutine dpptrs(uplo, n, nrhs, ap, b, ldb, info)
DPPTRS
Definition dpptrs.f:108
dpbcon
subroutine dpbcon(uplo, n, kd, ab, ldab, anorm, rcond, work, iwork, info)
DPBCON
Definition dpbcon.f:132
dpbrfs
subroutine dpbrfs(uplo, n, kd, nrhs, ab, ldab, afb, ldafb, b, ldb, x, ldx, ferr, berr, work, iwork, info)
DPBRFS
Definition dpbrfs.f:189
dpptri
subroutine dpptri(uplo, n, ap, info)
DPPTRI
Definition dpptri.f:93
dpptrf
subroutine dpptrf(uplo, n, ap, info)
DPPTRF
Definition dpptrf.f:119
dpprfs
subroutine dpprfs(uplo, n, nrhs, ap, afp, b, ldb, x, ldx, ferr, berr, work, iwork, info)
DPPRFS
Definition dpprfs.f:171
dppequ
subroutine dppequ(uplo, n, ap, s, scond, amax, info)
DPPEQU
Definition dppequ.f:116
dpbtrf
subroutine dpbtrf(uplo, n, kd, ab, ldab, info)
DPBTRF
Definition dpbtrf.f:142
dporfsx
subroutine dporfsx(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, iwork, info)
DPORFSX
Definition dporfsx.f:394
dpoequb
subroutine dpoequb(n, a, lda, s, scond, amax, info)
DPOEQUB
Definition dpoequb.f:118
dpoequ
subroutine dpoequ(n, a, lda, s, scond, amax, info)
DPOEQU
Definition dpoequ.f:112
dpotrs
subroutine dpotrs(uplo, n, nrhs, a, lda, b, ldb, info)
DPOTRS
Definition dpotrs.f:110
dpotrf
subroutine dpotrf(uplo, n, a, lda, info)
DPOTRF
Definition dpotrf.f:107
dpocon
subroutine dpocon(uplo, n, a, lda, anorm, rcond, work, iwork, info)
DPOCON
Definition dpocon.f:121
dpotf2
subroutine dpotf2(uplo, n, a, lda, info)
DPOTF2 computes the Cholesky factorization of a symmetric/Hermitian positive definite matrix (unblock...
Definition dpotf2.f:109
dpotri
subroutine dpotri(uplo, n, a, lda, info)
DPOTRI
Definition dpotri.f:95
dporfs
subroutine dporfs(uplo, n, nrhs, a, lda, af, ldaf, b, ldb, x, ldx, ferr, berr, work, iwork, info)
DPORFS
Definition dporfs.f:183
derrpo
subroutine derrpo(path, nunit)
DERRPO
Definition derrpo.f:55