1 SUBROUTINE clatb4( PATH, IMAT, M, N, TYPE, KL, KU, ANORM, MODE,
11 INTEGER imat, kl, ku, m, mode, n
68 parameter( shrink = 0.25e0, tenth = 0.1e+0 )
70 parameter( one = 1.0e+0 )
72 parameter( two = 2.0e+0 )
78 REAL badc1, badc2, eps, large, small
86 INTRINSIC abs,
max, sqrt
92 SAVE eps, small, large, badc1, badc2, first
103 eps = slamch(
'Precision' )
105 badc1 = sqrt( badc2 )
106 small = slamch(
'Safe minimum' )
112 CALL
slabad( small, large )
113 small = shrink*( small / eps )
127 IF(
lsamen( 2, c2,
'QR' ) .OR.
lsamen( 2, c2,
'LQ' ) .OR.
128 $
lsamen( 2, c2,
'QL' ) .OR.
lsamen( 2, c2,
'RQ' ) )
THEN
139 ELSE IF( imat.EQ.2 )
THEN
142 ELSE IF( imat.EQ.3 )
THEN
154 ELSE IF( imat.EQ.6 )
THEN
162 ELSE IF( imat.EQ.8 )
THEN
168 ELSE IF(
lsamen( 2, c2,
'GE' ) )
THEN
181 ELSE IF( imat.EQ.2 )
THEN
184 ELSE IF( imat.EQ.3 )
THEN
196 ELSE IF( imat.EQ.9 )
THEN
202 IF( imat.EQ.10 )
THEN
204 ELSE IF( imat.EQ.11 )
THEN
210 ELSE IF(
lsamen( 2, c2,
'GB' ) )
THEN
222 ELSE IF( imat.EQ.6 )
THEN
230 ELSE IF( imat.EQ.8 )
THEN
236 ELSE IF(
lsamen( 2, c2,
'GT' ) )
THEN
257 ELSE IF( imat.EQ.4 )
THEN
263 IF( imat.EQ.5 .OR. imat.EQ.11 )
THEN
265 ELSE IF( imat.EQ.6 .OR. imat.EQ.12 )
THEN
271 ELSE IF(
lsamen( 2, c2,
'PO' ) .OR.
lsamen( 2, c2,
'PP' ) .OR.
273 $
lsamen( 2, c2,
'SY' ) .OR.
lsamen( 2, c2,
'SP' ) )
THEN
295 ELSE IF( imat.EQ.7 )
THEN
303 ELSE IF( imat.EQ.9 )
THEN
309 ELSE IF(
lsamen( 2, c2,
'PB' ) )
THEN
321 ELSE IF( imat.EQ.6 )
THEN
329 ELSE IF( imat.EQ.8 )
THEN
335 ELSE IF(
lsamen( 2, c2,
'PT' ) )
THEN
352 ELSE IF( imat.EQ.4 )
THEN
358 IF( imat.EQ.5 .OR. imat.EQ.11 )
THEN
360 ELSE IF( imat.EQ.6 .OR. imat.EQ.12 )
THEN
366 ELSE IF(
lsamen( 2, c2,
'TR' ) .OR.
lsamen( 2, c2,
'TP' ) )
THEN
377 IF( mat.EQ.1 .OR. mat.EQ.7 )
THEN
380 ELSE IF( imat.LT.0 )
THEN
390 IF( mat.EQ.3 .OR. mat.EQ.9 )
THEN
392 ELSE IF( mat.EQ.4 .OR. mat.EQ.10 )
THEN
400 ELSE IF( mat.EQ.6 )
THEN
406 ELSE IF(
lsamen( 2, c2,
'TB' ) )
THEN
416 IF( imat.EQ.2 .OR. imat.EQ.8 )
THEN
418 ELSE IF( imat.EQ.3 .OR. imat.EQ.9 )
THEN
426 ELSE IF( imat.EQ.5 )
THEN