Contents


NAME

     saxpy - compute y := alpha * x + y

SYNOPSIS

     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);

PURPOSE

     saxpy compute y := alpha * x + y where alpha is a scalar and
     x and y are n-vectors.

ARGUMENTS

     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.  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.  Unchanged on exit.