185 SUBROUTINE ctbmv(UPLO,TRANS,DIAG,N,K,A,LDA,X,INCX)
193 CHARACTER DIAG,TRANS,UPLO
196 COMPLEX A(LDA,*),X(*)
203 parameter(zero= (0.0e+0,0.0e+0))
207 INTEGER I,INFO,IX,J,JX,KPLUS1,KX,L
208 LOGICAL NOCONJ,NOUNIT
218 INTRINSIC conjg,max,min
224 IF (.NOT.lsame(uplo,
'U') .AND. .NOT.lsame(uplo,
'L'))
THEN
226 ELSE IF (.NOT.lsame(trans,
'N') .AND. .NOT.lsame(trans,
'T') .AND.
227 + .NOT.lsame(trans,
'C'))
THEN
229 ELSE IF (.NOT.lsame(diag,
'U') .AND. .NOT.lsame(diag,
'N'))
THEN
231 ELSE IF (n.LT.0)
THEN
233 ELSE IF (k.LT.0)
THEN
235 ELSE IF (lda.LT. (k+1))
THEN
237 ELSE IF (incx.EQ.0)
THEN
241 CALL xerbla(
'CTBMV ',info)
249 noconj = lsame(trans,
'T')
250 nounit = lsame(diag,
'N')
257 ELSE IF (incx.NE.1)
THEN
264 IF (lsame(trans,
'N'))
THEN
268 IF (lsame(uplo,
'U'))
THEN
272 IF (x(j).NE.zero)
THEN
275 DO 10 i = max(1,j-k),j - 1
276 x(i) = x(i) + temp*a(l+i,j)
278 IF (nounit) x(j) = x(j)*a(kplus1,j)
284 IF (x(jx).NE.zero)
THEN
288 DO 30 i = max(1,j-k),j - 1
289 x(ix) = x(ix) + temp*a(l+i,j)
292 IF (nounit) x(jx) = x(jx)*a(kplus1,j)
295 IF (j.GT.k) kx = kx + incx
301 IF (x(j).NE.zero)
THEN
304 DO 50 i = min(n,j+k),j + 1,-1
305 x(i) = x(i) + temp*a(l+i,j)
307 IF (nounit) x(j) = x(j)*a(1,j)
314 IF (x(jx).NE.zero)
THEN
318 DO 70 i = min(n,j+k),j + 1,-1
319 x(ix) = x(ix) + temp*a(l+i,j)
322 IF (nounit) x(jx) = x(jx)*a(1,j)
325 IF ((n-j).GE.k) kx = kx - incx
333 IF (lsame(uplo,
'U'))
THEN
340 IF (nounit) temp = temp*a(kplus1,j)
341 DO 90 i = j - 1,max(1,j-k),-1
342 temp = temp + a(l+i,j)*x(i)
345 IF (nounit) temp = temp*conjg(a(kplus1,j))
346 DO 100 i = j - 1,max(1,j-k),-1
347 temp = temp + conjg(a(l+i,j))*x(i)
361 IF (nounit) temp = temp*a(kplus1,j)
362 DO 120 i = j - 1,max(1,j-k),-1
363 temp = temp + a(l+i,j)*x(ix)
367 IF (nounit) temp = temp*conjg(a(kplus1,j))
368 DO 130 i = j - 1,max(1,j-k),-1
369 temp = temp + conjg(a(l+i,j))*x(ix)
383 IF (nounit) temp = temp*a(1,j)
384 DO 150 i = j + 1,min(n,j+k)
385 temp = temp + a(l+i,j)*x(i)
388 IF (nounit) temp = temp*conjg(a(1,j))
389 DO 160 i = j + 1,min(n,j+k)
390 temp = temp + conjg(a(l+i,j))*x(i)
403 IF (nounit) temp = temp*a(1,j)
404 DO 180 i = j + 1,min(n,j+k)
405 temp = temp + a(l+i,j)*x(ix)
409 IF (nounit) temp = temp*conjg(a(1,j))
410 DO 190 i = j + 1,min(n,j+k)
411 temp = temp + conjg(a(l+i,j))*x(ix)
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine ctbmv(UPLO, TRANS, DIAG, N, K, A, LDA, X, INCX)
CTBMV