#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 */
typedef double coord;

/* Variables and functions */
double cos(int);
double fabs(double);
double floor(double);
double sin(double);
void make_lookup_sinc(double coord, int q, int dim, int *d1, double *table,
                      double **ptpend) {
  register int d2, d, fcoord;
  register double *tp, *tpend, dtmp, sm;
  static int oq = 0, k0, k1;
  coord--;
  if (fabs(coord - floor((coord) + 0.5)) < 0.00001) {
    *d1 = floor((coord) + 0.5);
    if (*d1 < 0 || *d1 >= dim)
      *ptpend = table - 1;
    else {
      table[0] = 1.0;
      *ptpend = table;
    }
  } else {
    if (q != oq) {
      oq = q;
      if (q % 2) {
        k0 = -(q - 1) / 2;
        k1 = (q - 1) / 2;
      } else {
        k0 = -(q - 2) / 2;
        k1 = q / 2;
      }
    }
    fcoord = floor(q % 2 ? coord + 0.5 : coord);
    *d1 = fcoord + k0;
    if (*d1 < 0)
      *d1 = 0;
    d2 = fcoord + k1;
    if (d2 >= dim)
      d2 = dim - 1;
    *ptpend = tpend = table + (d2 - *d1);
    d = *d1, tp = table;
    sm = 0;
    while (tp <= tpend) {
      dtmp = 3.14159265358979323846 * (coord - (d++));
      sm += (*(tp++) = sin(dtmp) / dtmp * 0.5 * (1.0 + cos(2 * dtmp / q)));
    }
    tp = table;
    while (tp <= tpend)
      *(tp++) /= sm;
  }
}


