1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
|
# MVMUL -- Multply a matrix (left-hand side) by a one dimensional vector
# (right-hand side) and return the resultant vector.
procedure mvmul (matrix, maxdim, dim, vector, result)
real matrix [maxdim, maxdim] # input matrix
int maxdim # maximum size of input matrix
int dim # dimension of matrix and vectors
real vector[maxdim] # input vector
real result[maxdim] # iutput vector
double sum
int i, j
begin
do i = 1, dim {
sum = 0.0
do j = 1, dim
sum = sum + double (matrix[j,i]) * double(vector[j])
result[i] = sum
}
end
# DMVMUL -- Multply a matrix (left-hand side) by a one dimensional vector
# (right-hand side) and return the resultant vector.
procedure dmvmul (matrix, maxdim, dim, vector, result)
double matrix [maxdim, maxdim] # input matrix
int maxdim # maximum size of input matrix
int dim # dimension of matrix and vectors
double vector[maxdim] # input vector
double result[maxdim] # iutput vector
double sum
int i, j
begin
do i = 1, dim {
sum = 0.0d0
do j = 1, dim
sum = sum + (matrix[j,i] * vector[j])
result[i] = sum
}
end
|
