double fir_shuffle(double input, int ntaps, double h[], double z[]) {
  int ii;
  double accum;
  z[0] = input;
  accum = h[ntaps - 1] * z[ntaps - 1];
  for (ii = ntaps - 2; ii >= 0; ii--) {
    accum += h[ii] * z[ii];
    z[ii + 1] = z[ii];
  }
  return accum;
}


x86 assembly (clang, -O0):
<fir_shuffle>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x8(%rbp)
mov    %edi,-0xc(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
movsd  -0x8(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x18(%rbp),%rax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0xc(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x0,-0x24(%rbp)
jl     11b6 <fir_shuffle+0xb6>
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x24(%rbp)
jmp    1153 <fir_shuffle+0x53>
movsd  -0x30(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<fir_shuffle>:
movsd  %xmm0,(%rdx)
movslq %edi,%rax
movsd  -0x8(%rsi,%rax,8),%xmm1
mulsd  -0x8(%rdx,%rax,8),%xmm1
cmp    $0x2,%eax
jl     1148 <fir_shuffle+0x48>
add    $0xfffffffe,%edi
test   $0x1,%dil
jne    114d <fir_shuffle+0x4d>
movsd  (%rdx,%rdi,8),%xmm2
movsd  (%rsi,%rdi,8),%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm2,0x8(%rdx,%rdi,8)
lea    -0x1(%rdi),%rax
movapd %xmm0,%xmm1
test   %rdi,%rdi
jne    1155 <fir_shuffle+0x55>
jmp    119b <fir_shuffle+0x9b>
movapd %xmm1,%xmm0
ret
mov    %rdi,%rax
test   %rdi,%rdi
je     119b <fir_shuffle+0x9b>
movapd %xmm1,%xmm0
nopl   0x0(%rax)
movsd  -0x8(%rdx,%rax,8),%xmm1
movsd  (%rdx,%rax,8),%xmm2
movsd  (%rsi,%rax,8),%xmm3
mulsd  %xmm2,%xmm3
addsd  %xmm0,%xmm3
movsd  %xmm2,0x8(%rdx,%rax,8)
movsd  -0x8(%rsi,%rax,8),%xmm0
mulsd  %xmm1,%xmm0
addsd  %xmm3,%xmm0
movsd  %xmm1,(%rdx,%rax,8)
add    $0xfffffffffffffffe,%rax
cmp    $0xffffffffffffffff,%rax
jne    1160 <fir_shuffle+0x60>
ret

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

LLVM IR (-O2):
define double @fir_shuffle(double %0, i32 %1, double* nocapture %2, double* nocapture %3) {
store double %0, double* %3
%5 = add nsw i32 %1, -1
%6 = sext i32 %5 to i64
%7 = getelementptr inbounds double, double* %2, i64 %6
%8 = load double, double* %7
%9 = getelementptr inbounds double, double* %3, i64 %6
%10 = load double, double* %9
%11 = fmul double %8, %10
%12 = icmp sgt i32 %1, 1
br i1 %12, label %13, label %51
13: ; preds = %4
%14 = add nsw i32 %1, -2
%15 = zext i32 %14 to i64
%16 = and i64 %15, 1
%17 = icmp eq i64 %16, 0
br i1 %17, label %18, label %27
18: ; preds = %13
%19 = getelementptr inbounds double, double* %2, i64 %15
%20 = load double, double* %19
%21 = getelementptr inbounds double, double* %3, i64 %15
%22 = load double, double* %21
%23 = tail call double @llvm.fmuladd.f64(double %20, double %22, double %11)
%24 = add nuw nsw i64 %15, 1
%25 = getelementptr inbounds double, double* %3, i64 %24
store double %22, double* %25
%26 = add nsw i64 %15, -1
br label %27
27: ; preds = %18, %13
%28 = phi double [ undef, %13 ], [ %23, %18 ]
%29 = phi i64 [ %15, %13 ], [ %26, %18 ]
%30 = phi double [ %11, %13 ], [ %23, %18 ]
%31 = icmp eq i32 %14, 0
br i1 %31, label %51, label %32
32: ; preds = %27, %32
%33 = phi i64 [ %49, %32 ], [ %29, %27 ]
%34 = phi double [ %47, %32 ], [ %30, %27 ]
%35 = getelementptr inbounds double, double* %2, i64 %33
%36 = load double, double* %35
%37 = getelementptr inbounds double, double* %3, i64 %33
%38 = load double, double* %37
%39 = tail call double @llvm.fmuladd.f64(double %36, double %38, double %34)
%40 = add nuw nsw i64 %33, 1
%41 = getelementptr inbounds double, double* %3, i64 %40
store double %38, double* %41
%42 = add nsw i64 %33, -1
%43 = getelementptr inbounds double, double* %2, i64 %42
%44 = load double, double* %43
%45 = getelementptr inbounds double, double* %3, i64 %42
%46 = load double, double* %45
%47 = tail call double @llvm.fmuladd.f64(double %44, double %46, double %39)
%48 = getelementptr inbounds double, double* %3, i64 %33
store double %46, double* %48
%49 = add nsw i64 %33, -2
%50 = icmp eq i64 %42, 0
br i1 %50, label %51, label %32
51: ; preds = %27, %32, %4
%52 = phi double [ %11, %4 ], [ %28, %27 ], [ %47, %32 ]
ret double %52
}