5 SUBROUTINE SSPR(UPLO,N,ALPHA,X,INCX,AP)
11 PARAMETER (ZERO=0.0E+0)
13 INTEGER I,INFO,IX,J,JX,K,KK,KX
19 IF (.NOT.LSAME(UPLO,'U') .AND. .NOT.LSAME(UPLO,'L')) THEN
23 ELSE IF (INCX.EQ.0) THEN
27 CALL XERBLA('SSPR ',INFO)
30 IF ((N.EQ.0) .OR. (ALPHA.EQ.ZERO)) RETURN
33 ELSE IF (INCX.NE.1) THEN
37 IF (LSAME(UPLO,'U')) THEN
40 IF (X(J).NE.ZERO) THEN
44 AP(K) = AP(K) + X(I)*TEMP
53 IF (X(JX).NE.ZERO) THEN
56 DO 30 K = KK,KK + J - 1
57 AP(K) = AP(K) + X(IX)*TEMP
68 IF (X(J).NE.ZERO) THEN
72 AP(K) = AP(K) + X(I)*TEMP
81 IF (X(JX).NE.ZERO) THEN
84 DO 70 K = KK,KK + N - J
85 AP(K) = AP(K) + X(IX)*TEMP