void split_block(float A[32][32], float L[32][32], float U[32][32]) {
  int i, j, k;
  for (i = 0; i < 32; i++)
    for (j = 0; j < 32; j++) {
      if (i == j) {
        L[i][j] = 1.0;
        U[i][j] = A[i][j];
      } else if (i > j) {
        L[i][j] = A[i][j];
        U[i][j] = 0.0;
      } else {
        L[i][j] = 0.0;
        U[i][j] = A[i][j];
      }
    }
}


x86 assembly (clang, -O0):
<split_block>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
cmpl   $0x20,-0x1c(%rbp)
jge    1265 <split_block+0x165>
movl   $0x0,-0x20(%rbp)
cmpl   $0x20,-0x20(%rbp)
jge    1252 <split_block+0x152>
mov    -0x1c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jne    1193 <split_block+0x93>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
shl    $0x7,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
movss  0xea7(%rip),%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
shl    $0x7,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
shl    $0x7,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    123f <split_block+0x13f>
mov    -0x1c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jle    11ef <split_block+0xef>
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
shl    $0x7,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
shl    $0x7,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
shl    $0x7,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    123a <split_block+0x13a>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
shl    $0x7,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
shl    $0x7,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
shl    $0x7,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    123f <split_block+0x13f>
jmp    1244 <split_block+0x144>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1128 <split_block+0x28>
jmp    1257 <split_block+0x157>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1117 <split_block+0x17>
pop    %rbp
ret

x86 assembly (clang, -O2):
<split_block>:
push   %r14
push   %rbx
xor    %eax,%eax
mov    %rsi,%r8
mov    %rdi,%r9
mov    %rdx,%r10
jmp    112a <split_block+0x2a>
add    $0x1,%rax
sub    $0xffffffffffffff80,%r10
sub    $0xffffffffffffff80,%r9
sub    $0xffffffffffffff80,%r8
cmp    $0x20,%rax
je     11a6 <split_block+0xa6>
mov    %rax,%rcx
shl    $0x7,%rcx
lea    (%rsi,%rcx,1),%rbx
lea    (%rbx,%rax,4),%r11
lea    (%rdi,%rcx,1),%rbx
lea    (%rbx,%rax,4),%r14
add    %rdx,%rcx
lea    (%rcx,%rax,4),%rbx
xor    %ecx,%ecx
jmp    116a <split_block+0x6a>
nopl   0x0(%rax)
movl   $0x3f800000,(%r11)
movss  (%r14),%xmm0
movss  %xmm0,(%rbx)
add    $0x1,%rcx
cmp    $0x20,%rcx
je     1110 <split_block+0x10>
cmp    %rcx,%rax
je     1150 <split_block+0x50>
jbe    1190 <split_block+0x90>
movss  (%r9,%rcx,4),%xmm0
movss  %xmm0,(%r8,%rcx,4)
movl   $0x0,(%r10,%rcx,4)
jmp    1160 <split_block+0x60>
nopw   0x0(%rax,%rax,1)
movl   $0x0,(%r8,%rcx,4)
movss  (%r9,%rcx,4),%xmm0
movss  %xmm0,(%r10,%rcx,4)
jmp    1160 <split_block+0x60>
pop    %rbx
pop    %r14
ret

