void multtensor(float (*c)[4][4], float (*a)[4], float (*b)[4][4]) {
  int i, j, k;
  float x1, x2, x3, x4;
  for (i = 0; i < 4; i++)
    for (j = 0; j < 4; j++) {
      x1 = x2 = x3 = x4 = 0.0;
      for (k = 0; k < 4; k++) {
        x1 += b[0][i][k] * a[k][j];
        x2 += b[1][i][k] * a[k][j];
        x3 += b[2][i][k] * a[k][j];
        x4 += b[3][i][k] * a[k][j];
      }
      c[0][i][j] = x1;
      c[1][i][j] = x2;
      c[2][i][j] = x3;
      c[3][i][j] = x4;
    }
}


x86 assembly (clang, -O0):
<multtensor>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
cmpl   $0x4,-0x1c(%rbp)
jge    132e <multtensor+0x22e>
movl   $0x0,-0x20(%rbp)
cmpl   $0x4,-0x20(%rbp)
jge    131b <multtensor+0x21b>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x34(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x30(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x2c(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x28(%rbp)
movl   $0x0,-0x24(%rbp)
cmpl   $0x4,-0x24(%rbp)
jge    1289 <multtensor+0x189>
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
shl    $0x4,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
shl    $0x4,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x28(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x28(%rbp)
mov    -0x18(%rbp),%rax
add    $0x40,%rax
movslq -0x1c(%rbp),%rcx
shl    $0x4,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
shl    $0x4,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x2c(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x2c(%rbp)
mov    -0x18(%rbp),%rax
add    $0x80,%rax
movslq -0x1c(%rbp),%rcx
shl    $0x4,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
shl    $0x4,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x30(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x30(%rbp)
mov    -0x18(%rbp),%rax
add    $0xc0,%rax
movslq -0x1c(%rbp),%rcx
shl    $0x4,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
shl    $0x4,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x34(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x34(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1159 <multtensor+0x59>
movss  -0x28(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
shl    $0x4,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x2c(%rbp),%xmm0
mov    -0x8(%rbp),%rax
add    $0x40,%rax
movslq -0x1c(%rbp),%rcx
shl    $0x4,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x30(%rbp),%xmm0
mov    -0x8(%rbp),%rax
add    $0x80,%rax
movslq -0x1c(%rbp),%rcx
shl    $0x4,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x34(%rbp),%xmm0
mov    -0x8(%rbp),%rax
add    $0xc0,%rax
movslq -0x1c(%rbp),%rcx
shl    $0x4,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1128 <multtensor+0x28>
jmp    1320 <multtensor+0x220>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1117 <multtensor+0x17>
pop    %rbp
ret

x86 assembly (clang, -O2):
<multtensor>:
add    $0xc0,%rdi
xor    %r8d,%r8d
xorps  %xmm0,%xmm0
nopl   (%rax)
mov    %r8,%rax
shl    $0x4,%rax
lea    (%rdx,%rax,1),%r9
add    %rdx,%rax
add    $0x40,%rax
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movss  (%rsi,%rcx,4),%xmm2
movss  0x10(%rsi,%rcx,4),%xmm1
movss  (%r9),%xmm3
mulss  %xmm2,%xmm3
addss  %xmm0,%xmm3
movss  (%rax),%xmm4
mulss  %xmm2,%xmm4
addss  %xmm0,%xmm4
movss  0x40(%rax),%xmm5
mulss  %xmm2,%xmm5
addss  %xmm0,%xmm5
mulss  0x80(%rax),%xmm2
addss  %xmm0,%xmm2
movss  -0x3c(%rax),%xmm6
mulss  %xmm1,%xmm6
addss  %xmm3,%xmm6
movss  0x4(%rax),%xmm3
mulss  %xmm1,%xmm3
addss  %xmm4,%xmm3
movss  0x44(%rax),%xmm4
mulss  %xmm1,%xmm4
addss  %xmm5,%xmm4
mulss  0x84(%rax),%xmm1
addss  %xmm2,%xmm1
movss  0x20(%rsi,%rcx,4),%xmm2
movss  -0x38(%rax),%xmm5
mulss  %xmm2,%xmm5
addss  %xmm6,%xmm5
movss  0x8(%rax),%xmm6
mulss  %xmm2,%xmm6
addss  %xmm3,%xmm6
movss  0x48(%rax),%xmm3
mulss  %xmm2,%xmm3
addss  %xmm4,%xmm3
mulss  0x88(%rax),%xmm2
addss  %xmm1,%xmm2
movss  0x30(%rsi,%rcx,4),%xmm1
movss  -0x34(%rax),%xmm4
mulss  %xmm1,%xmm4
addss  %xmm5,%xmm4
movss  0xc(%rax),%xmm5
mulss  %xmm1,%xmm5
addss  %xmm6,%xmm5
movss  0x4c(%rax),%xmm6
mulss  %xmm1,%xmm6
mulss  0x8c(%rax),%xmm1
addss  %xmm3,%xmm6
addss  %xmm2,%xmm1
movss  %xmm4,-0xc0(%rdi,%rcx,4)
ff ff
movss  %xmm5,-0x80(%rdi,%rcx,4)
movss  %xmm6,-0x40(%rdi,%rcx,4)
movss  %xmm1,(%rdi,%rcx,4)
add    $0x1,%rcx
cmp    $0x4,%rcx
jne    1130 <multtensor+0x30>
add    $0x1,%r8
add    $0x10,%rdi
cmp    $0x4,%r8
jne    1110 <multtensor+0x10>
ret

LLVM IR (-O0):
define void @multtensor([4 x [4 x float]]* %0, [4 x float]* %1, [4 x [4 x float]]* %2) {
%4 = alloca [4 x [4 x float]]*
%5 = alloca [4 x float]*
%6 = alloca [4 x [4 x float]]*
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca float
%11 = alloca float
%12 = alloca float
%13 = alloca float
store [4 x [4 x float]]* %0, [4 x [4 x float]]** %4
store [4 x float]* %1, [4 x float]** %5
store [4 x [4 x float]]* %2, [4 x [4 x float]]** %6
store i32 0, i32* %7
br label %14
14: ; preds = %146, %3
%15 = load i32, i32* %7
%16 = icmp slt i32 %15, 4
br i1 %16, label %17, label %149
17: ; preds = %14
store i32 0, i32* %8
br label %18
18: ; preds = %142, %17
%19 = load i32, i32* %8
%20 = icmp slt i32 %19, 4
br i1 %20, label %21, label %145
21: ; preds = %18
store float 0.000000e+00, float* %13
store float 0.000000e+00, float* %12
store float 0.000000e+00, float* %11
store float 0.000000e+00, float* %10
store i32 0, i32* %9
br label %22
22: ; preds = %102, %21
%23 = load i32, i32* %9
%24 = icmp slt i32 %23, 4
br i1 %24, label %25, label %105
25: ; preds = %22
%26 = load [4 x [4 x float]]*, [4 x [4 x float]]** %6
%27 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %26, i64 0
%28 = load i32, i32* %7
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %27, i64 0, i64 %29
%31 = load i32, i32* %9
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds [4 x float], [4 x float]* %30, i64 0, i64 %32
%34 = load float, float* %33
%35 = load [4 x float]*, [4 x float]** %5
%36 = load i32, i32* %9
%37 = sext i32 %36 to i64
%38 = getelementptr inbounds [4 x float], [4 x float]* %35, i64 %37
%39 = load i32, i32* %8
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds [4 x float], [4 x float]* %38, i64 0, i64 %40
%42 = load float, float* %41
%43 = load float, float* %10
%44 = call float @llvm.fmuladd.f32(float %34, float %42, float %43)
store float %44, float* %10
%45 = load [4 x [4 x float]]*, [4 x [4 x float]]** %6
%46 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %45, i64 1
%47 = load i32, i32* %7
%48 = sext i32 %47 to i64
%49 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %46, i64 0, i64 %48
%50 = load i32, i32* %9
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds [4 x float], [4 x float]* %49, i64 0, i64 %51
%53 = load float, float* %52
%54 = load [4 x float]*, [4 x float]** %5
%55 = load i32, i32* %9
%56 = sext i32 %55 to i64
%57 = getelementptr inbounds [4 x float], [4 x float]* %54, i64 %56
%58 = load i32, i32* %8
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds [4 x float], [4 x float]* %57, i64 0, i64 %59
%61 = load float, float* %60
%62 = load float, float* %11
%63 = call float @llvm.fmuladd.f32(float %53, float %61, float %62)
store float %63, float* %11
%64 = load [4 x [4 x float]]*, [4 x [4 x float]]** %6
%65 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %64, i64 2
%66 = load i32, i32* %7
%67 = sext i32 %66 to i64
%68 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %65, i64 0, i64 %67
%69 = load i32, i32* %9
%70 = sext i32 %69 to i64
%71 = getelementptr inbounds [4 x float], [4 x float]* %68, i64 0, i64 %70
%72 = load float, float* %71
%73 = load [4 x float]*, [4 x float]** %5
%74 = load i32, i32* %9
%75 = sext i32 %74 to i64
%76 = getelementptr inbounds [4 x float], [4 x float]* %73, i64 %75
%77 = load i32, i32* %8
%78 = sext i32 %77 to i64
%79 = getelementptr inbounds [4 x float], [4 x float]* %76, i64 0, i64 %78
%80 = load float, float* %79
%81 = load float, float* %12
%82 = call float @llvm.fmuladd.f32(float %72, float %80, float %81)
store float %82, float* %12
%83 = load [4 x [4 x float]]*, [4 x [4 x float]]** %6
%84 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %83, i64 3
%85 = load i32, i32* %7
%86 = sext i32 %85 to i64
%87 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %84, i64 0, i64 %86
%88 = load i32, i32* %9
%89 = sext i32 %88 to i64
%90 = getelementptr inbounds [4 x float], [4 x float]* %87, i64 0, i64 %89
%91 = load float, float* %90
%92 = load [4 x float]*, [4 x float]** %5
%93 = load i32, i32* %9
%94 = sext i32 %93 to i64
%95 = getelementptr inbounds [4 x float], [4 x float]* %92, i64 %94
%96 = load i32, i32* %8
%97 = sext i32 %96 to i64
%98 = getelementptr inbounds [4 x float], [4 x float]* %95, i64 0, i64 %97
%99 = load float, float* %98
%100 = load float, float* %13
%101 = call float @llvm.fmuladd.f32(float %91, float %99, float %100)
store float %101, float* %13
br label %102
102: ; preds = %25
%103 = load i32, i32* %9
%104 = add nsw i32 %103, 1
store i32 %104, i32* %9
br label %22
105: ; preds = %22
%106 = load float, float* %10
%107 = load [4 x [4 x float]]*, [4 x [4 x float]]** %4
%108 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %107, i64 0
%109 = load i32, i32* %7
%110 = sext i32 %109 to i64
%111 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %108, i64 0, i64 %110
%112 = load i32, i32* %8
%113 = sext i32 %112 to i64
%114 = getelementptr inbounds [4 x float], [4 x float]* %111, i64 0, i64 %113
store float %106, float* %114
%115 = load float, float* %11
%116 = load [4 x [4 x float]]*, [4 x [4 x float]]** %4
%117 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %116, i64 1
%118 = load i32, i32* %7
%119 = sext i32 %118 to i64
%120 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %117, i64 0, i64 %119
%121 = load i32, i32* %8
%122 = sext i32 %121 to i64
%123 = getelementptr inbounds [4 x float], [4 x float]* %120, i64 0, i64 %122
store float %115, float* %123
%124 = load float, float* %12
%125 = load [4 x [4 x float]]*, [4 x [4 x float]]** %4
%126 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %125, i64 2
%127 = load i32, i32* %7
%128 = sext i32 %127 to i64
%129 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %126, i64 0, i64 %128
%130 = load i32, i32* %8
%131 = sext i32 %130 to i64
%132 = getelementptr inbounds [4 x float], [4 x float]* %129, i64 0, i64 %131
store float %124, float* %132
%133 = load float, float* %13
%134 = load [4 x [4 x float]]*, [4 x [4 x float]]** %4
%135 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %134, i64 3
%136 = load i32, i32* %7
%137 = sext i32 %136 to i64
%138 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %135, i64 0, i64 %137
%139 = load i32, i32* %8
%140 = sext i32 %139 to i64
%141 = getelementptr inbounds [4 x float], [4 x float]* %138, i64 0, i64 %140
store float %133, float* %141
br label %142
142: ; preds = %105
%143 = load i32, i32* %8
%144 = add nsw i32 %143, 1
store i32 %144, i32* %8
br label %18
145: ; preds = %18
br label %146
146: ; preds = %145
%147 = load i32, i32* %7
%148 = add nsw i32 %147, 1
store i32 %148, i32* %7
br label %14
149: ; preds = %14
ret void
}

LLVM IR (-O2):
define void @multtensor([4 x [4 x float]]* nocapture writeonly %0, [4 x float]* nocapture %1, [4 x [4 x float]]* nocapture %2) {
br label %4
4: ; preds = %3, %70
%5 = phi i64 [ 0, %3 ], [ %71, %70 ]
%6 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %2, i64 0, i64 %5, i64 0
%7 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %2, i64 1, i64 %5, i64 0
%8 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %2, i64 2, i64 %5, i64 0
%9 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %2, i64 3, i64 %5, i64 0
%10 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %2, i64 0, i64 %5, i64 1
%11 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %2, i64 1, i64 %5, i64 1
%12 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %2, i64 2, i64 %5, i64 1
%13 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %2, i64 3, i64 %5, i64 1
%14 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %2, i64 0, i64 %5, i64 2
%15 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %2, i64 1, i64 %5, i64 2
%16 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %2, i64 2, i64 %5, i64 2
%17 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %2, i64 3, i64 %5, i64 2
%18 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %2, i64 0, i64 %5, i64 3
%19 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %2, i64 1, i64 %5, i64 3
%20 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %2, i64 2, i64 %5, i64 3
%21 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %2, i64 3, i64 %5, i64 3
br label %22
22: ; preds = %4, %22
%23 = phi i64 [ 0, %4 ], [ %68, %22 ]
%24 = load float, float* %6
%25 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 0, i64 %23
%26 = load float, float* %25
%27 = tail call float @llvm.fmuladd.f32(float %24, float %26, float 0.000000e+00)
%28 = load float, float* %7
%29 = tail call float @llvm.fmuladd.f32(float %28, float %26, float 0.000000e+00)
%30 = load float, float* %8
%31 = tail call float @llvm.fmuladd.f32(float %30, float %26, float 0.000000e+00)
%32 = load float, float* %9
%33 = tail call float @llvm.fmuladd.f32(float %32, float %26, float 0.000000e+00)
%34 = load float, float* %10
%35 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 1, i64 %23
%36 = load float, float* %35
%37 = tail call float @llvm.fmuladd.f32(float %34, float %36, float %27)
%38 = load float, float* %11
%39 = tail call float @llvm.fmuladd.f32(float %38, float %36, float %29)
%40 = load float, float* %12
%41 = tail call float @llvm.fmuladd.f32(float %40, float %36, float %31)
%42 = load float, float* %13
%43 = tail call float @llvm.fmuladd.f32(float %42, float %36, float %33)
%44 = load float, float* %14
%45 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 2, i64 %23
%46 = load float, float* %45
%47 = tail call float @llvm.fmuladd.f32(float %44, float %46, float %37)
%48 = load float, float* %15
%49 = tail call float @llvm.fmuladd.f32(float %48, float %46, float %39)
%50 = load float, float* %16
%51 = tail call float @llvm.fmuladd.f32(float %50, float %46, float %41)
%52 = load float, float* %17
%53 = tail call float @llvm.fmuladd.f32(float %52, float %46, float %43)
%54 = load float, float* %18
%55 = getelementptr inbounds [4 x float], [4 x float]* %1, i64 3, i64 %23
%56 = load float, float* %55
%57 = tail call float @llvm.fmuladd.f32(float %54, float %56, float %47)
%58 = load float, float* %19
%59 = tail call float @llvm.fmuladd.f32(float %58, float %56, float %49)
%60 = load float, float* %20
%61 = tail call float @llvm.fmuladd.f32(float %60, float %56, float %51)
%62 = load float, float* %21
%63 = tail call float @llvm.fmuladd.f32(float %62, float %56, float %53)
%64 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %0, i64 0, i64 %5, i64 %23
store float %57, float* %64
%65 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %0, i64 1, i64 %5, i64 %23
store float %59, float* %65
%66 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %0, i64 2, i64 %5, i64 %23
store float %61, float* %66
%67 = getelementptr inbounds [4 x [4 x float]], [4 x [4 x float]]* %0, i64 3, i64 %5, i64 %23
store float %63, float* %67
%68 = add nuw nsw i64 %23, 1
%69 = icmp eq i64 %68, 4
br i1 %69, label %70, label %22
70: ; preds = %22
%71 = add nuw nsw i64 %5, 1
%72 = icmp eq i64 %71, 4
br i1 %72, label %73, label %4
73: ; preds = %70
ret void
}