------------------------------------------------------------------------
IMSL Name:  DAXPY (Double precision version)
 
Revised:    August 9, 1986
 
Purpose:    Compute the scalar times a vector plus a vector,
            y = ax + y, all double precision.
 
Usage:      CALL DAXPY (N, DA, DX, INCX, DY, INCY)
 
Arguments:
   N      - Length of vectors X and Y.  (Input)
   DA     - Double precision scalar.  (Input)
   DX     - Double precision vector of length MAX(N*IABS(INCX),1).
               (Input)
   INCX   - Displacement between elements of DX.  (Input)
            X(I) is defined to be
               DX(1+(I-1)*INCX) if INCX.GE.0  or
               DX(1+(I-N)*INCX) if INCX.LT.0.
   DY     - Double precision vector of length MAX(N*IABS(INCY),1).
               (Input/Output)
            DAXPY replaces Y(I) with DA*X(I) + Y(I) for I=1,...,N.
            X(I) and Y(I) refer to specific elements of DX and DY.
   INCY   - Displacement between elements of DY.  (Input)
            Y(I) is defined to be
               DY(1+(I-1)*INCY) if INCY.GE.0  or
               DY(1+(I-N)*INCY) if INCY.LT.0.
 
Keywords:   Level 1 BLAS; DAXPY
 
GAMS:       D1a7
 
Chapters:   MATH/LIBRARY Basic Matrix/Vector Operations
            STAT/LIBRARY Mathematical Support
 
Page No.:   MATH/LIBRARY User's Manual page 1143
            STAT/LIBRARY User's Manual page 1505