1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
|
# MW_MMUL -- Matrix multiply.
procedure mw_mmuld (a, b, c, ndim)
double a[ndim,ndim] #I left input matrix
double b[ndim,ndim] #I right input matrix
double c[ndim,ndim] #O output matrix
int ndim #I dimensionality of system
int i, j, k
double 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
|
