#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
double sin(double);
void four1(double *data, unsigned long nn, int isign) {
  unsigned long n, mmax, m, j, istep, i;
  double wtemp, wr, wpr, wpi, wi, theta;
  double 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
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x80(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  0x8(%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%rcx
movsd  0x8(%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  %xmm0,0x8(%rax,%rcx,8)
movsd  -0x80(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%rcx
movsd  %xmm0,0x8(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
shr    $0x1,%rax
mov    %rax,-0x30(%rbp)
xor    %eax,%eax
cmpq   $0x2,-0x30(%rbp)
mov    %al,-0x89(%rbp)
jb     120e <four1+0xfe>
mov    -0x38(%rbp),%rax
cmp    -0x30(%rbp),%rax
seta   %al
mov    %al,-0x89(%rbp)
mov    -0x89(%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    14d5 <four1+0x3c5>
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    14c8 <four1+0x3b8>
mov    -0x30(%rbp),%rax
mov    %rax,-0x48(%rbp)
mov    -0x48(%rbp),%rax
cmp    -0x20(%rbp),%rax
ja     144f <four1+0x33f>
mov    -0x48(%rbp),%rax
add    -0x28(%rbp),%rax
mov    %rax,-0x38(%rbp)
movsd  -0x58(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x70(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
mulsd  0x8(%rax,%rcx,8),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x58(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  0x8(%rax,%rcx,8),%xmm2
movsd  -0x70(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x80(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%rcx
movsd  0x8(%rax,%rcx,8),%xmm0
subsd  -0x88(%rbp),%xmm0
ff
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  %xmm0,0x8(%rax,%rcx,8)
movsd  -0x80(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x88(%rbp),%xmm0
ff
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%rcx
addsd  0x8(%rax,%rcx,8),%xmm0
movsd  %xmm0,0x8(%rax,%rcx,8)
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     1363 <four1+0x253>
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>
movsd  (%rbx,%rax,8),%xmm0
movsd  (%rbx,%r8,8),%xmm1
movsd  %xmm1,(%rbx,%rax,8)
movsd  %xmm0,(%rbx,%r8,8)
movsd  0x8(%rbx,%rax,8),%xmm0
movsd  0x8(%rbx,%r8,8),%xmm1
movsd  %xmm1,0x8(%rbx,%rax,8)
movsd  %xmm0,0x8(%rbx,%r8,8)
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     1363 <four1+0x253>
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    1363 <four1+0x253>
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),%xmm10
movapd 0xda8(%rip),%xmm9
test   %r12,%r12
je     11f0 <four1+0xe0>
mulsd  0x10(%rsp),%xmm10
unpcklpd %xmm0,%xmm0
unpcklpd %xmm10,%xmm10
lea    (%rbx,%r12,8),%rax
mov    $0x1,%ecx
movapd 0x20(%rsp),%xmm1
jmp    12dc <four1+0x1cc>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movapd %xmm1,%xmm2
shufpd $0x1,%xmm8,%xmm2
mulpd  %xmm0,%xmm2
movapd %xmm1,%xmm3
mulpd  %xmm10,%xmm3
addpd  %xmm2,%xmm3
addpd  %xmm3,%xmm1
add    $0x2,%rcx
cmp    %r12,%rcx
jae    11f0 <four1+0xe0>
movapd %xmm1,%xmm8
xorpd  %xmm9,%xmm8
cmp    %r14,%rcx
ja     12b0 <four1+0x1a0>
movapd %xmm1,%xmm3
unpcklpd %xmm8,%xmm3
movapd %xmm1,%xmm4
unpckhpd %xmm1,%xmm4
mov    %rcx,%rdx
nop
movupd (%rax,%rdx,8),%xmm5
movapd %xmm5,%xmm6
mulpd  %xmm3,%xmm6
shufpd $0x1,%xmm6,%xmm6
mulpd  %xmm4,%xmm5
addpd  %xmm6,%xmm5
movsd  (%rbx,%rdx,8),%xmm6
subsd  %xmm5,%xmm6
movsd  %xmm6,(%rax,%rdx,8)
movsd  0x8(%rbx,%rdx,8),%xmm2
movapd %xmm5,%xmm7
unpckhpd %xmm5,%xmm7
movapd %xmm2,%xmm6
subsd  %xmm7,%xmm6
movsd  %xmm6,0x8(%rax,%rdx,8)
movsd  (%rbx,%rdx,8),%xmm6
unpcklpd %xmm2,%xmm6
addpd  %xmm5,%xmm6
movupd %xmm6,(%rbx,%rdx,8)
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(double* %0, i64 %1, i32 %2) {
%4 = alloca double*
%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 double
%20 = alloca double
store double* %0, double** %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 double*, double** %4
%33 = load i64, i64* %10
%34 = getelementptr inbounds double, double* %32, i64 %33
%35 = load double, double* %34
store double %35, double* %19
%36 = load double*, double** %4
%37 = load i64, i64* %12
%38 = getelementptr inbounds double, double* %36, i64 %37
%39 = load double, double* %38
%40 = load double*, double** %4
%41 = load i64, i64* %10
%42 = getelementptr inbounds double, double* %40, i64 %41
store double %39, double* %42
%43 = load double, double* %19
%44 = load double*, double** %4
%45 = load i64, i64* %12
%46 = getelementptr inbounds double, double* %44, i64 %45
store double %43, double* %46
%47 = load double*, double** %4
%48 = load i64, i64* %10
%49 = add i64 %48, 1
%50 = getelementptr inbounds double, double* %47, i64 %49
%51 = load double, double* %50
store double %51, double* %19
%52 = load double*, double** %4
%53 = load i64, i64* %12
%54 = add i64 %53, 1
%55 = getelementptr inbounds double, double* %52, i64 %54
%56 = load double, double* %55
%57 = load double*, double** %4
%58 = load i64, i64* %10
%59 = add i64 %58, 1
%60 = getelementptr inbounds double, double* %57, i64 %59
store double %56, double* %60
%61 = load double, double* %19
%62 = load double*, double** %4
%63 = load i64, i64* %12
%64 = add i64 %63, 1
%65 = getelementptr inbounds double, double* %62, i64 %64
store double %61, double* %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 = %213, %91
%93 = load i64, i64* %7
%94 = load i64, i64* %8
%95 = icmp ugt i64 %93, %94
br i1 %95, label %96, label %215
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 = %210, %96
%115 = load i64, i64* %9
%116 = load i64, i64* %8
%117 = icmp ult i64 %115, %116
br i1 %117, label %118, label %213
118: ; preds = %114
%119 = load i64, i64* %9
store i64 %119, i64* %12
br label %120
120: ; preds = %188, %118
%121 = load i64, i64* %12
%122 = load i64, i64* %7
%123 = icmp ule i64 %121, %122
br i1 %123, label %124, label %192
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 double*, double** %4
%130 = load i64, i64* %10
%131 = getelementptr inbounds double, double* %129, i64 %130
%132 = load double, double* %131
%133 = load double, double* %17
%134 = load double*, double** %4
%135 = load i64, i64* %10
%136 = add i64 %135, 1
%137 = getelementptr inbounds double, double* %134, i64 %136
%138 = load double, double* %137
%139 = fmul double %133, %138
%140 = fneg double %139
%141 = call double @llvm.fmuladd.f64(double %128, double %132, double %140)
store double %141, double* %19
%142 = load double, double* %14
%143 = load double*, double** %4
%144 = load i64, i64* %10
%145 = add i64 %144, 1
%146 = getelementptr inbounds double, double* %143, i64 %145
%147 = load double, double* %146
%148 = load double, double* %17
%149 = load double*, double** %4
%150 = load i64, i64* %10
%151 = getelementptr inbounds double, double* %149, i64 %150
%152 = load double, double* %151
%153 = fmul double %148, %152
%154 = call double @llvm.fmuladd.f64(double %142, double %147, double %153)
store double %154, double* %20
%155 = load double*, double** %4
%156 = load i64, i64* %12
%157 = getelementptr inbounds double, double* %155, i64 %156
%158 = load double, double* %157
%159 = load double, double* %19
%160 = fsub double %158, %159
%161 = load double*, double** %4
%162 = load i64, i64* %10
%163 = getelementptr inbounds double, double* %161, i64 %162
store double %160, double* %163
%164 = load double*, double** %4
%165 = load i64, i64* %12
%166 = add i64 %165, 1
%167 = getelementptr inbounds double, double* %164, i64 %166
%168 = load double, double* %167
%169 = load double, double* %20
%170 = fsub double %168, %169
%171 = load double*, double** %4
%172 = load i64, i64* %10
%173 = add i64 %172, 1
%174 = getelementptr inbounds double, double* %171, i64 %173
store double %170, double* %174
%175 = load double, double* %19
%176 = load double*, double** %4
%177 = load i64, i64* %12
%178 = getelementptr inbounds double, double* %176, i64 %177
%179 = load double, double* %178
%180 = fadd double %179, %175
store double %180, double* %178
%181 = load double, double* %20
%182 = load double*, double** %4
%183 = load i64, i64* %12
%184 = add i64 %183, 1
%185 = getelementptr inbounds double, double* %182, i64 %184
%186 = load double, double* %185
%187 = fadd double %186, %181
store double %187, double* %185
br label %188
188: ; preds = %124
%189 = load i64, i64* %11
%190 = load i64, i64* %12
%191 = add i64 %190, %189
store i64 %191, i64* %12
br label %120
192: ; preds = %120
%193 = load double, double* %14
store double %193, double* %13
%194 = load double, double* %15
%195 = load double, double* %17
%196 = load double, double* %16
%197 = fmul double %195, %196
%198 = fneg double %197
%199 = call double @llvm.fmuladd.f64(double %193, double %194, double %198)
%200 = load double, double* %14
%201 = fadd double %199, %200
store double %201, double* %14
%202 = load double, double* %17
%203 = load double, double* %15
%204 = load double, double* %13
%205 = load double, double* %16
%206 = fmul double %204, %205
%207 = call double @llvm.fmuladd.f64(double %202, double %203, double %206)
%208 = load double, double* %17
%209 = fadd double %207, %208
store double %209, double* %17
br label %210
210: ; preds = %192
%211 = load i64, i64* %9
%212 = add i64 %211, 2
store i64 %212, i64* %9
br label %114
213: ; preds = %114
%214 = load i64, i64* %11
store i64 %214, i64* %8
br label %92
215: ; preds = %92
ret void
}

LLVM IR (-O2):
define void @four1(double* 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 double, double* %0, i64 %14
%19 = load double, double* %18
%20 = getelementptr inbounds double, double* %0, i64 %15
%21 = load double, double* %20
store double %21, double* %18
store double %19, double* %20
%22 = add i64 %14, 1
%23 = getelementptr inbounds double, double* %0, i64 %22
%24 = load double, double* %23
%25 = add nuw i64 %15, 1
%26 = getelementptr inbounds double, double* %0, i64 %25
%27 = load double, double* %26
store double %27, double* %23
store double %24, double* %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 = insertelement <2 x double> %66, double %69, i64 1
%72 = shufflevector <2 x double> %66, <2 x double> undef, <2 x i32> <i32 1, i32 1>
br label %73
73: ; preds = %70, %73
%74 = phi i64 [ %98, %73 ], [ %65, %70 ]
%75 = add i64 %74, %48
%76 = getelementptr inbounds double, double* %0, i64 %75
%77 = add i64 %75, 1
%78 = getelementptr inbounds double, double* %0, i64 %77
%79 = bitcast double* %76 to <2 x double>*
%80 = load <2 x double>, <2 x double>* %79
%81 = fmul <2 x double> %80, %71
%82 = shufflevector <2 x double> %81, <2 x double> poison, <2 x i32> <i32 1, i32 0>
%83 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %72, <2 x double> %80, <2 x double> %82)
%84 = getelementptr inbounds double, double* %0, i64 %74
%85 = load double, double* %84
%86 = extractelement <2 x double> %83, i64 0
%87 = fsub double %85, %86
store double %87, double* %76
%88 = add i64 %74, 1
%89 = getelementptr inbounds double, double* %0, i64 %88
%90 = load double, double* %89
%91 = extractelement <2 x double> %83, i64 1
%92 = fsub double %90, %91
store double %92, double* %78
%93 = load double, double* %84
%94 = insertelement <2 x double> poison, double %93, i64 0
%95 = insertelement <2 x double> %94, double %90, i64 1
%96 = fadd <2 x double> %83, %95
%97 = bitcast double* %84 to <2 x double>*
store <2 x double> %96, <2 x double>* %97
%98 = add i64 %74, %49
%99 = icmp ugt i64 %98, %4
br i1 %99, label %100, label %73
100: ; preds = %73, %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
}