void twofft(float data1[], float data2[], float fft1[], float fft2[],
            unsigned long n) {
  void four1(float data[], unsigned long nn, int isign);
  unsigned long nn3, nn2, jj, j;
  float rep, rem, aip, aim;
  nn3 = 1 + (nn2 = 2 + n + n);
  for (j = 1, jj = 2; j <= n; j++, jj += 2) {
    fft1[jj - 1] = data1[j];
    fft1[jj] = data2[j];
  }
  four1(fft1, n, 1);
  fft2[1] = fft1[2];
  fft1[2] = fft2[2] = 0.0;
  for (j = 3; j <= n + 1; j += 2) {
    rep = 0.5 * (fft1[j] + fft1[nn2 - j]);
    rem = 0.5 * (fft1[j] - fft1[nn2 - j]);
    aip = 0.5 * (fft1[j + 1] + fft1[nn3 - j]);
    aim = 0.5 * (fft1[j + 1] - fft1[nn3 - j]);
    fft1[j] = rep;
    fft1[j + 1] = aim;
    fft1[nn2 - j] = rep;
    fft1[nn3 - j] = -aim;
    fft2[j] = aip;
    fft2[j + 1] = -rem;
    fft2[nn2 - j] = aip;
    fft2[nn3 - j] = rem;
  }
}


