#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 */
float catrom(float);
float fabsf(float);
void calc_coeffs(float *coeffs, float tx, int taps, int ltrim, int rtrim) {
  int i;
  float tmp, tap_mult, fudge;
  tap_mult = (float)taps / 4;
  tx = 1 - tx - taps / 2 + ltrim;
  fudge = 0.0f;
  for (i = ltrim; i < taps - rtrim; i++) {
    tmp = catrom(fabsf(tx) / tap_mult) / tap_mult;
    fudge += tmp;
    coeffs[i] = tmp;
    tx += 1;
  }
  fudge = 1 / fudge;
  for (i = ltrim; i < taps - rtrim; i++) {
    coeffs[i] *= fudge;
  }
}


x86 assembly (clang, -O0):
<calc_coeffs>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
movss  %xmm0,-0xc(%rbp)
mov    %esi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
cvtsi2ssl -0x10(%rbp),%xmm0
movss  0xecd(%rip),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x24(%rbp)
movss  0xeb8(%rip),%xmm0
subss  -0xc(%rbp),%xmm0
mov    -0x10(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cvtsi2ss %eax,%xmm1
subss  %xmm1,%xmm0
cvtsi2ssl -0x14(%rbp),%xmm1
addss  %xmm1,%xmm0
movss  %xmm0,-0xc(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x28(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    -0x10(%rbp),%ecx
sub    -0x18(%rbp),%ecx
cmp    %ecx,%eax
jge    11f2 <calc_coeffs+0xe2>
movss  -0xc(%rbp),%xmm0
movaps 0xe77(%rip),%xmm1
pand   %xmm1,%xmm0
divss  -0x24(%rbp),%xmm0
call   1030 <catrom@plt>
divss  -0x24(%rbp),%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0x20(%rbp),%xmm0
addss  -0x28(%rbp),%xmm0
movss  %xmm0,-0x28(%rbp)
movss  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  0xe26(%rip),%xmm0
addss  -0xc(%rbp),%xmm0
movss  %xmm0,-0xc(%rbp)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    117c <calc_coeffs+0x6c>
movss  0xe06(%rip),%xmm0
divss  -0x28(%rbp),%xmm0
movss  %xmm0,-0x28(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    -0x10(%rbp),%ecx
sub    -0x18(%rbp),%ecx
cmp    %ecx,%eax
jge    1240 <calc_coeffs+0x130>
movss  -0x28(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    120a <calc_coeffs+0xfa>
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<calc_coeffs>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %esi,%r12d
sub    %ecx,%r12d
cmp    %edx,%r12d
jle    12cb <calc_coeffs+0x1bb>
mov    %edx,%r14d
mov    %rdi,%r13
cvtsi2ss %esi,%xmm1
mulss  0xec1(%rip),%xmm1
movss  %xmm1,0x8(%rsp)
movss  0xeb7(%rip),%xmm1
subss  %xmm0,%xmm1
mov    %esi,%eax
shr    $0x1f,%eax
add    %esi,%eax
sar    %eax
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
subss  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsi2ss %edx,%xmm0
addss  %xmm1,%xmm0
movslq %edx,%rbx
movslq %r12d,%r15
xorps  %xmm1,%xmm1
mov    %rbx,%rbp
nopl   0x0(%rax)
movaps %xmm0,0x10(%rsp)
movss  %xmm1,0xc(%rsp)
movaps 0x10(%rsp),%xmm0
andps  0xe79(%rip),%xmm0
divss  0x8(%rsp),%xmm0
call   1030 <catrom@plt>
movss  0xc(%rsp),%xmm1
movss  0xe54(%rip),%xmm5
divss  0x8(%rsp),%xmm0
addss  %xmm0,%xmm1
movss  %xmm0,0x0(%r13,%rbp,4)
movaps 0x10(%rsp),%xmm0
addss  %xmm5,%xmm0
add    $0x1,%rbp
cmp    %r15,%rbp
jl     1180 <calc_coeffs+0x70>
cmp    %r14d,%r12d
jle    12cb <calc_coeffs+0x1bb>
divss  %xmm1,%xmm5
lea    0x1(%rbx),%rax
cmp    %rax,%r15
cmovg  %r15,%rax
sub    %rbx,%rax
cmp    $0x8,%rax
jb     12b0 <calc_coeffs+0x1a0>
mov    %rax,%r8
and    $0xfffffffffffffff8,%r8
movaps %xmm5,%xmm0
shufps $0x0,%xmm5,%xmm0
lea    -0x8(%r8),%rcx
mov    %rcx,%rsi
shr    $0x3,%rsi
add    $0x1,%rsi
test   %rcx,%rcx
je     12da <calc_coeffs+0x1ca>
mov    %rsi,%rdi
and    $0xfffffffffffffffe,%rdi
lea    0x30(,%rbx,4),%rcx
add    %r13,%rcx
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movups -0x30(%rcx,%rdx,4),%xmm1
movups -0x20(%rcx,%rdx,4),%xmm2
movups -0x10(%rcx,%rdx,4),%xmm3
movups (%rcx,%rdx,4),%xmm4
mulps  %xmm0,%xmm1
mulps  %xmm0,%xmm2
movups %xmm1,-0x30(%rcx,%rdx,4)
movups %xmm2,-0x20(%rcx,%rdx,4)
mulps  %xmm0,%xmm3
mulps  %xmm0,%xmm4
movups %xmm3,-0x10(%rcx,%rdx,4)
movups %xmm4,(%rcx,%rdx,4)
add    $0x10,%rdx
add    $0xfffffffffffffffe,%rdi
jne    1240 <calc_coeffs+0x130>
test   $0x1,%sil
je     12a3 <calc_coeffs+0x193>
add    %rbx,%rdx
movups 0x0(%r13,%rdx,4),%xmm1
movups 0x10(%r13,%rdx,4),%xmm2
mulps  %xmm0,%xmm1
mulps  %xmm0,%xmm2
movups %xmm1,0x0(%r13,%rdx,4)
movups %xmm2,0x10(%r13,%rdx,4)
cmp    %r8,%rax
je     12cb <calc_coeffs+0x1bb>
add    %r8,%rbx
nopl   0x0(%rax,%rax,1)
movss  0x0(%r13,%rbx,4),%xmm0
mulss  %xmm5,%xmm0
movss  %xmm0,0x0(%r13,%rbx,4)
add    $0x1,%rbx
cmp    %r15,%rbx
jl     12b0 <calc_coeffs+0x1a0>
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %edx,%edx
test   $0x1,%sil
jne    1282 <calc_coeffs+0x172>
jmp    12a3 <calc_coeffs+0x193>

LLVM IR (-O0):
define void @calc_coeffs(float* %0, float %1, i32 %2, i32 %3, i32 %4) {
%6 = alloca float*
%7 = alloca float
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca float
%13 = alloca float
%14 = alloca float
store float* %0, float** %6
store float %1, float* %7
store i32 %2, i32* %8
store i32 %3, i32* %9
store i32 %4, i32* %10
%15 = load i32, i32* %8
%16 = sitofp i32 %15 to float
%17 = fdiv float %16, 4.000000e+00
store float %17, float* %13
%18 = load float, float* %7
%19 = fsub float 1.000000e+00, %18
%20 = load i32, i32* %8
%21 = sdiv i32 %20, 2
%22 = sitofp i32 %21 to float
%23 = fsub float %19, %22
%24 = load i32, i32* %9
%25 = sitofp i32 %24 to float
%26 = fadd float %23, %25
store float %26, float* %7
store float 0.000000e+00, float* %14
%27 = load i32, i32* %9
store i32 %27, i32* %11
br label %28
28: ; preds = %52, %5
%29 = load i32, i32* %11
%30 = load i32, i32* %8
%31 = load i32, i32* %10
%32 = sub nsw i32 %30, %31
%33 = icmp slt i32 %29, %32
br i1 %33, label %34, label %55
34: ; preds = %28
%35 = load float, float* %7
%36 = call float @llvm.fabs.f32(float %35)
%37 = load float, float* %13
%38 = fdiv float %36, %37
%39 = call float @catrom(float %38)
%40 = load float, float* %13
%41 = fdiv float %39, %40
store float %41, float* %12
%42 = load float, float* %12
%43 = load float, float* %14
%44 = fadd float %43, %42
store float %44, float* %14
%45 = load float, float* %12
%46 = load float*, float** %6
%47 = load i32, i32* %11
%48 = sext i32 %47 to i64
%49 = getelementptr inbounds float, float* %46, i64 %48
store float %45, float* %49
%50 = load float, float* %7
%51 = fadd float %50, 1.000000e+00
store float %51, float* %7
br label %52
52: ; preds = %34
%53 = load i32, i32* %11
%54 = add nsw i32 %53, 1
store i32 %54, i32* %11
br label %28
55: ; preds = %28
%56 = load float, float* %14
%57 = fdiv float 1.000000e+00, %56
store float %57, float* %14
%58 = load i32, i32* %9
store i32 %58, i32* %11
br label %59
59: ; preds = %73, %55
%60 = load i32, i32* %11
%61 = load i32, i32* %8
%62 = load i32, i32* %10
%63 = sub nsw i32 %61, %62
%64 = icmp slt i32 %60, %63
br i1 %64, label %65, label %76
65: ; preds = %59
%66 = load float, float* %14
%67 = load float*, float** %6
%68 = load i32, i32* %11
%69 = sext i32 %68 to i64
%70 = getelementptr inbounds float, float* %67, i64 %69
%71 = load float, float* %70
%72 = fmul float %71, %66
store float %72, float* %70
br label %73
73: ; preds = %65
%74 = load i32, i32* %11
%75 = add nsw i32 %74, 1
store i32 %75, i32* %11
br label %59
76: ; preds = %59
ret void
}

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