188 SUBROUTINE cgbmv(TRANS,M,N,KL,KU,ALPHA,A,LDA,X,INCX,BETA,Y,INCY)
197 INTEGER incx,incy,kl,ku,lda,m,n
201 COMPLEX a(lda,*),
x(*),y(*)
208 parameter(one= (1.0e+0,0.0e+0))
210 parameter(zero= (0.0e+0,0.0e+0))
214 INTEGER i,info,ix,iy,
j,jx,jy,k,kup1,kx,ky,lenx,leny
225 INTRINSIC conjg,max,min
231 IF (.NOT.
lsame(trans,
'N') .AND. .NOT.
lsame(trans,
'T') .AND.
232 + .NOT.
lsame(trans,
'C'))
THEN
234 ELSE IF (m.LT.0)
THEN
236 ELSE IF (n.LT.0)
THEN
238 ELSE IF (kl.LT.0)
THEN
240 ELSE IF (ku.LT.0)
THEN
242 ELSE IF (lda.LT. (kl+ku+1))
THEN
244 ELSE IF (incx.EQ.0)
THEN
246 ELSE IF (incy.EQ.0)
THEN
250 CALL
xerbla(
'CGBMV ',info)
256 IF ((m.EQ.0) .OR. (n.EQ.0) .OR.
257 + ((alpha.EQ.zero).AND. (beta.EQ.one)))
RETURN
259 noconj =
lsame(trans,
'T')
264 IF (
lsame(trans,
'N'))
THEN
274 kx = 1 - (lenx-1)*incx
279 ky = 1 - (leny-1)*incy
287 IF (beta.NE.one)
THEN
289 IF (beta.EQ.zero)
THEN
300 IF (beta.EQ.zero)
THEN
313 IF (alpha.EQ.zero)
RETURN
315 IF (
lsame(trans,
'N'))
THEN
322 IF (
x(jx).NE.zero)
THEN
325 DO 50 i = max(1,
j-ku),min(m,
j+kl)
326 y(i) = y(i) + temp*a(k+i,
j)
333 IF (
x(jx).NE.zero)
THEN
337 DO 70 i = max(1,
j-ku),min(m,
j+kl)
338 y(iy) = y(iy) + temp*a(k+i,
j)
343 IF (
j.GT.ku) ky = ky + incy
356 DO 90 i = max(1,
j-ku),min(m,
j+kl)
357 temp = temp + a(k+i,
j)*
x(i)
360 DO 100 i = max(1,
j-ku),min(m,
j+kl)
361 temp = temp + conjg(a(k+i,
j))*
x(i)
364 y(jy) = y(jy) + alpha*temp
373 DO 120 i = max(1,
j-ku),min(m,
j+kl)
374 temp = temp + a(k+i,
j)*
x(ix)
378 DO 130 i = max(1,
j-ku),min(m,
j+kl)
379 temp = temp + conjg(a(k+i,
j))*
x(ix)
383 y(jy) = y(jy) + alpha*temp
385 IF (
j.GT.ku) kx = kx + incx
LOGICAL function lsame(CA, CB)
LSAME
subroutine cgbmv(TRANS, M, N, KL, KU, ALPHA, A, LDA, X, INCX, BETA, Y, INCY)
CGBMV
subroutine xerbla(SRNAME, INFO)
XERBLA
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j