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

include	"../mwcs.h"

# MW_LTRAN -- Perform a general N-dimensional linear transformation, i.e.,
# matrix multiply and translation.

procedure mw_ltran$t (p1, p2, ltm, ltv, ndim)

PIXEL	p1[ndim]			#I input point
PIXEL	p2[ndim]			#O transformed output point
PIXEL	ltm[ndim,ndim]			#I linear transformation matrix
PIXEL	ltv[ndim]			#I linear translation vector
int	ndim				#I dimension of system

int	i, j
PIXEL	p3[MAX_DIM]

begin
	call amov$t (p1, p3, ndim)
	do j = 1, ndim {
	    p2[j] = ltv[j]
	    do i = 1, ndim
		p2[j] = p2[j] + ltm[i,j] * p3[i]
	}
end