#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 */
scalar_t__ cos(float);
scalar_t__ sin(float);
void dctsub(int n, float *a) {
  int i, i0, j, k, m;
  float ec, w1r, w1i, wkr, wki, wdr, wdi, ss, xr, xi, yr, yi;
  ec = ((float)1.570796326794896619231321691639751442098584699687) / n;
  wkr = 0.5;
  wki = 0.5;
  w1r = (float)cos(ec);
  w1i = (float)sin(ec);
  wdr = ((float)0.5) * (w1r - w1i);
  wdi = ((float)0.5) * (w1r + w1i);
  ss = 2 * w1i;
  m = n >> 1;
  i = 0;
  for (;;) {
    i0 = i + 2 * 64;
    if (i0 > m - 2) {
      i0 = m - 2;
    }
    for (j = i + 2; j <= i0; j += 2) {
      k = n - j;
      xr = wdi * a[j - 1] - wdr * a[k + 1];
      xi = wdr * a[j - 1] + wdi * a[k + 1];
      wkr -= ss * wdi;
      wki += ss * wdr;
      yr = wki * a[j] - wkr * a[k];
      yi = wkr * a[j] + wki * a[k];
      wdr -= ss * wki;
      wdi += ss * wkr;
      a[k + 1] = xr;
      a[k] = yr;
      a[j - 1] = xi;
      a[j] = yi;
    }
    if (i0 == m - 2) {
      break;
    }
    wdr = (float)cos(ec * i0);
    wdi = (float)sin(ec * i0);
    wkr = ((float)0.5) * (wdr - wdi);
    wki = ((float)0.5) * (wdr + wdi);
    wdr = wkr * w1r - wki * w1i;
    wdi = wkr * w1i + wki * w1r;
    i = i0;
  }
  xr = wdi * a[m - 1] - wdr * a[m + 1];
  a[m - 1] = wdr * a[m - 1] + wdi * a[m + 1];
  a[m + 1] = xr;
  a[m] *= ((float)0.707106781186547524400844362104849039284835937688);
}


