Contents
     saxpy - compute y := alpha * x + y
     SUBROUTINE SAXPY(N, ALPHA, X, INCX, Y, INCY)
     INTEGER N, INCX, INCY
     REAL ALPHA
     REAL X(*), Y(*)
     SUBROUTINE SAXPY_64(N, ALPHA, X, INCX, Y, INCY)
     INTEGER*8 N, INCX, INCY
     REAL ALPHA
     REAL X(*), Y(*)
  F95 INTERFACE
     SUBROUTINE AXPY([N], ALPHA, X, [INCX], Y, [INCY])
     INTEGER :: N, INCX, INCY
     REAL :: ALPHA
     REAL, DIMENSION(:) :: X, Y
     SUBROUTINE AXPY_64([N], ALPHA, X, [INCX], Y, [INCY])
     INTEGER(8) :: N, INCX, INCY
     REAL :: ALPHA
     REAL, DIMENSION(:) :: X, Y
  C INTERFACE
     #include <sunperf.h>
     void saxpy(int n, float alpha, float *x, int incx, float *y,
               int incy);
     void saxpy_64(long n, float  alpha,  float  *x,  long  incx,
               float *y, long incy);
     saxpy compute y := alpha * x + y where alpha is a scalar and
     x and y are n-vectors.
     N (input)
               On entry, N specifies the number  of  elements  in
               the  vector.   N must be at least one for the sub-
               routine to have any visible effect.  Unchanged  on
               exit.
     ALPHA (input)
               On  entry,  ALPHA  specifies  the  scalar   alpha.
               Unchanged on exit.
     X (input)
               ( 1 + ( n - 1 )*abs( INCX ) ).  Before entry,  the
               incremented  array  X  must  contain the vector x.
               Unchanged on exit.
     INCX (input)
               On entry, INCX specifies  the  increment  for  the
               elements  of  X. INCX must not be zero.  Unchanged
               on exit.
     Y (input/output)
               ( 1 + ( n - 1 )*abs( INCY  )  ).   On  entry,  the
               incremented  array Y must contain the vector y. On
               exit, Y is overwritten by the updated vector y.
     INCY (input)
               On entry, INCY specifies  the  increment  for  the
               elements  of  Y. INCY must not be zero.  Unchanged
               on exit.