126 SUBROUTINE csyt03( UPLO, N, A, LDA, AINV, LDAINV, WORK, LDWORK,
127 $ rwork, rcond, resid )
136 INTEGER lda, ldainv, ldwork, n
141 COMPLEX a( lda, * ), ainv( ldainv, * ),
150 parameter( zero = 0.0e+0, one = 1.0e+0 )
152 parameter( czero = ( 0.0e+0, 0.0e+0 ),
153 $ cone = ( 1.0e+0, 0.0e+0 ) )
157 REAL ainvnm, anorm, eps
183 anorm =
clansy(
'1', uplo, n, a, lda, rwork )
184 ainvnm =
clansy(
'1', uplo, n, ainv, ldainv, rwork )
185 IF( anorm.LE.zero .OR. ainvnm.LE.zero )
THEN
190 rcond = ( one/anorm ) / ainvnm
195 IF(
lsame( uplo,
'U' ) )
THEN
198 ainv(
j, i ) = ainv( i,
j )
204 ainv(
j, i ) = ainv( i,
j )
208 CALL
csymm(
'Left', uplo, n, n, -cone, a, lda, ainv, ldainv,
209 $ czero, work, ldwork )
214 work( i, i ) = work( i, i ) + cone
219 resid =
clange(
'1', n, n, work, ldwork, rwork )
221 resid = ( ( resid*rcond )/eps ) /
REAL( n )
LOGICAL function lsame(CA, CB)
LSAME
REAL function slamch(CMACH)
SLAMCH
subroutine csyt03(UPLO, N, A, LDA, AINV, LDAINV, WORK, LDWORK, RWORK, RCOND, RESID)
CSYT03
REAL function clange(NORM, M, N, A, LDA, WORK)
CLANGE returns the value of the 1-norm, Frobenius norm, infinity-norm, or the largest absolute value ...
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
REAL function clansy(NORM, UPLO, N, A, LDA, WORK)
CLANSY returns the value of the 1-norm, or the Frobenius norm, or the infinity norm, or the element of largest absolute value of a complex symmetric matrix.
subroutine csymm(SIDE, UPLO, M, N, ALPHA, A, LDA, B, LDB, BETA, C, LDC)
CSYMM