1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
# 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
|
