void fillin(double *t_sim, double *t_obs, int length, int nodata) {
  int t, t2;
  int flag = 0;
  for (t = 1; t < length - 1; t++) {
    if ((t_obs[t] - t_obs[t + 1]) > 3500 && (t_obs[t - 1] - t_obs[t]) < -3500) {
      t_obs[t] = 0.5 * (t_obs[t - 1] + t_obs[t + 1]);
    }
  }
  if (t_obs[0] == nodata) {
    for (t = 1; t < length; t++) {
      if (flag == 0 && t_obs[t] != nodata) {
        t_sim[0] = t_obs[t];
        flag = 1;
      }
    }
  }
  flag = 0;
  if (t_obs[length - 1] == nodata) {
    for (t = length - 2; t > 0; t--) {
      if (flag == 0 && t_obs[t] != nodata) {
        t_sim[length - 1] = t_obs[t];
        flag = 1;
      }
    }
  }
  flag = 0;
  double previous, next;
  for (t = 1; t < length - 1; t++) {
    if (t_obs[t] == nodata) {
      for (t2 = t - 1; t2 > 0; t2--) {
        if (flag == 0 && t_obs[t2] != nodata) {
          previous = t_obs[t2];
          flag = 1;
        }
      }
      flag = 0;
      for (t2 = t + 1; t2 < length; t2++) {
        if (flag == 0 && t_obs[t2] != nodata) {
          next = t_obs[t2];
          flag = 1;
        }
      }
      flag = 0;
      t_sim[t] = 0.5 * (previous + next);
    }
  }
}


