void fwd(float diag[32][32], float col[32][32]) {
  int i, j, k;
  for (j = 0; j < 32; j++)
    for (k = 0; k < 32; k++)
      for (i = k + 1; i < 32; i++)
        col[i][j] -= diag[i][k] * col[k][j];
}


x86 assembly (clang, -O0):
<fwd>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movl   $0x0,-0x18(%rbp)
cmpl   $0x20,-0x18(%rbp)
jge    11d8 <fwd+0xd8>
movl   $0x0,-0x1c(%rbp)
cmpl   $0x20,-0x1c(%rbp)
jge    11c5 <fwd+0xc5>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
cmpl   $0x20,-0x14(%rbp)
jge    11b2 <fwd+0xb2>
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
shl    $0x7,%rcx
add    %rcx,%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
shl    $0x7,%rcx
add    %rcx,%rax
movslq -0x18(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
shl    $0x7,%rcx
add    %rcx,%rax
movslq -0x18(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
movd   %xmm0,%edx
xor    $0x80000000,%edx
movd   %edx,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1137 <fwd+0x37>
jmp    11b7 <fwd+0xb7>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1124 <fwd+0x24>
jmp    11ca <fwd+0xca>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1113 <fwd+0x13>
pop    %rbp
ret

x86 assembly (clang, -O2):
<fwd>:
push   %r14
push   %rbx
xor    %r9d,%r9d
mov    %rsi,%r8
jmp    111e <fwd+0x1e>
nopl   0x0(%rax,%rax,1)
add    $0x1,%r9
add    $0x4,%r8
cmp    $0x20,%r9
je     1194 <fwd+0x94>
mov    $0x1f,%r11d
mov    %r8,%rcx
mov    %rdi,%r14
xor    %r10d,%r10d
jmp    1145 <fwd+0x45>
nop
add    $0x84,%r14
sub    $0xffffffffffffff80,%rcx
add    $0xffffffffffffffff,%r11
cmp    $0x20,%r10
je     1110 <fwd+0x10>
mov    %r10,%rax
add    $0x1,%r10
cmp    $0x1e,%rax
ja     1130 <fwd+0x30>
shl    $0x7,%rax
add    %rsi,%rax
lea    (%rax,%r9,4),%rdx
mov    $0x80,%ebx
mov    %r11,%rax
cs nopw 0x0(%rax,%rax,1)
nop
movss  (%r14,%rbx,1),%xmm0
movss  (%rcx,%rbx,1),%xmm1
mulss  (%rdx),%xmm0
subss  %xmm0,%xmm1
movss  %xmm1,(%rcx,%rbx,1)
sub    $0xffffffffffffff80,%rbx
add    $0xffffffffffffffff,%rax
jne    1170 <fwd+0x70>
jmp    1130 <fwd+0x30>
pop    %rbx
pop    %r14
ret

LLVM IR (-O0):
define void @fwd([32 x float]* %0, [32 x float]* %1) {
%3 = alloca [32 x float]*
%4 = alloca [32 x float]*
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32
store [32 x float]* %0, [32 x float]** %3
store [32 x float]* %1, [32 x float]** %4
store i32 0, i32* %6
br label %8
8: ; preds = %56, %2
%9 = load i32, i32* %6
%10 = icmp slt i32 %9, 32
br i1 %10, label %11, label %59
11: ; preds = %8
store i32 0, i32* %7
br label %12
12: ; preds = %52, %11
%13 = load i32, i32* %7
%14 = icmp slt i32 %13, 32
br i1 %14, label %15, label %55
15: ; preds = %12
%16 = load i32, i32* %7
%17 = add nsw i32 %16, 1
store i32 %17, i32* %5
br label %18
18: ; preds = %48, %15
%19 = load i32, i32* %5
%20 = icmp slt i32 %19, 32
br i1 %20, label %21, label %51
21: ; preds = %18
%22 = load [32 x float]*, [32 x float]** %3
%23 = load i32, i32* %5
%24 = sext i32 %23 to i64
%25 = getelementptr inbounds [32 x float], [32 x float]* %22, i64 %24
%26 = load i32, i32* %7
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds [32 x float], [32 x float]* %25, i64 0, i64 %27
%29 = load float, float* %28
%30 = load [32 x float]*, [32 x float]** %4
%31 = load i32, i32* %7
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds [32 x float], [32 x float]* %30, i64 %32
%34 = load i32, i32* %6
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds [32 x float], [32 x float]* %33, i64 0, i64 %35
%37 = load float, float* %36
%38 = load [32 x float]*, [32 x float]** %4
%39 = load i32, i32* %5
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds [32 x float], [32 x float]* %38, i64 %40
%42 = load i32, i32* %6
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds [32 x float], [32 x float]* %41, i64 0, i64 %43
%45 = load float, float* %44
%46 = fneg float %29
%47 = call float @llvm.fmuladd.f32(float %46, float %37, float %45)
store float %47, float* %44
br label %48
48: ; preds = %21
%49 = load i32, i32* %5
%50 = add nsw i32 %49, 1
store i32 %50, i32* %5
br label %18
51: ; preds = %18
br label %52
52: ; preds = %51
%53 = load i32, i32* %7
%54 = add nsw i32 %53, 1
store i32 %54, i32* %7
br label %12
55: ; preds = %12
br label %56
56: ; preds = %55
%57 = load i32, i32* %6
%58 = add nsw i32 %57, 1
store i32 %58, i32* %6
br label %8
59: ; preds = %8
ret void
}

LLVM IR (-O2):
define void @fwd([32 x float]* nocapture %0, [32 x float]* nocapture %1) {
br label %3
3: ; preds = %2, %26
%4 = phi i64 [ 0, %2 ], [ %27, %26 ]
br label %8
5: ; preds = %15, %8
%6 = add nuw nsw i64 %10, 1
%7 = icmp eq i64 %11, 32
br i1 %7, label %26, label %8
8: ; preds = %3, %5
%9 = phi i64 [ 0, %3 ], [ %11, %5 ]
%10 = phi i64 [ 1, %3 ], [ %6, %5 ]
%11 = add nuw nsw i64 %9, 1
%12 = icmp ult i64 %9, 31
br i1 %12, label %13, label %5
13: ; preds = %8
%14 = getelementptr inbounds [32 x float], [32 x float]* %1, i64 %9, i64 %4
br label %15
15: ; preds = %13, %15
%16 = phi i64 [ %10, %13 ], [ %24, %15 ]
%17 = getelementptr inbounds [32 x float], [32 x float]* %0, i64 %16, i64 %9
%18 = load float, float* %17
%19 = load float, float* %14
%20 = getelementptr inbounds [32 x float], [32 x float]* %1, i64 %16, i64 %4
%21 = load float, float* %20
%22 = fneg float %18
%23 = tail call float @llvm.fmuladd.f32(float %22, float %19, float %21)
store float %23, float* %20
%24 = add nuw nsw i64 %16, 1
%25 = icmp eq i64 %24, 32
br i1 %25, label %5, label %15
26: ; preds = %5
%27 = add nuw nsw i64 %4, 1
%28 = icmp eq i64 %27, 32
br i1 %28, label %29, label %3
29: ; preds = %26
ret void
}