x86 assembly (clang, -O0):
<dctsub>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
cvtsi2ssl -0x4(%rbp),%xmm1
movss  0xecc(%rip),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,-0x28(%rbp)
movss  0xeb7(%rip),%xmm0
movss  %xmm0,-0x34(%rbp)
movss  0xeaa(%rip),%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x28(%rbp),%xmm0
call   1040 <cos@plt>
cvtsi2ss %rax,%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  -0x28(%rbp),%xmm0
call   1030 <sin@plt>
cvtsi2ss %rax,%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x2c(%rbp),%xmm1
subss  -0x30(%rbp),%xmm1
movss  0xe6b(%rip),%xmm0
mulss  %xmm1,%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x2c(%rbp),%xmm1
addss  -0x30(%rbp),%xmm1
movss  0xe50(%rip),%xmm0
mulss  %xmm1,%xmm0
movss  %xmm0,-0x40(%rbp)
movss  0xe3b(%rip),%xmm0
mulss  -0x30(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x4(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x24(%rbp)
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
add    $0x80,%eax
mov    %eax,-0x18(%rbp)
mov    -0x18(%rbp),%eax
mov    -0x24(%rbp),%ecx
sub    $0x2,%ecx
cmp    %ecx,%eax
jle    1204 <dctsub+0xe4>
mov    -0x24(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x18(%rbp)
mov    -0x14(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     13fe <dctsub+0x2de>
mov    -0x4(%rbp),%eax
sub    -0x1c(%rbp),%eax
mov    %eax,-0x20(%rbp)
movss  -0x40(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x3c(%rbp),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mulss  (%rax,%rcx,4),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x48(%rbp)
movss  -0x3c(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x40(%rbp),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mulss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x4c(%rbp)
movss  -0x44(%rbp),%xmm0
movss  -0x40(%rbp),%xmm2
movss  -0x34(%rbp),%xmm1
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x44(%rbp),%xmm0
movss  -0x3c(%rbp),%xmm2
movss  -0x38(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x38(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x34(%rbp),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x50(%rbp)
movss  -0x34(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x38(%rbp),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x54(%rbp)
movss  -0x44(%rbp),%xmm0
movss  -0x38(%rbp),%xmm2
movss  -0x3c(%rbp),%xmm1
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x44(%rbp),%xmm0
movss  -0x34(%rbp),%xmm2
movss  -0x40(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x40(%rbp)
movss  -0x48(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x50(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x4c(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x54(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x1c(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x1c(%rbp)
jmp    120d <dctsub+0xed>
mov    -0x18(%rbp),%eax
mov    -0x24(%rbp),%ecx
sub    $0x2,%ecx
cmp    %ecx,%eax
jne    1414 <dctsub+0x2f4>
jmp    14de <dctsub+0x3be>
movss  -0x28(%rbp),%xmm0
cvtsi2ssl -0x18(%rbp),%xmm1
mulss  %xmm1,%xmm0
call   1040 <cos@plt>
cvtsi2ss %rax,%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x28(%rbp),%xmm0
cvtsi2ssl -0x18(%rbp),%xmm1
mulss  %xmm1,%xmm0
call   1030 <sin@plt>
cvtsi2ss %rax,%xmm0
movss  %xmm0,-0x40(%rbp)
movss  -0x3c(%rbp),%xmm1
subss  -0x40(%rbp),%xmm1
movss  0xba4(%rip),%xmm0
mulss  %xmm1,%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x3c(%rbp),%xmm1
addss  -0x40(%rbp),%xmm1
movss  0xb89(%rip),%xmm0
mulss  %xmm1,%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x34(%rbp),%xmm0
movss  -0x2c(%rbp),%xmm2
movss  -0x38(%rbp),%xmm1
mulss  -0x30(%rbp),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x34(%rbp),%xmm0
movss  -0x30(%rbp),%xmm2
movss  -0x38(%rbp),%xmm1
mulss  -0x2c(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x40(%rbp)
mov    -0x18(%rbp),%eax
mov    %eax,-0x14(%rbp)
jmp    11df <dctsub+0xbf>
movss  -0x40(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x3c(%rbp),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mulss  (%rax,%rcx,4),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x48(%rbp)
movss  -0x3c(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x40(%rbp),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mulss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x48(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  0xa77(%rip),%xmm0
mulss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<dctsub>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x40,%rsp
mov    %edi,%r14d
cvtsi2ss %edi,%xmm1
mov    %rsi,%rbx
movss  0xec2(%rip),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,0x10(%rsp)
call   1040 <cos@plt>
xorps  %xmm0,%xmm0
cvtsi2ss %rax,%xmm0
movss  %xmm0,0xc(%rsp)
movss  0x10(%rsp),%xmm0
call   1030 <sin@plt>
xorps  %xmm0,%xmm0
cvtsi2ss %rax,%xmm0
movss  0xc(%rsp),%xmm1
movaps %xmm1,%xmm8
subss  %xmm0,%xmm8
movss  0xe7f(%rip),%xmm2
mulss  %xmm2,%xmm8
addss  %xmm0,%xmm1
mulss  %xmm2,%xmm1
movss  %xmm0,0x14(%rsp)
movaps %xmm0,%xmm9
addss  %xmm0,%xmm9
movslq %r14d,%rax
sar    %r14d
lea    -0x2(%r14),%r12d
movaps 0xe5d(%rip),%xmm10
xorps  %xmm9,%xmm10
lea    (%rbx,%rax,4),%r15
add    $0xfffffffffffffffc,%r15
xor    %ecx,%ecx
movaps %xmm2,%xmm11
movaps %xmm9,0x30(%rsp)
movaps %xmm10,0x20(%rsp)
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
lea    0x80(%rcx),%eax
cmp    %r12d,%eax
mov    %eax,%ebp
cmovg  %r12d,%ebp
lea    0x2(%rcx),%edx
cmp    %ebp,%edx
jg     12bd <dctsub+0x19d>
movslq %ecx,%rdi
lea    0x2(%rdi),%rcx
movslq %ebp,%rdx
shl    $0x2,%rdi
mov    %r15,%rsi
sub    %rdi,%rsi
xchg   %ax,%ax
movss  -0x4(%rsi),%xmm6
movss  (%rsi),%xmm7
movaps %xmm7,%xmm0
mulss  %xmm8,%xmm0
movss  -0x4(%rbx,%rcx,4),%xmm3
movss  (%rbx,%rcx,4),%xmm4
movaps %xmm1,%xmm5
mulss  %xmm3,%xmm5
subss  %xmm0,%xmm5
mulss  %xmm1,%xmm7
mulss  %xmm8,%xmm3
addss  %xmm7,%xmm3
movaps %xmm10,%xmm0
mulss  %xmm1,%xmm0
addss  %xmm0,%xmm2
movaps %xmm9,%xmm0
mulss  %xmm8,%xmm0
addss  %xmm0,%xmm11
movaps %xmm11,%xmm0
mulss  %xmm6,%xmm0
mulss  %xmm2,%xmm6
movaps %xmm11,%xmm7
mulss  %xmm4,%xmm7
subss  %xmm6,%xmm7
mulss  %xmm2,%xmm4
addss  %xmm0,%xmm4
movaps %xmm10,%xmm0
mulss  %xmm11,%xmm0
addss  %xmm0,%xmm8
movaps %xmm9,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm0,%xmm1
movss  %xmm5,(%rsi)
movss  %xmm7,-0x4(%rsi)
movss  %xmm3,-0x4(%rbx,%rcx,4)
movss  %xmm4,(%rbx,%rcx,4)
add    $0x2,%rcx
add    $0xfffffffffffffff8,%rsi
cmp    %rdx,%rcx
jle    1210 <dctsub+0xf0>
cmp    %eax,%r12d
jle    136c <dctsub+0x24c>
xorps  %xmm0,%xmm0
cvtsi2ss %ebp,%xmm0
mulss  0x10(%rsp),%xmm0
movss  %xmm0,0x1c(%rsp)
call   1040 <cos@plt>
xorps  %xmm0,%xmm0
cvtsi2ss %rax,%xmm0
movss  %xmm0,0x18(%rsp)
movss  0x1c(%rsp),%xmm0
call   1030 <sin@plt>
movss  0x18(%rsp),%xmm11
movaps 0x20(%rsp),%xmm10
movaps 0x30(%rsp),%xmm9
xorps  %xmm0,%xmm0
cvtsi2ss %rax,%xmm0
movaps %xmm11,%xmm2
subss  %xmm0,%xmm2
movss  0xce2(%rip),%xmm1
mulss  %xmm1,%xmm2
addss  %xmm0,%xmm11
mulss  %xmm1,%xmm11
movss  0x14(%rsp),%xmm3
movaps %xmm3,%xmm0
mulss  %xmm11,%xmm0
movaps %xmm2,%xmm8
movss  0xc(%rsp),%xmm1
mulss  %xmm1,%xmm8
subss  %xmm0,%xmm8
movaps %xmm11,%xmm0
mulss  %xmm1,%xmm0
movaps %xmm2,%xmm1
mulss  %xmm3,%xmm1
addss  %xmm0,%xmm1
mov    %ebp,%ecx
jmp    11e0 <dctsub+0xc0>
movslq %r14d,%rax
movss  0x4(%rbx,%rax,4),%xmm0
movaps %xmm1,%xmm2
mulss  %xmm0,%xmm1
mulss  %xmm8,%xmm0
movss  -0x4(%rbx,%rax,4),%xmm3
movss  (%rbx,%rax,4),%xmm4
mulss  %xmm3,%xmm2
subss  %xmm0,%xmm2
mulss  %xmm3,%xmm8
addss  %xmm1,%xmm8
movss  %xmm8,-0x4(%rbx,%rax,4)
movss  %xmm2,0x4(%rbx,%rax,4)
mulss  0xc55(%rip),%xmm4
movss  %xmm4,(%rbx,%rax,4)
add    $0x40,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @dctsub(i32 %0, float* nocapture %1) {
%3 = sitofp i32 %0 to float
%4 = fdiv float 0x3FF921FB60000000, %3
%5 = tail call i64 @cos(float %4)
%6 = sitofp i64 %5 to float
%7 = tail call i64 @sin(float %4)
%8 = sitofp i64 %7 to float
%9 = fsub float %6, %8
%10 = fmul float %9, 5.000000e-01
%11 = fadd float %6, %8
%12 = fmul float %11, 5.000000e-01
%13 = fmul float %8, 2.000000e+00
%14 = ashr i32 %0, 1
%15 = add nsw i32 %14, -2
%16 = fneg float %13
%17 = sext i32 %0 to i64
br label %18
18: ; preds = %72, %2
%19 = phi float [ 5.000000e-01, %2 ], [ %82, %72 ]
%20 = phi float [ %10, %2 ], [ %85, %72 ]
%21 = phi float [ %12, %2 ], [ %87, %72 ]
%22 = phi float [ 5.000000e-01, %2 ], [ %80, %72 ]
%23 = phi i32 [ 0, %2 ], [ %26, %72 ]
%24 = add nsw i32 %23, 128
%25 = icmp sgt i32 %24, %15
%26 = select i1 %25, i32 %15, i32 %24
%27 = add nsw i32 %23, 2
%28 = icmp sgt i32 %27, %26
br i1 %28, label %68, label %29
29: ; preds = %18
%30 = sext i32 %23 to i64
%31 = add nsw i64 %30, 2
%32 = sext i32 %26 to i64
br label %33
33: ; preds = %29, %33
%34 = phi i64 [ %30, %29 ], [ %40, %33 ]
%35 = phi i64 [ %31, %29 ], [ %66, %33 ]
%36 = phi float [ %22, %29 ], [ %53, %33 ]
%37 = phi float [ %21, %29 ], [ %65, %33 ]
%38 = phi float [ %20, %29 ], [ %64, %33 ]
%39 = phi float [ %19, %29 ], [ %54, %33 ]
%40 = add nsw i64 %34, 2
%41 = sub nsw i64 %17, %35
%42 = add nsw i64 %34, 1
%43 = getelementptr inbounds float, float* %1, i64 %42
%44 = load float, float* %43
%45 = add nsw i64 %41, 1
%46 = getelementptr inbounds float, float* %1, i64 %45
%47 = load float, float* %46
%48 = fneg float %38
%49 = fmul float %47, %48
%50 = tail call float @llvm.fmuladd.f32(float %37, float %44, float %49)
%51 = fmul float %37, %47
%52 = tail call float @llvm.fmuladd.f32(float %38, float %44, float %51)
%53 = tail call float @llvm.fmuladd.f32(float %16, float %37, float %36)
%54 = tail call float @llvm.fmuladd.f32(float %13, float %38, float %39)
%55 = getelementptr inbounds float, float* %1, i64 %35
%56 = load float, float* %55
%57 = getelementptr inbounds float, float* %1, i64 %41
%58 = load float, float* %57
%59 = fneg float %53
%60 = fmul float %58, %59
%61 = tail call float @llvm.fmuladd.f32(float %54, float %56, float %60)
%62 = fmul float %54, %58
%63 = tail call float @llvm.fmuladd.f32(float %53, float %56, float %62)
%64 = tail call float @llvm.fmuladd.f32(float %16, float %54, float %38)
%65 = tail call float @llvm.fmuladd.f32(float %13, float %53, float %37)
store float %50, float* %46
store float %61, float* %57
store float %52, float* %43
store float %63, float* %55
%66 = add nsw i64 %35, 2
%67 = icmp sgt i64 %66, %32
br i1 %67, label %68, label %33
68: ; preds = %33, %18
%69 = phi float [ %20, %18 ], [ %64, %33 ]
%70 = phi float [ %21, %18 ], [ %65, %33 ]
%71 = icmp sgt i32 %15, %24
br i1 %71, label %72, label %88
72: ; preds = %68
%73 = sitofp i32 %26 to float
%74 = fmul float %4, %73
%75 = tail call i64 @cos(float %74)
%76 = sitofp i64 %75 to float
%77 = tail call i64 @sin(float %74)
%78 = sitofp i64 %77 to float
%79 = fsub float %76, %78
%80 = fmul float %79, 5.000000e-01
%81 = fadd float %76, %78
%82 = fmul float %81, 5.000000e-01
%83 = fneg float %82
%84 = fmul float %8, %83
%85 = tail call float @llvm.fmuladd.f32(float %80, float %6, float %84)
%86 = fmul float %82, %6
%87 = tail call float @llvm.fmuladd.f32(float %80, float %8, float %86)
br label %18
88: ; preds = %68
%89 = add nsw i32 %14, -1
%90 = sext i32 %89 to i64
%91 = getelementptr inbounds float, float* %1, i64 %90
%92 = load float, float* %91
%93 = add nsw i32 %14, 1
%94 = sext i32 %93 to i64
%95 = getelementptr inbounds float, float* %1, i64 %94
%96 = load float, float* %95
%97 = fneg float %69
%98 = fmul float %96, %97
%99 = tail call float @llvm.fmuladd.f32(float %70, float %92, float %98)
%100 = fmul float %70, %96
%101 = tail call float @llvm.fmuladd.f32(float %69, float %92, float %100)
store float %101, float* %91
store float %99, float* %95
%102 = sext i32 %14 to i64
%103 = getelementptr inbounds float, float* %1, i64 %102
%104 = load float, float* %103
%105 = fmul float %104, 0x3FE6A09E60000000
store float %105, float* %103
ret void
}