1      subroutine yaxpy(n,da,dx,incx,dy,incy)
2*
3* $Id$
4*
5      implicit none
6#include "y64.fh"
7      double precision dx(*),dy(*),da
8      integer incx,incy,n
9c
10      INTGR4 incx4,incy4,n4
11c
12      if(n.le.0)return
13      if (da .eq. 0.0d0) return
14#ifdef USE_YFLOP
15      flops_ycount = flops_ycount + n*2
16#endif
17      incx4=incx
18      incy4=incy
19      n4=n
20c
21      call daxpy(n4,da,dx,incx4,dy,incy4)
22      return
23      end
24