zgthrz - Gather and zero.
SUBROUTINE ZGTHRZ(NZ, Y, X, INDX)
DOUBLE COMPLEX Y(*), X(*) INTEGER NZ INTEGER INDX(*)
SUBROUTINE ZGTHRZ_64(NZ, Y, X, INDX)
DOUBLE COMPLEX Y(*), X(*) INTEGER*8 NZ INTEGER*8 INDX(*)
F95 INTERFACE SUBROUTINE GTHRZ([NZ], Y, X, INDX)
COMPLEX(8), DIMENSION(:) :: Y, X INTEGER :: NZ INTEGER, DIMENSION(:) :: INDX
SUBROUTINE GTHRZ_64([NZ], Y, X, INDX)
COMPLEX(8), DIMENSION(:) :: Y, X INTEGER(8) :: NZ INTEGER(8), DIMENSION(:) :: INDX
ZGTHRZ - Gathers the specified elements from a vector y in full storage form into a vector x in compressed form. The gathered elements of y are set to zero. Only the elements of y whose indices are listed in indx are referenced.
do i = 1, n x(i) = y(indx(i)) y(indx(i)) = 0 enddo