void zerflt(float input[], float coeff[], float output[], int order, int npts) {
  int i, j;
  float accum;
  for (i = npts - 1; i >= 0; i--) {
    accum = 0.0;
    for (j = 0; j <= order; j++)
      accum += input[i - j] * coeff[j];
    output[i] = accum;
  }
}


x86 assembly (clang, -O0):
<zerflt>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
mov    -0x20(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x0,-0x24(%rbp)
jl     11a4 <zerflt+0xa4>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x2c(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jg     1184 <zerflt+0x84>
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x2c(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x2c(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1139 <zerflt+0x39>
movss  -0x2c(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x24(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x24(%rbp)
jmp    1120 <zerflt+0x20>
pop    %rbp
ret

x86 assembly (clang, -O2):
<zerflt>:
push   %r15
push   %r14
push   %r12
push   %rbx
test   %r8d,%r8d
jle    121f <zerflt+0x11f>
lea    0x1(%rcx),%r10d
mov    %r8d,%r8d
lea    -0x1(%r10),%r9
mov    %r10d,%r15d
and    $0x3,%r15d
lea    (%rdi,%r8,4),%r12
add    $0xfffffffffffffffc,%r12
lea    0xc(%rsi),%r11
and    $0xfffffffc,%r10d
neg    %r10
jmp    1154 <zerflt+0x54>
nopw   0x0(%rax,%rax,1)
movss  %xmm0,(%rdx,%r8,4)
add    $0xfffffffffffffffc,%r12
cmp    $0x1,%r14
jle    121f <zerflt+0x11f>
mov    %r8,%r14
add    $0xffffffffffffffff,%r8
xorps  %xmm0,%xmm0
test   %ecx,%ecx
js     1140 <zerflt+0x40>
cmp    $0x3,%r9
jae    1180 <zerflt+0x80>
xor    %ebx,%ebx
test   %r15,%r15
jne    11e7 <zerflt+0xe7>
jmp    1140 <zerflt+0x40>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %r11,%rax
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nop
movss  (%r12,%rbx,4),%xmm1
mulss  -0xc(%rax),%xmm1
addss  %xmm0,%xmm1
movss  -0x4(%r12,%rbx,4),%xmm2
mulss  -0x8(%rax),%xmm2
addss  %xmm1,%xmm2
movss  -0x8(%r12,%rbx,4),%xmm1
mulss  -0x4(%rax),%xmm1
movss  -0xc(%r12,%rbx,4),%xmm0
addss  %xmm2,%xmm1
mulss  (%rax),%xmm0
addss  %xmm1,%xmm0
add    $0xfffffffffffffffc,%rbx
add    $0x10,%rax
cmp    %rbx,%r10
jne    1190 <zerflt+0x90>
neg    %rbx
test   %r15,%r15
je     1140 <zerflt+0x40>
lea    (%rsi,%rbx,4),%rax
shl    $0x2,%rbx
mov    %r12,%rdi
sub    %rbx,%rdi
xor    %ebx,%ebx
nopw   0x0(%rax,%rax,1)
movss  (%rdi),%xmm1
mulss  (%rax,%rbx,4),%xmm1
addss  %xmm1,%xmm0
add    $0x1,%rbx
add    $0xfffffffffffffffc,%rdi
cmp    %rbx,%r15
jne    1200 <zerflt+0x100>
jmp    1140 <zerflt+0x40>
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define void @zerflt(float* %0, float* %1, float* %2, i32 %3, i32 %4) {
%6 = alloca float*
%7 = alloca float*
%8 = alloca float*
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca float
store float* %0, float** %6
store float* %1, float** %7
store float* %2, float** %8
store i32 %3, i32* %9
store i32 %4, i32* %10
%14 = load i32, i32* %10
%15 = sub nsw i32 %14, 1
store i32 %15, i32* %11
br label %16
16: ; preds = %48, %5
%17 = load i32, i32* %11
%18 = icmp sge i32 %17, 0
br i1 %18, label %19, label %51
19: ; preds = %16
store float 0.000000e+00, float* %13
store i32 0, i32* %12
br label %20
20: ; preds = %39, %19
%21 = load i32, i32* %12
%22 = load i32, i32* %9
%23 = icmp sle i32 %21, %22
br i1 %23, label %24, label %42
24: ; preds = %20
%25 = load float*, float** %6
%26 = load i32, i32* %11
%27 = load i32, i32* %12
%28 = sub nsw i32 %26, %27
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds float, float* %25, i64 %29
%31 = load float, float* %30
%32 = load float*, float** %7
%33 = load i32, i32* %12
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds float, float* %32, i64 %34
%36 = load float, float* %35
%37 = load float, float* %13
%38 = call float @llvm.fmuladd.f32(float %31, float %36, float %37)
store float %38, float* %13
br label %39
39: ; preds = %24
%40 = load i32, i32* %12
%41 = add nsw i32 %40, 1
store i32 %41, i32* %12
br label %20
42: ; preds = %20
%43 = load float, float* %13
%44 = load float*, float** %8
%45 = load i32, i32* %11
%46 = sext i32 %45 to i64
%47 = getelementptr inbounds float, float* %44, i64 %46
store float %43, float* %47
br label %48
48: ; preds = %42
%49 = load i32, i32* %11
%50 = add nsw i32 %49, -1
store i32 %50, i32* %11
br label %16
51: ; preds = %16
ret void
}

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