void ExpandStandard_1D(double In[], long NxIn, double Out[], double h[],
                       long nh) {
  long k, j, i, i1, i2;
  long kn, nexp, n;
  nexp = NxIn * 2L;
  n = NxIn;
  kn = n - 1;
  if (nh < 2L) {
    for (i = 0L; i < NxIn; i++) {
      j = i * 2L;
      Out[j] = In[i];
      Out[j + 1] = In[i];
    }
  } else {
    for (i = 0L; i < nexp; i++) {
      Out[i] = 0.0;
      for (k = (i % 2L); k < nh; k += 2L) {
        i1 = (i - k) / 2L;
        if (i1 < 0L) {
          i1 = (-i1) % kn;
          if (i1 > kn)
            i1 = kn - i1;
        }
        Out[i] = Out[i] + h[k] * In[i1];
      }
      for (k = 2L - (i % 2L); k < nh; k += 2L) {
        i2 = (i + k) / 2L;
        if (i2 > kn) {
          i2 = i2 % kn;
          i2 = kn - i2;
          if (i2 > kn)
            i2 = kn - i2;
        }
        Out[i] = Out[i] + h[k] * In[i2];
      }
    }
  }
}


x86 assembly (clang, -O0):
<ExpandStandard_1D>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    -0x10(%rbp),%rax
shl    $0x1,%rax
mov    %rax,-0x60(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x68(%rbp)
mov    -0x68(%rbp),%rax
sub    $0x1,%rax
mov    %rax,-0x58(%rbp)
cmpq   $0x2,-0x28(%rbp)
jge    11b0 <ExpandStandard_1D+0xb0>
movq   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    11ab <ExpandStandard_1D+0xab>
mov    -0x40(%rbp),%rax
shl    $0x1,%rax
mov    %rax,-0x38(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  %xmm0,0x8(%rax,%rcx,8)
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    114b <ExpandStandard_1D+0x4b>
jmp    1380 <ExpandStandard_1D+0x280>
movq   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%rax
cmp    -0x60(%rbp),%rax
jge    137b <ExpandStandard_1D+0x27b>
mov    -0x18(%rbp),%rax
mov    -0x40(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%rax
mov    $0x2,%ecx
cqto
idiv   %rcx
mov    %rdx,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x28(%rbp),%rax
jge    1293 <ExpandStandard_1D+0x193>
mov    -0x40(%rbp),%rax
sub    -0x30(%rbp),%rax
mov    $0x2,%ecx
cqto
idiv   %rcx
mov    %rax,-0x48(%rbp)
cmpq   $0x0,-0x48(%rbp)
jge    1246 <ExpandStandard_1D+0x146>
xor    %eax,%eax
sub    -0x48(%rbp),%rax
cqto
idivq  -0x58(%rbp)
mov    %rdx,-0x48(%rbp)
mov    -0x48(%rbp),%rax
cmp    -0x58(%rbp),%rax
jle    1241 <ExpandStandard_1D+0x141>
mov    -0x58(%rbp),%rax
sub    -0x48(%rbp),%rax
mov    %rax,-0x48(%rbp)
jmp    1246 <ExpandStandard_1D+0x146>
mov    -0x18(%rbp),%rax
mov    -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
mov    -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%rax
add    $0x2,%rax
mov    %rax,-0x30(%rbp)
jmp    11e8 <ExpandStandard_1D+0xe8>
mov    -0x40(%rbp),%rax
mov    $0x2,%ecx
cqto
idiv   %rcx
mov    $0x2,%eax
sub    %rdx,%rax
mov    %rax,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x28(%rbp),%rax
jge    1365 <ExpandStandard_1D+0x265>
mov    -0x40(%rbp),%rax
add    -0x30(%rbp),%rax
mov    $0x2,%ecx
cqto
idiv   %rcx
mov    %rax,-0x50(%rbp)
mov    -0x50(%rbp),%rax
cmp    -0x58(%rbp),%rax
jle    1318 <ExpandStandard_1D+0x218>
mov    -0x50(%rbp),%rax
cqto
idivq  -0x58(%rbp)
mov    %rdx,-0x50(%rbp)
mov    -0x58(%rbp),%rax
sub    -0x50(%rbp),%rax
mov    %rax,-0x50(%rbp)
mov    -0x50(%rbp),%rax
cmp    -0x58(%rbp),%rax
jle    1313 <ExpandStandard_1D+0x213>
mov    -0x58(%rbp),%rax
sub    -0x50(%rbp),%rax
mov    %rax,-0x50(%rbp)
jmp    1318 <ExpandStandard_1D+0x218>
mov    -0x18(%rbp),%rax
mov    -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
mov    -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%rax
add    $0x2,%rax
mov    %rax,-0x30(%rbp)
jmp    12ad <ExpandStandard_1D+0x1ad>
jmp    136a <ExpandStandard_1D+0x26a>
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    11b8 <ExpandStandard_1D+0xb8>
jmp    1380 <ExpandStandard_1D+0x280>
pop    %rbp
ret

x86 assembly (clang, -O2):
<ExpandStandard_1D>:
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rdx,%r9
lea    (%rsi,%rsi,1),%rax
cmp    $0x2,%r8
jge    11cb <ExpandStandard_1D+0xcb>
test   %rsi,%rsi
jle    13cf <ExpandStandard_1D+0x2cf>
cmp    $0x1,%rsi
je     1143 <ExpandStandard_1D+0x43>
lea    (%rdi,%rsi,8),%rcx
cmp    %r9,%rcx
jbe    1328 <ExpandStandard_1D+0x228>
lea    (%r9,%rax,8),%rax
cmp    %rdi,%rax
jbe    1328 <ExpandStandard_1D+0x228>
xor    %eax,%eax
mov    %rax,%rcx
not    %rcx
test   $0x1,%sil
je     1173 <ExpandStandard_1D+0x73>
movsd  (%rdi,%rax,8),%xmm0
mov    %rax,%rdx
shl    $0x4,%rdx
movsd  %xmm0,(%r9,%rdx,1)
movsd  (%rdi,%rax,8),%xmm0
movsd  %xmm0,0x8(%r9,%rdx,1)
or     $0x1,%rax
add    %rsi,%rcx
je     13cf <ExpandStandard_1D+0x2cf>
mov    %rax,%rcx
shl    $0x4,%rcx
add    %r9,%rcx
add    $0x18,%rcx
nopw   0x0(%rax,%rax,1)
movsd  (%rdi,%rax,8),%xmm0
movsd  %xmm0,-0x18(%rcx)
movsd  (%rdi,%rax,8),%xmm0
movsd  %xmm0,-0x10(%rcx)
movsd  0x8(%rdi,%rax,8),%xmm0
movsd  %xmm0,-0x8(%rcx)
movsd  0x8(%rdi,%rax,8),%xmm0
movsd  %xmm0,(%rcx)
add    $0x2,%rax
add    $0x20,%rcx
cmp    %rax,%rsi
jne    1190 <ExpandStandard_1D+0x90>
jmp    13cf <ExpandStandard_1D+0x2cf>
test   %rsi,%rsi
jle    13cf <ExpandStandard_1D+0x2cf>
lea    -0x1(%rsi),%r11
cmp    $0x2,%rax
mov    $0x1,%r14d
cmovge %rax,%r14
mov    $0x2,%r15d
xor    %r13d,%r13d
jmp    1211 <ExpandStandard_1D+0x111>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
add    $0x1,%r13
add    $0x1,%r15
cmp    %r14,%r13
je     13cf <ExpandStandard_1D+0x2cf>
mov    %r13d,%r12d
and    $0x1,%r12d
movq   $0x0,(%r9,%r13,8)
xorpd  %xmm0,%xmm0
cmp    %r8,%r12
jge    12b0 <ExpandStandard_1D+0x1b0>
mov    %r13,%rbx
and    $0xfffffffffffffffe,%rbx
xorpd  %xmm0,%xmm0
mov    %r12,%r10
jmp    127a <ExpandStandard_1D+0x17a>
nopl   (%rax)
cqto
idiv   %r11
mov    %rdx,%rax
mov    %rax,%rdx
neg    %rdx
add    %r11,%rax
cmp    %rsi,%rdx
cmovl  %rdx,%rax
movsd  (%rcx,%r10,8),%xmm1
mulsd  (%rdi,%rax,8),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,(%r9,%r13,8)
add    $0x2,%r10
add    $0xfffffffffffffffe,%rbx
cmp    %r8,%r10
jge    12b0 <ExpandStandard_1D+0x1b0>
mov    %rbx,%rax
shr    $0x3f,%rax
add    %rbx,%rax
sar    %rax
cmp    $0xfffffffffffffffe,%rbx
jg     1258 <ExpandStandard_1D+0x158>
mov    %rax,%rdx
or     %r11,%rdx
shr    $0x20,%rdx
jne    1240 <ExpandStandard_1D+0x140>
xor    %edx,%edx
div    %r11d
mov    %edx,%eax
jmp    1248 <ExpandStandard_1D+0x148>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    $0x2,%ebx
sub    %r12,%rbx
cmp    %r8,%rbx
jge    1200 <ExpandStandard_1D+0x100>
mov    %r15,%r10
sub    %r12,%r10
shr    %r10
jmp    130a <ExpandStandard_1D+0x20a>
nopl   0x0(%rax)
mov    %r10,%rax
cqto
idiv   %r11
mov    %r11,%rax
sub    %rdx,%rax
cmp    %rsi,%rax
cmovge %rdx,%rax
movsd  (%rcx,%rbx,8),%xmm1
mulsd  (%rdi,%rax,8),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,(%r9,%r13,8)
add    $0x2,%rbx
add    $0x1,%r10
cmp    %r8,%rbx
jge    1200 <ExpandStandard_1D+0x100>
mov    %r10,%rax
cmp    %rsi,%r10
jl     12e5 <ExpandStandard_1D+0x1e5>
mov    %r10,%rax
or     %r11,%rax
shr    $0x20,%rax
jne    12d0 <ExpandStandard_1D+0x1d0>
mov    %r10d,%eax
xor    %edx,%edx
div    %r11d
jmp    12d8 <ExpandStandard_1D+0x1d8>
mov    %rsi,%rax
and    $0xfffffffffffffffe,%rax
lea    -0x2(%rax),%rcx
mov    %rcx,%r8
shr    %r8
add    $0x1,%r8
test   %rcx,%rcx
je     13d9 <ExpandStandard_1D+0x2d9>
mov    %r8,%rdx
and    $0xfffffffffffffffe,%rdx
lea    0x20(%r9),%rbx
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movups (%rdi,%rcx,8),%xmm0
movaps %xmm0,%xmm1
movlhps %xmm0,%xmm1
movhlps %xmm0,%xmm0
movups %xmm0,-0x10(%rbx)
movups %xmm1,-0x20(%rbx)
movupd 0x10(%rdi,%rcx,8),%xmm0
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
unpckhpd %xmm0,%xmm0
movupd %xmm0,0x10(%rbx)
movupd %xmm1,(%rbx)
add    $0x4,%rcx
add    $0x40,%rbx
add    $0xfffffffffffffffe,%rdx
jne    1360 <ExpandStandard_1D+0x260>
test   $0x1,%r8b
je     13c6 <ExpandStandard_1D+0x2c6>
movupd (%rdi,%rcx,8),%xmm0
shl    $0x4,%rcx
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
unpckhpd %xmm0,%xmm0
movupd %xmm0,0x10(%r9,%rcx,1)
movupd %xmm1,(%r9,%rcx,1)
cmp    %rsi,%rax
jne    1145 <ExpandStandard_1D+0x45>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret
xor    %ecx,%ecx
test   $0x1,%r8b
jne    13a4 <ExpandStandard_1D+0x2a4>
jmp    13c6 <ExpandStandard_1D+0x2c6>

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

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