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


x86 assembly (clang, -O0):
<calc_coeffs>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
mov    %rdi,-0x8(%rbp)
movss  %xmm0,-0xc(%rbp)
mov    %esi,-0x10(%rbp)
cvtsi2ssl -0x10(%rbp),%xmm0
movss  0xed3(%rip),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x1c(%rbp)
movss  0xebe(%rip),%xmm0
subss  -0xc(%rbp),%xmm0
mov    -0x10(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cvtsi2ss %eax,%xmm1
subss  %xmm1,%xmm0
movss  %xmm0,-0xc(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x20(%rbp)
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    11df <calc_coeffs+0xcf>
movss  -0xc(%rbp),%xmm0
movaps 0xe8a(%rip),%xmm1
pand   %xmm1,%xmm0
divss  -0x1c(%rbp),%xmm0
call   1030 <catrom@plt>
divss  -0x1c(%rbp),%xmm0
movss  %xmm0,-0x18(%rbp)
movss  -0x18(%rbp),%xmm0
addss  -0x20(%rbp),%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0x18(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  0xe39(%rip),%xmm0
addss  -0xc(%rbp),%xmm0
movss  %xmm0,-0xc(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    116e <calc_coeffs+0x5e>
movss  0xe19(%rip),%xmm0
divss  -0x20(%rbp),%xmm0
movss  %xmm0,-0x20(%rbp)
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1229 <calc_coeffs+0x119>
movss  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    11f8 <calc_coeffs+0xe8>
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<calc_coeffs>:
push   %rbp
push   %r15
push   %r14
push   %rbx
sub    $0x28,%rsp
test   %esi,%esi
jle    1277 <calc_coeffs+0x167>
mov    %esi,%r14d
mov    %rdi,%rbx
cvtsi2ss %esi,%xmm1
mulss  0xecc(%rip),%xmm1
movss  %xmm1,0x8(%rsp)
movss  0xec2(%rip),%xmm1
subss  %xmm0,%xmm1
mov    %esi,%eax
shr    %eax
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
subss  %xmm0,%xmm1
mov    %esi,%r15d
xorps  %xmm2,%xmm2
xor    %ebp,%ebp
nopl   (%rax)
movss  %xmm2,0xc(%rsp)
movaps %xmm1,0x10(%rsp)
movaps %xmm1,%xmm0
andps  0xe9b(%rip),%xmm0
divss  0x8(%rsp),%xmm0
call   1030 <catrom@plt>
movss  0xc(%rsp),%xmm2
movaps 0x10(%rsp),%xmm1
movss  0xe71(%rip),%xmm5
divss  0x8(%rsp),%xmm0
addss  %xmm0,%xmm2
movss  %xmm0,(%rbx,%rbp,4)
addss  %xmm5,%xmm1
add    $0x1,%rbp
cmp    %rbp,%r15
jne    1160 <calc_coeffs+0x50>
test   %r14d,%r14d
jle    1277 <calc_coeffs+0x167>
divss  %xmm2,%xmm5
cmp    $0x8,%r14d
jae    11c9 <calc_coeffs+0xb9>
xor    %eax,%eax
jmp    1260 <calc_coeffs+0x150>
mov    %r15d,%eax
and    $0xfffffff8,%eax
movaps %xmm5,%xmm0
shufps $0x0,%xmm5,%xmm0
lea    -0x8(%rax),%rcx
mov    %rcx,%rdx
shr    $0x3,%rdx
add    $0x1,%rdx
test   %rcx,%rcx
je     1282 <calc_coeffs+0x172>
mov    %rdx,%rsi
and    $0xfffffffffffffffe,%rsi
xor    %ecx,%ecx
nopw   0x0(%rax,%rax,1)
movups (%rbx,%rcx,4),%xmm1
movups 0x10(%rbx,%rcx,4),%xmm2
movups 0x20(%rbx,%rcx,4),%xmm3
movups 0x30(%rbx,%rcx,4),%xmm4
mulps  %xmm0,%xmm1
mulps  %xmm0,%xmm2
movups %xmm1,(%rbx,%rcx,4)
movups %xmm2,0x10(%rbx,%rcx,4)
mulps  %xmm0,%xmm3
mulps  %xmm0,%xmm4
movups %xmm3,0x20(%rbx,%rcx,4)
movups %xmm4,0x30(%rbx,%rcx,4)
add    $0x10,%rcx
add    $0xfffffffffffffffe,%rsi
jne    1200 <calc_coeffs+0xf0>
test   $0x1,%dl
je     1259 <calc_coeffs+0x149>
movups (%rbx,%rcx,4),%xmm1
movups 0x10(%rbx,%rcx,4),%xmm2
mulps  %xmm0,%xmm1
mulps  %xmm0,%xmm2
movups %xmm1,(%rbx,%rcx,4)
movups %xmm2,0x10(%rbx,%rcx,4)
cmp    %r15,%rax
je     1277 <calc_coeffs+0x167>
xchg   %ax,%ax
movss  (%rbx,%rax,4),%xmm0
mulss  %xmm5,%xmm0
movss  %xmm0,(%rbx,%rax,4)
add    $0x1,%rax
cmp    %rax,%r15
jne    1260 <calc_coeffs+0x150>
add    $0x28,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %ecx,%ecx
test   $0x1,%dl
jne    1241 <calc_coeffs+0x131>
jmp    1259 <calc_coeffs+0x149>

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

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