x86 assembly (clang, -O0):
<make_lookup_sinc>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
movsd  %xmm0,-0x8(%rbp)
mov    %edi,-0xc(%rbp)
mov    %esi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
movsd  0xeb9(%rip),%xmm0
addsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  0xe95(%rip),%xmm0
addsd  -0x8(%rbp),%xmm0
call   1050 <floor@plt>
movsd  -0x60(%rbp),%xmm1
subsd  %xmm0,%xmm1
movaps 0xea3(%rip),%xmm0
pand   %xmm0,%xmm1
movsd  0xe67(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    120f <make_lookup_sinc+0xdf>
movsd  0xe5d(%rip),%xmm0
addsd  -0x8(%rbp),%xmm0
call   1050 <floor@plt>
cvttsd2si %xmm0,%ecx
mov    -0x18(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x18(%rbp),%rax
cmpl   $0x0,(%rax)
jl     11db <make_lookup_sinc+0xab>
mov    -0x18(%rbp),%rax
mov    (%rax),%eax
cmp    -0x10(%rbp),%eax
jl     11ef <make_lookup_sinc+0xbf>
mov    -0x20(%rbp),%rcx
add    $0xfffffffffffffff8,%rcx
mov    -0x28(%rbp),%rax
mov    %rcx,(%rax)
jmp    120a <make_lookup_sinc+0xda>
mov    -0x20(%rbp),%rax
movsd  0xe1d(%rip),%xmm0
movsd  %xmm0,(%rax)
mov    -0x20(%rbp),%rcx
mov    -0x28(%rbp),%rax
mov    %rcx,(%rax)
jmp    1451 <make_lookup_sinc+0x321>
mov    -0xc(%rbp),%eax
cmp    0x2e24(%rip),%eax
je     129a <make_lookup_sinc+0x16a>
mov    -0xc(%rbp),%eax
mov    %eax,0x2e15(%rip)
mov    -0xc(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cmp    $0x0,%edx
je     126c <make_lookup_sinc+0x13c>
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
xor    %eax,%eax
sub    %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,0x2ded(%rip)
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,0x2ddd(%rip)
jmp    1295 <make_lookup_sinc+0x165>
mov    -0xc(%rbp),%ecx
sub    $0x2,%ecx
xor    %eax,%eax
sub    %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,0x2dbc(%rip)
mov    -0xc(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,0x2daf(%rip)
jmp    129a <make_lookup_sinc+0x16a>
mov    -0xc(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cmp    $0x0,%edx
je     12c5 <make_lookup_sinc+0x195>
movsd  0xd52(%rip),%xmm0
addsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
jmp    12cf <make_lookup_sinc+0x19f>
movsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x68(%rbp),%xmm0
call   1050 <floor@plt>
cvttsd2si %xmm0,%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%ecx
add    0x2d57(%rip),%ecx
mov    -0x18(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x18(%rbp),%rax
cmpl   $0x0,(%rax)
jge    1306 <make_lookup_sinc+0x1d6>
mov    -0x18(%rbp),%rax
movl   $0x0,(%rax)
mov    -0x34(%rbp),%eax
add    0x2d35(%rip),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jl     1327 <make_lookup_sinc+0x1f7>
mov    -0x10(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x20(%rbp),%rcx
mov    -0x2c(%rbp),%eax
mov    -0x18(%rbp),%rdx
sub    (%rdx),%eax
cltq
shl    $0x3,%rax
add    %rax,%rcx
mov    %rcx,-0x48(%rbp)
mov    -0x28(%rbp),%rax
mov    %rcx,(%rax)
mov    -0x18(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x30(%rbp)
mov    -0x20(%rbp),%rax
mov    %rax,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x40(%rbp),%rax
cmp    -0x48(%rbp),%rax
ja     1411 <make_lookup_sinc+0x2e1>
movsd  -0x8(%rbp),%xmm1
mov    -0x30(%rbp),%eax
mov    %eax,%ecx
add    $0x1,%ecx
mov    %ecx,-0x30(%rbp)
cvtsi2sd %eax,%xmm0
subsd  %xmm0,%xmm1
movsd  0xc99(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
call   1030 <sin@plt>
divsd  -0x50(%rbp),%xmm0
movsd  0xc59(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  0xc60(%rip),%xmm0
mulsd  -0x50(%rbp),%xmm0
cvtsi2sdl -0xc(%rbp),%xmm1
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1040 <cos@plt>
movaps %xmm0,%xmm2
movsd  -0x70(%rbp),%xmm0
movsd  0xc31(%rip),%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  %xmm0,(%rax)
addsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
jmp    1361 <make_lookup_sinc+0x231>
mov    -0x20(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x40(%rbp),%rax
cmp    -0x48(%rbp),%rax
ja     144c <make_lookup_sinc+0x31c>
movsd  -0x58(%rbp),%xmm1
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
jmp    1419 <make_lookup_sinc+0x2e9>
jmp    1451 <make_lookup_sinc+0x321>
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<make_lookup_sinc>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
mov    %r8,%rbp
mov    %rcx,%r15
mov    %rdx,%r13
mov    %esi,%r12d
mov    %edi,%ebx
movapd %xmm0,%xmm1
addsd  0xea8(%rip),%xmm1
movsd  %xmm1,0x18(%rsp)
movsd  0xea2(%rip),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x20(%rsp)
call   1050 <floor@plt>
movsd  0x18(%rsp),%xmm3
movapd %xmm3,%xmm1
subsd  %xmm0,%xmm1
andpd  0xea5(%rip),%xmm1
movsd  0xe7d(%rip),%xmm2
ucomisd %xmm1,%xmm2
jbe    11b9 <make_lookup_sinc+0x89>
cvttsd2si %xmm0,%eax
mov    %eax,0x0(%r13)
test   %eax,%eax
js     11e4 <make_lookup_sinc+0xb4>
cmp    %r12d,%eax
jge    11e4 <make_lookup_sinc+0xb4>
movabs $0x3ff0000000000000,%rax
mov    %rax,(%r15)
jmp    11e8 <make_lookup_sinc+0xb8>
cmp    %ebx,0x2e7d(%rip)
mov    %rbp,0x8(%rsp)
jne    11f1 <make_lookup_sinc+0xc1>
movapd %xmm3,%xmm0
test   $0x1,%bl
je     11d5 <make_lookup_sinc+0xa5>
movsd  0x20(%rsp),%xmm0
mov    0x2e64(%rip),%r14d
mov    0x2e62(%rip),%ebp
jmp    1261 <make_lookup_sinc+0x131>
add    $0xfffffffffffffff8,%r15
mov    %r15,0x0(%rbp)
jmp    142a <make_lookup_sinc+0x2fa>
mov    %ebx,0x2e45(%rip)
test   $0x1,%bl
jne    122c <make_lookup_sinc+0xfc>
mov    $0x2,%eax
sub    %ebx,%eax
mov    %eax,%r14d
shr    $0x1f,%r14d
add    %eax,%r14d
sar    %r14d
mov    %r14d,0x2e29(%rip)
mov    %ebx,%ebp
shr    $0x1f,%ebp
add    %ebx,%ebp
sar    %ebp
mov    %ebp,0x2e1e(%rip)
movapd %xmm3,%xmm0
jmp    1261 <make_lookup_sinc+0x131>
lea    -0x1(%rbx),%eax
mov    $0x1,%ecx
sub    %ebx,%ecx
mov    %ecx,%r14d
shr    $0x1f,%r14d
add    %ecx,%r14d
sar    %r14d
mov    %r14d,0x2df6(%rip)
shr    $0x1f,%eax
lea    (%rbx,%rax,1),%ebp
add    $0xffffffff,%ebp
sar    %ebp
mov    %ebp,0x2de9(%rip)
movsd  0x20(%rsp),%xmm0
call   1050 <floor@plt>
cvttsd2si %xmm0,%eax
xor    %ecx,%ecx
add    %eax,%r14d
cmovs  %ecx,%r14d
mov    %r14d,0x0(%r13)
add    %eax,%ebp
lea    -0x1(%r12),%edx
cmp    %r12d,%ebp
cmovl  %ebp,%edx
mov    %edx,%r13d
sub    %r14d,%r13d
movslq %r13d,%rax
lea    (%r15,%rax,8),%r12
mov    0x8(%rsp),%rcx
mov    %r12,(%rcx)
test   %eax,%eax
js     142a <make_lookup_sinc+0x2fa>
mov    %edx,0x14(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sd %ebx,%xmm0
movsd  %xmm0,0x38(%rsp)
xorpd  %xmm3,%xmm3
mov    %r15,%rbx
mov    %r14d,%ebp
nopl   0x0(%rax)
movapd %xmm3,0x20(%rsp)
xorps  %xmm1,%xmm1
cvtsi2sd %ebp,%xmm1
add    $0x1,%ebp
movsd  0x18(%rsp),%xmm0
subsd  %xmm1,%xmm0
mulsd  0xd36(%rip),%xmm0
movsd  %xmm0,0x40(%rsp)
call   1030 <sin@plt>
movsd  0x40(%rsp),%xmm1
divsd  %xmm1,%xmm0
mulsd  0xd09(%rip),%xmm0
movsd  %xmm0,0x8(%rsp)
addsd  %xmm1,%xmm1
divsd  0x38(%rsp),%xmm1
cvttsd2si %xmm1,%edi
call   1040 <cos@plt>
movapd 0x20(%rsp),%xmm3
addsd  0xcfa(%rip),%xmm0
mulsd  0x8(%rsp),%xmm0
movsd  %xmm0,(%rbx)
add    $0x8,%rbx
addsd  %xmm0,%xmm3
cmp    %r12,%rbx
jbe    12c0 <make_lookup_sinc+0x190>
test   %r13d,%r13d
mov    0x14(%rsp),%eax
js     142a <make_lookup_sinc+0x2fa>
cltq
lea    (%r15,%rax,8),%rcx
mov    %r14d,%eax
shl    $0x3,%rax
sub    %rax,%rcx
add    $0x1,%rcx
lea    0x8(%r15),%rax
cmp    %rax,%rcx
cmovbe %rax,%rcx
mov    %r15,%rax
not    %rax
add    %rcx,%rax
cmp    $0x8,%rax
jb     1410 <make_lookup_sinc+0x2e0>
shr    $0x3,%rax
add    $0x1,%rax
mov    %rax,%rcx
and    $0xfffffffffffffffe,%rcx
movapd %xmm3,%xmm0
unpcklpd %xmm3,%xmm0
lea    -0x2(%rcx),%rsi
mov    %rsi,%rdx
shr    %rdx
add    $0x1,%rdx
test   %rsi,%rsi
je     1439 <make_lookup_sinc+0x309>
mov    %rdx,%rdi
and    $0xfffffffffffffffe,%rdi
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movupd (%r15,%rsi,8),%xmm1
movupd 0x10(%r15,%rsi,8),%xmm2
divpd  %xmm0,%xmm1
movupd %xmm1,(%r15,%rsi,8)
divpd  %xmm0,%xmm2
movupd %xmm2,0x10(%r15,%rsi,8)
add    $0x4,%rsi
add    $0xfffffffffffffffe,%rdi
jne    13c0 <make_lookup_sinc+0x290>
test   $0x1,%dl
je     1401 <make_lookup_sinc+0x2d1>
movupd (%r15,%rsi,8),%xmm1
divpd  %xmm0,%xmm1
movupd %xmm1,(%r15,%rsi,8)
cmp    %rcx,%rax
je     142a <make_lookup_sinc+0x2fa>
lea    (%r15,%rcx,8),%r15
nopw   0x0(%rax,%rax,1)
movsd  (%r15),%xmm0
divsd  %xmm3,%xmm0
movsd  %xmm0,(%r15)
lea    0x8(%r15),%rax
mov    %rax,%r15
cmp    %r12,%rax
jbe    1410 <make_lookup_sinc+0x2e0>
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %esi,%esi
test   $0x1,%dl
jne    13f1 <make_lookup_sinc+0x2c1>
jmp    1401 <make_lookup_sinc+0x2d1>

LLVM IR (-O0):
@make_lookup_sinc.oq = internal global i32 0
@make_lookup_sinc.k0 = internal global i32 0
@make_lookup_sinc.k1 = internal global i32 0
define void @make_lookup_sinc(double %0, i32 %1, i32 %2, i32* %3, double* %4, double** %5) {
%7 = alloca double
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32*
%11 = alloca double*
%12 = alloca double**
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca double*
%17 = alloca double*
%18 = alloca double
%19 = alloca double
store double %0, double* %7
store i32 %1, i32* %8
store i32 %2, i32* %9
store i32* %3, i32** %10
store double* %4, double** %11
store double** %5, double*** %12
%20 = load double, double* %7
%21 = fadd double %20, -1.000000e+00
store double %21, double* %7
%22 = load double, double* %7
%23 = load double, double* %7
%24 = fadd double %23, 5.000000e-01
%25 = call double @llvm.floor.f64(double %24)
%26 = fsub double %22, %25
%27 = call double @llvm.fabs.f64(double %26)
%28 = fcmp olt double %27, 1.000000e-05
br i1 %28, label %29, label %53
29: ; preds = %6
%30 = load double, double* %7
%31 = fadd double %30, 5.000000e-01
%32 = call double @llvm.floor.f64(double %31)
%33 = fptosi double %32 to i32
%34 = load i32*, i32** %10
store i32 %33, i32* %34
%35 = load i32*, i32** %10
%36 = load i32, i32* %35
%37 = icmp slt i32 %36, 0
br i1 %37, label %43, label %38
38: ; preds = %29
%39 = load i32*, i32** %10
%40 = load i32, i32* %39
%41 = load i32, i32* %9
%42 = icmp sge i32 %40, %41
br i1 %42, label %43, label %47
43: ; preds = %38, %29
%44 = load double*, double** %11
%45 = getelementptr inbounds double, double* %44, i64 -1
%46 = load double**, double*** %12
store double* %45, double** %46
br label %52
47: ; preds = %38
%48 = load double*, double** %11
%49 = getelementptr inbounds double, double* %48, i64 0
store double 1.000000e+00, double* %49
%50 = load double*, double** %11
%51 = load double**, double*** %12
store double* %50, double** %51
br label %52
52: ; preds = %47, %43
br label %164
53: ; preds = %6
%54 = load i32, i32* %8
%55 = load i32, i32* @make_lookup_sinc.oq
%56 = icmp ne i32 %54, %55
br i1 %56, label %57, label %78
57: ; preds = %53
%58 = load i32, i32* %8
store i32 %58, i32* @make_lookup_sinc.oq
%59 = load i32, i32* %8
%60 = srem i32 %59, 2
%61 = icmp ne i32 %60, 0
br i1 %61, label %62, label %70
62: ; preds = %57
%63 = load i32, i32* %8
%64 = sub nsw i32 %63, 1
%65 = sub nsw i32 0, %64
%66 = sdiv i32 %65, 2
store i32 %66, i32* @make_lookup_sinc.k0
%67 = load i32, i32* %8
%68 = sub nsw i32 %67, 1
%69 = sdiv i32 %68, 2
store i32 %69, i32* @make_lookup_sinc.k1
br label %77
70: ; preds = %57
%71 = load i32, i32* %8
%72 = sub nsw i32 %71, 2
%73 = sub nsw i32 0, %72
%74 = sdiv i32 %73, 2
store i32 %74, i32* @make_lookup_sinc.k0
%75 = load i32, i32* %8
%76 = sdiv i32 %75, 2
store i32 %76, i32* @make_lookup_sinc.k1
br label %77
77: ; preds = %70, %62
br label %78
78: ; preds = %77, %53
%79 = load i32, i32* %8
%80 = srem i32 %79, 2
%81 = icmp ne i32 %80, 0
br i1 %81, label %82, label %85
82: ; preds = %78
%83 = load double, double* %7
%84 = fadd double %83, 5.000000e-01
br label %87
85: ; preds = %78
%86 = load double, double* %7
br label %87
87: ; preds = %85, %82
%88 = phi double [ %84, %82 ], [ %86, %85 ]
%89 = call double @llvm.floor.f64(double %88)
%90 = fptosi double %89 to i32
store i32 %90, i32* %15
%91 = load i32, i32* %15
%92 = load i32, i32* @make_lookup_sinc.k0
%93 = add nsw i32 %91, %92
%94 = load i32*, i32** %10
store i32 %93, i32* %94
%95 = load i32*, i32** %10
%96 = load i32, i32* %95
%97 = icmp slt i32 %96, 0
br i1 %97, label %98, label %100
98: ; preds = %87
%99 = load i32*, i32** %10
store i32 0, i32* %99
br label %100
100: ; preds = %98, %87
%101 = load i32, i32* %15
%102 = load i32, i32* @make_lookup_sinc.k1
%103 = add nsw i32 %101, %102
store i32 %103, i32* %13
%104 = load i32, i32* %13
%105 = load i32, i32* %9
%106 = icmp sge i32 %104, %105
br i1 %106, label %107, label %110
107: ; preds = %100
%108 = load i32, i32* %9
%109 = sub nsw i32 %108, 1
store i32 %109, i32* %13
br label %110
110: ; preds = %107, %100
%111 = load double*, double** %11
%112 = load i32, i32* %13
%113 = load i32*, i32** %10
%114 = load i32, i32* %113
%115 = sub nsw i32 %112, %114
%116 = sext i32 %115 to i64
%117 = getelementptr inbounds double, double* %111, i64 %116
store double* %117, double** %17
%118 = load double**, double*** %12
store double* %117, double** %118
%119 = load i32*, i32** %10
%120 = load i32, i32* %119
store i32 %120, i32* %14
%121 = load double*, double** %11
store double* %121, double** %16
store double 0.000000e+00, double* %19
br label %122
122: ; preds = %126, %110
%123 = load double*, double** %16
%124 = load double*, double** %17
%125 = icmp ule double* %123, %124
br i1 %125, label %126, label %151
126: ; preds = %122
%127 = load double, double* %7
%128 = load i32, i32* %14
%129 = add nsw i32 %128, 1
store i32 %129, i32* %14
%130 = sitofp i32 %128 to double
%131 = fsub double %127, %130
%132 = fmul double 0x400921FB54442D18, %131
store double %132, double* %18
%133 = load double, double* %18
%134 = call double @sin(double %133)
%135 = load double, double* %18
%136 = fdiv double %134, %135
%137 = fmul double %136, 5.000000e-01
%138 = load double, double* %18
%139 = fmul double 2.000000e+00, %138
%140 = load i32, i32* %8
%141 = sitofp i32 %140 to double
%142 = fdiv double %139, %141
%143 = fptosi double %142 to i32
%144 = call double @cos(i32 %143)
%145 = fadd double 1.000000e+00, %144
%146 = fmul double %137, %145
%147 = load double*, double** %16
%148 = getelementptr inbounds double, double* %147, i32 1
store double* %148, double** %16
store double %146, double* %147
%149 = load double, double* %19
%150 = fadd double %149, %146
store double %150, double* %19
br label %122
151: ; preds = %122
%152 = load double*, double** %11
store double* %152, double** %16
br label %153
153: ; preds = %157, %151
%154 = load double*, double** %16
%155 = load double*, double** %17
%156 = icmp ule double* %154, %155
br i1 %156, label %157, label %163
157: ; preds = %153
%158 = load double, double* %19
%159 = load double*, double** %16
%160 = getelementptr inbounds double, double* %159, i32 1
store double* %160, double** %16
%161 = load double, double* %159
%162 = fdiv double %161, %158
store double %162, double* %159
br label %153
163: ; preds = %153
br label %164
164: ; preds = %163, %52
ret void
}

LLVM IR (-O2):
@make_lookup_sinc.oq = internal unnamed_addr global i32 0
@make_lookup_sinc.k0 = internal unnamed_addr global i32 0
@make_lookup_sinc.k1 = internal unnamed_addr global i32 0
define void @make_lookup_sinc(double %0, i32 %1, i32 %2, i32* nocapture %3, double* %4, double** nocapture writeonly %5) {
%7 = ptrtoint double* %4 to i64
%8 = fadd double %0, -1.000000e+00
%9 = fadd double %8, 5.000000e-01
%10 = tail call double @llvm.floor.f64(double %9)
%11 = fsub double %8, %10
%12 = tail call double @llvm.fabs.f64(double %11)
%13 = fcmp olt double %12, 1.000000e-05
br i1 %13, label %14, label %22
14: ; preds = %6
%15 = fptosi double %10 to i32
store i32 %15, i32* %3
%16 = icmp sgt i32 %15, -1
%17 = icmp slt i32 %15, %2
%18 = and i1 %16, %17
br i1 %18, label %21, label %19
19: ; preds = %14
%20 = getelementptr inbounds double, double* %4, i64 -1
br label %146
21: ; preds = %14
store double 1.000000e+00, double* %4
br label %146
22: ; preds = %6
%23 = load i32, i32* @make_lookup_sinc.oq
%24 = icmp eq i32 %23, %1
br i1 %24, label %25, label %31
25: ; preds = %22
%26 = load i32, i32* @make_lookup_sinc.k0
%27 = load i32, i32* @make_lookup_sinc.k1
%28 = and i32 %1, 1
%29 = icmp eq i32 %28, 0
%30 = select i1 %29, double %8, double %9
br label %43
31: ; preds = %22
store i32 %1, i32* @make_lookup_sinc.oq
%32 = and i32 %1, 1
%33 = icmp eq i32 %32, 0
br i1 %33, label %39, label %34
34: ; preds = %31
%35 = add nsw i32 %1, -1
%36 = sub i32 1, %1
%37 = sdiv i32 %36, 2
store i32 %37, i32* @make_lookup_sinc.k0
%38 = sdiv i32 %35, 2
store i32 %38, i32* @make_lookup_sinc.k1
br label %43
39: ; preds = %31
%40 = sub i32 2, %1
%41 = sdiv i32 %40, 2
store i32 %41, i32* @make_lookup_sinc.k0
%42 = sdiv i32 %1, 2
store i32 %42, i32* @make_lookup_sinc.k1
br label %43
43: ; preds = %25, %34, %39
%44 = phi double [ %30, %25 ], [ %9, %34 ], [ %8, %39 ]
%45 = phi i32 [ %27, %25 ], [ %38, %34 ], [ %42, %39 ]
%46 = phi i32 [ %26, %25 ], [ %37, %34 ], [ %41, %39 ]
%47 = tail call double @llvm.floor.f64(double %44)
%48 = fptosi double %47 to i32
%49 = add nsw i32 %46, %48
%50 = icmp sgt i32 %49, 0
%51 = select i1 %50, i32 %49, i32 0
store i32 %51, i32* %3
%52 = add nsw i32 %45, %48
%53 = icmp slt i32 %52, %2
%54 = add nsw i32 %2, -1
%55 = select i1 %53, i32 %52, i32 %54
%56 = sub nsw i32 %55, %51
%57 = sext i32 %56 to i64
%58 = getelementptr inbounds double, double* %4, i64 %57
store double* %58, double** %5
%59 = icmp slt i32 %56, 0
br i1 %59, label %148, label %60
60: ; preds = %43
%61 = sitofp i32 %1 to double
br label %120
62: ; preds = %120
br i1 %59, label %148, label %63
63: ; preds = %62
%64 = sext i32 %55 to i64
%65 = shl nsw i64 %64, 3
%66 = add i64 %65, %7
%67 = add i64 %66, 1
%68 = zext i32 %51 to i64
%69 = mul nsw i64 %68, -8
%70 = add i64 %69, %67
%71 = add i64 %7, 8
%72 = call i64 @llvm.umax.i64(i64 %70, i64 %71)
%73 = xor i64 %7, -1
%74 = add i64 %72, %73
%75 = lshr i64 %74, 3
%76 = add nuw nsw i64 %75, 1
%77 = icmp ult i64 %74, 8
br i1 %77, label %118, label %78
78: ; preds = %63
%79 = and i64 %76, 4611686018427387902
%80 = getelementptr double, double* %4, i64 %79
%81 = insertelement <2 x double> poison, double %138, i64 0
%82 = shufflevector <2 x double> %81, <2 x double> poison, <2 x i32> zeroinitializer
%83 = add nsw i64 %79, -2
%84 = lshr exact i64 %83, 1
%85 = add nuw i64 %84, 1
%86 = and i64 %85, 1
%87 = icmp eq i64 %83, 0
br i1 %87, label %107, label %88
88: ; preds = %78
%89 = and i64 %85, -2
br label %90
90: ; preds = %90, %88
%91 = phi i64 [ 0, %88 ], [ %104, %90 ]
%92 = phi i64 [ 0, %88 ], [ %105, %90 ]
%93 = getelementptr double, double* %4, i64 %91
%94 = bitcast double* %93 to <2 x double>*
%95 = load <2 x double>, <2 x double>* %94
%96 = fdiv <2 x double> %95, %82
%97 = bitcast double* %93 to <2 x double>*
store <2 x double> %96, <2 x double>* %97
%98 = or i64 %91, 2
%99 = getelementptr double, double* %4, i64 %98
%100 = bitcast double* %99 to <2 x double>*
%101 = load <2 x double>, <2 x double>* %100
%102 = fdiv <2 x double> %101, %82
%103 = bitcast double* %99 to <2 x double>*
store <2 x double> %102, <2 x double>* %103
%104 = add nuw i64 %91, 4
%105 = add nuw i64 %92, 2
%106 = icmp eq i64 %105, %89
br i1 %106, label %107, label %90
107: ; preds = %90, %78
%108 = phi i64 [ 0, %78 ], [ %104, %90 ]
%109 = icmp eq i64 %86, 0
br i1 %109, label %116, label %110
110: ; preds = %107
%111 = getelementptr double, double* %4, i64 %108
%112 = bitcast double* %111 to <2 x double>*
%113 = load <2 x double>, <2 x double>* %112
%114 = fdiv <2 x double> %113, %82
%115 = bitcast double* %111 to <2 x double>*
store <2 x double> %114, <2 x double>* %115
br label %116
116: ; preds = %107, %110
%117 = icmp eq i64 %76, %79
br i1 %117, label %148, label %118
118: ; preds = %63, %116
%119 = phi double* [ %4, %63 ], [ %80, %116 ]
br label %140
120: ; preds = %60, %120
%121 = phi double [ 0.000000e+00, %60 ], [ %138, %120 ]
%122 = phi double* [ %4, %60 ], [ %137, %120 ]
%123 = phi i32 [ %51, %60 ], [ %124, %120 ]
%124 = add nuw nsw i32 %123, 1
%125 = sitofp i32 %123 to double
%126 = fsub double %8, %125
%127 = fmul double %126, 0x400921FB54442D18
%128 = tail call double @sin(double %127)
%129 = fdiv double %128, %127
%130 = fmul double %129, 5.000000e-01
%131 = fmul double %127, 2.000000e+00
%132 = fdiv double %131, %61
%133 = fptosi double %132 to i32
%134 = tail call double @cos(i32 %133)
%135 = fadd double %134, 1.000000e+00
%136 = fmul double %130, %135
%137 = getelementptr inbounds double, double* %122, i64 1
store double %136, double* %122
%138 = fadd double %121, %136
%139 = icmp ugt double* %137, %58
br i1 %139, label %62, label %120
140: ; preds = %118, %140
%141 = phi double* [ %142, %140 ], [ %119, %118 ]
%142 = getelementptr inbounds double, double* %141, i64 1
%143 = load double, double* %141
%144 = fdiv double %143, %138
store double %144, double* %141
%145 = icmp ugt double* %142, %58
br i1 %145, label %148, label %140
146: ; preds = %21, %19
%147 = phi double* [ %20, %19 ], [ %4, %21 ]
store double* %147, double** %5
br label %148
148: ; preds = %140, %116, %146, %43, %62
ret void
}