00001 SUBROUTINE CHBT21( UPLO, N, KA, KS, A, LDA, D, E, U, LDU, WORK,
00002 $ RWORK, RESULT )
00003
00004
00005
00006
00007
00008
00009 CHARACTER UPLO
00010 INTEGER KA, KS, LDA, LDU, N
00011
00012
00013 REAL D( * ), E( * ), RESULT( 2 ), RWORK( * )
00014 COMPLEX A( LDA, * ), U( LDU, * ), WORK( * )
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044
00045
00046
00047
00048
00049
00050
00051
00052
00053
00054
00055
00056
00057
00058
00059
00060
00061
00062
00063
00064
00065
00066
00067
00068
00069
00070
00071
00072
00073
00074
00075
00076
00077
00078
00079
00080
00081
00082
00083
00084
00085
00086
00087
00088
00089
00090
00091
00092
00093 COMPLEX CZERO, CONE
00094 PARAMETER ( CZERO = ( 0.0E+0, 0.0E+0 ),
00095 $ CONE = ( 1.0E+0, 0.0E+0 ) )
00096 REAL ZERO, ONE
00097 PARAMETER ( ZERO = 0.0E+0, ONE = 1.0E+0 )
00098
00099
00100 LOGICAL LOWER
00101 CHARACTER CUPLO
00102 INTEGER IKA, J, JC, JR
00103 REAL ANORM, ULP, UNFL, WNORM
00104
00105
00106 LOGICAL LSAME
00107 REAL CLANGE, CLANHB, CLANHP, SLAMCH
00108 EXTERNAL LSAME, CLANGE, CLANHB, CLANHP, SLAMCH
00109
00110
00111 EXTERNAL CGEMM, CHPR, CHPR2
00112
00113
00114 INTRINSIC CMPLX, MAX, MIN, REAL
00115
00116
00117
00118
00119
00120 RESULT( 1 ) = ZERO
00121 RESULT( 2 ) = ZERO
00122 IF( N.LE.0 )
00123 $ RETURN
00124
00125 IKA = MAX( 0, MIN( N-1, KA ) )
00126
00127 IF( LSAME( UPLO, 'U' ) ) THEN
00128 LOWER = .FALSE.
00129 CUPLO = 'U'
00130 ELSE
00131 LOWER = .TRUE.
00132 CUPLO = 'L'
00133 END IF
00134
00135 UNFL = SLAMCH( 'Safe minimum' )
00136 ULP = SLAMCH( 'Epsilon' )*SLAMCH( 'Base' )
00137
00138
00139
00140
00141
00142
00143
00144 ANORM = MAX( CLANHB( '1', CUPLO, N, IKA, A, LDA, RWORK ), UNFL )
00145
00146
00147
00148
00149
00150 J = 0
00151 DO 50 JC = 1, N
00152 IF( LOWER ) THEN
00153 DO 10 JR = 1, MIN( IKA+1, N+1-JC )
00154 J = J + 1
00155 WORK( J ) = A( JR, JC )
00156 10 CONTINUE
00157 DO 20 JR = IKA + 2, N + 1 - JC
00158 J = J + 1
00159 WORK( J ) = ZERO
00160 20 CONTINUE
00161 ELSE
00162 DO 30 JR = IKA + 2, JC
00163 J = J + 1
00164 WORK( J ) = ZERO
00165 30 CONTINUE
00166 DO 40 JR = MIN( IKA, JC-1 ), 0, -1
00167 J = J + 1
00168 WORK( J ) = A( IKA+1-JR, JC )
00169 40 CONTINUE
00170 END IF
00171 50 CONTINUE
00172
00173 DO 60 J = 1, N
00174 CALL CHPR( CUPLO, N, -D( J ), U( 1, J ), 1, WORK )
00175 60 CONTINUE
00176
00177 IF( N.GT.1 .AND. KS.EQ.1 ) THEN
00178 DO 70 J = 1, N - 1
00179 CALL CHPR2( CUPLO, N, -CMPLX( E( J ) ), U( 1, J ), 1,
00180 $ U( 1, J+1 ), 1, WORK )
00181 70 CONTINUE
00182 END IF
00183 WNORM = CLANHP( '1', CUPLO, N, WORK, RWORK )
00184
00185 IF( ANORM.GT.WNORM ) THEN
00186 RESULT( 1 ) = ( WNORM / ANORM ) / ( N*ULP )
00187 ELSE
00188 IF( ANORM.LT.ONE ) THEN
00189 RESULT( 1 ) = ( MIN( WNORM, N*ANORM ) / ANORM ) / ( N*ULP )
00190 ELSE
00191 RESULT( 1 ) = MIN( WNORM / ANORM, REAL( N ) ) / ( N*ULP )
00192 END IF
00193 END IF
00194
00195
00196
00197
00198
00199 CALL CGEMM( 'N', 'C', N, N, N, CONE, U, LDU, U, LDU, CZERO, WORK,
00200 $ N )
00201
00202 DO 80 J = 1, N
00203 WORK( ( N+1 )*( J-1 )+1 ) = WORK( ( N+1 )*( J-1 )+1 ) - CONE
00204 80 CONTINUE
00205
00206 RESULT( 2 ) = MIN( CLANGE( '1', N, N, WORK, N, RWORK ),
00207 $ REAL( N ) ) / ( N*ULP )
00208
00209 RETURN
00210
00211
00212
00213 END