float allpole(float x, int *jcount, int npoles, float *past, float *c) {
  register int j, nfint;
  for (j = *jcount, nfint = 0; nfint < npoles; nfint++, j++)
    x += (*(c + nfint) * *(past + j));
  *(past + *jcount) = *(past + *jcount + npoles) = x;
  *jcount = (*jcount + 1) % npoles;
  return (x);
}


x86 assembly (clang, -O0):
<allpole>:
push   %rbp
mov    %rsp,%rbp
movss  %xmm0,-0x4(%rbp)
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
mov    -0x10(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x2c(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1177 <allpole+0x77>
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x4(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x4(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1128 <allpole+0x28>
movss  -0x4(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x10(%rbp),%rcx
movslq (%rcx),%rcx
shl    $0x2,%rcx
add    %rcx,%rax
movslq -0x14(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x20(%rbp),%rax
mov    -0x10(%rbp),%rcx
movslq (%rcx),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
mov    (%rax),%eax
add    $0x1,%eax
cltd
idivl  -0x14(%rbp)
mov    -0x10(%rbp),%rax
mov    %edx,(%rax)
movss  -0x4(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<allpole>:
movslq (%rdi),%r8
test   %esi,%esi
jle    11b8 <allpole+0xb8>
mov    %esi,%r10d
lea    -0x1(%r10),%rax
mov    %r10d,%r9d
and    $0x3,%r9d
cmp    $0x3,%rax
jae    112e <allpole+0x2e>
xor    %eax,%eax
mov    %r8,%r10
test   %r9,%r9
jne    1194 <allpole+0x94>
jmp    11b8 <allpole+0xb8>
and    $0xfffffffc,%r10d
lea    (%rdx,%r8,4),%r11
add    $0xc,%r11
xor    %eax,%eax
nopl   0x0(%rax)
movss  (%rcx,%rax,4),%xmm1
movss  0x4(%rcx,%rax,4),%xmm2
mulss  -0xc(%r11,%rax,4),%xmm1
mulss  -0x8(%r11,%rax,4),%xmm2
addss  %xmm0,%xmm1
movss  0x8(%rcx,%rax,4),%xmm3
mulss  -0x4(%r11,%rax,4),%xmm3
addss  %xmm1,%xmm2
movss  0xc(%rcx,%rax,4),%xmm0
mulss  (%r11,%rax,4),%xmm0
addss  %xmm2,%xmm3
addss  %xmm3,%xmm0
add    $0x4,%rax
cmp    %rax,%r10
jne    1140 <allpole+0x40>
lea    (%r8,%rax,1),%r10
test   %r9,%r9
je     11b8 <allpole+0xb8>
lea    (%rdx,%r10,4),%r10
lea    (%rcx,%rax,4),%rax
xor    %ecx,%ecx
xchg   %ax,%ax
movss  (%rax,%rcx,4),%xmm1
mulss  (%r10,%rcx,4),%xmm1
addss  %xmm1,%xmm0
add    $0x1,%rcx
cmp    %rcx,%r9
jne    11a0 <allpole+0xa0>
lea    (%rdx,%r8,4),%rax
movslq %esi,%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  %xmm0,(%rdx,%r8,4)
lea    0x1(%r8),%eax
cltd
idiv   %esi
mov    %edx,(%rdi)
ret

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

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