x86 assembly (clang, -O0):
<twofft>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    -0x28(%rbp),%rax
add    $0x2,%rax
add    -0x28(%rbp),%rax
mov    %rax,-0x38(%rbp)
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
movq   $0x1,-0x48(%rbp)
movq   $0x2,-0x40(%rbp)
mov    -0x48(%rbp),%rax
cmp    -0x28(%rbp),%rax
ja     11b7 <twofft+0xa7>
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x40(%rbp),%rcx
sub    $0x1,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
mov    -0x48(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x40(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x48(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x48(%rbp)
mov    -0x40(%rbp),%rax
add    $0x2,%rax
mov    %rax,-0x40(%rbp)
jmp    1154 <twofft+0x44>
mov    -0x18(%rbp),%rdi
mov    -0x28(%rbp),%rsi
mov    $0x1,%edx
call   1030 <four1@plt>
mov    -0x18(%rbp),%rax
movss  0x8(%rax),%xmm0
mov    -0x20(%rbp),%rax
movss  %xmm0,0x4(%rax)
mov    -0x20(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,0x8(%rax)
mov    -0x18(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,0x8(%rax)
movq   $0x3,-0x48(%rbp)
mov    -0x48(%rbp),%rax
mov    -0x28(%rbp),%rcx
add    $0x1,%rcx
cmp    %rcx,%rax
ja     13bb <twofft+0x2ab>
mov    -0x18(%rbp),%rax
mov    -0x48(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%rcx
sub    -0x48(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xdc6(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x4c(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x48(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%rcx
sub    -0x48(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xd8f(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x50(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x48(%rbp),%rcx
movss  0x4(%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%rcx
sub    -0x48(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xd57(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x54(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x48(%rbp),%rcx
movss  0x4(%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%rcx
sub    -0x48(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xd1f(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x58(%rbp)
movss  -0x4c(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x48(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x58(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x48(%rbp),%rcx
movss  %xmm0,0x4(%rax,%rcx,4)
movss  -0x4c(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%rcx
sub    -0x48(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x58(%rbp),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%rcx
sub    -0x48(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x54(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x48(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x50(%rbp),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mov    -0x20(%rbp),%rax
mov    -0x48(%rbp),%rcx
movss  %xmm0,0x4(%rax,%rcx,4)
movss  -0x54(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%rcx
sub    -0x48(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x50(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x30(%rbp),%rcx
sub    -0x48(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x48(%rbp),%rax
add    $0x2,%rax
mov    %rax,-0x48(%rbp)
jmp    11fb <twofft+0xeb>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<twofft>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
mov    %r8,%rbx
mov    %rcx,%r14
mov    %rdx,%r15
test   %r8,%r8
je     1346 <twofft+0x236>
mov    $0x1,%r8d
mov    $0x2,%r10d
cmp    $0x4,%rbx
jae    1156 <twofft+0x46>
test   $0x1,%bl
jne    12c0 <twofft+0x1b0>
mov    %r8,%r9
cmp    %rbx,%r8
jne    12e6 <twofft+0x1d6>
jmp    1346 <twofft+0x236>
lea    -0x1(%rbx),%rax
mov    $0x8,%ecx
mul    %rcx
lea    0x4(%r15),%r9
seto   %cl
lea    (%r15,%rax,1),%rdx
add    $0x4,%rdx
cmp    %r9,%rdx
jb     113c <twofft+0x2c>
test   %cl,%cl
jne    113c <twofft+0x2c>
lea    0x8(%r15),%rdx
add    %r15,%rax
add    $0x8,%rax
cmp    %rdx,%rax
jb     113c <twofft+0x2c>
test   %cl,%cl
jne    113c <twofft+0x2c>
lea    (%r15,%rbx,8),%rcx
add    $0x4,%rcx
lea    0x4(%rdi),%r12
lea    (%rdi,%rbx,4),%rdx
add    $0x4,%rdx
lea    0x4(%rsi),%r11
lea    (%rsi,%rbx,4),%rax
add    $0x4,%rax
cmp    %rdx,%r9
setb   %bpl
cmp    %rcx,%r12
setb   %dl
cmp    %rax,%r9
setb   %al
cmp    %rcx,%r11
setb   %cl
test   %dl,%bpl
jne    113c <twofft+0x2c>
and    %cl,%al
jne    113c <twofft+0x2c>
mov    %rbx,%r8
and    $0xfffffffffffffffc,%r8
lea    -0x4(%r8),%rax
mov    %rax,%rcx
shr    $0x2,%rcx
add    $0x1,%rcx
test   %rax,%rax
je     1437 <twofft+0x327>
mov    %rcx,%rdx
and    $0xfffffffffffffffe,%rdx
xor    %eax,%eax
movsd  0x4(%rdi,%rax,4),%xmm0
movsd  0xc(%rdi,%rax,4),%xmm1
movsd  0x4(%rsi,%rax,4),%xmm2
unpcklps %xmm2,%xmm0
movsd  0xc(%rsi,%rax,4),%xmm2
unpcklps %xmm2,%xmm1
movups %xmm0,0x4(%r15,%rax,8)
movups %xmm1,0x14(%r15,%rax,8)
movsd  0x14(%rdi,%rax,4),%xmm0
movsd  0x1c(%rdi,%rax,4),%xmm1
movsd  0x14(%rsi,%rax,4),%xmm2
unpcklps %xmm2,%xmm0
movsd  0x1c(%rsi,%rax,4),%xmm2
unpcklps %xmm2,%xmm1
movups %xmm0,0x24(%r15,%rax,8)
movups %xmm1,0x34(%r15,%rax,8)
add    $0x8,%rax
add    $0xfffffffffffffffe,%rdx
jne    1200 <twofft+0xf0>
test   $0x1,%cl
je     12a6 <twofft+0x196>
lea    0x0(,%rax,8),%rcx
shl    $0x2,%rax
or     $0x4,%rax
movsd  (%rdi,%rax,1),%xmm0
movsd  0x8(%rdi,%rax,1),%xmm1
movsd  (%rsi,%rax,1),%xmm2
unpcklps %xmm2,%xmm0
movsd  0x8(%rsi,%rax,1),%xmm2
unpcklps %xmm2,%xmm1
mov    %rcx,%rax
or     $0x8,%rax
or     $0x18,%rcx
movups %xmm0,-0x4(%r15,%rax,1)
movups %xmm1,-0x4(%r15,%rcx,1)
cmp    %rbx,%r8
je     1346 <twofft+0x236>
lea    (%r8,%r8,1),%r10
add    $0x2,%r10
or     $0x1,%r8
jmp    113c <twofft+0x2c>
movss  (%rdi,%r8,4),%xmm0
movss  %xmm0,-0x4(%r15,%r10,4)
movss  (%rsi,%r8,4),%xmm0
movss  %xmm0,(%r15,%r10,4)
lea    0x1(%r8),%r9
add    $0x2,%r10
cmp    %rbx,%r8
je     1346 <twofft+0x236>
lea    (%r15,%r10,4),%rax
mov    %rbx,%rcx
sub    %r9,%rcx
add    $0x1,%rcx
lea    (%rsi,%r9,4),%rdx
add    $0x4,%rdx
lea    (%rdi,%r9,4),%rsi
add    $0x4,%rsi
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
movss  -0x4(%rsi,%rdi,4),%xmm0
movss  %xmm0,-0x4(%rax,%rdi,8)
movss  -0x4(%rdx,%rdi,4),%xmm0
movss  %xmm0,(%rax,%rdi,8)
movss  (%rsi,%rdi,4),%xmm0
movss  %xmm0,0x4(%rax,%rdi,8)
movss  (%rdx,%rdi,4),%xmm0
movss  %xmm0,0x8(%rax,%rdi,8)
add    $0x2,%rdi
cmp    %rdi,%rcx
jne    1310 <twofft+0x200>
lea    0x1(%rbx),%r12
mov    %r15,%rdi
mov    %rbx,%rsi
mov    $0x1,%edx
call   1030 <four1@plt>
movss  0x8(%r15),%xmm0
movss  %xmm0,0x4(%r14)
movl   $0x0,0x8(%r14)
movl   $0x0,0x8(%r15)
cmp    $0x3,%r12
jb     142e <twofft+0x31e>
shl    $0x3,%rbx
mov    $0x3,%eax
movss  0xc6f(%rip),%xmm0
movaps 0xc78(%rip),%xmm1
nopl   0x0(%rax,%rax,1)
movss  (%r15,%rax,4),%xmm2
movss  0x4(%r15,%rax,4),%xmm3
movss  -0x4(%r15,%rbx,1),%xmm4
movss  (%r15,%rbx,1),%xmm5
movaps %xmm2,%xmm6
addss  %xmm4,%xmm6
mulss  %xmm0,%xmm6
subss  %xmm4,%xmm2
mulss  %xmm0,%xmm2
movaps %xmm3,%xmm4
addss  %xmm5,%xmm4
mulss  %xmm0,%xmm4
subss  %xmm5,%xmm3
mulss  %xmm0,%xmm3
movss  %xmm6,(%r15,%rax,4)
movss  %xmm3,0x4(%r15,%rax,4)
movss  %xmm6,-0x4(%r15,%rbx,1)
xorps  %xmm1,%xmm3
movss  %xmm3,(%r15,%rbx,1)
movss  %xmm4,(%r14,%rax,4)
movaps %xmm2,%xmm3
xorps  %xmm1,%xmm3
movss  %xmm3,0x4(%r14,%rax,4)
movss  %xmm4,-0x4(%r14,%rbx,1)
movss  %xmm2,(%r14,%rbx,1)
add    $0x2,%rax
add    $0xfffffffffffffff8,%rbx
cmp    %r12,%rax
jbe    13a0 <twofft+0x290>
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %eax,%eax
test   $0x1,%cl
jne    1263 <twofft+0x153>
jmp    12a6 <twofft+0x196>

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

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