Gather and zero sparse vector

SGTHRZ/DGTHRZ/IGTHRZ/CGTHRZ/ZGTHRZ

INTEGER*8

m, indx(m)

REAL*4

y(n), x(m)

CALL SGTHRZ(m, y, x, indx)

INTEGER*8

m, indx(m)

REAL*8

y(n), x(m)

CALL DGTHRZ(m, y, x, indx)

INTEGER*8

m, indx(m), y(n), x(m)

CALL IGTHRZ(m, y, x, indx)

INTEGER*8

m, indx(m)

COMPLEX*8

y(n), x(m)

CALL CGTHRZ(m, y, x, indx)

INTEGER*8

m, indx(m)

COMPLEX*16

y(n), x(m)

CALL ZGTHRZ(m, y, x, indx)

Chapter 2 Basic Vector Operations 97