blob: 985ecaa5e5e4d3d937d48845724f559edf8d4978 (
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
|
// this source is derived from CHILL AST originally from file 'mm.c' as parsed by frontend compiler clang
void mm(float **A, float **B, float **C, int ambn, int an, int bm) {
int t6;
int t4;
int t2;
int over1;
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 = 0 < ambn - over1 ? ambn - over1 : 0; t6 <= ambn - 1; t6 += 1)
C[t2][t4] += A[t2][t6] * B[t6][t4];
}
}
|