void uphi(x, n, hpf, m, y) double x[], hpf[], y[];
int n, m;
{
  int meven, modd, i, h, j, mmin;
  double s;
  meven = (m + 1) / 2;
  modd = m / 2;
  for (i = 0; i + meven < n; i++) {
    s = 0.;
    for (h = 0; h < meven; h++)
      s += hpf[2 * h] * x[i + h];
    y[2 * i + 1] = s;
    s = 0.;
    for (h = 0; h < modd; h++)
      s += hpf[2 * h + 1] * x[i + h];
    y[2 * i] = s;
  }
  mmin = n - meven;
  if (mmin < 0)
    mmin = 0;
  for (i = mmin; i < n; i++) {
    s = 0.;
    j = i;
    for (h = 0; h < meven; h++) {
      if (j >= n)
        j -= n;
      s += hpf[2 * h] * x[j];
      j++;
    }
    y[2 * i + 1] = s;
    s = 0.;
    j = i;
    for (h = 0; h < modd; h++) {
      if (j >= n)
        j -= n;
      s += hpf[2 * h + 1] * x[j];
      j++;
    }
    y[2 * i] = s;
  }
}


x86 assembly (clang, -O0):
<uphi>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8,-0x28(%rbp)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x2c(%rbp)
mov    -0x1c(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x30(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
add    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    124b <uphi+0x14b>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    11aa <uphi+0xaa>
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    115a <uphi+0x5a>
movsd  -0x48(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x34(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x30(%rbp),%eax
jge    1226 <uphi+0x126>
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    11d3 <uphi+0xd3>
movsd  -0x48(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x34(%rbp),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    113c <uphi+0x3c>
mov    -0xc(%rbp),%eax
sub    -0x2c(%rbp),%eax
mov    %eax,-0x40(%rbp)
cmpl   $0x0,-0x40(%rbp)
jge    1265 <uphi+0x165>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    13b9 <uphi+0x2b9>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x34(%rbp),%eax
mov    %eax,-0x3c(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    12f7 <uphi+0x1f7>
mov    -0x3c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jl     12af <uphi+0x1af>
mov    -0xc(%rbp),%ecx
mov    -0x3c(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    128c <uphi+0x18c>
movsd  -0x48(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x34(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x34(%rbp),%eax
mov    %eax,-0x3c(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x30(%rbp),%eax
jge    1394 <uphi+0x294>
mov    -0x3c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jl     1349 <uphi+0x249>
mov    -0xc(%rbp),%ecx
mov    -0x3c(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1326 <uphi+0x226>
movsd  -0x48(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x34(%rbp),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    126b <uphi+0x16b>
pop    %rbp
ret

x86 assembly (clang, -O2):
<uphi>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
lea    0x1(%rcx),%eax
shr    $0x1f,%eax
lea    (%rcx,%rax,1),%ebp
add    $0x1,%ebp
sar    %ebp
mov    %ecx,%eax
shr    $0x1f,%eax
add    %ecx,%eax
sar    %eax
mov    %eax,%r15d
cmp    %esi,%ebp
mov    %rdx,-0x30(%rsp)
mov    %r15,-0x10(%rsp)
jge    12ac <uphi+0x1ac>
cmp    $0x2,%ebp
mov    $0x1,%r9d
cmovge %ebp,%r9d
mov    %esi,%r12d
mov    %ebp,-0x28(%rsp)
sub    %ebp,%r12d
mov    %r9d,%r14d
and    $0x7ffffffe,%r14d
mov    %ecx,%r10d
and    $0xfffffffe,%r10d
and    $0xfffffffe,%r15d
lea    0x10(%rdx),%rax
mov    %rax,-0x18(%rsp)
lea    0x8(%rdi),%rbx
lea    0x18(%rdx),%rax
mov    %rax,-0x20(%rsp)
xor    %r13d,%r13d
jmp    1197 <uphi+0x97>
nopl   (%rax)
movsd  %xmm0,(%r8,%r11,8)
add    $0x1,%r13
add    $0x8,%rbx
cmp    %r12,%r13
je     12b4 <uphi+0x1b4>
xorpd  %xmm0,%xmm0
xorpd  %xmm1,%xmm1
test   %ecx,%ecx
jle    1212 <uphi+0x112>
cmp    $0x3,%ecx
jge    11b0 <uphi+0xb0>
xorpd  %xmm1,%xmm1
xor    %eax,%eax
jmp    11ee <uphi+0xee>
xorpd  %xmm2,%xmm2
mov    -0x18(%rsp),%rbp
xor    %eax,%eax
nopl   0x0(%rax,%rax,1)
movsd  -0x10(%rbp),%xmm3
mulsd  -0x8(%rbx,%rax,8),%xmm3
movsd  0x0(%rbp),%xmm1
addsd  %xmm2,%xmm3
mulsd  (%rbx,%rax,8),%xmm1
addsd  %xmm3,%xmm1
add    $0x2,%rax
add    $0x20,%rbp
movapd %xmm1,%xmm2
cmp    %rax,%r14
jne    11c0 <uphi+0xc0>
test   $0x1,%r9b
je     1212 <uphi+0x112>
mov    %rax,%rdx
shl    $0x4,%rdx
mov    -0x30(%rsp),%rbp
movsd  0x0(%rbp,%rdx,1),%xmm2
add    %r13,%rax
mulsd  (%rdi,%rax,8),%xmm2
addsd  %xmm2,%xmm1
mov    %r13,%r11
add    %r13,%r11
mov    %r13,%rax
shl    $0x4,%rax
movsd  %xmm1,0x8(%r8,%rax,1)
cmp    $0x2,%ecx
jl     1180 <uphi+0x80>
cmp    $0x2,%r10d
jne    1240 <uphi+0x140>
xorpd  %xmm0,%xmm0
xor    %edx,%edx
jmp    127e <uphi+0x17e>
nopl   (%rax)
xorpd  %xmm1,%xmm1
mov    -0x20(%rsp),%rbp
xor    %edx,%edx
nopl   0x0(%rax,%rax,1)
movsd  -0x10(%rbp),%xmm2
mulsd  -0x8(%rbx,%rdx,8),%xmm2
movsd  0x0(%rbp),%xmm0
addsd  %xmm1,%xmm2
mulsd  (%rbx,%rdx,8),%xmm0
addsd  %xmm2,%xmm0
add    $0x2,%rdx
add    $0x20,%rbp
movapd %xmm0,%xmm1
cmp    %rdx,%r15
jne    1250 <uphi+0x150>
testb  $0x1,-0x10(%rsp)
je     1180 <uphi+0x80>
mov    %rdx,%rax
shl    $0x4,%rax
mov    -0x30(%rsp),%rbp
movsd  0x8(%rbp,%rax,1),%xmm1
add    %r13,%rdx
mulsd  (%rdi,%rdx,8),%xmm1
addsd  %xmm1,%xmm0
jmp    1180 <uphi+0x80>
mov    %esi,%r12d
sub    %ebp,%r12d
jmp    12b8 <uphi+0x1b8>
mov    -0x28(%rsp),%ebp
xor    %r15d,%r15d
test   %r12d,%r12d
cmovs  %r15d,%r12d
cmp    %esi,%r12d
jge    14ba <uphi+0x3ba>
cmp    $0x2,%ebp
mov    $0x1,%r10d
cmovge %ebp,%r10d
mov    %r12d,%r11d
mov    %esi,%r14d
mov    %r10,-0x18(%rsp)
and    $0x7ffffffe,%r10d
mov    %ecx,%eax
and    $0xfffffffe,%eax
mov    %eax,-0x20(%rsp)
mov    -0x10(%rsp),%r9
and    $0xfffffffe,%r9d
mov    -0x30(%rsp),%rax
lea    0x10(%rax),%rdx
mov    %rdx,-0x28(%rsp)
add    $0x18,%rax
mov    %rax,-0x8(%rsp)
jmp    1333 <uphi+0x233>
cs nopw 0x0(%rax,%rax,1)
nop
movsd  %xmm0,(%r8,%r13,8)
add    $0x1,%r11
cmp    %r14,%r11
je     14ba <uphi+0x3ba>
xorpd  %xmm0,%xmm0
xorpd  %xmm1,%xmm1
test   %ecx,%ecx
jle    13ee <uphi+0x2ee>
cmp    $0x3,%ecx
jge    1360 <uphi+0x260>
xorpd  %xmm1,%xmm1
xor    %r13d,%r13d
mov    %r11d,%ebx
jmp    13c2 <uphi+0x2c2>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
xorpd  %xmm2,%xmm2
mov    -0x28(%rsp),%rax
xor    %r13d,%r13d
mov    %r11d,%ebx
nop
cmp    %esi,%ebx
mov    %esi,%edx
cmovl  %r15d,%edx
mov    %ebx,%ebp
sub    %edx,%ebp
movslq %ebp,%rbp
movsd  -0x10(%rax),%xmm3
mulsd  (%rdi,%rbp,8),%xmm3
movsd  (%rax),%xmm1
addsd  %xmm2,%xmm3
neg    %edx
add    %ebx,%edx
add    $0x1,%edx
cmp    %esi,%edx
mov    %esi,%ebx
cmovl  %r15d,%ebx
sub    %ebx,%edx
movslq %edx,%rbx
mulsd  (%rdi,%rbx,8),%xmm1
addsd  %xmm3,%xmm1
add    $0x1,%ebx
add    $0x2,%r13
add    $0x20,%rax
movapd %xmm1,%xmm2
cmp    %r13,%r10
jne    1370 <uphi+0x270>
testb  $0x1,-0x18(%rsp)
je     13ee <uphi+0x2ee>
cmp    %esi,%ebx
mov    %esi,%eax
cmovl  %r15d,%eax
sub    %eax,%ebx
shl    $0x4,%r13
mov    -0x30(%rsp),%rax
movsd  (%rax,%r13,1),%xmm2
movslq %ebx,%rax
mulsd  (%rdi,%rax,8),%xmm2
addsd  %xmm2,%xmm1
lea    (%r11,%r11,1),%r13
mov    %r11,%rax
shl    $0x4,%rax
movsd  %xmm1,0x8(%r8,%rax,1)
cmp    $0x2,%ecx
jl     1320 <uphi+0x220>
cmpl   $0x2,-0x20(%rsp)
jne    1420 <uphi+0x320>
xorpd  %xmm0,%xmm0
xor    %ebp,%ebp
mov    %r11d,%ebx
jmp    1485 <uphi+0x385>
nopl   0x0(%rax,%rax,1)
xorpd  %xmm1,%xmm1
mov    -0x8(%rsp),%rax
xor    %ebp,%ebp
mov    %r11d,%ebx
xchg   %ax,%ax
cmp    %esi,%ebx
mov    %esi,%edx
cmovl  %r15d,%edx
mov    %ebx,%r12d
sub    %edx,%r12d
movslq %r12d,%r12
movsd  -0x10(%rax),%xmm2
mulsd  (%rdi,%r12,8),%xmm2
movsd  (%rax),%xmm0
addsd  %xmm1,%xmm2
neg    %edx
add    %ebx,%edx
add    $0x1,%edx
cmp    %esi,%edx
mov    %esi,%ebx
cmovl  %r15d,%ebx
sub    %ebx,%edx
movslq %edx,%rbx
mulsd  (%rdi,%rbx,8),%xmm0
addsd  %xmm2,%xmm0
add    $0x1,%ebx
add    $0x2,%rbp
add    $0x20,%rax
movapd %xmm0,%xmm1
cmp    %rbp,%r9
jne    1430 <uphi+0x330>
testb  $0x1,-0x10(%rsp)
je     1320 <uphi+0x220>
cmp    %esi,%ebx
mov    %esi,%eax
cmovl  %r15d,%eax
sub    %eax,%ebx
shl    $0x4,%rbp
mov    -0x30(%rsp),%rax
movsd  0x8(%rax,%rbp,1),%xmm1
movslq %ebx,%rax
mulsd  (%rdi,%rax,8),%xmm1
addsd  %xmm1,%xmm0
jmp    1320 <uphi+0x220>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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