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