int linear_expansion_sum(elen, e, flen, f, h)
int elen;
double *e;
int flen;
double *f;
double *h;
{
  double Q, q;
  double Qnew;
  double R;
  double bvirt;
  double avirt, bround, around;
  int eindex, findex, hindex;
  double enow, fnow;
  double g0;
  enow = e[0];
  fnow = f[0];
  eindex = findex = 0;
  if ((fnow > enow) == (fnow > -enow)) {
    g0 = enow;
    enow = e[++eindex];
  } else {
    g0 = fnow;
    fnow = f[++findex];
  }
  if ((eindex < elen) &&
      ((findex >= flen) || ((fnow > enow) == (fnow > -enow)))) {
    Qnew = (double)(enow + g0);
    bvirt = Qnew - enow;
    q = g0 - bvirt;
    enow = e[++eindex];
  } else {
    Qnew = (double)(fnow + g0);
    bvirt = Qnew - fnow;
    q = g0 - bvirt;
    fnow = f[++findex];
  }
  Q = Qnew;
  for (hindex = 0; hindex < elen + flen - 2; hindex++) {
    if ((eindex < elen) &&
        ((findex >= flen) || ((fnow > enow) == (fnow > -enow)))) {
      R = (double)(enow + q);
      bvirt = R - enow;
      h[hindex] = q - bvirt;
      enow = e[++eindex];
    } else {
      R = (double)(fnow + q);
      bvirt = R - fnow;
      h[hindex] = q - bvirt;
      fnow = f[++findex];
    }
    Qnew = (double)(Q + R);
    bvirt = (double)(Qnew - Q);
    avirt = Qnew - bvirt;
    bround = R - bvirt;
    around = Q - avirt;
    q = around + bround;
    Q = Qnew;
  }
  h[hindex] = q;
  h[hindex + 1] = Q;
  return hindex + 2;
}


