daxpyi - Compute y := alpha * x + y
SUBROUTINE DAXPYI(NZ, A, X, INDX, Y) DOUBLE PRECISION A DOUBLE PRECISION X(*), Y(*) INTEGER NZ INTEGER INDX(*) SUBROUTINE DAXPYI_64(NZ, A, X, INDX, Y) DOUBLE PRECISION A DOUBLE PRECISION X(*), Y(*) INTEGER*8 NZ INTEGER*8 INDX(*) F95 INTERFACE SUBROUTINE AXPYI(NZ, A, X, INDX, Y) REAL(8) :: A REAL(8), DIMENSION(:) :: X, Y INTEGER :: NZ INTEGER, DIMENSION(:) :: INDX SUBROUTINE AXPYI_64(NZ, A, X, INDX, Y) REAL(8) :: A REAL(8), DIMENSION(:) :: X, Y INTEGER(8) :: NZ INTEGER(8), DIMENSION(:) :: INDX C INTERFACE #include <sunperf.h> void daxpyi (const int nz, const double a, const double* x, const int* indx, double* y); void daxpyi_64 (const long nz, const double a, const double* x, const long* indx, double* y);
Oracle Solaris Studio Performance Library daxpyi(3P) NAME daxpyi - Compute y := alpha * x + y SYNOPSIS SUBROUTINE DAXPYI(NZ, A, X, INDX, Y) DOUBLE PRECISION A DOUBLE PRECISION X(*), Y(*) INTEGER NZ INTEGER INDX(*) SUBROUTINE DAXPYI_64(NZ, A, X, INDX, Y) DOUBLE PRECISION A DOUBLE PRECISION X(*), Y(*) INTEGER*8 NZ INTEGER*8 INDX(*) F95 INTERFACE SUBROUTINE AXPYI(NZ, A, X, INDX, Y) REAL(8) :: A REAL(8), DIMENSION(:) :: X, Y INTEGER :: NZ INTEGER, DIMENSION(:) :: INDX SUBROUTINE AXPYI_64(NZ, A, X, INDX, Y) REAL(8) :: A REAL(8), DIMENSION(:) :: X, Y INTEGER(8) :: NZ INTEGER(8), DIMENSION(:) :: INDX C INTERFACE #include <sunperf.h> void daxpyi (const int nz, const double a, const double* x, const int* indx, double* y); void daxpyi_64 (const long nz, const double a, const double* x, const long* indx, double* y); PURPOSE DAXPYI Compute y := alpha * x + y where alpha is a scalar, x is a sparse vector, and y is a vector in full storage form do i = 1, n y(indx(i)) = alpha * x(i) + y(indx(i)) enddo ARGUMENTS NZ (input) - INTEGER Number of elements in the compressed form. Unchanged on exit. A (input) On entry, A(LPHA) specifies the scaling value. Unchanged on exit. X (input) Vector containing the values of the compressed form. Unchanged on exit. INDX (input) - INTEGER Vector containing the indices of the compressed form. It is assumed that the elements in INDX are distinct and greater than zero. Unchanged on exit. Y (output) Vector on input which contains the vector Y in full storage form. On exit, only the elements corresponding to the indices in INDX have been modified. 3rd Berkeley Distribution 7 Nov 2015 daxpyi(3P)