#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 */
int /*<<< orphan*/ __ogg_fdrfftf(int, double *, double *, int *);
void dcsqf1(int n, double *x, double *w, double *xh, int *ifac) {
  int modn, i, k, kc;
  int np2, ns2;
  double xim1;
  ns2 = (n + 1) >> 1;
  np2 = n;
  kc = np2;
  for (k = 1; k < ns2; k++) {
    kc--;
    xh[k] = x[k] + x[kc];
    xh[kc] = x[k] - x[kc];
  }
  modn = n % 2;
  if (modn == 0)
    xh[ns2] = x[ns2] + x[ns2];
  for (k = 1; k < ns2; k++) {
    kc = np2 - k;
    x[k] = w[k - 1] * xh[kc] + w[kc - 1] * xh[k];
    x[kc] = w[k - 1] * xh[k] - w[kc - 1] * xh[kc];
  }
  if (modn == 0)
    x[ns2] = w[ns2 - 1] * xh[ns2];
  __ogg_fdrfftf(n, x, xh, ifac);
  for (i = 2; i < n; i += 2) {
    xim1 = x[i - 1] - x[i];
    x[i] = x[i - 1] + x[i];
    x[i - 1] = xim1;
  }
}


x86 assembly (clang, -O0):
<dcsqf1>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    -0x4(%rbp),%eax
add    $0x1,%eax
sar    $0x1,%eax
mov    %eax,-0x40(%rbp)
mov    -0x4(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
mov    %eax,-0x38(%rbp)
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x40(%rbp),%eax
jge    11bb <dcsqf1+0xab>
mov    -0x38(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x38(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    114a <dcsqf1+0x3a>
mov    -0x4(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %edx,-0x2c(%rbp)
cmpl   $0x0,-0x2c(%rbp)
jne    11fa <dcsqf1+0xea>
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x40(%rbp),%eax
jge    12e1 <dcsqf1+0x1d1>
mov    -0x3c(%rbp),%eax
sub    -0x34(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1201 <dcsqf1+0xf1>
cmpl   $0x0,-0x2c(%rbp)
jne    1317 <dcsqf1+0x207>
mov    -0x18(%rbp),%rax
mov    -0x40(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x40(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4(%rbp),%edi
mov    -0x10(%rbp),%rsi
mov    -0x20(%rbp),%rdx
mov    -0x28(%rbp),%rcx
call   1030 <__ogg_fdrfftf@plt>
movl   $0x2,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    13b3 <dcsqf1+0x2a3>
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x48(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x30(%rbp)
jmp    1332 <dcsqf1+0x222>
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<dcsqf1>:
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rsi,%r12
mov    %edi,%r14d
lea    0x1(%r14),%r10d
sar    %r10d
mov    %edi,%r15d
cmp    $0x3,%edi
jl     121a <dcsqf1+0x10a>
cmp    $0x3,%r10d
mov    $0x2,%r9d
cmovge %r10d,%r9d
add    $0xffffffffffffffff,%r9
mov    $0x1,%edi
cmp    $0x5,%r14d
jge    1161 <dcsqf1+0x51>
mov    %r15,%rsi
test   $0x1,%r9b
jne    11f5 <dcsqf1+0xe5>
jmp    121a <dcsqf1+0x10a>
lea    (%rcx,%r15,8),%rdi
add    $0xfffffffffffffff8,%rdi
lea    (%r12,%r15,8),%rsi
add    $0xfffffffffffffff8,%rsi
mov    %r9,%r11
and    $0xfffffffffffffffe,%r11
neg    %r11
mov    $0x10,%ebx
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  -0x8(%r12,%rbx,1),%xmm0
addsd  (%rsi,%rax,8),%xmm0
movsd  %xmm0,-0x8(%rcx,%rbx,1)
movsd  -0x8(%r12,%rbx,1),%xmm0
subsd  (%rsi,%rax,8),%xmm0
movsd  %xmm0,(%rdi,%rax,8)
movsd  (%r12,%rbx,1),%xmm0
addsd  -0x8(%rsi,%rax,8),%xmm0
movsd  %xmm0,(%rcx,%rbx,1)
movsd  (%r12,%rbx,1),%xmm0
subsd  -0x8(%rsi,%rax,8),%xmm0
movsd  %xmm0,-0x8(%rdi,%rax,8)
add    $0xfffffffffffffffe,%rax
add    $0x10,%rbx
cmp    %rax,%r11
jne    1190 <dcsqf1+0x80>
lea    (%r15,%rax,1),%rsi
mov    $0x1,%edi
sub    %rax,%rdi
test   $0x1,%r9b
je     121a <dcsqf1+0x10a>
movsd  (%r12,%rdi,8),%xmm0
addsd  -0x8(%r12,%rsi,8),%xmm0
movsd  %xmm0,(%rcx,%rdi,8)
movsd  (%r12,%rdi,8),%xmm0
subsd  -0x8(%r12,%rsi,8),%xmm0
movsd  %xmm0,-0x8(%rcx,%rsi,8)
movslq %r10d,%r9
test   $0x1,%r14b
jne    1233 <dcsqf1+0x123>
movsd  (%r12,%r9,8),%xmm0
addsd  %xmm0,%xmm0
movsd  %xmm0,(%rcx,%r9,8)
cmp    $0x3,%r14d
jl     12c2 <dcsqf1+0x1b2>
cmp    $0x3,%r10d
mov    $0x2,%esi
cmovge %r10d,%esi
lea    -0x10(%rdx,%r15,8),%r13
lea    -0x8(%rcx,%r15,8),%rdi
lea    (%r12,%r15,8),%r10
add    $0xfffffffffffffff8,%r10
lea    -0x8(,%rsi,8),%r11
ff
mov    $0x1,%esi
xor    %ebx,%ebx
nopl   0x0(%rax,%rax,1)
movsd  0x0(%r13,%rbx,1),%xmm0
mulsd  (%rcx,%rsi,8),%xmm0
movsd  -0x8(%rdx,%rsi,8),%xmm1
mulsd  (%rdi,%rbx,1),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%r12,%rsi,8)
movsd  (%rdi,%rbx,1),%xmm0
mulsd  0x0(%r13,%rbx,1),%xmm0
movsd  -0x8(%rdx,%rsi,8),%xmm1
mulsd  (%rcx,%rsi,8),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,(%r10,%rbx,1)
add    $0x1,%rsi
add    $0xfffffffffffffff8,%rbx
mov    %r11,%rax
add    %rbx,%rax
jne    1270 <dcsqf1+0x160>
test   $0x1,%r14b
jne    12db <dcsqf1+0x1cb>
movsd  -0x8(%rdx,%r9,8),%xmm0
mulsd  (%rcx,%r9,8),%xmm0
movsd  %xmm0,(%r12,%r9,8)
mov    %r14d,%edi
mov    %r12,%rsi
mov    %rcx,%rdx
mov    %r8,%rcx
call   1030 <__ogg_fdrfftf@plt>
cmp    $0x3,%r14d
jl     145f <dcsqf1+0x34f>
cmp    $0x5,%r15
mov    $0x4,%eax
cmovae %r15,%rax
add    $0xfffffffffffffffd,%rax
mov    $0x2,%ecx
cmp    $0x2,%rax
jb     1430 <dcsqf1+0x320>
shr    %rax
add    $0x1,%rax
movabs $0x7ffffffffffffffe,%rcx
ff ff 7f
and    %rax,%rcx
lea    -0x2(%rcx),%rdx
mov    %rdx,%rsi
shr    %rsi
add    $0x1,%rsi
test   %rdx,%rdx
je     1469 <dcsqf1+0x359>
mov    %rsi,%rdi
and    $0xfffffffffffffffe,%rdi
mov    $0x1,%edx
nopl   (%rax)
movupd (%r12,%rdx,8),%xmm0
movupd 0x10(%r12,%rdx,8),%xmm1
movupd 0x20(%r12,%rdx,8),%xmm2
movupd 0x30(%r12,%rdx,8),%xmm3
movapd %xmm0,%xmm4
unpcklpd %xmm1,%xmm4
unpckhpd %xmm1,%xmm0
movapd %xmm4,%xmm1
subpd  %xmm0,%xmm1
addpd  %xmm4,%xmm0
movapd %xmm1,%xmm4
unpckhpd %xmm0,%xmm4
unpcklpd %xmm0,%xmm1
movupd %xmm1,(%r12,%rdx,8)
movupd %xmm4,0x10(%r12,%rdx,8)
movapd %xmm2,%xmm0
unpcklpd %xmm3,%xmm0
unpckhpd %xmm3,%xmm2
movapd %xmm0,%xmm1
subpd  %xmm2,%xmm1
addpd  %xmm0,%xmm2
movapd %xmm1,%xmm0
unpckhpd %xmm2,%xmm0
unpcklpd %xmm2,%xmm1
movupd %xmm1,0x20(%r12,%rdx,8)
movupd %xmm0,0x30(%r12,%rdx,8)
add    $0x8,%rdx
add    $0xfffffffffffffffe,%rdi
jne    1350 <dcsqf1+0x240>
test   $0x1,%sil
je     1420 <dcsqf1+0x310>
movupd (%r12,%rdx,8),%xmm0
movupd 0x10(%r12,%rdx,8),%xmm1
movapd %xmm0,%xmm2
unpcklpd %xmm1,%xmm2
unpckhpd %xmm1,%xmm0
movapd %xmm2,%xmm1
subpd  %xmm0,%xmm1
addpd  %xmm2,%xmm0
movapd %xmm1,%xmm2
unpckhpd %xmm0,%xmm2
unpcklpd %xmm0,%xmm1
movupd %xmm1,(%r12,%rdx,8)
movupd %xmm2,0x10(%r12,%rdx,8)
cmp    %rcx,%rax
je     145f <dcsqf1+0x34f>
add    %rax,%rax
or     $0x2,%rax
mov    %rax,%rcx
nop
movsd  -0x8(%r12,%rcx,8),%xmm0
movsd  (%r12,%rcx,8),%xmm1
movapd %xmm0,%xmm2
subsd  %xmm1,%xmm2
addsd  %xmm1,%xmm0
movsd  %xmm0,(%r12,%rcx,8)
movsd  %xmm2,-0x8(%r12,%rcx,8)
add    $0x2,%rcx
cmp    %r15,%rcx
jb     1430 <dcsqf1+0x320>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret
mov    $0x1,%edx
test   $0x1,%sil
jne    13e2 <dcsqf1+0x2d2>
jmp    1420 <dcsqf1+0x310>

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

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