112 SUBROUTINE dlatb5( PATH, IMAT, N, TYPE, KL, KU, ANORM, MODE,
120 DOUBLE PRECISION ANORM, CNDNUM
121 INTEGER IMAT, KL, KU, MODE, N
129 DOUBLE PRECISION SHRINK, TENTH
130 parameter( shrink = 0.25d0, tenth = 0.1d+0 )
132 parameter( one = 1.0d+0 )
134 parameter( two = 2.0d+0 )
137 DOUBLE PRECISION BADC1, BADC2, EPS, LARGE, SMALL
142 DOUBLE PRECISION DLAMCH
152 SAVE eps, small, large, badc1, badc2, first
155 DATA first / .true. /
163 eps = dlamch(
'Precision' )
165 badc1 = sqrt( badc2 )
166 small = dlamch(
'Safe minimum' )
172 CALL dlabad( small, large )
173 small = shrink*( small / eps )
202 ELSE IF( imat.EQ.4 )
THEN
205 ELSE IF( imat.EQ.5 )
THEN
208 ELSE IF( imat.EQ.6 )
THEN
210 ELSE IF( imat.EQ.7 )
THEN
218 ELSE IF( imat.EQ.9 )
THEN
subroutine dlabad(SMALL, LARGE)
DLABAD
subroutine dlatb5(PATH, IMAT, N, TYPE, KL, KU, ANORM, MODE, CNDNUM, DIST)
DLATB5