NAG Library Routine Document

F06GTF (ZAXPYI)

+− Contents

1 Purpose

9 Example

F06GTF (ZAXPYI) adds a scaled sparse complex vector to an unscaled complex vector.

SUBROUTINE F06GTF (

INTEGER	NZ, INDX(*)
COMPLEX (KIND=nag_wp)	A, X(), Y()

The routine may be called by its BLAS name zaxpyi.

F06GTF (ZAXPYI) performs the operation

y \leftarrow α x + y

where

x

is a sparse complex vector stored in compressed form, and

y

is a complex vector in full storage form.

Dodson D S, Grimes R G and Lewis J G (1991) Sparse extensions to the Fortran basic linear algebra subprograms ACM Trans. Math. Software 17 253–263

1: NZ – INTEGERInput: On entry: the number of nonzeros in the sparse vector $x$ .
2: A – COMPLEX (KIND=nag_wp)Input: On entry: the scalar $α$ .
3: X( $*$ ) – COMPLEX (KIND=nag_wp) arrayInput: Note: the dimension of the array X must be at least $\max (1, NZ)$ .
On entry: the compressed vector $x$ . X contains $x_{i}$ for $i \in J$ .
4: INDX( $*$ ) – INTEGER arrayInput: Note: the dimension of the array INDX must be at least $\max (1, NZ)$ .
On entry: the indices of the elements in the compressed vector $x$ .
Constraint: the indices must be distinct.
5: Y( $*$ ) – COMPLEX (KIND=nag_wp) arrayInput/Output: Note: the dimension of the array Y must be at least $\max_{k} \{INDX (k)\}$ .
On entry: the vector $y$ . Only elements corresponding to indices in INDX are accessed.

On exit: the updated vector $y$ .