HP MLIB User's Guide Vol. 1 7th Ed.

Chapter 2 Basic Vector Operations 97
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)