float allpole(x, jcount, npoles, past, c)
float x, *past, *c;
long *jcount, npoles;
{
  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
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x4(%rbp)
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
mov    -0x10(%rbp),%rax
mov    (%rax),%rax
mov    %eax,-0x2c(%rbp)
movl   $0x0,-0x30(%rbp)
movslq -0x30(%rbp),%rax
cmp    -0x18(%rbp),%rax
jge    117f <allpole+0x7f>
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    112e <allpole+0x2e>
movss  -0x4(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x10(%rbp),%rcx
mov    (%rcx),%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    -0x18(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x20(%rbp),%rax
mov    -0x10(%rbp),%rcx
mov    (%rcx),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
mov    (%rax),%rax
add    $0x1,%rax
cqto
idivq  -0x18(%rbp)
mov    -0x10(%rbp),%rax
mov    %rdx,(%rax)
movss  -0x4(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<allpole>:
push   %rbx
cvtsd2ss %xmm0,%xmm0
mov    (%rdi),%rax
test   %rsi,%rsi
jle    11d8 <allpole+0xd8>
movslq %eax,%r9
lea    -0x1(%rsi),%rbx
mov    %esi,%r8d
and    $0x3,%r8d
cmp    $0x3,%rbx
jae    1132 <allpole+0x32>
xor    %r10d,%r10d
test   %r8,%r8
jne    11a7 <allpole+0xa7>
jmp    11d8 <allpole+0xd8>
mov    %rsi,%r11
and    $0xfffffffffffffffc,%r11
lea    (%rdx,%r9,4),%rbx
add    $0xc,%rbx
xor    %r10d,%r10d
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movss  (%rcx,%r10,4),%xmm1
movss  0x4(%rcx,%r10,4),%xmm2
mulss  -0xc(%rbx,%r10,4),%xmm1
mulss  -0x8(%rbx,%r10,4),%xmm2
addss  %xmm0,%xmm1
movss  0x8(%rcx,%r10,4),%xmm3
mulss  -0x4(%rbx,%r10,4),%xmm3
addss  %xmm1,%xmm2
movss  0xc(%rcx,%r10,4),%xmm0
mulss  (%rbx,%r10,4),%xmm0
addss  %xmm2,%xmm3
addss  %xmm3,%xmm0
add    $0x4,%r10
cmp    %r10,%r11
jne    1150 <allpole+0x50>
add    %r10,%r9
test   %r8,%r8
je     11d8 <allpole+0xd8>
lea    (%rdx,%r9,4),%r9
lea    (%rcx,%r10,4),%rcx
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movss  (%rcx,%rbx,4),%xmm1
mulss  (%r9,%rbx,4),%xmm1
addss  %xmm1,%xmm0
add    $0x1,%rbx
cmp    %rbx,%r8
jne    11c0 <allpole+0xc0>
lea    (%rdx,%rax,4),%rcx
movss  %xmm0,(%rcx,%rsi,4)
movss  %xmm0,(%rdx,%rax,4)
add    $0x1,%rax
mov    %rax,%rcx
or     %rsi,%rcx
shr    $0x20,%rcx
je     11fd <allpole+0xfd>
cqto
idiv   %rsi
jmp    1201 <allpole+0x101>
xor    %edx,%edx
div    %esi
mov    %rdx,(%rdi)
pop    %rbx
ret

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

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