#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 */
void dddsp_filter_iir_fwd_float(float out_samples[], const float in_samples[],
                                size_t count, const float a[], const float b[],
                                int poles) {
  if (count < poles)
    return;
  size_t step = 0;
  int i;
  for (; count; step++, count--, out_samples++, in_samples++) {
    *out_samples = 0;
    for (i = 0; i <= poles; i++) {
      if (step < i)
        continue;
      *out_samples += a[i] * in_samples[-i];
      *out_samples += b[i] * out_samples[-i];
    }
  }
}


x86 assembly (clang, -O0):
<dddsp_filter_iir_fwd_float>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9d,-0x2c(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
cmp    %rcx,%rax
jae    1132 <dddsp_filter_iir_fwd_float+0x32>
jmp    1225 <dddsp_filter_iir_fwd_float+0x125>
movq   $0x0,-0x38(%rbp)
cmpq   $0x0,-0x18(%rbp)
je     1225 <dddsp_filter_iir_fwd_float+0x125>
mov    -0x8(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jg     11eb <dddsp_filter_iir_fwd_float+0xeb>
mov    -0x38(%rbp),%rax
movslq -0x3c(%rbp),%rcx
cmp    %rcx,%rax
jae    1179 <dddsp_filter_iir_fwd_float+0x79>
jmp    11dd <dddsp_filter_iir_fwd_float+0xdd>
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
xor    %ecx,%ecx
sub    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
xor    %ecx,%ecx
sub    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1157 <dddsp_filter_iir_fwd_float+0x57>
jmp    11f0 <dddsp_filter_iir_fwd_float+0xf0>
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
mov    -0x18(%rbp),%rax
add    $0xffffffffffffffff,%rax
mov    %rax,-0x18(%rbp)
mov    -0x8(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x8(%rbp)
mov    -0x10(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x10(%rbp)
jmp    113a <dddsp_filter_iir_fwd_float+0x3a>
pop    %rbp
ret

x86 assembly (clang, -O2):
<dddsp_filter_iir_fwd_float>:
push   %r15
push   %r14
push   %rbx
movslq %r9d,%rax
cmp    %rdx,%rax
ja     119b <dddsp_filter_iir_fwd_float+0x9b>
test   %rdx,%rdx
je     119b <dddsp_filter_iir_fwd_float+0x9b>
lea    0x1(%r9),%r11d
xor    %r10d,%r10d
xor    %r14d,%r14d
mov    %rdi,%r15
jmp    1142 <dddsp_filter_iir_fwd_float+0x42>
nopl   0x0(%rax)
add    $0x1,%r14
add    $0x4,%r15
add    $0x4,%r10
add    $0xffffffffffffffff,%rdx
je     119b <dddsp_filter_iir_fwd_float+0x9b>
movl   $0x0,(%r15)
test   %r9d,%r9d
js     1130 <dddsp_filter_iir_fwd_float+0x30>
xorps  %xmm0,%xmm0
mov    %r10,%rax
xor    %ebx,%ebx
jmp    116d <dddsp_filter_iir_fwd_float+0x6d>
nopl   0x0(%rax,%rax,1)
add    $0x1,%rbx
add    $0xfffffffffffffffc,%rax
cmp    %rbx,%r11
je     1130 <dddsp_filter_iir_fwd_float+0x30>
cmp    %rbx,%r14
jb     1160 <dddsp_filter_iir_fwd_float+0x60>
movss  (%rcx,%rbx,4),%xmm1
mulss  (%rsi,%rax,1),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,(%r15)
movss  (%r8,%rbx,4),%xmm0
mulss  (%rdi,%rax,1),%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%r15)
jmp    1160 <dddsp_filter_iir_fwd_float+0x60>
pop    %rbx
pop    %r14
pop    %r15
ret

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

LLVM IR (-O2):
define void @dddsp_filter_iir_fwd_float(float* nocapture %0, float* nocapture %1, i64 %2, float* nocapture %3, float* nocapture %4, i32 %5) {
%7 = sext i32 %5 to i64
%8 = icmp ugt i64 %7, %2
%9 = icmp eq i64 %2, 0
%10 = or i1 %8, %9
br i1 %10, label %46, label %11
11: ; preds = %6
%12 = icmp slt i32 %5, 0
%13 = add i32 %5, 1
%14 = zext i32 %13 to i64
br label %15
15: ; preds = %11, %40
%16 = phi i64 [ 0, %11 ], [ %41, %40 ]
%17 = phi float* [ %0, %11 ], [ %43, %40 ]
%18 = phi float* [ %1, %11 ], [ %44, %40 ]
%19 = phi i64 [ %2, %11 ], [ %42, %40 ]
store float 0.000000e+00, float* %17
br i1 %12, label %40, label %20
20: ; preds = %15, %36
%21 = phi float [ %37, %36 ], [ 0.000000e+00, %15 ]
%22 = phi i64 [ %38, %36 ], [ 0, %15 ]
%23 = icmp ult i64 %16, %22
br i1 %23, label %36, label %24
24: ; preds = %20
%25 = getelementptr inbounds float, float* %3, i64 %22
%26 = load float, float* %25
%27 = sub nsw i64 0, %22
%28 = getelementptr inbounds float, float* %18, i64 %27
%29 = load float, float* %28
%30 = tail call float @llvm.fmuladd.f32(float %26, float %29, float %21)
store float %30, float* %17
%31 = getelementptr inbounds float, float* %4, i64 %22
%32 = load float, float* %31
%33 = getelementptr inbounds float, float* %17, i64 %27
%34 = load float, float* %33
%35 = tail call float @llvm.fmuladd.f32(float %32, float %34, float %30)
store float %35, float* %17
br label %36
36: ; preds = %20, %24
%37 = phi float [ %21, %20 ], [ %35, %24 ]
%38 = add nuw nsw i64 %22, 1
%39 = icmp eq i64 %38, %14
br i1 %39, label %40, label %20
40: ; preds = %36, %15
%41 = add nuw i64 %16, 1
%42 = add i64 %19, -1
%43 = getelementptr inbounds float, float* %17, i64 1
%44 = getelementptr inbounds float, float* %18, i64 1
%45 = icmp eq i64 %42, 0
br i1 %45, label %46, label %15
46: ; preds = %40, %6
ret void
}