161 SUBROUTINE sspgv( ITYPE, JOBZ, UPLO, N, AP, BP, W, Z, LDZ, WORK,
171 INTEGER info, itype, ldz, n
174 REAL ap( * ), bp( * ), w( * ), work( * ),
196 wantz =
lsame( jobz,
'V' )
197 upper =
lsame( uplo,
'U' )
200 IF( itype.LT.1 .OR. itype.GT.3 )
THEN
202 ELSE IF( .NOT.( wantz .OR.
lsame( jobz,
'N' ) ) )
THEN
204 ELSE IF( .NOT.( upper .OR.
lsame( uplo,
'L' ) ) )
THEN
206 ELSE IF( n.LT.0 )
THEN
208 ELSE IF( ldz.LT.1 .OR. ( wantz .AND. ldz.LT.n ) )
THEN
212 CALL
xerbla(
'SSPGV ', -info )
223 CALL
spptrf( uplo, n, bp, info )
231 CALL
sspgst( itype, uplo, n, ap, bp, info )
232 CALL
sspev( jobz, uplo, n, ap, w, z, ldz, work, info )
241 IF( itype.EQ.1 .OR. itype.EQ.2 )
THEN
253 CALL
stpsv( uplo, trans,
'Non-unit', n, bp, z( 1,
j ),
257 ELSE IF( itype.EQ.3 )
THEN
269 CALL
stpmv( uplo, trans,
'Non-unit', n, bp, z( 1,
j ),
subroutine sspgv(ITYPE, JOBZ, UPLO, N, AP, BP, W, Z, LDZ, WORK, INFO)
SSPGST
subroutine sspgst(ITYPE, UPLO, N, AP, BP, INFO)
SSPGST
subroutine xerbla(SRNAME, INFO)
XERBLA
logical function lsame(CA, CB)
LSAME
subroutine stpmv(UPLO, TRANS, DIAG, N, AP, X, INCX)
STPMV
subroutine stpsv(UPLO, TRANS, DIAG, N, AP, X, INCX)
STPSV
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
subroutine sspev(JOBZ, UPLO, N, AP, W, Z, LDZ, WORK, INFO)
SSPEV computes the eigenvalues and, optionally, the left and/or right eigenvectors for OTHER matrice...
subroutine spptrf(UPLO, N, AP, INFO)
SPPTRF