LAPACK 3.11.0
LAPACK: Linear Algebra PACKage
Loading...
Searching...
No Matches
alahd.f
Go to the documentation of this file.
1*> \brief \b ALAHD
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 ALAHD( IOUNIT, PATH )
12*
13* .. Scalar Arguments ..
14* CHARACTER*3 PATH
15* INTEGER IOUNIT
16* ..
17*
18*
19*> \par Purpose:
20* =============
21*>
22*> \verbatim
23*>
24*> ALAHD prints header information for the different test paths.
25*> \endverbatim
26*
27* Arguments:
28* ==========
29*
30*> \param[in] IOUNIT
31*> \verbatim
32*> IOUNIT is INTEGER
33*> The unit number to which the header information should be
34*> printed.
35*> \endverbatim
36*>
37*> \param[in] PATH
38*> \verbatim
39*> PATH is CHARACTER*3
40*> The name of the path for which the header information is to
41*> be printed. Current paths are
42*> _GE: General matrices
43*> _GB: General band
44*> _GT: General Tridiagonal
45*> _PO: Symmetric or Hermitian positive definite
46*> _PS: Symmetric or Hermitian positive semi-definite
47*> _PP: Symmetric or Hermitian positive definite packed
48*> _PB: Symmetric or Hermitian positive definite band
49*> _PT: Symmetric or Hermitian positive definite tridiagonal
50*> _SY: Symmetric indefinite,
51*> with partial (Bunch-Kaufman) pivoting
52*> _SR: Symmetric indefinite,
53*> with rook (bounded Bunch-Kaufman) pivoting
54*> _SK: Symmetric indefinite,
55*> with rook (bounded Bunch-Kaufman) pivoting
56*> ( new storage format for factors:
57*> L and diagonal of D is stored in A,
58*> subdiagonal of D is stored in E )
59*> _SP: Symmetric indefinite packed,
60*> with partial (Bunch-Kaufman) pivoting
61*> _HA: (complex) Hermitian ,
62*> with Aasen Algorithm
63*> _HE: (complex) Hermitian indefinite,
64*> with partial (Bunch-Kaufman) pivoting
65*> _HR: (complex) Hermitian indefinite,
66*> with rook (bounded Bunch-Kaufman) pivoting
67*> _HK: (complex) Hermitian indefinite,
68*> with rook (bounded Bunch-Kaufman) pivoting
69*> ( new storage format for factors:
70*> L and diagonal of D is stored in A,
71*> subdiagonal of D is stored in E )
72*> _HP: (complex) Hermitian indefinite packed,
73*> with partial (Bunch-Kaufman) pivoting
74*> _TR: Triangular
75*> _TP: Triangular packed
76*> _TB: Triangular band
77*> _QR: QR (general matrices)
78*> _LQ: LQ (general matrices)
79*> _QL: QL (general matrices)
80*> _RQ: RQ (general matrices)
81*> _QP: QR with column pivoting
82*> _TZ: Trapezoidal
83*> _LS: Least Squares driver routines
84*> _LU: LU variants
85*> _CH: Cholesky variants
86*> _QS: QR variants
87*> _QT: QRT (general matrices)
88*> _QX: QRT (triangular-pentagonal matrices)
89*> _TS: QR routines for tall-skinny and short-wide matrices
90*> _HH: Householder reconstruction for tall-skinny matrices
91*> The first character must be one of S, D, C, or Z (C or Z only
92*> if complex).
93*> \endverbatim
94*
95* Authors:
96* ========
97*
98*> \author Univ. of Tennessee
99*> \author Univ. of California Berkeley
100*> \author Univ. of Colorado Denver
101*> \author NAG Ltd.
102*
103*> \ingroup aux_lin
104*
105* =====================================================================
106 SUBROUTINE alahd( IOUNIT, PATH )
107*
108* -- LAPACK test routine --
109* -- LAPACK is a software package provided by Univ. of Tennessee, --
110* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
111*
112* .. Scalar Arguments ..
113 CHARACTER*3 PATH
114 INTEGER IOUNIT
115* ..
116*
117* =====================================================================
118*
119* .. Local Scalars ..
120 LOGICAL CORZ, SORD
121 CHARACTER C1, C3
122 CHARACTER*2 P2
123 CHARACTER*4 EIGCNM
124 CHARACTER*32 SUBNAM
125 CHARACTER*9 SYM
126* ..
127* .. External Functions ..
128 LOGICAL LSAME, LSAMEN
129 EXTERNAL lsame, lsamen
130* ..
131* .. Intrinsic Functions ..
132 INTRINSIC len_trim
133* ..
134* .. Executable Statements ..
135*
136 IF( iounit.LE.0 )
137 $ RETURN
138 c1 = path( 1: 1 )
139 c3 = path( 3: 3 )
140 p2 = path( 2: 3 )
141 sord = lsame( c1, 'S' ) .OR. lsame( c1, 'D' )
142 corz = lsame( c1, 'C' ) .OR. lsame( c1, 'Z' )
143 IF( .NOT.( sord .OR. corz ) )
144 $ RETURN
145*
146 IF( lsamen( 2, p2, 'GE' ) ) THEN
147*
148* GE: General dense
149*
150 WRITE( iounit, fmt = 9999 )path
151 WRITE( iounit, fmt = '( '' Matrix types:'' )' )
152 WRITE( iounit, fmt = 9979 )
153 WRITE( iounit, fmt = '( '' Test ratios:'' )' )
154 WRITE( iounit, fmt = 9962 )1
155 WRITE( iounit, fmt = 9961 )2
156 WRITE( iounit, fmt = 9960 )3
157 WRITE( iounit, fmt = 9959 )4
158 WRITE( iounit, fmt = 9958 )5
159 WRITE( iounit, fmt = 9957 )6
160 WRITE( iounit, fmt = 9956 )7
161 WRITE( iounit, fmt = 9955 )8
162 WRITE( iounit, fmt = '( '' Messages:'' )' )
163*
164 ELSE IF( lsamen( 2, p2, 'GB' ) ) THEN
165*
166* GB: General band
167*
168 WRITE( iounit, fmt = 9998 )path
169 WRITE( iounit, fmt = '( '' Matrix types:'' )' )
170 WRITE( iounit, fmt = 9978 )
171 WRITE( iounit, fmt = '( '' Test ratios:'' )' )
172 WRITE( iounit, fmt = 9962 )1
173 WRITE( iounit, fmt = 9960 )2
174 WRITE( iounit, fmt = 9959 )3
175 WRITE( iounit, fmt = 9958 )4
176 WRITE( iounit, fmt = 9957 )5
177 WRITE( iounit, fmt = 9956 )6
178 WRITE( iounit, fmt = 9955 )7
179 WRITE( iounit, fmt = '( '' Messages:'' )' )
180*
181 ELSE IF( lsamen( 2, p2, 'GT' ) ) THEN
182*
183* GT: General tridiagonal
184*
185 WRITE( iounit, fmt = 9997 )path
186 WRITE( iounit, fmt = 9977 )
187 WRITE( iounit, fmt = '( '' Test ratios:'' )' )
188 WRITE( iounit, fmt = 9962 )1
189 WRITE( iounit, fmt = 9960 )2
190 WRITE( iounit, fmt = 9959 )3
191 WRITE( iounit, fmt = 9958 )4
192 WRITE( iounit, fmt = 9957 )5
193 WRITE( iounit, fmt = 9956 )6
194 WRITE( iounit, fmt = 9955 )7
195 WRITE( iounit, fmt = '( '' Messages:'' )' )
196*
197 ELSE IF( lsamen( 2, p2, 'PO' ) .OR. lsamen( 2, p2, 'PP' ) ) THEN
198*
199* PO: Positive definite full
200* PP: Positive definite packed
201*
202 IF( sord ) THEN
203 sym = 'Symmetric'
204 ELSE
205 sym = 'Hermitian'
206 END IF
207 IF( lsame( c3, 'O' ) ) THEN
208 WRITE( iounit, fmt = 9996 )path, sym
209 ELSE
210 WRITE( iounit, fmt = 9995 )path, sym
211 END IF
212 WRITE( iounit, fmt = '( '' Matrix types:'' )' )
213 WRITE( iounit, fmt = 9975 )path
214 WRITE( iounit, fmt = '( '' Test ratios:'' )' )
215 WRITE( iounit, fmt = 9954 )1
216 WRITE( iounit, fmt = 9961 )2
217 WRITE( iounit, fmt = 9960 )3
218 WRITE( iounit, fmt = 9959 )4
219 WRITE( iounit, fmt = 9958 )5
220 WRITE( iounit, fmt = 9957 )6
221 WRITE( iounit, fmt = 9956 )7
222 WRITE( iounit, fmt = 9955 )8
223 WRITE( iounit, fmt = '( '' Messages:'' )' )
224*
225 ELSE IF( lsamen( 2, p2, 'PS' ) ) THEN
226*
227* PS: Positive semi-definite full
228*
229 IF( sord ) THEN
230 sym = 'Symmetric'
231 ELSE
232 sym = 'Hermitian'
233 END IF
234 IF( lsame( c1, 'S' ) .OR. lsame( c1, 'C' ) ) THEN
235 eigcnm = '1E04'
236 ELSE
237 eigcnm = '1D12'
238 END IF
239 WRITE( iounit, fmt = 9995 )path, sym
240 WRITE( iounit, fmt = '( '' Matrix types:'' )' )
241 WRITE( iounit, fmt = 8973 )eigcnm, eigcnm, eigcnm
242 WRITE( iounit, fmt = '( '' Difference:'' )' )
243 WRITE( iounit, fmt = 8972 )c1
244 WRITE( iounit, fmt = '( '' Test ratio:'' )' )
245 WRITE( iounit, fmt = 8950 )
246 WRITE( iounit, fmt = '( '' Messages:'' )' )
247 ELSE IF( lsamen( 2, p2, 'PB' ) ) THEN
248*
249* PB: Positive definite band
250*
251 IF( sord ) THEN
252 WRITE( iounit, fmt = 9994 )path, 'Symmetric'
253 ELSE
254 WRITE( iounit, fmt = 9994 )path, 'Hermitian'
255 END IF
256 WRITE( iounit, fmt = '( '' Matrix types:'' )' )
257 WRITE( iounit, fmt = 9973 )path
258 WRITE( iounit, fmt = '( '' Test ratios:'' )' )
259 WRITE( iounit, fmt = 9954 )1
260 WRITE( iounit, fmt = 9960 )2
261 WRITE( iounit, fmt = 9959 )3
262 WRITE( iounit, fmt = 9958 )4
263 WRITE( iounit, fmt = 9957 )5
264 WRITE( iounit, fmt = 9956 )6
265 WRITE( iounit, fmt = 9955 )7
266 WRITE( iounit, fmt = '( '' Messages:'' )' )
267*
268 ELSE IF( lsamen( 2, p2, 'PT' ) ) THEN
269*
270* PT: Positive definite tridiagonal
271*
272 IF( sord ) THEN
273 WRITE( iounit, fmt = 9993 )path, 'Symmetric'
274 ELSE
275 WRITE( iounit, fmt = 9993 )path, 'Hermitian'
276 END IF
277 WRITE( iounit, fmt = 9976 )
278 WRITE( iounit, fmt = '( '' Test ratios:'' )' )
279 WRITE( iounit, fmt = 9952 )1
280 WRITE( iounit, fmt = 9960 )2
281 WRITE( iounit, fmt = 9959 )3
282 WRITE( iounit, fmt = 9958 )4
283 WRITE( iounit, fmt = 9957 )5
284 WRITE( iounit, fmt = 9956 )6
285 WRITE( iounit, fmt = 9955 )7
286 WRITE( iounit, fmt = '( '' Messages:'' )' )
287*
288 ELSE IF( lsamen( 2, p2, 'SY' ) ) THEN
289*
290* SY: Symmetric indefinite full,
291* with partial (Bunch-Kaufman) pivoting algorithm
292*
293 IF( lsame( c3, 'Y' ) ) THEN
294 WRITE( iounit, fmt = 9992 )path, 'Symmetric'
295 ELSE
296 WRITE( iounit, fmt = 9991 )path, 'Symmetric'
297 END IF
298 WRITE( iounit, fmt = '( '' Matrix types:'' )' )
299 IF( sord ) THEN
300 WRITE( iounit, fmt = 9972 )
301 ELSE
302 WRITE( iounit, fmt = 9971 )
303 END IF
304 WRITE( iounit, fmt = '( '' Test ratios:'' )' )
305 WRITE( iounit, fmt = 9953 )1
306 WRITE( iounit, fmt = 9961 )2
307 WRITE( iounit, fmt = 9960 )3
308 WRITE( iounit, fmt = 9960 )4
309 WRITE( iounit, fmt = 9959 )5
310 WRITE( iounit, fmt = 9958 )6
311 WRITE( iounit, fmt = 9956 )7
312 WRITE( iounit, fmt = 9957 )8
313 WRITE( iounit, fmt = 9955 )9
314 WRITE( iounit, fmt = '( '' Messages:'' )' )
315*
316 ELSE IF( lsamen( 2, p2, 'SR' ) .OR. lsamen( 2, p2, 'SK') ) THEN
317*
318* SR: Symmetric indefinite full,
319* with rook (bounded Bunch-Kaufman) pivoting algorithm
320*
321* SK: Symmetric indefinite full,
322* with rook (bounded Bunch-Kaufman) pivoting algorithm,
323* ( new storage format for factors:
324* L and diagonal of D is stored in A,
325* subdiagonal of D is stored in E )
326*
327 WRITE( iounit, fmt = 9892 )path, 'Symmetric'
328*
329 WRITE( iounit, fmt = '( '' Matrix types:'' )' )
330 IF( sord ) THEN
331 WRITE( iounit, fmt = 9972 )
332 ELSE
333 WRITE( iounit, fmt = 9971 )
334 END IF
335*
336 WRITE( iounit, fmt = '( '' Test ratios:'' )' )
337 WRITE( iounit, fmt = 9953 )1
338 WRITE( iounit, fmt = 9961 )2
339 WRITE( iounit, fmt = 9927 )3
340 WRITE( iounit, fmt = 9928 )
341 WRITE( iounit, fmt = 9926 )4
342 WRITE( iounit, fmt = 9928 )
343 WRITE( iounit, fmt = 9960 )5
344 WRITE( iounit, fmt = 9959 )6
345 WRITE( iounit, fmt = 9955 )7
346 WRITE( iounit, fmt = '( '' Messages:'' )' )
347*
348 ELSE IF( lsamen( 2, p2, 'SP' ) ) THEN
349*
350* SP: Symmetric indefinite packed,
351* with partial (Bunch-Kaufman) pivoting algorithm
352*
353 IF( lsame( c3, 'Y' ) ) THEN
354 WRITE( iounit, fmt = 9992 )path, 'Symmetric'
355 ELSE
356 WRITE( iounit, fmt = 9991 )path, 'Symmetric'
357 END IF
358 WRITE( iounit, fmt = '( '' Matrix types:'' )' )
359 IF( sord ) THEN
360 WRITE( iounit, fmt = 9972 )
361 ELSE
362 WRITE( iounit, fmt = 9971 )
363 END IF
364 WRITE( iounit, fmt = '( '' Test ratios:'' )' )
365 WRITE( iounit, fmt = 9953 )1
366 WRITE( iounit, fmt = 9961 )2
367 WRITE( iounit, fmt = 9960 )3
368 WRITE( iounit, fmt = 9959 )4
369 WRITE( iounit, fmt = 9958 )5
370 WRITE( iounit, fmt = 9956 )6
371 WRITE( iounit, fmt = 9957 )7
372 WRITE( iounit, fmt = 9955 )8
373 WRITE( iounit, fmt = '( '' Messages:'' )' )
374*
375 ELSE IF( lsamen( 2, p2, 'HA' ) ) THEN
376*
377* HA: Hermitian,
378* with Assen Algorithm
379*
380 WRITE( iounit, fmt = 9992 )path, 'Hermitian'
381*
382 WRITE( iounit, fmt = '( '' Matrix types:'' )' )
383 WRITE( iounit, fmt = 9972 )
384*
385 WRITE( iounit, fmt = '( '' Test ratios:'' )' )
386 WRITE( iounit, fmt = 9953 )1
387 WRITE( iounit, fmt = 9961 )2
388 WRITE( iounit, fmt = 9960 )3
389 WRITE( iounit, fmt = 9960 )4
390 WRITE( iounit, fmt = 9959 )5
391 WRITE( iounit, fmt = 9958 )6
392 WRITE( iounit, fmt = 9956 )7
393 WRITE( iounit, fmt = 9957 )8
394 WRITE( iounit, fmt = 9955 )9
395 WRITE( iounit, fmt = '( '' Messages:'' )' )
396*
397 ELSE IF( lsamen( 2, p2, 'HE' ) ) THEN
398*
399* HE: Hermitian indefinite full,
400* with partial (Bunch-Kaufman) pivoting algorithm
401*
402 WRITE( iounit, fmt = 9992 )path, 'Hermitian'
403*
404 WRITE( iounit, fmt = '( '' Matrix types:'' )' )
405 WRITE( iounit, fmt = 9972 )
406*
407 WRITE( iounit, fmt = '( '' Test ratios:'' )' )
408 WRITE( iounit, fmt = 9953 )1
409 WRITE( iounit, fmt = 9961 )2
410 WRITE( iounit, fmt = 9960 )3
411 WRITE( iounit, fmt = 9960 )4
412 WRITE( iounit, fmt = 9959 )5
413 WRITE( iounit, fmt = 9958 )6
414 WRITE( iounit, fmt = 9956 )7
415 WRITE( iounit, fmt = 9957 )8
416 WRITE( iounit, fmt = 9955 )9
417 WRITE( iounit, fmt = '( '' Messages:'' )' )
418*
419 ELSE IF( lsamen( 2, p2, 'HR' ) .OR. lsamen( 2, p2, 'HR' ) ) THEN
420*
421* HR: Hermitian indefinite full,
422* with rook (bounded Bunch-Kaufman) pivoting algorithm
423*
424* HK: Hermitian indefinite full,
425* with rook (bounded Bunch-Kaufman) pivoting algorithm,
426* ( new storage format for factors:
427* L and diagonal of D is stored in A,
428* subdiagonal of D is stored in E )
429*
430 WRITE( iounit, fmt = 9892 )path, 'Hermitian'
431*
432 WRITE( iounit, fmt = '( '' Matrix types:'' )' )
433 WRITE( iounit, fmt = 9972 )
434*
435 WRITE( iounit, fmt = '( '' Test ratios:'' )' )
436 WRITE( iounit, fmt = 9953 )1
437 WRITE( iounit, fmt = 9961 )2
438 WRITE( iounit, fmt = 9927 )3
439 WRITE( iounit, fmt = 9928 )
440 WRITE( iounit, fmt = 9926 )4
441 WRITE( iounit, fmt = 9928 )
442 WRITE( iounit, fmt = 9960 )5
443 WRITE( iounit, fmt = 9959 )6
444 WRITE( iounit, fmt = 9955 )7
445 WRITE( iounit, fmt = '( '' Messages:'' )' )
446*
447 ELSE IF( lsamen( 2, p2, 'HP' ) ) THEN
448*
449* HP: Hermitian indefinite packed,
450* with partial (Bunch-Kaufman) pivoting algorithm
451*
452 IF( lsame( c3, 'E' ) ) THEN
453 WRITE( iounit, fmt = 9992 )path, 'Hermitian'
454 ELSE
455 WRITE( iounit, fmt = 9991 )path, 'Hermitian'
456 END IF
457 WRITE( iounit, fmt = '( '' Matrix types:'' )' )
458 WRITE( iounit, fmt = 9972 )
459 WRITE( iounit, fmt = '( '' Test ratios:'' )' )
460 WRITE( iounit, fmt = 9953 )1
461 WRITE( iounit, fmt = 9961 )2
462 WRITE( iounit, fmt = 9960 )3
463 WRITE( iounit, fmt = 9959 )4
464 WRITE( iounit, fmt = 9958 )5
465 WRITE( iounit, fmt = 9956 )6
466 WRITE( iounit, fmt = 9957 )7
467 WRITE( iounit, fmt = 9955 )8
468 WRITE( iounit, fmt = '( '' Messages:'' )' )
469*
470 ELSE IF( lsamen( 2, p2, 'TR' ) .OR. lsamen( 2, p2, 'TP' ) ) THEN
471*
472* TR: Triangular full
473* TP: Triangular packed
474*
475 IF( lsame( c3, 'R' ) ) THEN
476 WRITE( iounit, fmt = 9990 )path
477 subnam = path( 1: 1 ) // 'LATRS'
478 ELSE
479 WRITE( iounit, fmt = 9989 )path
480 subnam = path( 1: 1 ) // 'LATPS'
481 END IF
482 WRITE( iounit, fmt = 9966 )path
483 WRITE( iounit, fmt = 9965 )subnam(1:len_trim( subnam ))
484 WRITE( iounit, fmt = '( '' Test ratios:'' )' )
485 WRITE( iounit, fmt = 9961 )1
486 WRITE( iounit, fmt = 9960 )2
487 WRITE( iounit, fmt = 9959 )3
488 WRITE( iounit, fmt = 9958 )4
489 WRITE( iounit, fmt = 9957 )5
490 WRITE( iounit, fmt = 9956 )6
491 WRITE( iounit, fmt = 9955 )7
492 WRITE( iounit, fmt = 9951 )subnam(1:len_trim( subnam )), 8
493 WRITE( iounit, fmt = '( '' Messages:'' )' )
494*
495 ELSE IF( lsamen( 2, p2, 'TB' ) ) THEN
496*
497* TB: Triangular band
498*
499 WRITE( iounit, fmt = 9988 )path
500 subnam = path( 1: 1 ) // 'LATBS'
501 WRITE( iounit, fmt = 9964 )path
502 WRITE( iounit, fmt = 9963 )subnam(1:len_trim( subnam ))
503 WRITE( iounit, fmt = '( '' Test ratios:'' )' )
504 WRITE( iounit, fmt = 9960 )1
505 WRITE( iounit, fmt = 9959 )2
506 WRITE( iounit, fmt = 9958 )3
507 WRITE( iounit, fmt = 9957 )4
508 WRITE( iounit, fmt = 9956 )5
509 WRITE( iounit, fmt = 9955 )6
510 WRITE( iounit, fmt = 9951 )subnam(1:len_trim( subnam )), 7
511 WRITE( iounit, fmt = '( '' Messages:'' )' )
512*
513 ELSE IF( lsamen( 2, p2, 'QR' ) ) THEN
514*
515* QR decomposition of rectangular matrices
516*
517 WRITE( iounit, fmt = 9987 )path, 'QR'
518 WRITE( iounit, fmt = '( '' Matrix types:'' )' )
519 WRITE( iounit, fmt = 9970 )
520 WRITE( iounit, fmt = '( '' Test ratios:'' )' )
521 WRITE( iounit, fmt = 9950 )1
522 WRITE( iounit, fmt = 6950 )8
523 WRITE( iounit, fmt = 9946 )2
524 WRITE( iounit, fmt = 9944 )3, 'M'
525 WRITE( iounit, fmt = 9943 )4, 'M'
526 WRITE( iounit, fmt = 9942 )5, 'M'
527 WRITE( iounit, fmt = 9941 )6, 'M'
528 WRITE( iounit, fmt = 9960 )7
529 WRITE( iounit, fmt = 6660 )9
530 WRITE( iounit, fmt = '( '' Messages:'' )' )
531*
532 ELSE IF( lsamen( 2, p2, 'LQ' ) ) THEN
533*
534* LQ decomposition of rectangular matrices
535*
536 WRITE( iounit, fmt = 9987 )path, 'LQ'
537 WRITE( iounit, fmt = '( '' Matrix types:'' )' )
538 WRITE( iounit, fmt = 9970 )
539 WRITE( iounit, fmt = '( '' Test ratios:'' )' )
540 WRITE( iounit, fmt = 9949 )1
541 WRITE( iounit, fmt = 9945 )2
542 WRITE( iounit, fmt = 9944 )3, 'N'
543 WRITE( iounit, fmt = 9943 )4, 'N'
544 WRITE( iounit, fmt = 9942 )5, 'N'
545 WRITE( iounit, fmt = 9941 )6, 'N'
546 WRITE( iounit, fmt = 9960 )7
547 WRITE( iounit, fmt = '( '' Messages:'' )' )
548*
549 ELSE IF( lsamen( 2, p2, 'QL' ) ) THEN
550*
551* QL decomposition of rectangular matrices
552*
553 WRITE( iounit, fmt = 9987 )path, 'QL'
554 WRITE( iounit, fmt = '( '' Matrix types:'' )' )
555 WRITE( iounit, fmt = 9970 )
556 WRITE( iounit, fmt = '( '' Test ratios:'' )' )
557 WRITE( iounit, fmt = 9948 )1
558 WRITE( iounit, fmt = 9946 )2
559 WRITE( iounit, fmt = 9944 )3, 'M'
560 WRITE( iounit, fmt = 9943 )4, 'M'
561 WRITE( iounit, fmt = 9942 )5, 'M'
562 WRITE( iounit, fmt = 9941 )6, 'M'
563 WRITE( iounit, fmt = 9960 )7
564 WRITE( iounit, fmt = '( '' Messages:'' )' )
565*
566 ELSE IF( lsamen( 2, p2, 'RQ' ) ) THEN
567*
568* RQ decomposition of rectangular matrices
569*
570 WRITE( iounit, fmt = 9987 )path, 'RQ'
571 WRITE( iounit, fmt = '( '' Matrix types:'' )' )
572 WRITE( iounit, fmt = 9970 )
573 WRITE( iounit, fmt = '( '' Test ratios:'' )' )
574 WRITE( iounit, fmt = 9947 )1
575 WRITE( iounit, fmt = 9945 )2
576 WRITE( iounit, fmt = 9944 )3, 'N'
577 WRITE( iounit, fmt = 9943 )4, 'N'
578 WRITE( iounit, fmt = 9942 )5, 'N'
579 WRITE( iounit, fmt = 9941 )6, 'N'
580 WRITE( iounit, fmt = 9960 )7
581 WRITE( iounit, fmt = '( '' Messages:'' )' )
582*
583 ELSE IF( lsamen( 2, p2, 'QP' ) ) THEN
584*
585* QR decomposition with column pivoting
586*
587 WRITE( iounit, fmt = 9986 )path
588 WRITE( iounit, fmt = 9969 )
589 WRITE( iounit, fmt = '( '' Test ratios:'' )' )
590 WRITE( iounit, fmt = 9940 )1
591 WRITE( iounit, fmt = 9939 )2
592 WRITE( iounit, fmt = 9938 )3
593 WRITE( iounit, fmt = '( '' Messages:'' )' )
594*
595 ELSE IF( lsamen( 2, p2, 'TZ' ) ) THEN
596*
597* TZ: Trapezoidal
598*
599 WRITE( iounit, fmt = 9985 )path
600 WRITE( iounit, fmt = 9968 )
601 WRITE( iounit, fmt = 9929 )c1
602 WRITE( iounit, fmt = '( '' Test ratios:'' )' )
603 WRITE( iounit, fmt = 9940 )1
604 WRITE( iounit, fmt = 9937 )2
605 WRITE( iounit, fmt = 9938 )3
606 WRITE( iounit, fmt = '( '' Messages:'' )' )
607*
608 ELSE IF( lsamen( 2, p2, 'LS' ) ) THEN
609*
610* LS: Least Squares driver routines for
611* LS, LST, TSLS, LSD, LSS, LSX and LSY.
612*
613 WRITE( iounit, fmt = 9984 )path
614 WRITE( iounit, fmt = 9967 )
615 WRITE( iounit, fmt = 9921 )c1, c1, c1, c1, c1, c1
616 WRITE( iounit, fmt = 9935 )1
617 WRITE( iounit, fmt = 9931 )2
618 WRITE( iounit, fmt = 9919 )
619 WRITE( iounit, fmt = 9933 )7
620 WRITE( iounit, fmt = 9935 )8
621 WRITE( iounit, fmt = 9934 )9
622 WRITE( iounit, fmt = 9932 )10
623 WRITE( iounit, fmt = 9920 )
624 WRITE( iounit, fmt = '( '' Messages:'' )' )
625*
626 ELSE IF( lsamen( 2, p2, 'LU' ) ) THEN
627*
628* LU factorization variants
629*
630 WRITE( iounit, fmt = 9983 )path
631 WRITE( iounit, fmt = '( '' Matrix types:'' )' )
632 WRITE( iounit, fmt = 9979 )
633 WRITE( iounit, fmt = '( '' Test ratio:'' )' )
634 WRITE( iounit, fmt = 9962 )1
635 WRITE( iounit, fmt = '( '' Messages:'' )' )
636*
637 ELSE IF( lsamen( 2, p2, 'CH' ) ) THEN
638*
639* Cholesky factorization variants
640*
641 WRITE( iounit, fmt = 9982 )path
642 WRITE( iounit, fmt = '( '' Matrix types:'' )' )
643 WRITE( iounit, fmt = 9974 )
644 WRITE( iounit, fmt = '( '' Test ratio:'' )' )
645 WRITE( iounit, fmt = 9954 )1
646 WRITE( iounit, fmt = '( '' Messages:'' )' )
647*
648 ELSE IF( lsamen( 2, p2, 'QS' ) ) THEN
649*
650* QR factorization variants
651*
652 WRITE( iounit, fmt = 9981 )path
653 WRITE( iounit, fmt = '( '' Matrix types:'' )' )
654 WRITE( iounit, fmt = 9970 )
655 WRITE( iounit, fmt = '( '' Test ratios:'' )' )
656*
657 ELSE IF( lsamen( 2, p2, 'QT' ) ) THEN
658*
659* QRT (general matrices)
660*
661 WRITE( iounit, fmt = 8000 ) path
662 WRITE( iounit, fmt = '( '' Test ratios:'' )' )
663 WRITE( iounit, fmt = 8011 ) 1
664 WRITE( iounit, fmt = 8012 ) 2
665 WRITE( iounit, fmt = 8013 ) 3
666 WRITE( iounit, fmt = 8014 ) 4
667 WRITE( iounit, fmt = 8015 ) 5
668 WRITE( iounit, fmt = 8016 ) 6
669*
670 ELSE IF( lsamen( 2, p2, 'QX' ) ) THEN
671*
672* QRT (triangular-pentagonal)
673*
674 WRITE( iounit, fmt = 8001 ) path
675 WRITE( iounit, fmt = '( '' Test ratios:'' )' )
676 WRITE( iounit, fmt = 8017 ) 1
677 WRITE( iounit, fmt = 8018 ) 2
678 WRITE( iounit, fmt = 8019 ) 3
679 WRITE( iounit, fmt = 8020 ) 4
680 WRITE( iounit, fmt = 8021 ) 5
681 WRITE( iounit, fmt = 8022 ) 6
682*
683 ELSE IF( lsamen( 2, p2, 'TQ' ) ) THEN
684*
685* QRT (triangular-pentagonal)
686*
687 WRITE( iounit, fmt = 8002 ) path
688 WRITE( iounit, fmt = '( '' Test ratios:'' )' )
689 WRITE( iounit, fmt = 8023 ) 1
690 WRITE( iounit, fmt = 8024 ) 2
691 WRITE( iounit, fmt = 8025 ) 3
692 WRITE( iounit, fmt = 8026 ) 4
693 WRITE( iounit, fmt = 8027 ) 5
694 WRITE( iounit, fmt = 8028 ) 6
695*
696 ELSE IF( lsamen( 2, p2, 'XQ' ) ) THEN
697*
698* QRT (triangular-pentagonal)
699*
700 WRITE( iounit, fmt = 8003 ) path
701 WRITE( iounit, fmt = '( '' Test ratios:'' )' )
702 WRITE( iounit, fmt = 8029 ) 1
703 WRITE( iounit, fmt = 8030 ) 2
704 WRITE( iounit, fmt = 8031 ) 3
705 WRITE( iounit, fmt = 8032 ) 4
706 WRITE( iounit, fmt = 8033 ) 5
707 WRITE( iounit, fmt = 8034 ) 6
708*
709 ELSE IF( lsamen( 2, p2, 'TS' ) ) THEN
710*
711* TS: QR routines for tall-skinny and short-wide matrices
712*
713 WRITE( iounit, fmt = 8004 ) path
714 WRITE( iounit, fmt = '( '' Test ratios:'' )' )
715 WRITE( iounit, fmt = 8035 ) 1
716 WRITE( iounit, fmt = 8036 ) 2
717 WRITE( iounit, fmt = 8037 ) 3
718 WRITE( iounit, fmt = 8038 ) 4
719 WRITE( iounit, fmt = 8039 ) 5
720 WRITE( iounit, fmt = 8040 ) 6
721*
722 ELSE IF( lsamen( 2, p2, 'HH' ) ) THEN
723*
724* HH: Householder reconstruction for tall-skinny matrices
725*
726 WRITE( iounit, fmt = 8005 ) path
727 WRITE( iounit, fmt = '( '' Test ratios:'' )' )
728 WRITE( iounit, fmt = 8050 ) 1
729 WRITE( iounit, fmt = 8051 ) 2
730 WRITE( iounit, fmt = 8052 ) 3
731 WRITE( iounit, fmt = 8053 ) 4
732 WRITE( iounit, fmt = 8054 ) 5
733 WRITE( iounit, fmt = 8055 ) 6
734*
735 ELSE
736*
737* Print error message if no header is available.
738*
739 WRITE( iounit, fmt = 9980 )path
740 END IF
741*
742* First line of header
743*
744 9999 FORMAT( / 1x, a3, ': General dense matrices' )
745 9998 FORMAT( / 1x, a3, ': General band matrices' )
746 9997 FORMAT( / 1x, a3, ': General tridiagonal' )
747 9996 FORMAT( / 1x, a3, ': ', a9, ' positive definite matrices' )
748 9995 FORMAT( / 1x, a3, ': ', a9, ' positive definite packed matrices'
749 $ )
750 9994 FORMAT( / 1x, a3, ': ', a9, ' positive definite band matrices' )
751 9993 FORMAT( / 1x, a3, ': ', a9, ' positive definite tridiagonal' )
752 9992 FORMAT( / 1x, a3, ': ', a9, ' indefinite matrices',
753 $ ', partial (Bunch-Kaufman) pivoting' )
754 9991 FORMAT( / 1x, a3, ': ', a9, ' indefinite packed matrices',
755 $ ', partial (Bunch-Kaufman) pivoting' )
756 9892 FORMAT( / 1x, a3, ': ', a9, ' indefinite matrices',
757 $ ', "rook" (bounded Bunch-Kaufman) pivoting' )
758 9891 FORMAT( / 1x, a3, ': ', a9, ' indefinite packed matrices',
759 $ ', "rook" (bounded Bunch-Kaufman) pivoting' )
760 9990 FORMAT( / 1x, a3, ': Triangular matrices' )
761 9989 FORMAT( / 1x, a3, ': Triangular packed matrices' )
762 9988 FORMAT( / 1x, a3, ': Triangular band matrices' )
763 9987 FORMAT( / 1x, a3, ': ', a2, ' factorization of general matrices'
764 $ )
765 9986 FORMAT( / 1x, a3, ': QR factorization with column pivoting' )
766 9985 FORMAT( / 1x, a3, ': RQ factorization of trapezoidal matrix' )
767 9984 FORMAT( / 1x, a3, ': Least squares driver routines' )
768 9983 FORMAT( / 1x, a3, ': LU factorization variants' )
769 9982 FORMAT( / 1x, a3, ': Cholesky factorization variants' )
770 9981 FORMAT( / 1x, a3, ': QR factorization variants' )
771 9980 FORMAT( / 1x, a3, ': No header available' )
772 8000 FORMAT( / 1x, a3, ': QRT factorization for general matrices' )
773 8001 FORMAT( / 1x, a3, ': QRT factorization for ',
774 $ 'triangular-pentagonal matrices' )
775 8002 FORMAT( / 1x, a3, ': LQT factorization for general matrices' )
776 8003 FORMAT( / 1x, a3, ': LQT factorization for ',
777 $ 'triangular-pentagonal matrices' )
778 8004 FORMAT( / 1x, a3, ': TS factorization for ',
779 $ 'tall-skinny or short-wide matrices' )
780 8005 FORMAT( / 1x, a3, ': Householder recostruction from TSQR',
781 $ ' factorization output ', /,' for tall-skinny matrices.' )
782*
783* GE matrix types
784*
785 9979 FORMAT( 4x, '1. Diagonal', 24x, '7. Last n/2 columns zero', / 4x,
786 $ '2. Upper triangular', 16x,
787 $ '8. Random, CNDNUM = sqrt(0.1/EPS)', / 4x,
788 $ '3. Lower triangular', 16x, '9. Random, CNDNUM = 0.1/EPS',
789 $ / 4x, '4. Random, CNDNUM = 2', 13x,
790 $ '10. Scaled near underflow', / 4x, '5. First column zero',
791 $ 14x, '11. Scaled near overflow', / 4x,
792 $ '6. Last column zero' )
793*
794* GB matrix types
795*
796 9978 FORMAT( 4x, '1. Random, CNDNUM = 2', 14x,
797 $ '5. Random, CNDNUM = sqrt(0.1/EPS)', / 4x,
798 $ '2. First column zero', 15x, '6. Random, CNDNUM = .01/EPS',
799 $ / 4x, '3. Last column zero', 16x,
800 $ '7. Scaled near underflow', / 4x,
801 $ '4. Last n/2 columns zero', 11x, '8. Scaled near overflow' )
802*
803* GT matrix types
804*
805 9977 FORMAT( ' Matrix types (1-6 have specified condition numbers):',
806 $ / 4x, '1. Diagonal', 24x, '7. Random, unspecified CNDNUM',
807 $ / 4x, '2. Random, CNDNUM = 2', 14x, '8. First column zero',
808 $ / 4x, '3. Random, CNDNUM = sqrt(0.1/EPS)', 2x,
809 $ '9. Last column zero', / 4x, '4. Random, CNDNUM = 0.1/EPS',
810 $ 7x, '10. Last n/2 columns zero', / 4x,
811 $ '5. Scaled near underflow', 10x,
812 $ '11. Scaled near underflow', / 4x,
813 $ '6. Scaled near overflow', 11x, '12. Scaled near overflow' )
814*
815* PT matrix types
816*
817 9976 FORMAT( ' Matrix types (1-6 have specified condition numbers):',
818 $ / 4x, '1. Diagonal', 24x, '7. Random, unspecified CNDNUM',
819 $ / 4x, '2. Random, CNDNUM = 2', 14x,
820 $ '8. First row and column zero', / 4x,
821 $ '3. Random, CNDNUM = sqrt(0.1/EPS)', 2x,
822 $ '9. Last row and column zero', / 4x,
823 $ '4. Random, CNDNUM = 0.1/EPS', 7x,
824 $ '10. Middle row and column zero', / 4x,
825 $ '5. Scaled near underflow', 10x,
826 $ '11. Scaled near underflow', / 4x,
827 $ '6. Scaled near overflow', 11x, '12. Scaled near overflow' )
828*
829* PO, PP matrix types
830*
831 9975 FORMAT( 4x, '1. Diagonal', 24x,
832 $ '6. Random, CNDNUM = sqrt(0.1/EPS)', / 4x,
833 $ '2. Random, CNDNUM = 2', 14x, '7. Random, CNDNUM = 0.1/EPS',
834 $ / 3x, '*3. First row and column zero', 7x,
835 $ '8. Scaled near underflow', / 3x,
836 $ '*4. Last row and column zero', 8x,
837 $ '9. Scaled near overflow', / 3x,
838 $ '*5. Middle row and column zero', / 3x,
839 $ '(* - tests error exits from ', a3,
840 $ 'TRF, no test ratios are computed)' )
841*
842* CH matrix types
843*
844 9974 FORMAT( 4x, '1. Diagonal', 24x,
845 $ '6. Random, CNDNUM = sqrt(0.1/EPS)', / 4x,
846 $ '2. Random, CNDNUM = 2', 14x, '7. Random, CNDNUM = 0.1/EPS',
847 $ / 3x, '*3. First row and column zero', 7x,
848 $ '8. Scaled near underflow', / 3x,
849 $ '*4. Last row and column zero', 8x,
850 $ '9. Scaled near overflow', / 3x,
851 $ '*5. Middle row and column zero', / 3x,
852 $ '(* - tests error exits, no test ratios are computed)' )
853*
854* PS matrix types
855*
856 8973 FORMAT( 4x, '1. Diagonal', / 4x, '2. Random, CNDNUM = 2', 14x,
857 $ / 3x, '*3. Nonzero eigenvalues of: D(1:RANK-1)=1 and ',
858 $ 'D(RANK) = 1.0/', a4, / 3x,
859 $ '*4. Nonzero eigenvalues of: D(1)=1 and ',
860 $ ' D(2:RANK) = 1.0/', a4, / 3x,
861 $ '*5. Nonzero eigenvalues of: D(I) = ', a4,
862 $ '**(-(I-1)/(RANK-1)) ', ' I=1:RANK', / 4x,
863 $ '6. Random, CNDNUM = sqrt(0.1/EPS)', / 4x,
864 $ '7. Random, CNDNUM = 0.1/EPS', / 4x,
865 $ '8. Scaled near underflow', / 4x, '9. Scaled near overflow',
866 $ / 3x, '(* - Semi-definite tests )' )
867 8972 FORMAT( 3x, 'RANK minus computed rank, returned by ', a, 'PSTRF' )
868*
869* PB matrix types
870*
871 9973 FORMAT( 4x, '1. Random, CNDNUM = 2', 14x,
872 $ '5. Random, CNDNUM = sqrt(0.1/EPS)', / 3x,
873 $ '*2. First row and column zero', 7x,
874 $ '6. Random, CNDNUM = 0.1/EPS', / 3x,
875 $ '*3. Last row and column zero', 8x,
876 $ '7. Scaled near underflow', / 3x,
877 $ '*4. Middle row and column zero', 6x,
878 $ '8. Scaled near overflow', / 3x,
879 $ '(* - tests error exits from ', a3,
880 $ 'TRF, no test ratios are computed)' )
881*
882* SSY, SSR, SSP, CHE, CHR, CHP matrix types
883*
884 9972 FORMAT( 4x, '1. Diagonal', 24x,
885 $ '6. Last n/2 rows and columns zero', / 4x,
886 $ '2. Random, CNDNUM = 2', 14x,
887 $ '7. Random, CNDNUM = sqrt(0.1/EPS)', / 4x,
888 $ '3. First row and column zero', 7x,
889 $ '8. Random, CNDNUM = 0.1/EPS', / 4x,
890 $ '4. Last row and column zero', 8x,
891 $ '9. Scaled near underflow', / 4x,
892 $ '5. Middle row and column zero', 5x,
893 $ '10. Scaled near overflow' )
894*
895* CSY, CSR, CSP matrix types
896*
897 9971 FORMAT( 4x, '1. Diagonal', 24x,
898 $ '7. Random, CNDNUM = sqrt(0.1/EPS)', / 4x,
899 $ '2. Random, CNDNUM = 2', 14x, '8. Random, CNDNUM = 0.1/EPS',
900 $ / 4x, '3. First row and column zero', 7x,
901 $ '9. Scaled near underflow', / 4x,
902 $ '4. Last row and column zero', 7x,
903 $ '10. Scaled near overflow', / 4x,
904 $ '5. Middle row and column zero', 5x,
905 $ '11. Block diagonal matrix', / 4x,
906 $ '6. Last n/2 rows and columns zero' )
907*
908* QR matrix types
909*
910 9970 FORMAT( 4x, '1. Diagonal', 24x,
911 $ '5. Random, CNDNUM = sqrt(0.1/EPS)', / 4x,
912 $ '2. Upper triangular', 16x, '6. Random, CNDNUM = 0.1/EPS',
913 $ / 4x, '3. Lower triangular', 16x,
914 $ '7. Scaled near underflow', / 4x, '4. Random, CNDNUM = 2',
915 $ 14x, '8. Scaled near overflow' )
916*
917* QP matrix types
918*
919 9969 FORMAT( ' Matrix types (2-6 have condition 1/EPS):', / 4x,
920 $ '1. Zero matrix', 21x, '4. First n/2 columns fixed', / 4x,
921 $ '2. One small eigenvalue', 12x, '5. Last n/2 columns fixed',
922 $ / 4x, '3. Geometric distribution', 10x,
923 $ '6. Every second column fixed' )
924*
925* TZ matrix types
926*
927 9968 FORMAT( ' Matrix types (2-3 have condition 1/EPS):', / 4x,
928 $ '1. Zero matrix', / 4x, '2. One small eigenvalue', / 4x,
929 $ '3. Geometric distribution' )
930*
931* LS matrix types
932*
933 9967 FORMAT( ' Matrix types (1-3: full rank, 4-6: rank deficient):',
934 $ / 4x, '1 and 4. Normal scaling', / 4x,
935 $ '2 and 5. Scaled near overflow', / 4x,
936 $ '3 and 6. Scaled near underflow' )
937*
938* TR, TP matrix types
939*
940 9966 FORMAT( ' Matrix types for ', a3, ' routines:', / 4x,
941 $ '1. Diagonal', 24x, '6. Scaled near overflow', / 4x,
942 $ '2. Random, CNDNUM = 2', 14x, '7. Identity', / 4x,
943 $ '3. Random, CNDNUM = sqrt(0.1/EPS) ',
944 $ '8. Unit triangular, CNDNUM = 2', / 4x,
945 $ '4. Random, CNDNUM = 0.1/EPS', 8x,
946 $ '9. Unit, CNDNUM = sqrt(0.1/EPS)', / 4x,
947 $ '5. Scaled near underflow', 10x,
948 $ '10. Unit, CNDNUM = 0.1/EPS' )
949 9965 FORMAT( ' Special types for testing ', a, ':', / 3x,
950 $ '11. Matrix elements are O(1), large right hand side', / 3x,
951 $ '12. First diagonal causes overflow,',
952 $ ' offdiagonal column norms < 1', / 3x,
953 $ '13. First diagonal causes overflow,',
954 $ ' offdiagonal column norms > 1', / 3x,
955 $ '14. Growth factor underflows, solution does not overflow',
956 $ / 3x, '15. Small diagonal causes gradual overflow', / 3x,
957 $ '16. One zero diagonal element', / 3x,
958 $ '17. Large offdiagonals cause overflow when adding a column'
959 $ , / 3x, '18. Unit triangular with large right hand side' )
960*
961* TB matrix types
962*
963 9964 FORMAT( ' Matrix types for ', a3, ' routines:', / 4x,
964 $ '1. Random, CNDNUM = 2', 14x, '6. Identity', / 4x,
965 $ '2. Random, CNDNUM = sqrt(0.1/EPS) ',
966 $ '7. Unit triangular, CNDNUM = 2', / 4x,
967 $ '3. Random, CNDNUM = 0.1/EPS', 8x,
968 $ '8. Unit, CNDNUM = sqrt(0.1/EPS)', / 4x,
969 $ '4. Scaled near underflow', 11x,
970 $ '9. Unit, CNDNUM = 0.1/EPS', / 4x,
971 $ '5. Scaled near overflow' )
972 9963 FORMAT( ' Special types for testing ', a, ':', / 3x,
973 $ '10. Matrix elements are O(1), large right hand side', / 3x,
974 $ '11. First diagonal causes overflow,',
975 $ ' offdiagonal column norms < 1', / 3x,
976 $ '12. First diagonal causes overflow,',
977 $ ' offdiagonal column norms > 1', / 3x,
978 $ '13. Growth factor underflows, solution does not overflow',
979 $ / 3x, '14. Small diagonal causes gradual overflow', / 3x,
980 $ '15. One zero diagonal element', / 3x,
981 $ '16. Large offdiagonals cause overflow when adding a column'
982 $ , / 3x, '17. Unit triangular with large right hand side' )
983*
984* Test ratios
985*
986 9962 FORMAT( 3x, i2, ': norm( L * U - A ) / ( N * norm(A) * EPS )' )
987 9961 FORMAT( 3x, i2, ': norm( I - A*AINV ) / ',
988 $ '( N * norm(A) * norm(AINV) * EPS )' )
989 9960 FORMAT( 3x, i2, ': norm( B - A * X ) / ',
990 $ '( norm(A) * norm(X) * EPS )' )
991 6660 FORMAT( 3x, i2, ': diagonal is not non-negative')
992 9959 FORMAT( 3x, i2, ': norm( X - XACT ) / ',
993 $ '( norm(XACT) * CNDNUM * EPS )' )
994 9958 FORMAT( 3x, i2, ': norm( X - XACT ) / ',
995 $ '( norm(XACT) * CNDNUM * EPS ), refined' )
996 9957 FORMAT( 3x, i2, ': norm( X - XACT ) / ',
997 $ '( norm(XACT) * (error bound) )' )
998 9956 FORMAT( 3x, i2, ': (backward error) / EPS' )
999 9955 FORMAT( 3x, i2, ': RCOND * CNDNUM - 1.0' )
1000 9954 FORMAT( 3x, i2, ': norm( U'' * U - A ) / ( N * norm(A) * EPS )',
1001 $ ', or', / 7x, 'norm( L * L'' - A ) / ( N * norm(A) * EPS )'
1002 $ )
1003 8950 FORMAT( 3x,
1004 $ 'norm( P * U'' * U * P'' - A ) / ( N * norm(A) * EPS )',
1005 $ ', or', / 3x,
1006 $ 'norm( P * L * L'' * P'' - A ) / ( N * norm(A) * EPS )' )
1007 9953 FORMAT( 3x, i2, ': norm( U*D*U'' - A ) / ( N * norm(A) * EPS )',
1008 $ ', or', / 7x, 'norm( L*D*L'' - A ) / ( N * norm(A) * EPS )'
1009 $ )
1010 9952 FORMAT( 3x, i2, ': norm( U''*D*U - A ) / ( N * norm(A) * EPS )',
1011 $ ', or', / 7x, 'norm( L*D*L'' - A ) / ( N * norm(A) * EPS )'
1012 $ )
1013 9951 FORMAT( ' Test ratio for ', a, ':', / 3x, i2,
1014 $ ': norm( s*b - A*x ) / ( norm(A) * norm(x) * EPS )' )
1015 9950 FORMAT( 3x, i2, ': norm( R - Q'' * A ) / ( M * norm(A) * EPS )' )
1016 6950 FORMAT( 3x, i2, ': norm( R - Q'' * A ) / ( M * norm(A) * EPS )
1017 $ [RFPG]' )
1018 9949 FORMAT( 3x, i2, ': norm( L - A * Q'' ) / ( N * norm(A) * EPS )' )
1019 9948 FORMAT( 3x, i2, ': norm( L - Q'' * A ) / ( M * norm(A) * EPS )' )
1020 9947 FORMAT( 3x, i2, ': norm( R - A * Q'' ) / ( N * norm(A) * EPS )' )
1021 9946 FORMAT( 3x, i2, ': norm( I - Q''*Q ) / ( M * EPS )' )
1022 9945 FORMAT( 3x, i2, ': norm( I - Q*Q'' ) / ( N * EPS )' )
1023 9944 FORMAT( 3x, i2, ': norm( Q*C - Q*C ) / ', '( ', a1,
1024 $ ' * norm(C) * EPS )' )
1025 9943 FORMAT( 3x, i2, ': norm( C*Q - C*Q ) / ', '( ', a1,
1026 $ ' * norm(C) * EPS )' )
1027 9942 FORMAT( 3x, i2, ': norm( Q''*C - Q''*C )/ ', '( ', a1,
1028 $ ' * norm(C) * EPS )' )
1029 9941 FORMAT( 3x, i2, ': norm( C*Q'' - C*Q'' )/ ', '( ', a1,
1030 $ ' * norm(C) * EPS )' )
1031 9940 FORMAT( 3x, i2, ': norm(svd(A) - svd(R)) / ',
1032 $ '( M * norm(svd(R)) * EPS )' )
1033 9939 FORMAT( 3x, i2, ': norm( A*P - Q*R ) / ( M * norm(A) * EPS )'
1034 $ )
1035 9938 FORMAT( 3x, i2, ': norm( I - Q''*Q ) / ( M * EPS )' )
1036 9937 FORMAT( 3x, i2, ': norm( A - R*Q ) / ( M * norm(A) * EPS )'
1037 $ )
1038 9935 FORMAT( 3x, i2, ': norm( B - A * X ) / ',
1039 $ '( max(M,N) * norm(A) * norm(X) * EPS )' )
1040 9934 FORMAT( 3x, i2, ': norm( (A*X-B)'' *A ) / ',
1041 $ '( max(M,N,NRHS) * norm(A) * norm(B) * EPS )' )
1042 9933 FORMAT( 3x, i2, ': norm(svd(A)-svd(R)) / ',
1043 $ '( min(M,N) * norm(svd(R)) * EPS )' )
1044 9932 FORMAT( 3x, i2, ': Check if X is in the row space of A or A''' )
1045 9931 FORMAT( 3x, i2, ': norm( (A*X-B)'' *A ) / ',
1046 $ '( max(M,N,NRHS) * norm(A) * norm(B) * EPS )', / 7x,
1047 $ 'if TRANS=''N'.GE.' and MN or TRANS=''T'.LT.' and MN, ',
1048 $ 'otherwise', / 7x,
1049 $ 'check if X is in the row space of A or A'' ',
1050 $ '(overdetermined case)' )
1051 9929 FORMAT( ' Test ratios (1-3: ', a1, 'TZRZF):' )
1052 9919 FORMAT( 3x, ' 3-4: same as 1-2', 3x, ' 5-6: same as 1-2' )
1053 9920 FORMAT( 3x, ' 11-14: same as 7-10', 3x, ' 15-18: same as 7-10' )
1054 9921 FORMAT( ' Test ratios:', / ' (1-2: ', a1, 'GELS, 3-4: ', a1,
1055 $ 'GELST, 5-6: ', a1, 'GETSLS, 7-10: ', a1, 'GELSY, 11-14: ',
1056 $ a1, 'GETSS, 15-18: ', a1, 'GELSD)' )
1057 9928 FORMAT( 7x, 'where ALPHA = ( 1 + SQRT( 17 ) ) / 8' )
1058 9927 FORMAT( 3x, i2, ': ABS( Largest element in L )', / 12x,
1059 $ ' - ( 1 / ( 1 - ALPHA ) ) + THRESH' )
1060 9926 FORMAT( 3x, i2, ': Largest 2-Norm of 2-by-2 pivots', / 12x,
1061 $ ' - ( ( 1 + ALPHA ) / ( 1 - ALPHA ) ) + THRESH' )
1062 8011 FORMAT(3x,i2,': norm( R - Q''*A ) / ( M * norm(A) * EPS )' )
1063 8012 FORMAT(3x,i2,': norm( I - Q''*Q ) / ( M * EPS )' )
1064 8013 FORMAT(3x,i2,': norm( Q*C - Q*C ) / ( M * norm(C) * EPS )' )
1065 8014 FORMAT(3x,i2,': norm( Q''*C - Q''*C ) / ( M * norm(C) * EPS )')
1066 8015 FORMAT(3x,i2,': norm( C*Q - C*Q ) / ( M * norm(C) * EPS )' )
1067 8016 FORMAT(3x,i2,': norm( C*Q'' - C*Q'' ) / ( M * norm(C) * EPS )')
1068 8017 FORMAT(3x,i2,': norm( R - Q''*A ) / ( (M+N) * norm(A) * EPS )' )
1069 8018 FORMAT(3x,i2,': norm( I - Q''*Q ) / ( (M+N) * EPS )' )
1070 8019 FORMAT(3x,i2,': norm( Q*C - Q*C ) / ( (M+N) * norm(C) * EPS )' )
1071 8020 FORMAT(3x,i2,
1072 $ ': norm( Q''*C - Q''*C ) / ( (M+N) * norm(C) * EPS )')
1073 8021 FORMAT(3x,i2,': norm( C*Q - C*Q ) / ( (M+N) * norm(C) * EPS )' )
1074 8022 FORMAT(3x,i2,
1075 $ ': norm( C*Q'' - C*Q'' ) / ( (M+N) * norm(C) * EPS )')
1076 8023 FORMAT(3x,i2,': norm( L - A*Q'' ) / ( (M+N) * norm(A) * EPS )' )
1077 8024 FORMAT(3x,i2,': norm( I - Q*Q'' ) / ( (M+N) * EPS )' )
1078 8025 FORMAT(3x,i2,': norm( Q*C - Q*C ) / ( (M+N) * norm(C) * EPS )' )
1079 8026 FORMAT(3x,i2,
1080 $ ': norm( Q''*C - Q''*C ) / ( (M+N) * norm(C) * EPS )')
1081 8027 FORMAT(3x,i2,': norm( C*Q - C*Q ) / ( (M+N) * norm(C) * EPS )' )
1082 8028 FORMAT(3x,i2,
1083 $ ': norm( C*Q'' - C*Q'' ) / ( (M+N) * norm(C) * EPS )')
1084 8029 FORMAT(3x,i2,': norm( L - A*Q'' ) / ( (M+N) * norm(A) * EPS )' )
1085 8030 FORMAT(3x,i2,': norm( I - Q*Q'' ) / ( (M+N) * EPS )' )
1086 8031 FORMAT(3x,i2,': norm( Q*C - Q*C ) / ( (M+N) * norm(C) * EPS )' )
1087 8032 FORMAT(3x,i2,
1088 $ ': norm( Q''*C - Q''*C ) / ( (M+N) * norm(C) * EPS )')
1089 8033 FORMAT(3x,i2,': norm( C*Q - C*Q ) / ( (M+N) * norm(C) * EPS )' )
1090 8034 FORMAT(3x,i2,
1091 $ ': norm( C*Q'' - C*Q'' ) / ( (M+N) * norm(C) * EPS )')
1092 8035 FORMAT(3x,i2,': norm( R - Q''*A ) / ( (M+N) * norm(A) * EPS )' )
1093 8036 FORMAT(3x,i2,': norm( I - Q''*Q ) / ( (M+N) * EPS )' )
1094 8037 FORMAT(3x,i2,': norm( Q*C - Q*C ) / ( (M+N) * norm(C) * EPS )' )
1095 8038 FORMAT(3x,i2,
1096 $ ': norm( Q''*C - Q''*C ) / ( (M+N) * norm(C) * EPS )')
1097 8039 FORMAT(3x,i2,': norm( C*Q - C*Q ) / ( (M+N) * norm(C) * EPS )' )
1098 8040 FORMAT(3x,i2,
1099 $ ': norm( C*Q'' - C*Q'' ) / ( (M+N) * norm(C) * EPS )')
1100*
1101 8050 FORMAT(3x,i2,': norm( R - Q''*A ) / ( M * norm(A) * EPS )' )
1102 8051 FORMAT(3x,i2,': norm( I - Q''*Q ) / ( M * EPS )' )
1103 8052 FORMAT(3x,i2,': norm( Q*C - Q*C ) / ( M * norm(C) * EPS )' )
1104 8053 FORMAT(3x,i2,': norm( Q''*C - Q''*C ) / ( M * norm(C) * EPS )')
1105 8054 FORMAT(3x,i2,': norm( C*Q - C*Q ) / ( M * norm(C) * EPS )' )
1106 8055 FORMAT(3x,i2,': norm( C*Q'' - C*Q'' ) / ( M * norm(C) * EPS )')
1107
1108*
1109 RETURN
1110*
1111* End of ALAHD
1112*
1113 END
alahd
subroutine alahd(IOUNIT, PATH)
ALAHD
Definition: alahd.f:107