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