186 SUBROUTINE dgbmv(TRANS,M,N,KL,KU,ALPHA,A,LDA,X,INCX,BETA,Y,INCY)
194 DOUBLE PRECISION ALPHA,BETA
195 INTEGER INCX,INCY,KL,KU,LDA,M,N
199 DOUBLE PRECISION A(lda,*),X(*),Y(*)
205 DOUBLE PRECISION ONE,ZERO
206 parameter(one=1.0d+0,zero=0.0d+0)
209 DOUBLE PRECISION TEMP
210 INTEGER I,INFO,IX,IY,J,JX,JY,K,KUP1,KX,KY,LENX,LENY
226 IF (.NOT.lsame(trans,
'N') .AND. .NOT.lsame(trans,
'T') .AND.
227 + .NOT.lsame(trans,
'C'))
THEN
229 ELSE IF (m.LT.0)
THEN
231 ELSE IF (n.LT.0)
THEN
233 ELSE IF (kl.LT.0)
THEN
235 ELSE IF (ku.LT.0)
THEN
237 ELSE IF (lda.LT. (kl+ku+1))
THEN
239 ELSE IF (incx.EQ.0)
THEN
241 ELSE IF (incy.EQ.0)
THEN
245 CALL xerbla(
'DGBMV ',info)
251 IF ((m.EQ.0) .OR. (n.EQ.0) .OR.
252 + ((alpha.EQ.zero).AND. (beta.EQ.one)))
RETURN
257 IF (lsame(trans,
'N'))
THEN
267 kx = 1 - (lenx-1)*incx
272 ky = 1 - (leny-1)*incy
280 IF (beta.NE.one)
THEN
282 IF (beta.EQ.zero)
THEN
293 IF (beta.EQ.zero)
THEN
306 IF (alpha.EQ.zero)
RETURN
308 IF (lsame(trans,
'N'))
THEN
317 DO 50 i = max(1,j-ku),min(m,j+kl)
318 y(i) = y(i) + temp*a(k+i,j)
327 DO 70 i = max(1,j-ku),min(m,j+kl)
328 y(iy) = y(iy) + temp*a(k+i,j)
332 IF (j.GT.ku) ky = ky + incy
344 DO 90 i = max(1,j-ku),min(m,j+kl)
345 temp = temp + a(k+i,j)*x(i)
347 y(jy) = y(jy) + alpha*temp
355 DO 110 i = max(1,j-ku),min(m,j+kl)
356 temp = temp + a(k+i,j)*x(ix)
359 y(jy) = y(jy) + alpha*temp
361 IF (j.GT.ku) kx = kx + incx
subroutine dgbmv(TRANS, M, N, KL, KU, ALPHA, A, LDA, X, INCX, BETA, Y, INCY)
DGBMV
subroutine xerbla(SRNAME, INFO)
XERBLA