#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*/ assert(int);
long double fabsl(long double);
long double ref_catrom(long double);
void ref_calc_coeffs(long double *coeffs, long double offset, long taps,
                     long ltrim, long rtrim) {
  long i;
  long double tap_offset, tap_mult, fudge, total_check;
  assert(taps - ltrim - rtrim > 0);
  tap_mult = (long double)taps / 4;
  fudge = 0.0;
  for (i = 0; i < taps; i++) {
    if (i < ltrim || i >= taps - rtrim) {
      coeffs[i] = 0;
      continue;
    }
    tap_offset = 1 - offset - taps / 2 + i;
    coeffs[i] = ref_catrom(fabsl(tap_offset) / tap_mult) / tap_mult;
    fudge += coeffs[i];
  }
  total_check = 0.0;
  for (i = 0; i < taps; i++) {
    coeffs[i] /= fudge;
    total_check += coeffs[i];
  }
  assert(fabsl(total_check - 1.0) < 0.0000000001L);
}


x86 assembly (clang, -O0):
<ref_calc_coeffs>:
push   %rbp
mov    %rsp,%rbp
sub    $0xb0,%rsp
fldt   0x10(%rbp)
mov    %rdi,-0x8(%rbp)
fstpt  -0x20(%rbp)
mov    %rsi,-0x28(%rbp)
mov    %rdx,-0x30(%rbp)
mov    %rcx,-0x38(%rbp)
mov    -0x28(%rbp),%rax
mov    -0x30(%rbp),%rcx
sub    %rcx,%rax
mov    -0x38(%rbp),%rcx
sub    %rcx,%rax
setg   %al
movzbl %al,%edi
call   1040 <assert@plt>
mov    -0x28(%rbp),%rax
mov    %rax,-0x88(%rbp)
fildll -0x88(%rbp)
flds   0xe8b(%rip)
fdivrp %st,%st(1)
fstpt  -0x60(%rbp)
fldz
fstpt  -0x70(%rbp)
movq   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%rax
cmp    -0x28(%rbp),%rax
jge    125b <ref_calc_coeffs+0x13b>
mov    -0x40(%rbp),%rax
cmp    -0x30(%rbp),%rax
jl     11b8 <ref_calc_coeffs+0x98>
mov    -0x40(%rbp),%rax
mov    -0x28(%rbp),%rcx
sub    -0x38(%rbp),%rcx
cmp    %rcx,%rax
jl     11ce <ref_calc_coeffs+0xae>
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%rcx
shl    $0x4,%rcx
fldz
fstpt  (%rax,%rcx,1)
jmp    124a <ref_calc_coeffs+0x12a>
fldt   -0x20(%rbp)
fld1
fsubp  %st,%st(1)
mov    -0x28(%rbp),%rax
mov    %rax,%rcx
shr    $0x3f,%rcx
add    %rcx,%rax
sar    %rax
mov    %rax,-0x90(%rbp)
fildll -0x90(%rbp)
fsubrp %st,%st(1)
mov    -0x40(%rbp),%rax
mov    %rax,-0x98(%rbp)
fildll -0x98(%rbp)
faddp  %st,%st(1)
fstpt  -0x50(%rbp)
fldt   -0x50(%rbp)
fabs
fldt   -0x60(%rbp)
fdivrp %st,%st(1)
mov    %rsp,%rax
fstpt  (%rax)
call   1030 <ref_catrom@plt>
fldt   -0x60(%rbp)
fdivrp %st,%st(1)
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%rcx
shl    $0x4,%rcx
fstpt  (%rax,%rcx,1)
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%rcx
shl    $0x4,%rcx
fldt   (%rax,%rcx,1)
fldt   -0x70(%rbp)
faddp  %st,%st(1)
fstpt  -0x70(%rbp)
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    1187 <ref_calc_coeffs+0x67>
fldz
fstpt  -0x80(%rbp)
movq   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%rax
cmp    -0x28(%rbp),%rax
jge    12b5 <ref_calc_coeffs+0x195>
fldt   -0x70(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%rcx
shl    $0x4,%rcx
fldt   (%rax,%rcx,1)
fdivp  %st,%st(1)
fstpt  (%rax,%rcx,1)
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%rcx
shl    $0x4,%rcx
fldt   (%rax,%rcx,1)
fldt   -0x80(%rbp)
faddp  %st,%st(1)
fstpt  -0x80(%rbp)
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    1268 <ref_calc_coeffs+0x148>
fldt   -0x80(%rbp)
fld1
fchs
faddp  %st,%st(1)
fabs
fldt   0xd4a(%rip)
fucomip %st(1),%st
fstp   %st(0)
seta   %al
and    $0x1,%al
movzbl %al,%edi
call   1040 <assert@plt>
add    $0xb0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<ref_calc_coeffs>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x78,%rsp
mov    %rcx,%r12
mov    %rdx,%r14
mov    %rsi,%r15
mov    %rdi,%rbp
fldt   0xb0(%rsp)
fstpt  0x14(%rsp)
lea    (%rdx,%rcx,1),%rax
mov    %rsi,%rcx
sub    %rax,%rcx
xor    %edi,%edi
test   %rcx,%rcx
setg   %dil
call   1040 <assert@plt>
mov    %r15,0x40(%rsp)
fildll 0x40(%rsp)
fmuls  0xe94(%rip)
fstpt  0x2c(%rsp)
fld1
fld    %st(0)
fchs
test   %r15,%r15
jle    1241 <ref_calc_coeffs+0x121>
fstpt  0x48(%rsp)
mov    %r15,%r13
sub    %r12,%r13
fldt   0x14(%rsp)
fsubr  %st(1),%st
mov    %r15,%rax
shr    %rax
mov    %rax,0x38(%rsp)
fildll 0x38(%rsp)
fsubrp %st,%st(1)
fstpt  0x54(%rsp)
fldz
xor    %ebx,%ebx
mov    %rbp,%r12
fld    %st(0)
fld    %st(1)
fstpt  0x20(%rsp)
fld    %st(2)
fstpt  0x60(%rsp)
jmp    11e2 <ref_calc_coeffs+0xc2>
nopl   0x0(%rax,%rax,1)
fldt   0x20(%rsp)
fxch   %st(3)
fxch   %st(1)
fxch   %st(2)
fxch   %st(3)
fstpt  0x0(%rbp)
fadd   %st,%st(1)
add    $0x1,%rbx
add    $0x10,%rbp
cmp    %rbx,%r15
fxch   %st(2)
fxch   %st(1)
je     1224 <ref_calc_coeffs+0x104>
cmp    %r14,%rbx
jl     11c0 <ref_calc_coeffs+0xa0>
cmp    %r13,%rbx
jge    11c0 <ref_calc_coeffs+0xa0>
fstp   %st(2)
fldt   0x54(%rsp)
fadd   %st(2),%st
fabs
fldt   0x2c(%rsp)
fdivrp %st,%st(1)
fstpt  (%rsp)
fstpt  0x14(%rsp)
fstpt  0x6c(%rsp)
call   1030 <ref_catrom@plt>
fldt   0x6c(%rsp)
fldt   0x14(%rsp)
fldt   0x60(%rsp)
fldt   0x2c(%rsp)
fdivrp %st,%st(4)
fxch   %st(1)
fadd   %st(3),%st
jmp    11ca <ref_calc_coeffs+0xaa>
fstp   %st(0)
test   %r15,%r15
fldt   0x48(%rsp)
jle    12af <ref_calc_coeffs+0x18f>
fstp   %st(0)
fldz
cmp    $0x1,%r15
jne    1245 <ref_calc_coeffs+0x125>
xor    %eax,%eax
jmp    1283 <ref_calc_coeffs+0x163>
fstp   %st(1)
jmp    12a7 <ref_calc_coeffs+0x187>
fstp   %st(0)
mov    %r15,%rcx
and    $0xfffffffffffffffe,%rcx
lea    0x10(%r12),%rdx
fldz
xor    %eax,%eax
nopw   0x0(%rax,%rax,1)
fldt   -0x10(%rdx)
fdiv   %st(2),%st
fld    %st(0)
fstpt  -0x10(%rdx)
faddp  %st,%st(1)
fldt   (%rdx)
fdiv   %st(2),%st
fld    %st(0)
fstpt  (%rdx)
faddp  %st,%st(1)
add    $0x2,%rax
add    $0x20,%rdx
cmp    %rax,%rcx
jne    1260 <ref_calc_coeffs+0x140>
test   $0x1,%r15b
je     129f <ref_calc_coeffs+0x17f>
shl    $0x4,%rax
fldt   (%r12,%rax,1)
fdivp  %st,%st(2)
fld    %st(1)
fstpt  (%r12,%rax,1)
faddp  %st,%st(1)
fldz
fxch   %st(1)
fstp   %st(1)
fxch   %st(1)
fchs
faddp  %st,%st(1)
fldz
fldz
fxch   %st(1)
fxch   %st(2)
fstp   %st(2)
fstp   %st(0)
fabs
fldt   0xd55(%rip)
xor    %edi,%edi
fucomip %st(1),%st
fstp   %st(0)
seta   %dil
add    $0x78,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1040 <assert@plt>

LLVM IR (-O0):
define void @ref_calc_coeffs(x86_fp80* %0, x86_fp80 %1, i64 %2, i64 %3, i64 %4) {
%6 = alloca x86_fp80*
%7 = alloca x86_fp80
%8 = alloca i64
%9 = alloca i64
%10 = alloca i64
%11 = alloca i64
%12 = alloca x86_fp80
%13 = alloca x86_fp80
%14 = alloca x86_fp80
%15 = alloca x86_fp80
store x86_fp80* %0, x86_fp80** %6
store x86_fp80 %1, x86_fp80* %7
store i64 %2, i64* %8
store i64 %3, i64* %9
store i64 %4, i64* %10
%16 = load i64, i64* %8
%17 = load i64, i64* %9
%18 = sub nsw i64 %16, %17
%19 = load i64, i64* %10
%20 = sub nsw i64 %18, %19
%21 = icmp sgt i64 %20, 0
%22 = zext i1 %21 to i32
%23 = call i32 @assert(i32 %22)
%24 = load i64, i64* %8
%25 = sitofp i64 %24 to x86_fp80
%26 = fdiv x86_fp80 %25, 0xK40018000000000000000
store x86_fp80 %26, x86_fp80* %13
store x86_fp80 0xK00000000000000000000, x86_fp80* %14
store i64 0, i64* %11
br label %27
27: ; preds = %71, %5
%28 = load i64, i64* %11
%29 = load i64, i64* %8
%30 = icmp slt i64 %28, %29
br i1 %30, label %31, label %74
31: ; preds = %27
%32 = load i64, i64* %11
%33 = load i64, i64* %9
%34 = icmp slt i64 %32, %33
br i1 %34, label %41, label %35
35: ; preds = %31
%36 = load i64, i64* %11
%37 = load i64, i64* %8
%38 = load i64, i64* %10
%39 = sub nsw i64 %37, %38
%40 = icmp sge i64 %36, %39
br i1 %40, label %41, label %45
41: ; preds = %35, %31
%42 = load x86_fp80*, x86_fp80** %6
%43 = load i64, i64* %11
%44 = getelementptr inbounds x86_fp80, x86_fp80* %42, i64 %43
store x86_fp80 0xK00000000000000000000, x86_fp80* %44
br label %71
45: ; preds = %35
%46 = load x86_fp80, x86_fp80* %7
%47 = fsub x86_fp80 0xK3FFF8000000000000000, %46
%48 = load i64, i64* %8
%49 = sdiv i64 %48, 2
%50 = sitofp i64 %49 to x86_fp80
%51 = fsub x86_fp80 %47, %50
%52 = load i64, i64* %11
%53 = sitofp i64 %52 to x86_fp80
%54 = fadd x86_fp80 %51, %53
store x86_fp80 %54, x86_fp80* %12
%55 = load x86_fp80, x86_fp80* %12
%56 = call x86_fp80 @llvm.fabs.f80(x86_fp80 %55)
%57 = load x86_fp80, x86_fp80* %13
%58 = fdiv x86_fp80 %56, %57
%59 = call x86_fp80 @ref_catrom(x86_fp80 %58)
%60 = load x86_fp80, x86_fp80* %13
%61 = fdiv x86_fp80 %59, %60
%62 = load x86_fp80*, x86_fp80** %6
%63 = load i64, i64* %11
%64 = getelementptr inbounds x86_fp80, x86_fp80* %62, i64 %63
store x86_fp80 %61, x86_fp80* %64
%65 = load x86_fp80*, x86_fp80** %6
%66 = load i64, i64* %11
%67 = getelementptr inbounds x86_fp80, x86_fp80* %65, i64 %66
%68 = load x86_fp80, x86_fp80* %67
%69 = load x86_fp80, x86_fp80* %14
%70 = fadd x86_fp80 %69, %68
store x86_fp80 %70, x86_fp80* %14
br label %71
71: ; preds = %45, %41
%72 = load i64, i64* %11
%73 = add nsw i64 %72, 1
store i64 %73, i64* %11
br label %27
74: ; preds = %27
store x86_fp80 0xK00000000000000000000, x86_fp80* %15
store i64 0, i64* %11
br label %75
75: ; preds = %92, %74
%76 = load i64, i64* %11
%77 = load i64, i64* %8
%78 = icmp slt i64 %76, %77
br i1 %78, label %79, label %95
79: ; preds = %75
%80 = load x86_fp80, x86_fp80* %14
%81 = load x86_fp80*, x86_fp80** %6
%82 = load i64, i64* %11
%83 = getelementptr inbounds x86_fp80, x86_fp80* %81, i64 %82
%84 = load x86_fp80, x86_fp80* %83
%85 = fdiv x86_fp80 %84, %80
store x86_fp80 %85, x86_fp80* %83
%86 = load x86_fp80*, x86_fp80** %6
%87 = load i64, i64* %11
%88 = getelementptr inbounds x86_fp80, x86_fp80* %86, i64 %87
%89 = load x86_fp80, x86_fp80* %88
%90 = load x86_fp80, x86_fp80* %15
%91 = fadd x86_fp80 %90, %89
store x86_fp80 %91, x86_fp80* %15
br label %92
92: ; preds = %79
%93 = load i64, i64* %11
%94 = add nsw i64 %93, 1
store i64 %94, i64* %11
br label %75
95: ; preds = %75
%96 = load x86_fp80, x86_fp80* %15
%97 = fsub x86_fp80 %96, 0xK3FFF8000000000000000
%98 = call x86_fp80 @llvm.fabs.f80(x86_fp80 %97)
%99 = fcmp olt x86_fp80 %98, 0xK3FDDDBE6FECEBDEDD5BF
%100 = zext i1 %99 to i32
%101 = call i32 @assert(i32 %100)
ret void
}

LLVM IR (-O2):
define void @ref_calc_coeffs(x86_fp80* nocapture %0, x86_fp80 %1, i64 %2, i64 %3, i64 %4) {
%6 = add i64 %3, %4
%7 = sub i64 %2, %6
%8 = icmp sgt i64 %7, 0
%9 = zext i1 %8 to i32
%10 = tail call i32 @assert(i32 %9)
%11 = sitofp i64 %2 to x86_fp80
%12 = fmul x86_fp80 %11, 0xK3FFD8000000000000000
%13 = icmp sgt i64 %2, 0
br i1 %13, label %14, label %75
14: ; preds = %5
%15 = sub nsw i64 %2, %4
%16 = fsub x86_fp80 0xK3FFF8000000000000000, %1
%17 = lshr i64 %2, 1
%18 = sitofp i64 %17 to x86_fp80
%19 = fsub x86_fp80 %16, %18
br label %26
20: ; preds = %40
br i1 %13, label %21, label %75
21: ; preds = %20
%22 = and i64 %2, 1
%23 = icmp eq i64 %2, 1
br i1 %23, label %62, label %24
24: ; preds = %21
%25 = and i64 %2, -2
br label %46
26: ; preds = %14, %40
%27 = phi i64 [ 0, %14 ], [ %44, %40 ]
%28 = phi x86_fp80 [ 0xK00000000000000000000, %14 ], [ %42, %40 ]
%29 = icmp sge i64 %27, %3
%30 = icmp slt i64 %27, %15
%31 = select i1 %29, i1 %30, i1 false
br i1 %31, label %32, label %40
32: ; preds = %26
%33 = sitofp i64 %27 to x86_fp80
%34 = fadd x86_fp80 %19, %33
%35 = tail call x86_fp80 @llvm.fabs.f80(x86_fp80 %34)
%36 = fdiv x86_fp80 %35, %12
%37 = tail call x86_fp80 @ref_catrom(x86_fp80 %36)
%38 = fdiv x86_fp80 %37, %12
%39 = fadd x86_fp80 %28, %38
br label %40
40: ; preds = %26, %32
%41 = phi x86_fp80 [ %38, %32 ], [ 0xK00000000000000000000, %26 ]
%42 = phi x86_fp80 [ %39, %32 ], [ %28, %26 ]
%43 = getelementptr inbounds x86_fp80, x86_fp80* %0, i64 %27
store x86_fp80 %41, x86_fp80* %43
%44 = add nuw nsw i64 %27, 1
%45 = icmp eq i64 %44, %2
br i1 %45, label %20, label %26
46: ; preds = %46, %24
%47 = phi x86_fp80 [ 0xK00000000000000000000, %24 ], [ %58, %46 ]
%48 = phi i64 [ 0, %24 ], [ %59, %46 ]
%49 = phi i64 [ 0, %24 ], [ %60, %46 ]
%50 = getelementptr inbounds x86_fp80, x86_fp80* %0, i64 %48
%51 = load x86_fp80, x86_fp80* %50
%52 = fdiv x86_fp80 %51, %42
store x86_fp80 %52, x86_fp80* %50
%53 = fadd x86_fp80 %47, %52
%54 = or i64 %48, 1
%55 = getelementptr inbounds x86_fp80, x86_fp80* %0, i64 %54
%56 = load x86_fp80, x86_fp80* %55
%57 = fdiv x86_fp80 %56, %42
store x86_fp80 %57, x86_fp80* %55
%58 = fadd x86_fp80 %53, %57
%59 = add nuw nsw i64 %48, 2
%60 = add i64 %49, 2
%61 = icmp eq i64 %60, %25
br i1 %61, label %62, label %46
62: ; preds = %46, %21
%63 = phi x86_fp80 [ undef, %21 ], [ %58, %46 ]
%64 = phi x86_fp80 [ 0xK00000000000000000000, %21 ], [ %58, %46 ]
%65 = phi i64 [ 0, %21 ], [ %59, %46 ]
%66 = icmp eq i64 %22, 0
br i1 %66, label %72, label %67
67: ; preds = %62
%68 = getelementptr inbounds x86_fp80, x86_fp80* %0, i64 %65
%69 = load x86_fp80, x86_fp80* %68
%70 = fdiv x86_fp80 %69, %42
store x86_fp80 %70, x86_fp80* %68
%71 = fadd x86_fp80 %64, %70
br label %72
72: ; preds = %62, %67
%73 = phi x86_fp80 [ %63, %62 ], [ %71, %67 ]
%74 = fadd x86_fp80 %73, 0xKBFFF8000000000000000
br label %75
75: ; preds = %5, %72, %20
%76 = phi x86_fp80 [ 0xKBFFF8000000000000000, %20 ], [ %74, %72 ], [ 0xKBFFF8000000000000000, %5 ]
%77 = tail call x86_fp80 @llvm.fabs.f80(x86_fp80 %76)
%78 = fcmp olt x86_fp80 %77, 0xK3FDDDBE6FECEBDEDD5BF
%79 = zext i1 %78 to i32
%80 = tail call i32 @assert(i32 %79)
ret void
}