void rick_four1_nr(rdata, nn, isign) double *rdata;
int nn;
int isign;
{
  double tempr, tempi;
  double wr, wi, wpr, wpi, wtemp, theta, temp2;
  int n, m, i, j, mmax, istep;
  n = 2 * nn;
  j = 1;
  for (i = 1; i <= n; i += 2) {
    if (j > i) {
      tempr = rdata[j];
      tempi = rdata[j + 1];
      rdata[j] = rdata[i];
      rdata[j + 1] = rdata[i + 1];
      rdata[i] = tempr;
      rdata[i + 1] = tempi;
    }
    m = n / 2;
    while ((m >= 2) && (j > m)) {
      j = j - m;
      m /= 2;
    }
    j += m;
  }
  mmax = 2;
  while (n > mmax) {
    istep = 2 * mmax;
    theta = 6.28318530717959 / (isign * mmax);
    temp2 = sin(0.5 * theta);
    wpr = -2.0 * temp2 * temp2;
    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 * rdata[j] - wi * rdata[j + 1];
        tempi = wr * rdata[j + 1] + wi * rdata[j];
        rdata[j] = rdata[i] - tempr;
        rdata[j + 1] = rdata[i + 1] - tempi;
        rdata[i] = rdata[i] + tempr;
        rdata[i + 1] = rdata[i + 1] + tempi;
      }
      wtemp = wr;
      wr = wr * wpr - wi * wpi + wr;
      wi = wi * wpr + wtemp * wpi + wi;
    }
    mmax = istep;
  }
}


