subroutine daxpy0(n,da,dx,incx,dy,incy)
c
c     constant times a vector plus a vector.
! simple version of daxpy written by Greg
c
      real dx(n),dy(n),da
      integer i,incx,incy,ix,iy,m,mp1,n
c
      if(n.le.0)return
      if (da .eq. 0.0d0) return

      do 10 i = 1,n
        dy(i) = dy(i) + da*dx(i)
   10 continue
      return
      end

      subroutine daxpy2(n,da,dx,incx,dy,incy)
c
c     constant times a vector plus a vector.
! simple version of daxpy written by Greg, unrolled by 2
c
      real dx(n),dy(n),da
      integer i,incx,incy,ix,iy,m,mp1,n
c
      if(n.le.0)return
      if (da .eq. 0.0d0) return

      do 10 i = 1,n,2
        dy(i) = dy(i) + da*dx(i)
        dy(i+1) = dy(i+1) + da*dx(i+1)
   10 continue
      return
      end

      subroutine daxpy4(n,da,dx,incx,dy,incy)
c
c     constant times a vector plus a vector.
! simple version of daxpy written by Greg, unrolled by 4
c
      real dx(n),dy(n),da
      integer i,incx,incy,ix,iy,m,mp1,n
c
      if(n.le.0)return
      if (da .eq. 0.0d0) return

      do 10 i = 1,n,4
        dy(i) = dy(i) + da*dx(i)
        dy(i+1) = dy(i+1) + da*dx(i+1)
        dy(i+2) = dy(i+2) + da*dx(i+2)
        dy(i+3) = dy(i+3) + da*dx(i+3)
   10 continue
      return
      end

      subroutine daxpy10(n,da,dx,incx,dy,incy)
c
c     constant times a vector plus a vector.
! simple version of daxpy written by Greg, unrolled by 10
c
      real dx(n),dy(n),da
      integer i,incx,incy,ix,iy,m,mp1,n
c
      if(n.le.0)return
      if (da .eq. 0.0d0) return

      do 10 i = 1,n,10
        dy(i) = dy(i) + da*dx(i)
        dy(i+1) = dy(i+1) + da*dx(i+1)
        dy(i+2) = dy(i+2) + da*dx(i+2)
        dy(i+3) = dy(i+3) + da*dx(i+3)
        dy(i+4) = dy(i+4) + da*dx(i+4)
        dy(i+5) = dy(i+5) + da*dx(i+5)
        dy(i+6) = dy(i+6) + da*dx(i+6)
        dy(i+7) = dy(i+7) + da*dx(i+7)
        dy(i+8) = dy(i+8) + da*dx(i+8)
        dy(i+9) = dy(i+9) + da*dx(i+9)
   10 continue
      return
      end