LLVM IR (-O0):
define void @split_block([32 x float]* %0, [32 x float]* %1, [32 x float]* %2) {
%4 = alloca [32 x float]*
%5 = alloca [32 x float]*
%6 = alloca [32 x float]*
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
store [32 x float]* %0, [32 x float]** %4
store [32 x float]* %1, [32 x float]** %5
store [32 x float]* %2, [32 x float]** %6
store i32 0, i32* %7
br label %10
10: ; preds = %100, %3
%11 = load i32, i32* %7
%12 = icmp slt i32 %11, 32
br i1 %12, label %13, label %103
13: ; preds = %10
store i32 0, i32* %8
br label %14
14: ; preds = %96, %13
%15 = load i32, i32* %8
%16 = icmp slt i32 %15, 32
br i1 %16, label %17, label %99
17: ; preds = %14
%18 = load i32, i32* %7
%19 = load i32, i32* %8
%20 = icmp eq i32 %18, %19
br i1 %20, label %21, label %44
21: ; preds = %17
%22 = load [32 x float]*, [32 x float]** %5
%23 = load i32, i32* %7
%24 = sext i32 %23 to i64
%25 = getelementptr inbounds [32 x float], [32 x float]* %22, i64 %24
%26 = load i32, i32* %8
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds [32 x float], [32 x float]* %25, i64 0, i64 %27
store float 1.000000e+00, float* %28
%29 = load [32 x float]*, [32 x float]** %4
%30 = load i32, i32* %7
%31 = sext i32 %30 to i64
%32 = getelementptr inbounds [32 x float], [32 x float]* %29, i64 %31
%33 = load i32, i32* %8
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds [32 x float], [32 x float]* %32, i64 0, i64 %34
%36 = load float, float* %35
%37 = load [32 x float]*, [32 x float]** %6
%38 = load i32, i32* %7
%39 = sext i32 %38 to i64
%40 = getelementptr inbounds [32 x float], [32 x float]* %37, i64 %39
%41 = load i32, i32* %8
%42 = sext i32 %41 to i64
%43 = getelementptr inbounds [32 x float], [32 x float]* %40, i64 0, i64 %42
store float %36, float* %43
br label %95
44: ; preds = %17
%45 = load i32, i32* %7
%46 = load i32, i32* %8
%47 = icmp sgt i32 %45, %46
br i1 %47, label %48, label %71
48: ; preds = %44
%49 = load [32 x float]*, [32 x float]** %4
%50 = load i32, i32* %7
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds [32 x float], [32 x float]* %49, i64 %51
%53 = load i32, i32* %8
%54 = sext i32 %53 to i64
%55 = getelementptr inbounds [32 x float], [32 x float]* %52, i64 0, i64 %54
%56 = load float, float* %55
%57 = load [32 x float]*, [32 x float]** %5
%58 = load i32, i32* %7
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds [32 x float], [32 x float]* %57, i64 %59
%61 = load i32, i32* %8
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds [32 x float], [32 x float]* %60, i64 0, i64 %62
store float %56, float* %63
%64 = load [32 x float]*, [32 x float]** %6
%65 = load i32, i32* %7
%66 = sext i32 %65 to i64
%67 = getelementptr inbounds [32 x float], [32 x float]* %64, i64 %66
%68 = load i32, i32* %8
%69 = sext i32 %68 to i64
%70 = getelementptr inbounds [32 x float], [32 x float]* %67, i64 0, i64 %69
store float 0.000000e+00, float* %70
br label %94
71: ; preds = %44
%72 = load [32 x float]*, [32 x float]** %5
%73 = load i32, i32* %7
%74 = sext i32 %73 to i64
%75 = getelementptr inbounds [32 x float], [32 x float]* %72, i64 %74
%76 = load i32, i32* %8
%77 = sext i32 %76 to i64
%78 = getelementptr inbounds [32 x float], [32 x float]* %75, i64 0, i64 %77
store float 0.000000e+00, float* %78
%79 = load [32 x float]*, [32 x float]** %4
%80 = load i32, i32* %7
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds [32 x float], [32 x float]* %79, i64 %81
%83 = load i32, i32* %8
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds [32 x float], [32 x float]* %82, i64 0, i64 %84
%86 = load float, float* %85
%87 = load [32 x float]*, [32 x float]** %6
%88 = load i32, i32* %7
%89 = sext i32 %88 to i64
%90 = getelementptr inbounds [32 x float], [32 x float]* %87, i64 %89
%91 = load i32, i32* %8
%92 = sext i32 %91 to i64
%93 = getelementptr inbounds [32 x float], [32 x float]* %90, i64 0, i64 %92
store float %86, float* %93
br label %94
94: ; preds = %71, %48
br label %95
95: ; preds = %94, %21
br label %96
96: ; preds = %95
%97 = load i32, i32* %8
%98 = add nsw i32 %97, 1
store i32 %98, i32* %8
br label %14
99: ; preds = %14
br label %100
100: ; preds = %99
%101 = load i32, i32* %7
%102 = add nsw i32 %101, 1
store i32 %102, i32* %7
br label %10
103: ; preds = %10
ret void
}

LLVM IR (-O2):
define void @split_block([32 x float]* nocapture %0, [32 x float]* nocapture writeonly %1, [32 x float]* nocapture writeonly %2) {
br label %4
4: ; preds = %3, %29
%5 = phi i64 [ 0, %3 ], [ %30, %29 ]
%6 = getelementptr inbounds [32 x float], [32 x float]* %1, i64 %5, i64 %5
%7 = getelementptr inbounds [32 x float], [32 x float]* %0, i64 %5, i64 %5
%8 = getelementptr inbounds [32 x float], [32 x float]* %2, i64 %5, i64 %5
br label %9
9: ; preds = %4, %26
%10 = phi i64 [ 0, %4 ], [ %27, %26 ]
%11 = icmp eq i64 %5, %10
br i1 %11, label %12, label %14
12: ; preds = %9
store float 1.000000e+00, float* %6
%13 = load float, float* %7
store float %13, float* %8
br label %26
14: ; preds = %9
%15 = icmp ugt i64 %5, %10
br i1 %15, label %16, label %21
16: ; preds = %14
%17 = getelementptr inbounds [32 x float], [32 x float]* %0, i64 %5, i64 %10
%18 = load float, float* %17
%19 = getelementptr inbounds [32 x float], [32 x float]* %1, i64 %5, i64 %10
store float %18, float* %19
%20 = getelementptr inbounds [32 x float], [32 x float]* %2, i64 %5, i64 %10
store float 0.000000e+00, float* %20
br label %26
21: ; preds = %14
%22 = getelementptr inbounds [32 x float], [32 x float]* %1, i64 %5, i64 %10
store float 0.000000e+00, float* %22
%23 = getelementptr inbounds [32 x float], [32 x float]* %0, i64 %5, i64 %10
%24 = load float, float* %23
%25 = getelementptr inbounds [32 x float], [32 x float]* %2, i64 %5, i64 %10
store float %24, float* %25
br label %26
26: ; preds = %12, %21, %16
%27 = add nuw nsw i64 %10, 1
%28 = icmp eq i64 %27, 32
br i1 %28, label %29, label %9
29: ; preds = %26
%30 = add nuw nsw i64 %5, 1
%31 = icmp eq i64 %30, 32
br i1 %31, label %32, label %4
32: ; preds = %29
ret void
}