NAG Library Routine Document

F06GVF (ZGTHRZ)

+− Contents

1 Purpose

9 Example

F06GVF (ZGTHRZ) gathers specified (usually nonzero) elements of a complex vector

y

in full storage form into a sparse complex vector

x

in compressed form. The specified elements of

y

are set to zero.

SUBROUTINE F06GVF (

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

The routine may be called by its BLAS name zgthrz.

None.

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: 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 elements of $y$ corresponding to indices in INDX are set to zero.
3: X( $*$ ) – COMPLEX (KIND=nag_wp) arrayOutput: Note: the dimension of the array X must be at least $\max (1, NZ)$ .
On exit: the compressed vector $x$ .
4: INDX( $*$ ) – INTEGER arrayInput: Note: the dimension of the array INDX must be at least $\max (1, NZ)$ .
On entry: INDX must contain the set of indices $J$ .