void lu(int n, double *a, double *l, double *u) {
  int i, j, k;
  for (i = 0; i < n * n; i++) {
    u[i] = a[i];
    if ((i / n) == (i % n))
      l[i] = 1;
    else
      l[i] = 0;
  }
  for (j = 0; j < n - 1; j++) {
    for (i = j + 1; i < n; i++) {
      double m = -u[i * n + j] / u[j * n + j];
      for (k = j; k < n; k++)
        u[i * n + k] = (m * u[j * n + k]) + u[i * n + k];
      l[i * n + j] = -m;
    }
  }
}


x86 assembly (clang, -O0):
<lu>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x4(%rbp),%ecx
imul   -0x4(%rbp),%ecx
cmp    %ecx,%eax
jge    119f <lu+0x9f>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
cltd
idivl  -0x4(%rbp)
mov    %eax,-0x3c(%rbp)
mov    -0x24(%rbp),%eax
cltd
idivl  -0x4(%rbp)
mov    -0x3c(%rbp),%eax
cmp    %edx,%eax
jne    117c <lu+0x7c>
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  0xe8e(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    118c <lu+0x8c>
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1191 <lu+0x91>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    111a <lu+0x1a>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    12d6 <lu+0x1d6>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    12c3 <lu+0x1c3>
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x20(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x28(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1283 <lu+0x183>
movsd  -0x38(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    121a <lu+0x11a>
movsd  -0x38(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11c0 <lu+0xc0>
jmp    12c8 <lu+0x1c8>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    11a6 <lu+0xa6>
pop    %rbp
ret

x86 assembly (clang, -O2):
<lu>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %rdx,-0x48(%rsp)
mov    %edi,%eax
imul   %edi,%eax
test   %eax,%eax
je     1164 <lu+0x64>
mov    %eax,%ebp
xor    %r9d,%r9d
movsd  0xed7(%rip),%xmm0
jmp    1144 <lu+0x44>
nopl   0x0(%rax,%rax,1)
mov    -0x48(%rsp),%rax
movsd  %xmm1,(%rax,%r9,8)
add    $0x1,%r9
cmp    %r9,%rbp
je     1164 <lu+0x64>
movsd  (%rsi,%r9,8),%xmm1
movsd  %xmm1,(%rcx,%r9,8)
mov    %r9d,%eax
cltd
idiv   %edi
movapd %xmm0,%xmm1
cmp    %edx,%eax
je     1130 <lu+0x30>
xorpd  %xmm1,%xmm1
jmp    1130 <lu+0x30>
mov    %rdi,-0x78(%rsp)
cmp    $0x2,%edi
jl     144d <lu+0x34d>
mov    -0x78(%rsp),%rax
lea    -0x1(%rax),%edx
mov    %rdx,-0x30(%rsp)
mov    %eax,%edx
lea    0x0(,%rdx,8),%rdi
lea    0x8(,%rdx,8),%rax
mov    %rax,-0x50(%rsp)
lea    (%rdx,%rdx,1),%rax
mov    %rax,-0x38(%rsp)
mov    %rdx,%rax
neg    %rax
mov    %rax,0x18(%rsp)
lea    0x10(%rcx),%rbx
lea    (%rcx,%rdx,8),%rax
add    $0x10,%rax
mov    %rax,-0x58(%rsp)
lea    0x8(%rcx,%rdx,8),%rax
mov    %rax,-0x60(%rsp)
mov    $0x1,%eax
mov    %rax,-0x68(%rsp)
xor    %esi,%esi
movapd 0xe37(%rip),%xmm0
mov    %rcx,%rax
mov    %rdx,-0x70(%rsp)
jmp    1221 <lu+0x121>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
addq   $0x1,-0x68(%rsp)
addq   $0xffffffffffffffff,-0x70(%rsp)
mov    -0x50(%rsp),%rsi
add    %rsi,%rbx
add    %rsi,-0x58(%rsp)
add    %rdi,-0x60(%rsp)
add    %rdi,%rax
mov    -0x28(%rsp),%rsi
cmp    -0x30(%rsp),%rsi
je     144d <lu+0x34d>
mov    %rsi,-0x80(%rsp)
add    $0x1,%rsi
mov    %rsi,-0x28(%rsp)
cmp    %rdx,%rsi
jae    11f0 <lu+0xf0>
mov    -0x70(%rsp),%rbp
and    $0xfffffffffffffffc,%rbp
mov    %rdx,%r9
mov    -0x80(%rsp),%r8
sub    %r8,%r9
mov    -0x50(%rsp),%r10
imul   %r8,%r10
mov    %r8,%r11
imul   %rdx,%r11
mov    -0x38(%rsp),%rsi
add    %r11,%rsi
mov    %rsi,-0x8(%rsp)
mov    %r10,0x10(%rsp)
lea    (%rcx,%r10,1),%rsi
mov    %rsi,-0x10(%rsp)
mov    %r11,(%rsp)
lea    (%r11,%rdx,1),%rsi
lea    (%rcx,%rsi,8),%rsi
mov    %rsi,0x8(%rsp)
mov    %r8d,%esi
imul   -0x78(%rsp),%esi
movslq %esi,%rsi
add    %r8,%rsi
mov    %rsi,0x20(%rsp)
mov    %r9,-0x40(%rsp)
and    $0xfffffffffffffffc,%r9
mov    %r9,-0x18(%rsp)
lea    (%r8,%r9,1),%rsi
mov    %rsi,-0x20(%rsp)
mov    -0x60(%rsp),%r12
mov    -0x58(%rsp),%r11
xor    %r14d,%r14d
mov    -0x68(%rsp),%r8
jmp    12f6 <lu+0x1f6>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
xorpd  %xmm0,%xmm1
mov    -0x48(%rsp),%rsi
movlpd %xmm1,(%rsi,%r13,8)
add    $0x1,%r8
add    $0x1,%r14
add    %rdi,%r11
add    %rdi,%r12
cmp    %rdx,%r8
je     11f0 <lu+0xf0>
mov    %r8,%r9
imul   %rdx,%r9
mov    -0x80(%rsp),%r15
lea    (%r9,%r15,1),%r13
movsd  (%rcx,%r13,8),%xmm1
xorpd  %xmm0,%xmm1
mov    0x20(%rsp),%rsi
divsd  (%rcx,%rsi,8),%xmm1
cmpq   $0x4,-0x40(%rsp)
jb     13d0 <lu+0x2d0>
mov    %r14,%rsi
imul   %rdx,%rsi
mov    %rdi,%r10
lea    (%rsi,%rdx,1),%rdi
lea    (%rcx,%rdi,8),%rdi
add    0x10(%rsp),%rdi
cmp    0x8(%rsp),%rdi
mov    %r10,%rdi
jae    135c <lu+0x25c>
add    -0x8(%rsp),%rsi
lea    (%rcx,%rsi,8),%rsi
mov    -0x80(%rsp),%r15
cmp    %rsi,-0x10(%rsp)
jb     13d0 <lu+0x2d0>
movapd %xmm1,%xmm2
unpcklpd %xmm1,%xmm2
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
movupd -0x10(%rbx,%rsi,8),%xmm3
movupd (%rbx,%rsi,8),%xmm4
movupd -0x10(%r11,%rsi,8),%xmm5
movupd (%r11,%rsi,8),%xmm6
mulpd  %xmm2,%xmm3
addpd  %xmm5,%xmm3
mulpd  %xmm2,%xmm4
addpd  %xmm6,%xmm4
movupd %xmm3,-0x10(%r11,%rsi,8)
movupd %xmm4,(%r11,%rsi,8)
add    $0x4,%rsi
cmp    %rsi,%rbp
jne    1370 <lu+0x270>
mov    -0x20(%rsp),%r15
mov    -0x18(%rsp),%rsi
cmp    %rsi,-0x40(%rsp)
je     12d0 <lu+0x1d0>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %edx,%esi
sub    %r15d,%esi
mov    %r15,%r10
test   $0x1,%sil
je     1401 <lu+0x301>
mov    (%rsp),%rsi
add    %r15,%rsi
add    %r15,%r9
movsd  (%rcx,%rsi,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  (%rcx,%r9,8),%xmm2
movsd  %xmm2,(%rcx,%r9,8)
lea    0x1(%r15),%r10
not    %r15
cmp    0x18(%rsp),%r15
je     12d0 <lu+0x1d0>
nop
movsd  (%rax,%r10,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  -0x8(%r12,%r10,8),%xmm2
movsd  %xmm2,-0x8(%r12,%r10,8)
movsd  0x8(%rax,%r10,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  (%r12,%r10,8),%xmm2
movsd  %xmm2,(%r12,%r10,8)
add    $0x2,%r10
cmp    %r10,%rdx
jne    1410 <lu+0x310>
jmp    12d0 <lu+0x1d0>
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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