1 subroutine drot (n,dx,incx,dy,incy,c,s)
7 double precision dx(*),dy(*),dtemp,c,s
8 integer i,incx,incy,ix,iy,n
11 if(incx.eq.1.and.incy.eq.1)
go to 20
18 if(incx.lt.0)ix = (-n+1)*incx + 1
19 if(incy.lt.0)iy = (-n+1)*incy + 1
21 dtemp = c*dx(ix) + s*dy(iy)
22 dy(iy) = c*dy(iy) - s*dx(ix)
32 dtemp = c*dx(i) + s*dy(i)
33 dy(i) = c*dy(i) - s*dx(i)
subroutine drot(n, dx, incx, dy, incy, c, s)