x86 assembly (clang, -O0):
<rick_four1_nr>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    -0xc(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x5c(%rbp)
movl   $0x1,-0x68(%rbp)
movl   $0x1,-0x64(%rbp)
mov    -0x64(%rbp),%eax
cmp    -0x5c(%rbp),%eax
jg     1250 <rick_four1_nr+0x140>
mov    -0x68(%rbp),%eax
cmp    -0x64(%rbp),%eax
jle    11e4 <rick_four1_nr+0xd4>
mov    -0x8(%rbp),%rax
movslq -0x68(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x68(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x68(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    -0x64(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x68(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x18(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x64(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x5c(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x60(%rbp)
xor    %eax,%eax
cmpl   $0x2,-0x60(%rbp)
mov    %al,-0x71(%rbp)
jl     120d <rick_four1_nr+0xfd>
mov    -0x68(%rbp),%eax
cmp    -0x60(%rbp),%eax
setg   %al
mov    %al,-0x71(%rbp)
mov    -0x71(%rbp),%al
test   $0x1,%al
jne    121d <rick_four1_nr+0x10d>
jmp    1239 <rick_four1_nr+0x129>
mov    -0x68(%rbp),%eax
sub    -0x60(%rbp),%eax
mov    %eax,-0x68(%rbp)
mov    -0x60(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x60(%rbp)
jmp    11f2 <rick_four1_nr+0xe2>
mov    -0x60(%rbp),%eax
add    -0x68(%rbp),%eax
mov    %eax,-0x68(%rbp)
mov    -0x64(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x64(%rbp)
jmp    113c <rick_four1_nr+0x2c>
movl   $0x2,-0x6c(%rbp)
mov    -0x5c(%rbp),%eax
cmp    -0x6c(%rbp),%eax
jle    14af <rick_four1_nr+0x39f>
mov    -0x6c(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x70(%rbp)
mov    -0x10(%rbp),%eax
imul   -0x6c(%rbp),%eax
cvtsi2sd %eax,%xmm1
movsd  0xd99(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  0xd80(%rip),%xmm0
mulsd  -0x50(%rbp),%xmm0
call   1030 <sin@plt>
movsd  %xmm0,-0x58(%rbp)
movsd  0xd61(%rip),%xmm0
mulsd  -0x58(%rbp),%xmm0
mulsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x50(%rbp),%xmm0
call   1030 <sin@plt>
movsd  %xmm0,-0x40(%rbp)
movsd  0xd33(%rip),%xmm0
movsd  %xmm0,-0x28(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x1,-0x60(%rbp)
mov    -0x60(%rbp),%eax
cmp    -0x6c(%rbp),%eax
jg     14a4 <rick_four1_nr+0x394>
mov    -0x60(%rbp),%eax
mov    %eax,-0x64(%rbp)
mov    -0x64(%rbp),%eax
cmp    -0x5c(%rbp),%eax
jg     1429 <rick_four1_nr+0x319>
mov    -0x64(%rbp),%eax
add    -0x6c(%rbp),%eax
mov    %eax,-0x68(%rbp)
movsd  -0x28(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x68(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x30(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x68(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  -0x28(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x68(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x30(%rbp),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x68(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x18(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x68(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    -0x64(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x68(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x18(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    -0x64(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x64(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x70(%rbp),%eax
add    -0x64(%rbp),%eax
mov    %eax,-0x64(%rbp)
jmp    12f3 <rick_four1_nr+0x1e3>
movsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x28(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  -0x40(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
addsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x30(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x60(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x60(%rbp)
jmp    12e1 <rick_four1_nr+0x1d1>
mov    -0x70(%rbp),%eax
mov    %eax,-0x6c(%rbp)
jmp    1257 <rick_four1_nr+0x147>
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<rick_four1_nr>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
test   %esi,%esi
jle    1323 <rick_four1_nr+0x213>
mov    %edx,%r14d
mov    %rdi,%rbx
lea    (%rsi,%rsi,1),%r15d
movslq %r15d,%rbp
mov    $0x1,%eax
mov    $0x1,%ecx
jmp    114b <rick_four1_nr+0x3b>
nop
add    %edx,%eax
add    $0x2,%rcx
cmp    %rbp,%rcx
jg     11a1 <rick_four1_nr+0x91>
movslq %eax,%rdx
cmp    %rdx,%rcx
jge    117f <rick_four1_nr+0x6f>
movsd  (%rbx,%rdx,8),%xmm0
movsd  0x8(%rbx,%rdx,8),%xmm1
movsd  (%rbx,%rcx,8),%xmm2
movsd  %xmm2,(%rbx,%rdx,8)
movsd  0x8(%rbx,%rcx,8),%xmm2
movsd  %xmm2,0x8(%rbx,%rdx,8)
movsd  %xmm0,(%rbx,%rcx,8)
movsd  %xmm1,0x8(%rbx,%rcx,8)
mov    %esi,%edx
cmp    $0x1,%esi
je     1140 <rick_four1_nr+0x30>
cmp    %esi,%eax
jle    1140 <rick_four1_nr+0x30>
nopw   0x0(%rax,%rax,1)
mov    %edx,%edi
sub    %edx,%eax
shr    %edx
cmp    $0x4,%edi
jb     1140 <rick_four1_nr+0x30>
cmp    %edx,%eax
jg     1190 <rick_four1_nr+0x80>
jmp    1140 <rick_four1_nr+0x30>
cmp    $0x2,%esi
jl     1323 <rick_four1_nr+0x213>
mov    $0x2,%r12d
jmp    11c9 <rick_four1_nr+0xb9>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
cmp    %r12d,%r15d
jle    1323 <rick_four1_nr+0x213>
mov    %r12d,%r13d
lea    0x0(%r13,%r13,1),%r12d
mov    %r13d,%eax
imul   %r14d,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  0xe19(%rip),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,0x10(%rsp)
movapd %xmm1,%xmm0
mulsd  0xe0b(%rip),%xmm0
call   1030 <sin@plt>
movsd  %xmm0,0x8(%rsp)
movsd  0x10(%rsp),%xmm0
call   1030 <sin@plt>
movapd 0xe04(%rip),%xmm10
test   %r13d,%r13d
jle    11c0 <rick_four1_nr+0xb0>
movsd  0x8(%rsp),%xmm1
movapd %xmm1,%xmm8
mulsd  0xddb(%rip),%xmm8
mulsd  %xmm1,%xmm8
movslq %r12d,%rax
mov    %r13d,%ecx
unpcklpd %xmm8,%xmm8
unpcklpd %xmm0,%xmm0
lea    (%rbx,%rcx,8),%rcx
mov    $0x1,%edx
movsd  0xdbd(%rip),%xmm9
jmp    1289 <rick_four1_nr+0x179>
nopl   (%rax)
unpcklpd %xmm9,%xmm3
mulpd  %xmm0,%xmm3
movapd %xmm9,%xmm1
mulpd  %xmm8,%xmm1
addpd  %xmm3,%xmm1
addpd  %xmm1,%xmm9
add    $0x2,%rdx
cmp    %edx,%r13d
jl     11c0 <rick_four1_nr+0xb0>
movapd %xmm9,%xmm4
unpckhpd %xmm9,%xmm4
movapd %xmm4,%xmm3
xorpd  %xmm10,%xmm3
cmp    %edx,%r15d
jl     1260 <rick_four1_nr+0x150>
unpcklpd %xmm3,%xmm4
movapd %xmm9,%xmm5
unpcklpd %xmm9,%xmm5
mov    %rdx,%rsi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movupd (%rcx,%rsi,8),%xmm6
movapd %xmm6,%xmm7
mulpd  %xmm4,%xmm7
shufpd $0x1,%xmm7,%xmm7
mulpd  %xmm5,%xmm6
addpd  %xmm7,%xmm6
movsd  (%rbx,%rsi,8),%xmm7
subsd  %xmm6,%xmm7
movsd  %xmm7,(%rcx,%rsi,8)
movsd  0x8(%rbx,%rsi,8),%xmm2
movapd %xmm6,%xmm1
unpckhpd %xmm6,%xmm1
movapd %xmm2,%xmm7
subsd  %xmm1,%xmm7
movsd  %xmm7,0x8(%rcx,%rsi,8)
movsd  (%rbx,%rsi,8),%xmm1
unpcklpd %xmm2,%xmm1
addpd  %xmm6,%xmm1
movupd %xmm1,(%rbx,%rsi,8)
add    %rax,%rsi
cmp    %rbp,%rsi
jle    12c0 <rick_four1_nr+0x1b0>
jmp    1260 <rick_four1_nr+0x150>
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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