void render_trace(int steps, int timesteps, int nstates, double *x, double *y,
                  int *states, double *times, double *dwells,
                  double *emissions) {
  int i, j;
  double t0, t1;
  double tau = x[1] - x[0];
  int a = 0, b = 0, aflag = 1, bflag = 1;
  for (i = 0; i < steps; i++) {
    t1 = x[i];
    t0 = t1 - tau;
    aflag = 1;
    bflag = 1;
    for (j = a; j < timesteps + 1; j++) {
      if (times[j] > t0 && aflag) {
        a = j - 1;
        aflag = 0;
      }
      if (times[j] > t1 && bflag) {
        b = j;
        bflag = 0;
      }
      if (!aflag && !bflag) {
        y[i] = 0.;
        if (b - a == 1) {
          y[i] += (t1 - t0) / tau * emissions[states[a]];
        } else if (b - a > 1) {
          y[i] += (times[a + 1] - t0) / tau * emissions[states[a]];
          for (j = a + 1; j < b - 1; j++) {
            y[i] += dwells[j] / tau * emissions[states[j]];
          }
          y[i] += (t1 - times[b - 1]) / tau * emissions[states[b - 1]];
        }
        break;
      }
    }
  }
}


x86 assembly (clang, -O0):
<render_trace>:
push   %rbp
mov    %rsp,%rbp
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x18(%rbp),%rax
subsd  (%rax),%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x4c(%rbp)
movl   $0x0,-0x50(%rbp)
movl   $0x1,-0x54(%rbp)
movl   $0x1,-0x58(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    13ca <render_trace+0x2ca>
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
subsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x1,-0x54(%rbp)
movl   $0x1,-0x58(%rbp)
mov    -0x4c(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    -0x8(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    13b7 <render_trace+0x2b7>
mov    0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x38(%rbp),%xmm0
jbe    11e2 <render_trace+0xe2>
cmpl   $0x0,-0x54(%rbp)
je     11e2 <render_trace+0xe2>
mov    -0x30(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x4c(%rbp)
movl   $0x0,-0x54(%rbp)
mov    0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x40(%rbp),%xmm0
jbe    1211 <render_trace+0x111>
cmpl   $0x0,-0x58(%rbp)
je     1211 <render_trace+0x111>
mov    -0x30(%rbp),%eax
mov    %eax,-0x50(%rbp)
movl   $0x0,-0x58(%rbp)
cmpl   $0x0,-0x54(%rbp)
jne    13a4 <render_trace+0x2a4>
cmpl   $0x0,-0x58(%rbp)
jne    13a4 <render_trace+0x2a4>
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x50(%rbp),%eax
sub    -0x4c(%rbp),%eax
cmp    $0x1,%eax
jne    1287 <render_trace+0x187>
movsd  -0x40(%rbp),%xmm0
subsd  -0x38(%rbp),%xmm0
divsd  -0x48(%rbp),%xmm0
mov    0x20(%rbp),%rax
mov    -0x28(%rbp),%rcx
movslq -0x4c(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    139f <render_trace+0x29f>
mov    -0x50(%rbp),%eax
sub    -0x4c(%rbp),%eax
cmp    $0x1,%eax
jle    139a <render_trace+0x29a>
mov    0x10(%rbp),%rax
mov    -0x4c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x38(%rbp),%xmm0
divsd  -0x48(%rbp),%xmm0
mov    0x20(%rbp),%rax
mov    -0x28(%rbp),%rcx
movslq -0x4c(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    -0x50(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    134a <render_trace+0x24a>
mov    0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  -0x48(%rbp),%xmm0
mov    0x20(%rbp),%rax
mov    -0x28(%rbp),%rcx
movslq -0x30(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    12ea <render_trace+0x1ea>
movsd  -0x40(%rbp),%xmm0
mov    0x10(%rbp),%rax
mov    -0x50(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
divsd  -0x48(%rbp),%xmm0
mov    0x20(%rbp),%rax
mov    -0x28(%rbp),%rcx
mov    -0x50(%rbp),%edx
sub    $0x1,%edx
movslq %edx,%rdx
movslq (%rcx,%rdx,4),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    139f <render_trace+0x29f>
jmp    13b7 <render_trace+0x2b7>
jmp    13a9 <render_trace+0x2a9>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    119f <render_trace+0x9f>
jmp    13bc <render_trace+0x2bc>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    115e <render_trace+0x5e>
pop    %rbp
ret

x86 assembly (clang, -O2):
<render_trace>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rcx,-0x10(%rsp)
test   %edi,%edi
jle    1346 <render_trace+0x246>
mov    -0x10(%rsp),%rax
movsd  (%rax),%xmm2
movsd  0x8(%rax),%xmm0
subsd  %xmm2,%xmm0
mov    %edi,%eax
mov    %rax,-0x8(%rsp)
xor    %r13d,%r13d
xorpd  %xmm1,%xmm1
xor    %ecx,%ecx
xor    %eax,%eax
cmp    %esi,%eax
jle    1150 <render_trace+0x50>
nop
mov    %eax,%r12d
jmp    1320 <render_trace+0x220>
nopl   0x0(%rax,%rax,1)
movapd %xmm2,%xmm3
subsd  %xmm0,%xmm3
movslq %eax,%rdx
lea    -0x1(%rax),%edi
mov    0x38(%rsp),%rbp
lea    0x0(%rbp,%rdx,8),%rbp
mov    $0x1,%ebx
mov    $0x1,%r10d
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  0x0(%rbp),%xmm4
ucomisd %xmm3,%xmm4
seta   %r14b
test   %r10d,%r10d
setne  %r15b
mov    %edi,%r12d
mov    $0x0,%r11d
test   %r15b,%r14b
jne    11a8 <render_trace+0xa8>
mov    %eax,%r12d
mov    %r10d,%r11d
ucomisd %xmm2,%xmm4
seta   %al
test   %ebx,%ebx
setne  %dl
test   %dl,%al
je     11bd <render_trace+0xbd>
lea    0x1(%rdi),%ecx
xor    %ebx,%ebx
mov    %r11d,%eax
or     %ebx,%eax
je     11e0 <render_trace+0xe0>
add    $0x1,%edi
add    $0x8,%rbp
mov    %r11d,%r10d
mov    %r12d,%eax
cmp    %edi,%esi
jne    1180 <render_trace+0x80>
jmp    1320 <render_trace+0x220>
nopw   0x0(%rax,%rax,1)
movq   $0x0,(%r8,%r13,8)
mov    %ecx,%eax
sub    %r12d,%eax
cmp    $0x1,%eax
jne    1220 <render_trace+0x120>
subsd  %xmm3,%xmm2
divsd  %xmm0,%xmm2
movslq %r12d,%rax
movslq (%r9,%rax,4),%rax
mov    0x48(%rsp),%rdx
mulsd  (%rdx,%rax,8),%xmm2
addsd  %xmm1,%xmm2
jmp    1311 <render_trace+0x211>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
jle    1320 <render_trace+0x220>
movslq %r12d,%rdi
lea    0x1(%r12),%edx
mov    0x38(%rsp),%rbp
movsd  0x8(%rbp,%rdi,8),%xmm4
subsd  %xmm3,%xmm4
divsd  %xmm0,%xmm4
movslq (%r9,%rdi,4),%rbp
mov    0x48(%rsp),%rbx
mulsd  (%rbx,%rbp,8),%xmm4
addsd  %xmm1,%xmm4
movsd  %xmm4,(%r8,%r13,8)
lea    -0x1(%rcx),%r10d
cmp    %r10d,%edx
mov    0x40(%rsp),%rbp
jge    1279 <render_trace+0x179>
add    $0x1,%rdi
lea    -0x3(%rcx),%edx
test   $0x1,%al
jne    127f <render_trace+0x17f>
cmp    %r12d,%edx
jne    12a9 <render_trace+0x1a9>
jmp    12f3 <render_trace+0x1f3>
movapd %xmm4,%xmm3
jmp    12f3 <render_trace+0x1f3>
movsd  0x0(%rbp,%rdi,8),%xmm3
divsd  %xmm0,%xmm3
movslq (%r9,%rdi,4),%rax
mulsd  (%rbx,%rax,8),%xmm3
addsd  %xmm4,%xmm3
movsd  %xmm3,(%r8,%r13,8)
add    $0x1,%rdi
movapd %xmm3,%xmm4
cmp    %r12d,%edx
je     12f3 <render_trace+0x1f3>
mov    %r10d,%eax
movapd %xmm4,%xmm3
movsd  0x0(%rbp,%rdi,8),%xmm4
divsd  %xmm0,%xmm4
movslq (%r9,%rdi,4),%rdx
mulsd  (%rbx,%rdx,8),%xmm4
addsd  %xmm3,%xmm4
movsd  %xmm4,(%r8,%r13,8)
movsd  0x8(%rbp,%rdi,8),%xmm3
divsd  %xmm0,%xmm3
movslq 0x4(%r9,%rdi,4),%rdx
mulsd  (%rbx,%rdx,8),%xmm3
addsd  %xmm4,%xmm3
movsd  %xmm3,(%r8,%r13,8)
add    $0x2,%rdi
cmp    %edi,%eax
jne    12b0 <render_trace+0x1b0>
movslq %r10d,%rax
mov    0x38(%rsp),%rdx
subsd  (%rdx,%rax,8),%xmm2
divsd  %xmm0,%xmm2
movslq (%r9,%rax,4),%rax
mulsd  (%rbx,%rax,8),%xmm2
addsd  %xmm3,%xmm2
movsd  %xmm2,(%r8,%r13,8)
nopw   0x0(%rax,%rax,1)
add    $0x1,%r13
cmp    -0x8(%rsp),%r13
je     1346 <render_trace+0x246>
mov    -0x10(%rsp),%rax
movsd  (%rax,%r13,8),%xmm2
mov    %r12d,%eax
cmp    %esi,%eax
jle    1150 <render_trace+0x50>
jmp    1140 <render_trace+0x40>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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