x86 assembly (clang, -O0):
<fillin>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
movl   $0x0,-0x24(%rbp)
movl   $0x1,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    11e3 <fillin+0xe3>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  0xeb0(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    11d0 <fillin+0xd0>
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
movsd  0xe87(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    11d0 <fillin+0xd0>
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm1
movsd  0xe41(%rip),%xmm0
mulsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    11d5 <fillin+0xd5>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1120 <fillin+0x20>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
cvtsi2sdl -0x18(%rbp),%xmm1
ucomisd %xmm1,%xmm0
jne    1278 <fillin+0x178>
jp     1278 <fillin+0x178>
movl   $0x1,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1273 <fillin+0x173>
cmpl   $0x0,-0x24(%rbp)
jne    1260 <fillin+0x160>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
cvtsi2sdl -0x18(%rbp),%xmm1
ucomisd %xmm1,%xmm0
jne    1244 <fillin+0x144>
jp     1244 <fillin+0x144>
jmp    1260 <fillin+0x160>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x1,-0x24(%rbp)
jmp    1265 <fillin+0x165>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1207 <fillin+0x107>
jmp    1278 <fillin+0x178>
movl   $0x0,-0x24(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
cvtsi2sdl -0x18(%rbp),%xmm1
ucomisd %xmm1,%xmm0
jne    1328 <fillin+0x228>
jp     1328 <fillin+0x228>
mov    -0x14(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x1c(%rbp)
cmpl   $0x0,-0x1c(%rbp)
jle    1323 <fillin+0x223>
cmpl   $0x0,-0x24(%rbp)
jne    1310 <fillin+0x210>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
cvtsi2sdl -0x18(%rbp),%xmm1
ucomisd %xmm1,%xmm0
jne    12ea <fillin+0x1ea>
jp     12ea <fillin+0x1ea>
jmp    1310 <fillin+0x210>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x1,-0x24(%rbp)
jmp    1315 <fillin+0x215>
mov    -0x1c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x1c(%rbp)
jmp    12af <fillin+0x1af>
jmp    1328 <fillin+0x228>
movl   $0x0,-0x24(%rbp)
movl   $0x1,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    148f <fillin+0x38f>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
cvtsi2sdl -0x18(%rbp),%xmm1
ucomisd %xmm1,%xmm0
jne    147c <fillin+0x37c>
jp     147c <fillin+0x37c>
mov    -0x1c(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x20(%rbp)
cmpl   $0x0,-0x20(%rbp)
jle    13d9 <fillin+0x2d9>
cmpl   $0x0,-0x24(%rbp)
jne    13c6 <fillin+0x2c6>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
cvtsi2sdl -0x18(%rbp),%xmm1
ucomisd %xmm1,%xmm0
jne    13ad <fillin+0x2ad>
jp     13ad <fillin+0x2ad>
jmp    13c6 <fillin+0x2c6>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x1,-0x24(%rbp)
jmp    13cb <fillin+0x2cb>
mov    -0x20(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x20(%rbp)
jmp    1372 <fillin+0x272>
movl   $0x0,-0x24(%rbp)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1452 <fillin+0x352>
cmpl   $0x0,-0x24(%rbp)
jne    143f <fillin+0x33f>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
cvtsi2sdl -0x18(%rbp),%xmm1
ucomisd %xmm1,%xmm0
jne    1426 <fillin+0x326>
jp     1426 <fillin+0x326>
jmp    143f <fillin+0x33f>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x1,-0x24(%rbp)
jmp    1444 <fillin+0x344>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    13e9 <fillin+0x2e9>
movl   $0x0,-0x24(%rbp)
movsd  -0x30(%rbp),%xmm1
addsd  -0x38(%rbp),%xmm1
movsd  0xb95(%rip),%xmm0
mulsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1481 <fillin+0x381>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1336 <fillin+0x236>
pop    %rbp
ret

x86 assembly (clang, -O2):
<fillin>:
push   %r14
push   %rbx
lea    -0x1(%rdx),%r8d
cmp    $0x3,%edx
jl     1187 <fillin+0x87>
mov    %r8d,%r9d
movsd  0x8(%rsi),%xmm1
mov    $0x1,%eax
movsd  0xedf(%rip),%xmm0
movsd  0xedf(%rip),%xmm2
movsd  0xedf(%rip),%xmm3
jmp    1145 <fillin+0x45>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
cmp    %rax,%r9
je     1187 <fillin+0x87>
mov    %rax,%rbx
movapd %xmm1,%xmm4
add    $0x1,%rax
movsd  0x8(%rsi,%rbx,8),%xmm1
movapd %xmm4,%xmm5
subsd  %xmm1,%xmm5
ucomisd %xmm0,%xmm5
jbe    1140 <fillin+0x40>
movsd  -0x8(%rsi,%rbx,8),%xmm5
movapd %xmm5,%xmm6
subsd  %xmm4,%xmm6
ucomisd %xmm6,%xmm2
jbe    1140 <fillin+0x40>
addsd  %xmm1,%xmm5
mulsd  %xmm3,%xmm5
movsd  %xmm5,(%rsi,%rbx,8)
jmp    1140 <fillin+0x40>
movsd  (%rsi),%xmm1
xorps  %xmm0,%xmm0
cvtsi2sd %ecx,%xmm0
ucomisd %xmm0,%xmm1
jne    11d1 <fillin+0xd1>
jp     11d1 <fillin+0xd1>
cmp    $0x2,%edx
jl     11d1 <fillin+0xd1>
mov    %edx,%r9d
add    $0xffffffffffffffff,%r9
cmp    $0x2,%edx
jne    133d <fillin+0x23d>
xor    %ebx,%ebx
mov    $0x1,%eax
test   $0x1,%r9b
je     11d1 <fillin+0xd1>
test   %ebx,%ebx
jne    11d1 <fillin+0xd1>
movsd  (%rsi,%rax,8),%xmm1
ucomisd %xmm0,%xmm1
jne    11cd <fillin+0xcd>
jnp    11d1 <fillin+0xd1>
movsd  %xmm1,(%rdi)
movslq %r8d,%r8
movsd  (%rsi,%r8,8),%xmm1
ucomisd %xmm0,%xmm1
jne    122d <fillin+0x12d>
jp     122d <fillin+0x12d>
cmp    $0x3,%edx
jl     1339 <fillin+0x239>
lea    -0x2(%rdx),%ecx
add    $0x1,%rcx
xor    %eax,%eax
jmp    1210 <fillin+0x110>
cs nopw 0x0(%rax,%rax,1)
movsd  %xmm1,(%rdi,%r8,8)
add    $0xffffffffffffffff,%rcx
cmp    $0x1,%rcx
jle    122d <fillin+0x12d>
mov    %eax,%ebx
mov    $0x1,%eax
test   %ebx,%ebx
jne    1206 <fillin+0x106>
movsd  -0x8(%rsi,%rcx,8),%xmm1
ucomisd %xmm0,%xmm1
jne    1200 <fillin+0x100>
jp     1200 <fillin+0x100>
xor    %eax,%eax
jmp    1206 <fillin+0x106>
cmp    $0x3,%edx
jl     1339 <fillin+0x239>
mov    %edx,%r8d
mov    $0x1,%r9d
mov    $0x2,%r11d
movsd  0xdc3(%rip),%xmm1
mov    %r8,%rax
mov    %rsi,%r14
jmp    128a <fillin+0x18a>
cs nopw 0x0(%rax,%rax,1)
nop
movapd %xmm3,%xmm4
addsd  %xmm2,%xmm4
mulsd  %xmm1,%xmm4
movsd  %xmm4,(%rdi,%r9,8)
mov    %r10,%r9
add    $0x1,%r11
add    $0x8,%r14
add    $0xffffffffffffffff,%rax
cmp    %r8,%r11
je     1339 <fillin+0x239>
movsd  (%rsi,%r9,8),%xmm4
ucomisd %xmm0,%xmm4
jne    1330 <fillin+0x230>
jp     1330 <fillin+0x230>
cmp    $0x2,%r9
jb     12e0 <fillin+0x1e0>
lea    0x1(%r9),%rdx
xor    %ebx,%ebx
jmp    12be <fillin+0x1be>
xchg   %ax,%ax
movapd %xmm4,%xmm3
add    $0xffffffffffffffff,%rdx
cmp    $0x2,%rdx
jle    12e0 <fillin+0x1e0>
mov    %ebx,%ecx
mov    $0x1,%ebx
test   %ecx,%ecx
jne    12b4 <fillin+0x1b4>
lea    -0x2(%rdx),%rcx
mov    %ecx,%ecx
movsd  (%rsi,%rcx,8),%xmm4
ucomisd %xmm0,%xmm4
jne    12b0 <fillin+0x1b0>
jp     12b0 <fillin+0x1b0>
xor    %ebx,%ebx
jmp    12b4 <fillin+0x1b4>
lea    0x1(%r9),%r10
cmp    %r8,%r10
jge    1260 <fillin+0x160>
xor    %ecx,%ecx
mov    $0x2,%edx
jmp    1311 <fillin+0x211>
cs nopw 0x0(%rax,%rax,1)
movapd %xmm4,%xmm2
add    $0x1,%rdx
cmp    %rdx,%rax
je     1260 <fillin+0x160>
mov    %ecx,%ebx
mov    $0x1,%ecx
test   %ebx,%ebx
jne    1304 <fillin+0x204>
movsd  (%r14,%rdx,8),%xmm4
ucomisd %xmm0,%xmm4
jne    1300 <fillin+0x200>
jp     1300 <fillin+0x200>
xor    %ecx,%ecx
jmp    1304 <fillin+0x204>
xchg   %ax,%ax
add    $0x1,%r9
jmp    1275 <fillin+0x175>
pop    %rbx
pop    %r14
ret
mov    %r9,%r10
and    $0xfffffffffffffffe,%r10
neg    %r10
xor    %ebx,%ebx
mov    $0x1,%eax
jmp    136a <fillin+0x26a>
movsd  %xmm1,(%rdi)
lea    (%r10,%rax,1),%rcx
add    $0x2,%rcx
add    $0x2,%rax
cmp    $0x1,%rcx
je     11b6 <fillin+0xb6>
mov    %ebx,%ecx
mov    $0x1,%ebx
test   %ecx,%ecx
jne    1354 <fillin+0x254>
movsd  (%rsi,%rax,8),%xmm1
ucomisd %xmm0,%xmm1
jne    1350 <fillin+0x250>
jp     1350 <fillin+0x250>
movsd  0x8(%rsi,%rax,8),%xmm1
ucomisd %xmm0,%xmm1
jne    1350 <fillin+0x250>
jp     1350 <fillin+0x250>
xor    %ebx,%ebx
jmp    1354 <fillin+0x254>

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

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