zaxpyi(3P) Sun Performance Library zaxpyi(3P)NAMEzaxpyi - Compute y := alpha * x + y
SYNOPSIS
SUBROUTINE ZAXPYI(NZ, A, X, INDX, Y)
DOUBLE COMPLEX A
DOUBLE COMPLEX X(*), Y(*)
INTEGER NZ
INTEGER INDX(*)
SUBROUTINE ZAXPYI_64(NZ, A, X, INDX, Y)
DOUBLE COMPLEX A
DOUBLE COMPLEX X(*), Y(*)
INTEGER*8 NZ
INTEGER*8 INDX(*)
F95 INTERFACE
SUBROUTINE AXPYI([NZ], [A], X, INDX, Y)
COMPLEX(8) :: A
COMPLEX(8), DIMENSION(:) :: X, Y
INTEGER :: NZ
INTEGER, DIMENSION(:) :: INDX
SUBROUTINE AXPYI_64([NZ], [A], X, INDX, Y)
COMPLEX(8) :: A
COMPLEX(8), DIMENSION(:) :: X, Y
INTEGER(8) :: NZ
INTEGER(8), DIMENSION(:) :: INDX
C INTERFACE
#include <sunperf.h>
void zaxpyi (const int nz, const doublecomplex* a, const doublecomplex*
x, const int* indx, doublecomplex* y);
void zaxpyi_64 (const long nz, const doublecomplex* a, const doublecom‐
plex* x, const long* indx, doublecomplex* y);
PURPOSE
ZAXPYI 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. A is defaulted to (1.0D0,0.0D0) for F95 INTERFACE.
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 6 Mar 2009 zaxpyi(3P)
