blob: 7770c603a6e3890ba10617543f78bad1742198ea (
plain)
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
|
#define __rose_lt(x,y) ((x)<(y)?(x):(y))
#define __rose_gt(x,y) ((x)>(y)?(x):(y))
void mm(float **A,float **B,float **C,int ambn,int an,int bm)
{
int t6;
int t4;
int t2;
int chill_t9;
int chill_t8;
int chill_t7;
int chill_t6;
int chill_t5;
int chill_t4;
int over1;
int chill_t3;
int chill_t2;
int chill_t1;
int i;
int j;
int n;
over1 = 0;
for (t2 = 0; t2 <= an - 1; t2 += 1)
for (t4 = 0; t4 <= bm - 1; t4 += 1)
C[t2][t4] = 0.0f;
for (t2 = 0; t2 <= an - 1; t2 += 1)
for (t4 = 0; t4 <= bm - 1; t4 += 1) {
over1 = ambn % 4;
for (t6 = 0; t6 <= -over1 + ambn - 1; t6 += 4) {
C[t2][t4] += (A[t2][t6] * B[t6][t4]);
C[t2][t4] += (A[t2][t6 + 1] * B[t6 + 1][t4]);
C[t2][t4] += (A[t2][t6 + 2] * B[t6 + 2][t4]);
C[t2][t4] += (A[t2][t6 + 3] * B[t6 + 3][t4]);
}
for (t6 = __rose_gt(0,ambn - over1); t6 <= ambn - 1; t6 += 1)
C[t2][t4] += (A[t2][t6] * B[t6][t4]);
}
}
|