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

# MW_MMUL -- Matrix multiply.

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

PIXEL	a[ndim,ndim]		#I left input matrix
PIXEL	b[ndim,ndim]		#I right input matrix
PIXEL	c[ndim,ndim]		#O output matrix
int	ndim			#I dimensionality of system

int	i, j, k
PIXEL	v

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