101 DOUBLE PRECISION FUNCTION dlanst( NORM, N, D, E )
113 DOUBLE PRECISION d( * ), e( * )
119 DOUBLE PRECISION one, zero
120 parameter( one = 1.0d+0, zero = 0.0d+0 )
124 DOUBLE PRECISION anorm, scale, sum
134 INTRINSIC abs, max, sqrt
140 ELSE IF( lsame( norm,
'M' ) )
THEN 144 anorm = abs( d( n ) )
146 anorm = max( anorm, abs( d( i ) ) )
147 anorm = max( anorm, abs( e( i ) ) )
149 ELSE IF( lsame( norm,
'O' ) .OR. norm.EQ.
'1' .OR.
150 $ lsame( norm,
'I' ) )
THEN 155 anorm = abs( d( 1 ) )
157 anorm = max( abs( d( 1 ) )+abs( e( 1 ) ),
158 $ abs( e( n-1 ) )+abs( d( n ) ) )
160 anorm = max( anorm, abs( d( i ) )+abs( e( i ) )+
164 ELSE IF( ( lsame( norm,
'F' ) ) .OR. ( lsame( norm,
'E' ) ) )
THEN 171 CALL dlassq( n-1, e, 1, scale, sum )
174 CALL dlassq( n, d, 1, scale, sum )
175 anorm = scale*sqrt( sum )