void TriDiagSolver(double *d, double *ld, double *ud, double *b, double *res,
                   int size_d) {
  double den;
  double num;
  int i, j;
  int size_ld, size_ud;
  size_ld = size_d - 1;
  size_ud = size_ld;
  ud[0] = ud[0] / d[0];
  b[0] = b[0] / d[0];
  for (i = 1; i < size_ud; i++) {
    den = d[i] - ld[i - 1] * ud[i - 1];
    num = b[i] - ld[i - 1] * b[i - 1];
    ud[i] /= den;
    b[i] = num / den;
  }
  den = d[size_d - 1] - ld[size_ld - 1] * ud[size_ud - 1];
  num = b[size_d - 1] - ld[size_ld - 1] * b[size_d - 2];
  b[size_d - 1] = num / den;
  res[size_d - 1] = b[size_d - 1];
  for (j = size_d - 2; j >= 0; j--)
    res[j] = b[j] - ud[j] * res[j + 1];
}


x86 assembly (clang, -O0):
<TriDiagSolver>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9d,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
mov    %eax,-0x50(%rbp)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x8(%rbp),%rax
divsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x8(%rbp),%rax
divsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x1,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x50(%rbp),%eax
jge    1258 <TriDiagSolver+0x158>
mov    -0x8(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x44(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x38(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x40(%rbp),%xmm0
divsd  -0x38(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1162 <TriDiagSolver+0x62>
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x4c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x50(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x4c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
divsd  -0x38(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x48(%rbp)
cmpl   $0x0,-0x48(%rbp)
jl     13c5 <TriDiagSolver+0x2c5>
mov    -0x20(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x48(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x48(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x48(%rbp)
jmp    1355 <TriDiagSolver+0x255>
pop    %rbp
ret

x86 assembly (clang, -O2):
<TriDiagSolver>:
push   %rbx
lea    -0x1(%r9),%r10d
movsd  (%rdx),%xmm0
divsd  (%rdi),%xmm0
movsd  %xmm0,(%rdx)
movsd  (%rcx),%xmm0
divsd  (%rdi),%xmm0
movsd  %xmm0,(%rcx)
cmp    $0x3,%r9d
jl     11fa <TriDiagSolver+0xfa>
mov    %r10d,%r11d
lea    0x8(%rcx),%rax
lea    (%rcx,%r11,8),%rbx
cmp    %rdx,%rbx
jbe    119d <TriDiagSolver+0x9d>
lea    (%rdx,%r11,8),%rbx
cmp    %rbx,%rax
jae    119d <TriDiagSolver+0x9d>
add    $0xffffffffffffffff,%r11
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
movsd  0x8(%rdi,%rbx,8),%xmm1
movsd  (%rsi,%rbx,8),%xmm2
movsd  (%rdx,%rbx,8),%xmm3
mulsd  %xmm2,%xmm3
subsd  %xmm3,%xmm1
movsd  (%rax,%rbx,8),%xmm3
mulsd  %xmm2,%xmm0
subsd  %xmm0,%xmm3
movsd  0x8(%rdx,%rbx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,0x8(%rdx,%rbx,8)
divsd  %xmm1,%xmm3
movsd  %xmm3,(%rax,%rbx,8)
add    $0x1,%rbx
movapd %xmm3,%xmm0
cmp    %rbx,%r11
jne    1150 <TriDiagSolver+0x50>
jmp    11fa <TriDiagSolver+0xfa>
movhpd (%rdx),%xmm0
add    $0xffffffffffffffff,%r11
xor    %ebx,%ebx
nopw   0x0(%rax,%rax,1)
movsd  0x8(%rdi,%rbx,8),%xmm1
movsd  (%rsi,%rbx,8),%xmm2
movapd %xmm0,%xmm3
unpckhpd %xmm0,%xmm3
mulsd  %xmm2,%xmm3
subsd  %xmm3,%xmm1
mulsd  %xmm0,%xmm2
movsd  (%rax,%rbx,8),%xmm0
subsd  %xmm2,%xmm0
movhpd 0x8(%rdx,%rbx,8),%xmm0
unpcklpd %xmm1,%xmm1
divpd  %xmm1,%xmm0
movhpd %xmm0,0x8(%rdx,%rbx,8)
movlpd %xmm0,(%rax,%rbx,8)
add    $0x1,%rbx
cmp    %rbx,%r11
jne    11b0 <TriDiagSolver+0xb0>
movslq %r10d,%rax
movsd  (%rdi,%rax,8),%xmm0
movslq %r9d,%rdi
movsd  -0x10(%rsi,%rdi,8),%xmm1
movsd  -0x10(%rdx,%rdi,8),%xmm2
mulsd  %xmm1,%xmm2
subsd  %xmm2,%xmm0
movsd  (%rcx,%rax,8),%xmm2
mulsd  -0x10(%rcx,%rdi,8),%xmm1
subsd  %xmm1,%xmm2
divsd  %xmm0,%xmm2
movsd  %xmm2,(%rcx,%rax,8)
movsd  %xmm2,(%r8,%rax,8)
cmp    $0x2,%r9d
jl     12ac <TriDiagSolver+0x1ac>
add    $0xfffffffe,%r9d
mov    %r9d,%eax
test   $0x1,%r9b
jne    1269 <TriDiagSolver+0x169>
movsd  (%rdx,%rax,8),%xmm0
mulsd  0x8(%r8,%rax,8),%xmm0
movsd  (%rcx,%rax,8),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,(%r8,%rax,8)
add    $0xffffffffffffffff,%rax
test   %r9d,%r9d
je     12ac <TriDiagSolver+0x1ac>
xchg   %ax,%ax
movsd  (%rcx,%rax,8),%xmm0
movsd  (%rdx,%rax,8),%xmm1
mulsd  0x8(%r8,%rax,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,(%r8,%rax,8)
movsd  -0x8(%rcx,%rax,8),%xmm1
mulsd  -0x8(%rdx,%rax,8),%xmm0
subsd  %xmm0,%xmm1
movsd  %xmm1,-0x8(%r8,%rax,8)
add    $0xfffffffffffffffe,%rax
cmp    $0xffffffffffffffff,%rax
jne    1270 <TriDiagSolver+0x170>
pop    %rbx
ret

LLVM IR (-O0):
define void @TriDiagSolver(double* %0, double* %1, double* %2, double* %3, double* %4, i32 %5) {
%7 = alloca double*
%8 = alloca double*
%9 = alloca double*
%10 = alloca double*
%11 = alloca double*
%12 = alloca i32
%13 = alloca double
%14 = alloca double
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
store double* %0, double** %7
store double* %1, double** %8
store double* %2, double** %9
store double* %3, double** %10
store double* %4, double** %11
store i32 %5, i32* %12
%19 = load i32, i32* %12
%20 = sub nsw i32 %19, 1
store i32 %20, i32* %17
%21 = load i32, i32* %17
store i32 %21, i32* %18
%22 = load double*, double** %9
%23 = getelementptr inbounds double, double* %22, i64 0
%24 = load double, double* %23
%25 = load double*, double** %7
%26 = getelementptr inbounds double, double* %25, i64 0
%27 = load double, double* %26
%28 = fdiv double %24, %27
%29 = load double*, double** %9
%30 = getelementptr inbounds double, double* %29, i64 0
store double %28, double* %30
%31 = load double*, double** %10
%32 = getelementptr inbounds double, double* %31, i64 0
%33 = load double, double* %32
%34 = load double*, double** %7
%35 = getelementptr inbounds double, double* %34, i64 0
%36 = load double, double* %35
%37 = fdiv double %33, %36
%38 = load double*, double** %10
%39 = getelementptr inbounds double, double* %38, i64 0
store double %37, double* %39
store i32 1, i32* %15
br label %40
40: ; preds = %97, %6
%41 = load i32, i32* %15
%42 = load i32, i32* %18
%43 = icmp slt i32 %41, %42
br i1 %43, label %44, label %100
44: ; preds = %40
%45 = load double*, double** %7
%46 = load i32, i32* %15
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds double, double* %45, i64 %47
%49 = load double, double* %48
%50 = load double*, double** %8
%51 = load i32, i32* %15
%52 = sub nsw i32 %51, 1
%53 = sext i32 %52 to i64
%54 = getelementptr inbounds double, double* %50, i64 %53
%55 = load double, double* %54
%56 = load double*, double** %9
%57 = load i32, i32* %15
%58 = sub nsw i32 %57, 1
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds double, double* %56, i64 %59
%61 = load double, double* %60
%62 = fneg double %55
%63 = call double @llvm.fmuladd.f64(double %62, double %61, double %49)
store double %63, double* %13
%64 = load double*, double** %10
%65 = load i32, i32* %15
%66 = sext i32 %65 to i64
%67 = getelementptr inbounds double, double* %64, i64 %66
%68 = load double, double* %67
%69 = load double*, double** %8
%70 = load i32, i32* %15
%71 = sub nsw i32 %70, 1
%72 = sext i32 %71 to i64
%73 = getelementptr inbounds double, double* %69, i64 %72
%74 = load double, double* %73
%75 = load double*, double** %10
%76 = load i32, i32* %15
%77 = sub nsw i32 %76, 1
%78 = sext i32 %77 to i64
%79 = getelementptr inbounds double, double* %75, i64 %78
%80 = load double, double* %79
%81 = fneg double %74
%82 = call double @llvm.fmuladd.f64(double %81, double %80, double %68)
store double %82, double* %14
%83 = load double, double* %13
%84 = load double*, double** %9
%85 = load i32, i32* %15
%86 = sext i32 %85 to i64
%87 = getelementptr inbounds double, double* %84, i64 %86
%88 = load double, double* %87
%89 = fdiv double %88, %83
store double %89, double* %87
%90 = load double, double* %14
%91 = load double, double* %13
%92 = fdiv double %90, %91
%93 = load double*, double** %10
%94 = load i32, i32* %15
%95 = sext i32 %94 to i64
%96 = getelementptr inbounds double, double* %93, i64 %95
store double %92, double* %96
br label %97
97: ; preds = %44
%98 = load i32, i32* %15
%99 = add nsw i32 %98, 1
store i32 %99, i32* %15
br label %40
100: ; preds = %40
%101 = load double*, double** %7
%102 = load i32, i32* %12
%103 = sub nsw i32 %102, 1
%104 = sext i32 %103 to i64
%105 = getelementptr inbounds double, double* %101, i64 %104
%106 = load double, double* %105
%107 = load double*, double** %8
%108 = load i32, i32* %17
%109 = sub nsw i32 %108, 1
%110 = sext i32 %109 to i64
%111 = getelementptr inbounds double, double* %107, i64 %110
%112 = load double, double* %111
%113 = load double*, double** %9
%114 = load i32, i32* %18
%115 = sub nsw i32 %114, 1
%116 = sext i32 %115 to i64
%117 = getelementptr inbounds double, double* %113, i64 %116
%118 = load double, double* %117
%119 = fneg double %112
%120 = call double @llvm.fmuladd.f64(double %119, double %118, double %106)
store double %120, double* %13
%121 = load double*, double** %10
%122 = load i32, i32* %12
%123 = sub nsw i32 %122, 1
%124 = sext i32 %123 to i64
%125 = getelementptr inbounds double, double* %121, i64 %124
%126 = load double, double* %125
%127 = load double*, double** %8
%128 = load i32, i32* %17
%129 = sub nsw i32 %128, 1
%130 = sext i32 %129 to i64
%131 = getelementptr inbounds double, double* %127, i64 %130
%132 = load double, double* %131
%133 = load double*, double** %10
%134 = load i32, i32* %12
%135 = sub nsw i32 %134, 2
%136 = sext i32 %135 to i64
%137 = getelementptr inbounds double, double* %133, i64 %136
%138 = load double, double* %137
%139 = fneg double %132
%140 = call double @llvm.fmuladd.f64(double %139, double %138, double %126)
store double %140, double* %14
%141 = load double, double* %14
%142 = load double, double* %13
%143 = fdiv double %141, %142
%144 = load double*, double** %10
%145 = load i32, i32* %12
%146 = sub nsw i32 %145, 1
%147 = sext i32 %146 to i64
%148 = getelementptr inbounds double, double* %144, i64 %147
store double %143, double* %148
%149 = load double*, double** %10
%150 = load i32, i32* %12
%151 = sub nsw i32 %150, 1
%152 = sext i32 %151 to i64
%153 = getelementptr inbounds double, double* %149, i64 %152
%154 = load double, double* %153
%155 = load double*, double** %11
%156 = load i32, i32* %12
%157 = sub nsw i32 %156, 1
%158 = sext i32 %157 to i64
%159 = getelementptr inbounds double, double* %155, i64 %158
store double %154, double* %159
%160 = load i32, i32* %12
%161 = sub nsw i32 %160, 2
store i32 %161, i32* %16
br label %162
162: ; preds = %188, %100
%163 = load i32, i32* %16
%164 = icmp sge i32 %163, 0
br i1 %164, label %165, label %191
165: ; preds = %162
%166 = load double*, double** %10
%167 = load i32, i32* %16
%168 = sext i32 %167 to i64
%169 = getelementptr inbounds double, double* %166, i64 %168
%170 = load double, double* %169
%171 = load double*, double** %9
%172 = load i32, i32* %16
%173 = sext i32 %172 to i64
%174 = getelementptr inbounds double, double* %171, i64 %173
%175 = load double, double* %174
%176 = load double*, double** %11
%177 = load i32, i32* %16
%178 = add nsw i32 %177, 1
%179 = sext i32 %178 to i64
%180 = getelementptr inbounds double, double* %176, i64 %179
%181 = load double, double* %180
%182 = fneg double %175
%183 = call double @llvm.fmuladd.f64(double %182, double %181, double %170)
%184 = load double*, double** %11
%185 = load i32, i32* %16
%186 = sext i32 %185 to i64
%187 = getelementptr inbounds double, double* %184, i64 %186
store double %183, double* %187
br label %188
188: ; preds = %165
%189 = load i32, i32* %16
%190 = add nsw i32 %189, -1
store i32 %190, i32* %16
br label %162
191: ; preds = %162
ret void
}

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