141 SUBROUTINE stpmv(UPLO,TRANS,DIAG,N,AP,X,INCX)
149 CHARACTER DIAG,TRANS,UPLO
159 parameter(zero=0.0e+0)
163 INTEGER I,INFO,IX,J,JX,K,KK,KX
177 IF (.NOT.lsame(uplo,
'U') .AND. .NOT.lsame(uplo,
'L'))
THEN
179 ELSE IF (.NOT.lsame(trans,
'N') .AND. .NOT.lsame(trans,
'T') .AND.
180 + .NOT.lsame(trans,
'C'))
THEN
182 ELSE IF (.NOT.lsame(diag,
'U') .AND. .NOT.lsame(diag,
'N'))
THEN
184 ELSE IF (n.LT.0)
THEN
186 ELSE IF (incx.EQ.0)
THEN
190 CALL xerbla(
'STPMV ',info)
198 nounit = lsame(diag,
'N')
205 ELSE IF (incx.NE.1)
THEN
212 IF (lsame(trans,
'N'))
THEN
216 IF (lsame(uplo,
'U'))
THEN
220 IF (x(j).NE.zero)
THEN
224 x(i) = x(i) + temp*ap(k)
227 IF (nounit) x(j) = x(j)*ap(kk+j-1)
234 IF (x(jx).NE.zero)
THEN
237 DO 30 k = kk,kk + j - 2
238 x(ix) = x(ix) + temp*ap(k)
241 IF (nounit) x(jx) = x(jx)*ap(kk+j-1)
251 IF (x(j).NE.zero)
THEN
255 x(i) = x(i) + temp*ap(k)
258 IF (nounit) x(j) = x(j)*ap(kk-n+j)
266 IF (x(jx).NE.zero)
THEN
269 DO 70 k = kk,kk - (n- (j+1)),-1
270 x(ix) = x(ix) + temp*ap(k)
273 IF (nounit) x(jx) = x(jx)*ap(kk-n+j)
284 IF (lsame(uplo,
'U'))
THEN
289 IF (nounit) temp = temp*ap(kk)
292 temp = temp + ap(k)*x(i)
303 IF (nounit) temp = temp*ap(kk)
304 DO 110 k = kk - 1,kk - j + 1,-1
306 temp = temp + ap(k)*x(ix)
318 IF (nounit) temp = temp*ap(kk)
321 temp = temp + ap(k)*x(i)
332 IF (nounit) temp = temp*ap(kk)
333 DO 150 k = kk + 1,kk + n - j
335 temp = temp + ap(k)*x(ix)
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine stpmv(UPLO, TRANS, DIAG, N, AP, X, INCX)
STPMV