# Copyright(c) 1986 Association of Universities for Research in Astronomy Inc.

# ALUT -- Map an array using table lookup.  Note that an input value of zero
# indexes the first element of the lookup table.  No bounds checking is
# performed.

procedure alut$t (a, b, npix, lut)

$if (datatype == rd)
int	a[ARB]				# input array of indices
$else
PIXEL	a[ARB]
$endif

PIXEL	b[ARB]				# output data array
PIXEL	lut[ARB]			# lookup table
int	npix, i

begin
	do i = 1, npix
	    b[i] = lut[a[i]+1]
end