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

# MW_VMUL -- Vector multiply.

procedure mw_vmul$t (a, b, c, ndim)

PIXEL	a[ndim,ndim]		#I input matrix
PIXEL	b[ndim]			#I input vector
PIXEL	c[ndim]			#O output vector
int	ndim			#I system dimension

int	i, j
PIXEL	v

begin
	do j = 1, ndim {
	    v = 0
	    do i = 1, ndim
		v = v + a[i,j] * b[i]
	    c[j] = v
	}
end