void four1(data, nn, isign) float data[];
unsigned long nn;
int isign;
{
  unsigned long n, mmax, m, j, istep, i;
  double wtemp, wr, wpr, wpi, wi, theta;
  float tempr, tempi;
  n = nn << 1;
  j = 1;
  for (i = 1; i < n; i += 2) {
    if (j > i) {
      tempr = (data[j]);
      (data[j]) = (data[i]);
      (data[i]) = tempr;
      tempr = (data[j + 1]);
      (data[j + 1]) = (data[i + 1]);
      (data[i + 1]) = tempr;
    }
    m = n >> 1;
    while (m >= 2 && j > m) {
      j -= m;
      m >>= 1;
    }
    j += m;
  }
  mmax = 2;
  while (n > mmax) {
    istep = mmax << 1;
    theta = isign * (6.28318530717959 / mmax);
    wtemp = sin(0.5 * theta);
    wpr = -2.0 * wtemp * wtemp;
    wpi = sin(theta);
    wr = 1.0;
    wi = 0.0;
    for (m = 1; m < mmax; m += 2) {
      for (i = m; i <= n; i += istep) {
        j = i + mmax;
        tempr = wr * data[j] - wi * data[j + 1];
        tempi = wr * data[j + 1] + wi * data[j];
        data[j] = data[i] - tempr;
        data[j + 1] = data[i + 1] - tempi;
        data[i] += tempr;
        data[i + 1] += tempi;
      }
      wr = (wtemp = wr) * wpr - wi * wpi + wr;
      wi = wi * wpr + wtemp * wpi + wi;
    }
    mmax = istep;
  }
}