x86 assembly (clang, -O0):
<linear_expansion_sum>:
push   %rbp
mov    %rsp,%rbp
sub    $0x10,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movl   $0x0,-0x70(%rbp)
movl   $0x0,-0x6c(%rbp)
movsd  -0x88(%rbp),%xmm0
ff
ucomisd -0x80(%rbp),%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
movsd  -0x88(%rbp),%xmm0
ff
movsd  -0x80(%rbp),%xmm1
movq   %xmm1,%rcx
movabs $0x8000000000000000,%rdx
xor    %rdx,%rcx
movq   %rcx,%xmm1
ucomisd %xmm1,%xmm0
seta   %cl
and    $0x1,%cl
movzbl %cl,%ecx
cmp    %ecx,%eax
jne    11bf <linear_expansion_sum+0xbf>
movsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    -0x6c(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x6c(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
jmp    11ec <linear_expansion_sum+0xec>
movsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0x90(%rbp)
ff
mov    -0x20(%rbp),%rax
mov    -0x70(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x70(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
mov    -0x6c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    12a4 <linear_expansion_sum+0x1a4>
mov    -0x70(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1252 <linear_expansion_sum+0x152>
movsd  -0x88(%rbp),%xmm0
ff
ucomisd -0x80(%rbp),%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
movsd  -0x88(%rbp),%xmm0
ff
movsd  -0x80(%rbp),%xmm1
movq   %xmm1,%rcx
movabs $0x8000000000000000,%rdx
xor    %rdx,%rcx
movq   %rcx,%xmm1
ucomisd %xmm1,%xmm0
seta   %cl
and    $0x1,%cl
movzbl %cl,%ecx
cmp    %ecx,%eax
jne    12a4 <linear_expansion_sum+0x1a4>
movsd  -0x80(%rbp),%xmm0
addsd  -0x90(%rbp),%xmm0
ff
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
subsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x90(%rbp),%xmm0
ff
subsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x6c(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x6c(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
jmp    12fa <linear_expansion_sum+0x1fa>
movsd  -0x88(%rbp),%xmm0
ff
addsd  -0x90(%rbp),%xmm0
ff
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
subsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0x50(%rbp)
movsd  -0x90(%rbp),%xmm0
ff
subsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x70(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x70(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x0,-0x74(%rbp)
mov    -0x74(%rbp),%eax
mov    -0x4(%rbp),%ecx
add    -0x14(%rbp),%ecx
sub    $0x2,%ecx
cmp    %ecx,%eax
jge    14a3 <linear_expansion_sum+0x3a3>
mov    -0x6c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    13d9 <linear_expansion_sum+0x2d9>
mov    -0x70(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1385 <linear_expansion_sum+0x285>
movsd  -0x88(%rbp),%xmm0
ff
ucomisd -0x80(%rbp),%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
movsd  -0x88(%rbp),%xmm0
ff
movsd  -0x80(%rbp),%xmm1
movq   %xmm1,%rcx
movabs $0x8000000000000000,%rdx
xor    %rdx,%rcx
movq   %rcx,%xmm1
ucomisd %xmm1,%xmm0
seta   %cl
and    $0x1,%cl
movzbl %cl,%ecx
cmp    %ecx,%eax
jne    13d9 <linear_expansion_sum+0x2d9>
movsd  -0x80(%rbp),%xmm0
addsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
subsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x38(%rbp),%xmm0
subsd  -0x50(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x74(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
mov    -0x6c(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x6c(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
jmp    1431 <linear_expansion_sum+0x331>
movsd  -0x88(%rbp),%xmm0
ff
addsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
subsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0x50(%rbp)
movsd  -0x38(%rbp),%xmm0
subsd  -0x50(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x74(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
mov    -0x70(%rbp),%ecx
add    $0x1,%ecx
mov    %ecx,-0x70(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x30(%rbp),%xmm0
addsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
subsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x40(%rbp),%xmm0
subsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x48(%rbp),%xmm0
subsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x30(%rbp),%xmm0
subsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x68(%rbp),%xmm0
addsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x74(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x74(%rbp)
jmp    130b <linear_expansion_sum+0x20b>
movsd  -0x38(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x74(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x30(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x74(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x74(%rbp),%eax
add    $0x2,%eax
add    $0x10,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<linear_expansion_sum>:
push   %rbp
push   %r14
push   %rbx
movsd  (%rsi),%xmm0
movsd  (%rcx),%xmm3
ucomisd %xmm0,%xmm3
setbe  %al
movapd 0xee5(%rip),%xmm1
xorpd  %xmm0,%xmm1
ucomisd %xmm1,%xmm3
seta   %bl
cmp    %bl,%al
je     1143 <linear_expansion_sum+0x43>
movsd  0x8(%rsi),%xmm1
xor    %r9d,%r9d
mov    $0x1,%r10d
movapd %xmm3,%xmm2
cmp    %edi,%r10d
jl     115e <linear_expansion_sum+0x5e>
jmp    11a0 <linear_expansion_sum+0xa0>
movsd  0x8(%rcx),%xmm2
mov    $0x1,%r9d
xor    %r10d,%r10d
movapd %xmm0,%xmm1
movapd %xmm3,%xmm0
cmp    %edi,%r10d
jge    11a0 <linear_expansion_sum+0xa0>
cmp    %edx,%r9d
jge    1181 <linear_expansion_sum+0x81>
ucomisd %xmm1,%xmm2
setbe  %al
movapd 0xe8e(%rip),%xmm3
xorpd  %xmm1,%xmm3
ucomisd %xmm3,%xmm2
seta   %bl
cmp    %bl,%al
je     11a0 <linear_expansion_sum+0xa0>
movapd %xmm1,%xmm3
addsd  %xmm0,%xmm3
movapd %xmm3,%xmm4
subsd  %xmm1,%xmm4
mov    %r10d,%eax
add    $0x1,%r10d
movsd  0x8(%rsi,%rax,8),%xmm1
jmp    11bd <linear_expansion_sum+0xbd>
movapd %xmm2,%xmm3
addsd  %xmm0,%xmm3
movapd %xmm3,%xmm4
subsd  %xmm2,%xmm4
mov    %r9d,%eax
add    $0x1,%r9d
movsd  0x8(%rcx,%rax,8),%xmm2
subsd  %xmm4,%xmm0
lea    (%rdi,%rdx,1),%eax
add    $0xfffffffe,%eax
test   %eax,%eax
jle    1298 <linear_expansion_sum+0x198>
mov    %eax,%r11d
xor    %ebp,%ebp
movapd 0xe24(%rip),%xmm4
jmp    123c <linear_expansion_sum+0x13c>
xchg   %ax,%ax
movapd %xmm0,%xmm5
addsd  %xmm2,%xmm5
movapd %xmm5,%xmm6
subsd  %xmm2,%xmm6
subsd  %xmm6,%xmm0
movsd  %xmm0,(%r8,%rbp,8)
movslq %r9d,%rbx
add    $0x1,%r9d
movsd  0x8(%rcx,%rbx,8),%xmm2
movapd %xmm3,%xmm6
addsd  %xmm5,%xmm6
movapd %xmm6,%xmm0
subsd  %xmm3,%xmm0
movapd %xmm6,%xmm7
subsd  %xmm0,%xmm7
subsd  %xmm0,%xmm5
subsd  %xmm7,%xmm3
movapd %xmm3,%xmm0
addsd  %xmm5,%xmm0
add    $0x1,%rbp
movapd %xmm6,%xmm3
cmp    %rbp,%r11
je     1292 <linear_expansion_sum+0x192>
cmp    %edi,%r10d
jge    11e0 <linear_expansion_sum+0xe0>
cmp    %edx,%r9d
jge    1266 <linear_expansion_sum+0x166>
ucomisd %xmm1,%xmm2
setbe  %r14b
movapd %xmm1,%xmm5
xorpd  %xmm4,%xmm5
ucomisd %xmm5,%xmm2
seta   %bl
cmp    %bl,%r14b
je     11e0 <linear_expansion_sum+0xe0>
movapd %xmm0,%xmm5
addsd  %xmm1,%xmm5
movapd %xmm5,%xmm6
subsd  %xmm1,%xmm6
subsd  %xmm6,%xmm0
movsd  %xmm0,(%r8,%rbp,8)
movslq %r10d,%rbx
add    $0x1,%r10d
movsd  0x8(%rsi,%rbx,8),%xmm1
jmp    1207 <linear_expansion_sum+0x107>
movapd %xmm6,%xmm3
jmp    129a <linear_expansion_sum+0x19a>
xor    %eax,%eax
mov    %eax,%ecx
movsd  %xmm0,(%r8,%rcx,8)
movsd  %xmm3,0x8(%r8,%rcx,8)
add    $0x2,%eax
pop    %rbx
pop    %r14
pop    %rbp
ret

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

LLVM IR (-O2):
define i32 @linear_expansion_sum(i32 %0, double* nocapture %1, i32 %2, double* nocapture %3, double* nocapture writeonly %4) {
%6 = load double, double* %1
%7 = load double, double* %3
%8 = fcmp ule double %7, %6
%9 = fneg double %6
%10 = fcmp ogt double %7, %9
%11 = xor i1 %8, %10
br i1 %11, label %12, label %15
12: ; preds = %5
%13 = getelementptr inbounds double, double* %1, i64 1
%14 = load double, double* %13
br label %18
15: ; preds = %5
%16 = getelementptr inbounds double, double* %3, i64 1
%17 = load double, double* %16
br label %18
18: ; preds = %15, %12
%19 = phi i32 [ 1, %12 ], [ 0, %15 ]
%20 = phi i32 [ 0, %12 ], [ 1, %15 ]
%21 = phi double [ %14, %12 ], [ %6, %15 ]
%22 = phi double [ %7, %12 ], [ %17, %15 ]
%23 = phi double [ %6, %12 ], [ %7, %15 ]
%24 = icmp slt i32 %19, %0
br i1 %24, label %25, label %39
25: ; preds = %18
%26 = icmp slt i32 %20, %2
br i1 %26, label %27, label %32
27: ; preds = %25
%28 = fcmp ule double %22, %21
%29 = fneg double %21
%30 = fcmp ogt double %22, %29
%31 = xor i1 %28, %30
br i1 %31, label %32, label %39
32: ; preds = %27, %25
%33 = fadd double %21, %23
%34 = fsub double %33, %21
%35 = add nuw nsw i32 %19, 1
%36 = zext i32 %35 to i64
%37 = getelementptr inbounds double, double* %1, i64 %36
%38 = load double, double* %37
br label %46
39: ; preds = %27, %18
%40 = fadd double %22, %23
%41 = fsub double %40, %22
%42 = add nuw nsw i32 %20, 1
%43 = zext i32 %42 to i64
%44 = getelementptr inbounds double, double* %3, i64 %43
%45 = load double, double* %44
br label %46
46: ; preds = %39, %32
%47 = phi double [ %33, %32 ], [ %40, %39 ]
%48 = phi double [ %34, %32 ], [ %41, %39 ]
%49 = phi i32 [ %35, %32 ], [ %19, %39 ]
%50 = phi i32 [ %20, %32 ], [ %42, %39 ]
%51 = phi double [ %38, %32 ], [ %21, %39 ]
%52 = phi double [ %22, %32 ], [ %45, %39 ]
%53 = fsub double %23, %48
%54 = add i32 %0, -2
%55 = add i32 %54, %2
%56 = icmp sgt i32 %55, 0
br i1 %56, label %57, label %107
57: ; preds = %46
%58 = zext i32 %55 to i64
br label %59
59: ; preds = %57, %93
%60 = phi i64 [ 0, %57 ], [ %105, %93 ]
%61 = phi double [ %52, %57 ], [ %98, %93 ]
%62 = phi double [ %51, %57 ], [ %97, %93 ]
%63 = phi i32 [ %50, %57 ], [ %96, %93 ]
%64 = phi i32 [ %49, %57 ], [ %95, %93 ]
%65 = phi double [ %47, %57 ], [ %99, %93 ]
%66 = phi double [ %53, %57 ], [ %104, %93 ]
%67 = icmp slt i32 %64, %0
br i1 %67, label %68, label %84
68: ; preds = %59
%69 = icmp slt i32 %63, %2
br i1 %69, label %70, label %75
70: ; preds = %68
%71 = fcmp ule double %61, %62
%72 = fneg double %62
%73 = fcmp ogt double %61, %72
%74 = xor i1 %71, %73
br i1 %74, label %75, label %84
75: ; preds = %70, %68
%76 = fadd double %66, %62
%77 = fsub double %76, %62
%78 = fsub double %66, %77
%79 = getelementptr inbounds double, double* %4, i64 %60
store double %78, double* %79
%80 = add nsw i32 %64, 1
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds double, double* %1, i64 %81
%83 = load double, double* %82
br label %93
84: ; preds = %70, %59
%85 = fadd double %66, %61
%86 = fsub double %85, %61
%87 = fsub double %66, %86
%88 = getelementptr inbounds double, double* %4, i64 %60
store double %87, double* %88
%89 = add nsw i32 %63, 1
%90 = sext i32 %89 to i64
%91 = getelementptr inbounds double, double* %3, i64 %90
%92 = load double, double* %91
br label %93
93: ; preds = %84, %75
%94 = phi double [ %76, %75 ], [ %85, %84 ]
%95 = phi i32 [ %80, %75 ], [ %64, %84 ]
%96 = phi i32 [ %63, %75 ], [ %89, %84 ]
%97 = phi double [ %83, %75 ], [ %62, %84 ]
%98 = phi double [ %61, %75 ], [ %92, %84 ]
%99 = fadd double %65, %94
%100 = fsub double %99, %65
%101 = fsub double %99, %100
%102 = fsub double %94, %100
%103 = fsub double %65, %101
%104 = fadd double %102, %103
%105 = add nuw nsw i64 %60, 1
%106 = icmp eq i64 %105, %58
br i1 %106, label %107, label %59
107: ; preds = %93, %46
%108 = phi double [ %53, %46 ], [ %104, %93 ]
%109 = phi double [ %47, %46 ], [ %99, %93 ]
%110 = phi i32 [ 0, %46 ], [ %55, %93 ]
%111 = zext i32 %110 to i64
%112 = getelementptr inbounds double, double* %4, i64 %111
store double %108, double* %112
%113 = add nuw nsw i32 %110, 1
%114 = zext i32 %113 to i64
%115 = getelementptr inbounds double, double* %4, i64 %114
store double %109, double* %115
%116 = add nuw nsw i32 %110, 2
ret i32 %116
}