x86 assembly (clang, -O0):
<four1>:
push   %rbp
mov    %rsp,%rbp
sub    $0x90,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    -0x10(%rbp),%rax
shl    $0x1,%rax
mov    %rax,-0x20(%rbp)
movq   $0x1,-0x38(%rbp)
movq   $0x1,-0x48(%rbp)
mov    -0x48(%rbp),%rax
cmp    -0x20(%rbp),%rax
jae    125e <four1+0x14e>
mov    -0x38(%rbp),%rax
cmp    -0x48(%rbp),%rax
jbe    11de <four1+0xce>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x7c(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x7c(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  0x4(%rax,%rcx,4),%xmm0
movss  %xmm0,-0x7c(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%rcx
movss  0x4(%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  %xmm0,0x4(%rax,%rcx,4)
movss  -0x7c(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%rcx
movss  %xmm0,0x4(%rax,%rcx,4)
mov    -0x20(%rbp),%rax
shr    $0x1,%rax
mov    %rax,-0x30(%rbp)
xor    %eax,%eax
cmpq   $0x2,-0x30(%rbp)
mov    %al,-0x81(%rbp)
jb     120e <four1+0xfe>
mov    -0x38(%rbp),%rax
cmp    -0x30(%rbp),%rax
seta   %al
mov    %al,-0x81(%rbp)
mov    -0x81(%rbp),%al
test   $0x1,%al
jne    1221 <four1+0x111>
jmp    1241 <four1+0x131>
mov    -0x30(%rbp),%rcx
mov    -0x38(%rbp),%rax
sub    %rcx,%rax
mov    %rax,-0x38(%rbp)
mov    -0x30(%rbp),%rax
shr    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    11ea <four1+0xda>
mov    -0x30(%rbp),%rax
add    -0x38(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x48(%rbp),%rax
add    $0x2,%rax
mov    %rax,-0x48(%rbp)
jmp    1142 <four1+0x32>
movq   $0x2,-0x28(%rbp)
mov    -0x20(%rbp),%rax
cmp    -0x28(%rbp),%rax
jbe    14ec <four1+0x3dc>
mov    -0x28(%rbp),%rax
add    %rax,%rax
mov    %rax,-0x40(%rbp)
cvtsi2sdl -0x14(%rbp),%xmm0
movaps 0xd95(%rip),%xmm1
movq   -0x28(%rbp),%xmm2
punpckldq %xmm1,%xmm2
movapd 0xd94(%rip),%xmm1
subpd  %xmm1,%xmm2
movaps %xmm2,%xmm1
unpckhpd %xmm2,%xmm2
addsd  %xmm1,%xmm2
movsd  0xd65(%rip),%xmm1
divsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  0xd48(%rip),%xmm0
mulsd  -0x78(%rbp),%xmm0
call   1030 <sin@plt>
movsd  %xmm0,-0x50(%rbp)
movsd  0xd29(%rip),%xmm0
mulsd  -0x50(%rbp),%xmm0
mulsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x78(%rbp),%xmm0
call   1030 <sin@plt>
movsd  %xmm0,-0x68(%rbp)
movsd  0xcfb(%rip),%xmm0
movsd  %xmm0,-0x58(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x70(%rbp)
movq   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x28(%rbp),%rax
jae    14df <four1+0x3cf>
mov    -0x30(%rbp),%rax
mov    %rax,-0x48(%rbp)
mov    -0x48(%rbp),%rax
cmp    -0x20(%rbp),%rax
ja     1466 <four1+0x356>
mov    -0x48(%rbp),%rax
add    -0x28(%rbp),%rax
mov    %rax,-0x38(%rbp)
movsd  -0x58(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm2
movsd  -0x70(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  0x4(%rax,%rcx,4),%xmm3
cvtss2sd %xmm3,%xmm3
mulsd  %xmm3,%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x7c(%rbp)
movsd  -0x58(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  0x4(%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm2
movsd  -0x70(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm3
cvtss2sd %xmm3,%xmm3
mulsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x80(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x7c(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%rcx
movss  0x4(%rax,%rcx,4),%xmm0
subss  -0x80(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  %xmm0,0x4(%rax,%rcx,4)
movss  -0x7c(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
movss  -0x80(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%rcx
addss  0x4(%rax,%rcx,4),%xmm0
movss  %xmm0,0x4(%rax,%rcx,4)
mov    -0x40(%rbp),%rax
add    -0x48(%rbp),%rax
mov    %rax,-0x48(%rbp)
jmp    1330 <four1+0x220>
movsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x60(%rbp),%xmm2
movsd  -0x70(%rbp),%xmm1
mulsd  -0x68(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
addsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x70(%rbp),%xmm0
movsd  -0x60(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  -0x68(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
addsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x30(%rbp),%rax
add    $0x2,%rax
mov    %rax,-0x30(%rbp)
jmp    131a <four1+0x20a>
mov    -0x40(%rbp),%rax
mov    %rax,-0x28(%rbp)
jmp    1266 <four1+0x156>
add    $0x90,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<four1>:
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x48,%rsp
mov    %rsi,%r14
add    %rsi,%r14
je     1383 <four1+0x273>
mov    %rdi,%rbx
movabs $0x7fffffffffffffff,%rax
ff ff 7f
and    %rax,%rsi
mov    $0x1,%eax
mov    $0x1,%r8d
jmp    115c <four1+0x4c>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    %rdi,%rax
add    $0x2,%r8
cmp    %r14,%r8
jae    11b6 <four1+0xa6>
cmp    %r8,%rax
jbe    1191 <four1+0x81>
movss  (%rbx,%rax,4),%xmm0
movss  (%rbx,%r8,4),%xmm1
movss  %xmm1,(%rbx,%rax,4)
movss  %xmm0,(%rbx,%r8,4)
movss  0x4(%rbx,%rax,4),%xmm0
movss  0x4(%rbx,%r8,4),%xmm1
movss  %xmm1,0x4(%rbx,%rax,4)
movss  %xmm0,0x4(%rbx,%r8,4)
mov    %rsi,%rdi
cmp    $0x2,%rsi
jb     1150 <four1+0x40>
cmp    %rsi,%rax
jbe    1150 <four1+0x40>
nop
mov    %rdi,%rcx
sub    %rdi,%rax
shr    %rdi
cmp    $0x4,%rcx
jb     1150 <four1+0x40>
cmp    %rdi,%rax
ja     11a0 <four1+0x90>
jmp    1150 <four1+0x40>
cmp    $0x3,%r14
jb     1383 <four1+0x273>
xorps  %xmm0,%xmm0
cvtsi2sd %edx,%xmm0
movsd  %xmm0,0x8(%rsp)
mov    $0x2,%r15d
xorpd  %xmm0,%xmm0
movhpd 0xe69(%rip),%xmm0
movapd %xmm0,0x20(%rsp)
jmp    11f9 <four1+0xe9>
nopw   0x0(%rax,%rax,1)
cmp    %r15,%r14
jbe    1383 <four1+0x273>
mov    %r15,%r12
add    %r15,%r15
movq   %r12,%xmm0
punpckldq 0xdf4(%rip),%xmm0
subpd  0xdfc(%rip),%xmm0
movapd %xmm0,%xmm1
unpckhpd %xmm0,%xmm1
addsd  %xmm0,%xmm1
movsd  0xe08(%rip),%xmm0
divsd  %xmm1,%xmm0
mulsd  0x8(%rsp),%xmm0
movsd  %xmm0,0x18(%rsp)
mulsd  0xdf8(%rip),%xmm0
call   1030 <sin@plt>
movapd %xmm0,0x30(%rsp)
movapd %xmm0,%xmm1
mulsd  0xde9(%rip),%xmm1
movsd  %xmm1,0x10(%rsp)
movsd  0x18(%rsp),%xmm0
call   1030 <sin@plt>
movapd 0x30(%rsp),%xmm9
movapd 0xda8(%rip),%xmm8
test   %r12,%r12
je     11f0 <four1+0xe0>
mulsd  0x10(%rsp),%xmm9
unpcklpd %xmm0,%xmm0
unpcklpd %xmm9,%xmm9
lea    (%rbx,%r12,4),%rax
add    $0x4,%rax
mov    $0x1,%ecx
movapd 0x20(%rsp),%xmm1
jmp    12db <four1+0x1cb>
cs nopw 0x0(%rax,%rax,1)
movapd %xmm1,%xmm3
shufpd $0x1,%xmm2,%xmm3
mulpd  %xmm0,%xmm3
movapd %xmm1,%xmm2
mulpd  %xmm9,%xmm2
addpd  %xmm3,%xmm2
addpd  %xmm2,%xmm1
add    $0x2,%rcx
cmp    %r12,%rcx
jae    11f0 <four1+0xe0>
movapd %xmm1,%xmm2
xorpd  %xmm8,%xmm2
cmp    %r14,%rcx
ja     12b0 <four1+0x1a0>
movapd %xmm1,%xmm3
unpckhpd %xmm1,%xmm3
mov    %rcx,%rdx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movss  -0x4(%rax,%rdx,4),%xmm4
movss  (%rax,%rdx,4),%xmm5
cvtss2sd %xmm4,%xmm4
cvtss2sd %xmm5,%xmm5
movapd %xmm2,%xmm6
mulsd  %xmm5,%xmm6
movapd %xmm3,%xmm7
mulsd  %xmm4,%xmm7
addsd  %xmm6,%xmm7
xorps  %xmm6,%xmm6
cvtsd2ss %xmm7,%xmm6
mulsd  %xmm1,%xmm4
mulsd  %xmm3,%xmm5
addsd  %xmm4,%xmm5
xorps  %xmm4,%xmm4
cvtsd2ss %xmm5,%xmm4
movss  (%rbx,%rdx,4),%xmm5
subss  %xmm6,%xmm5
movss  %xmm5,-0x4(%rax,%rdx,4)
movss  0x4(%rbx,%rdx,4),%xmm5
movaps %xmm5,%xmm7
subss  %xmm4,%xmm7
movss  %xmm7,(%rax,%rdx,4)
addss  (%rbx,%rdx,4),%xmm6
movss  %xmm6,(%rbx,%rdx,4)
addss  %xmm4,%xmm5
movss  %xmm5,0x4(%rbx,%rdx,4)
add    %r15,%rdx
cmp    %r14,%rdx
jbe    1300 <four1+0x1f0>
jmp    12b0 <four1+0x1a0>
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
ret

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

LLVM IR (-O2):
define void @four1(float* nocapture %0, i64 %1, i32 %2) {
%4 = shl i64 %1, 1
%5 = icmp eq i64 %4, 0
br i1 %5, label %109, label %6
6: ; preds = %3
%7 = and i64 %1, 9223372036854775807
%8 = icmp ugt i64 %7, 1
br label %13
9: ; preds = %39
%10 = icmp ugt i64 %4, 2
br i1 %10, label %11, label %109
11: ; preds = %9
%12 = sitofp i32 %2 to double
br label %47
13: ; preds = %6, %39
%14 = phi i64 [ 1, %6 ], [ %42, %39 ]
%15 = phi i64 [ 1, %6 ], [ %43, %39 ]
%16 = icmp ugt i64 %14, %15
br i1 %16, label %17, label %28
17: ; preds = %13
%18 = getelementptr inbounds float, float* %0, i64 %14
%19 = load float, float* %18
%20 = getelementptr inbounds float, float* %0, i64 %15
%21 = load float, float* %20
store float %21, float* %18
store float %19, float* %20
%22 = add i64 %14, 1
%23 = getelementptr inbounds float, float* %0, i64 %22
%24 = load float, float* %23
%25 = add nuw i64 %15, 1
%26 = getelementptr inbounds float, float* %0, i64 %25
%27 = load float, float* %26
store float %27, float* %23
store float %24, float* %26
br label %28
28: ; preds = %17, %13
%29 = icmp ugt i64 %14, %7
%30 = select i1 %8, i1 %29, i1 false
br i1 %30, label %31, label %39
31: ; preds = %28, %31
%32 = phi i64 [ %35, %31 ], [ %7, %28 ]
%33 = phi i64 [ %34, %31 ], [ %14, %28 ]
%34 = sub i64 %33, %32
%35 = lshr i64 %32, 1
%36 = icmp ugt i64 %32, 3
%37 = icmp ugt i64 %34, %35
%38 = select i1 %36, i1 %37, i1 false
br i1 %38, label %31, label %39
39: ; preds = %31, %28
%40 = phi i64 [ %14, %28 ], [ %34, %31 ]
%41 = phi i64 [ %7, %28 ], [ %35, %31 ]
%42 = add i64 %41, %40
%43 = add nuw i64 %15, 2
%44 = icmp ult i64 %43, %4
br i1 %44, label %13, label %9
45: ; preds = %100, %47
%46 = icmp ugt i64 %4, %49
br i1 %46, label %47, label %109
47: ; preds = %11, %45
%48 = phi i64 [ 2, %11 ], [ %49, %45 ]
%49 = shl i64 %48, 1
%50 = uitofp i64 %48 to double
%51 = fdiv double 0x401921FB54442D1C, %50
%52 = fmul double %51, %12
%53 = fmul double %52, 5.000000e-01
%54 = tail call double @sin(double %53)
%55 = fmul double %54, -2.000000e+00
%56 = fmul double %54, %55
%57 = tail call double @sin(double %52)
%58 = icmp eq i64 %48, 0
br i1 %58, label %45, label %59
59: ; preds = %47
%60 = insertelement <2 x double> poison, double %57, i64 0
%61 = shufflevector <2 x double> %60, <2 x double> poison, <2 x i32> zeroinitializer
%62 = insertelement <2 x double> poison, double %56, i64 0
%63 = shufflevector <2 x double> %62, <2 x double> poison, <2 x i32> zeroinitializer
br label %64
64: ; preds = %59, %100
%65 = phi i64 [ %107, %100 ], [ 1, %59 ]
%66 = phi <2 x double> [ %106, %100 ], [ <double 0.000000e+00, double 1.000000e+00>, %59 ]
%67 = icmp ugt i64 %65, %4
%68 = extractelement <2 x double> %66, i64 0
%69 = fneg double %68
br i1 %67, label %100, label %70
70: ; preds = %64
%71 = extractelement <2 x double> %66, i64 1
br label %72
72: ; preds = %70, %72
%73 = phi i64 [ %98, %72 ], [ %65, %70 ]
%74 = add i64 %73, %48
%75 = getelementptr inbounds float, float* %0, i64 %74
%76 = load float, float* %75
%77 = fpext float %76 to double
%78 = add i64 %74, 1
%79 = getelementptr inbounds float, float* %0, i64 %78
%80 = load float, float* %79
%81 = fpext float %80 to double
%82 = fmul double %69, %81
%83 = tail call double @llvm.fmuladd.f64(double %71, double %77, double %82)
%84 = fptrunc double %83 to float
%85 = fmul double %68, %77
%86 = tail call double @llvm.fmuladd.f64(double %71, double %81, double %85)
%87 = fptrunc double %86 to float
%88 = getelementptr inbounds float, float* %0, i64 %73
%89 = load float, float* %88
%90 = fsub float %89, %84
store float %90, float* %75
%91 = add i64 %73, 1
%92 = getelementptr inbounds float, float* %0, i64 %91
%93 = load float, float* %92
%94 = fsub float %93, %87
store float %94, float* %79
%95 = load float, float* %88
%96 = fadd float %95, %84
store float %96, float* %88
%97 = fadd float %93, %87
store float %97, float* %92
%98 = add i64 %73, %49
%99 = icmp ugt i64 %98, %4
br i1 %99, label %100, label %72
100: ; preds = %72, %64
%101 = extractelement <2 x double> %66, i64 1
%102 = insertelement <2 x double> poison, double %101, i64 0
%103 = insertelement <2 x double> %102, double %69, i64 1
%104 = fmul <2 x double> %61, %103
%105 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %66, <2 x double> %63, <2 x double> %104)
%106 = fadd <2 x double> %66, %105
%107 = add nuw i64 %65, 2
%108 = icmp ult i64 %107, %48
br i1 %108, label %64, label %45
109: ; preds = %45, %3, %9
ret void
}