void UPDATE(int *taps, const short *history, const int nsuppr,
            const int ntaps) {
  int i;
  int correction;
  for (i = 0; i < ntaps; i++) {
    correction = history[i] * nsuppr;
    taps[i] += correction;
  }
}


x86 assembly (clang, -O0):
<UPDATE>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    1157 <UPDATE+0x57>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movswl (%rax,%rcx,2),%eax
imul   -0x14(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%edx
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
add    (%rax,%rcx,4),%edx
mov    %edx,(%rax,%rcx,4)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1119 <UPDATE+0x19>
pop    %rbp
ret

x86 assembly (clang, -O2):
<UPDATE>:
test   %ecx,%ecx
jle    12f3 <UPDATE+0x1f3>
mov    %ecx,%r8d
cmp    $0x8,%ecx
jae    1117 <UPDATE+0x17>
xor    %ecx,%ecx
jmp    12e0 <UPDATE+0x1e0>
mov    %r8d,%ecx
and    $0xfffffff8,%ecx
movd   %edx,%xmm0
pshufd $0x0,%xmm0,%xmm0
lea    -0x8(%rcx),%rax
mov    %rax,%r9
shr    $0x3,%r9
add    $0x1,%r9
pshufd $0xf5,%xmm0,%xmm1
test   %rax,%rax
je     12f4 <UPDATE+0x1f4>
mov    %r9,%r10
and    $0xfffffffffffffffe,%r10
xor    %eax,%eax
nopl   0x0(%rax)
movq   (%rsi,%rax,2),%xmm2
punpcklwd %xmm2,%xmm2
movq   0x8(%rsi,%rax,2),%xmm3
punpcklwd %xmm3,%xmm3
psrad  $0x10,%xmm2
psrad  $0x10,%xmm3
movdqa %xmm0,%xmm4
pmuludq %xmm2,%xmm4
pshufd $0xe8,%xmm4,%xmm4
pshufd $0xf5,%xmm2,%xmm2
pmuludq %xmm1,%xmm2
pshufd $0xe8,%xmm2,%xmm2
punpckldq %xmm2,%xmm4
movdqa %xmm0,%xmm2
pmuludq %xmm3,%xmm2
pshufd $0xe8,%xmm2,%xmm2
pshufd $0xf5,%xmm3,%xmm3
pmuludq %xmm1,%xmm3
pshufd $0xe8,%xmm3,%xmm3
punpckldq %xmm3,%xmm2
movdqu (%rdi,%rax,4),%xmm3
paddd  %xmm4,%xmm3
movdqu 0x10(%rdi,%rax,4),%xmm4
paddd  %xmm2,%xmm4
movdqu 0x20(%rdi,%rax,4),%xmm2
movdqu 0x30(%rdi,%rax,4),%xmm5
movdqu %xmm3,(%rdi,%rax,4)
movdqu %xmm4,0x10(%rdi,%rax,4)
movq   0x10(%rsi,%rax,2),%xmm3
punpcklwd %xmm3,%xmm3
movq   0x18(%rsi,%rax,2),%xmm4
punpcklwd %xmm4,%xmm4
psrad  $0x10,%xmm3
psrad  $0x10,%xmm4
movdqa %xmm0,%xmm6
pmuludq %xmm3,%xmm6
pshufd $0xe8,%xmm6,%xmm6
pshufd $0xf5,%xmm3,%xmm3
pmuludq %xmm1,%xmm3
pshufd $0xe8,%xmm3,%xmm3
punpckldq %xmm3,%xmm6
paddd  %xmm2,%xmm6
movdqa %xmm0,%xmm2
pmuludq %xmm4,%xmm2
pshufd $0xe8,%xmm2,%xmm2
pshufd $0xf5,%xmm4,%xmm3
pmuludq %xmm1,%xmm3
pshufd $0xe8,%xmm3,%xmm3
punpckldq %xmm3,%xmm2
paddd  %xmm5,%xmm2
movdqu %xmm6,0x20(%rdi,%rax,4)
movdqu %xmm2,0x30(%rdi,%rax,4)
add    $0x10,%rax
add    $0xfffffffffffffffe,%r10
jne    1150 <UPDATE+0x50>
test   $0x1,%r9b
je     12ce <UPDATE+0x1ce>
movq   (%rsi,%rax,2),%xmm2
movq   0x8(%rsi,%rax,2),%xmm3
punpcklwd %xmm2,%xmm2
psrad  $0x10,%xmm2
punpcklwd %xmm3,%xmm3
psrad  $0x10,%xmm3
movdqa %xmm0,%xmm4
pmuludq %xmm2,%xmm4
pshufd $0xe8,%xmm4,%xmm4
pshufd $0xf5,%xmm2,%xmm2
pmuludq %xmm1,%xmm2
pshufd $0xe8,%xmm2,%xmm2
punpckldq %xmm2,%xmm4
pmuludq %xmm3,%xmm0
pshufd $0xe8,%xmm0,%xmm0
pshufd $0xf5,%xmm3,%xmm2
pmuludq %xmm1,%xmm2
pshufd $0xe8,%xmm2,%xmm1
punpckldq %xmm1,%xmm0
movdqu (%rdi,%rax,4),%xmm1
paddd  %xmm4,%xmm1
movdqu 0x10(%rdi,%rax,4),%xmm2
paddd  %xmm0,%xmm2
movdqu %xmm1,(%rdi,%rax,4)
movdqu %xmm2,0x10(%rdi,%rax,4)
cmp    %r8,%rcx
je     12f3 <UPDATE+0x1f3>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movswl (%rsi,%rcx,2),%eax
imul   %edx,%eax
add    %eax,(%rdi,%rcx,4)
add    $0x1,%rcx
cmp    %rcx,%r8
jne    12e0 <UPDATE+0x1e0>
ret
xor    %eax,%eax
test   $0x1,%r9b
jne    1259 <UPDATE+0x159>
jmp    12ce <UPDATE+0x1ce>

LLVM IR (-O0):
define void @UPDATE(i32* %0, i16* %1, i32 %2, i32 %3) {
%5 = alloca i32*
%6 = alloca i16*
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
store i32* %0, i32** %5
store i16* %1, i16** %6
store i32 %2, i32* %7
store i32 %3, i32* %8
store i32 0, i32* %9
br label %11
11: ; preds = %31, %4
%12 = load i32, i32* %9
%13 = load i32, i32* %8
%14 = icmp slt i32 %12, %13
br i1 %14, label %15, label %34
15: ; preds = %11
%16 = load i16*, i16** %6
%17 = load i32, i32* %9
%18 = sext i32 %17 to i64
%19 = getelementptr inbounds i16, i16* %16, i64 %18
%20 = load i16, i16* %19
%21 = sext i16 %20 to i32
%22 = load i32, i32* %7
%23 = mul nsw i32 %21, %22
store i32 %23, i32* %10
%24 = load i32, i32* %10
%25 = load i32*, i32** %5
%26 = load i32, i32* %9
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds i32, i32* %25, i64 %27
%29 = load i32, i32* %28
%30 = add nsw i32 %29, %24
store i32 %30, i32* %28
br label %31
31: ; preds = %15
%32 = load i32, i32* %9
%33 = add nsw i32 %32, 1
store i32 %33, i32* %9
br label %11
34: ; preds = %11
ret void
}

LLVM IR (-O2):
define void @UPDATE(i32* nocapture %0, i16* nocapture %1, i32 %2, i32 %3) {
%5 = icmp sgt i32 %3, 0
br i1 %5, label %6, label %108
6: ; preds = %4
%7 = zext i32 %3 to i64
%8 = icmp ult i32 %3, 8
br i1 %8, label %95, label %9
9: ; preds = %6
%10 = and i64 %7, 4294967288
%11 = insertelement <4 x i32> poison, i32 %2, i64 0
%12 = shufflevector <4 x i32> %11, <4 x i32> poison, <4 x i32> zeroinitializer
%13 = insertelement <4 x i32> poison, i32 %2, i64 0
%14 = shufflevector <4 x i32> %13, <4 x i32> poison, <4 x i32> zeroinitializer
%15 = add nsw i64 %10, -8
%16 = lshr exact i64 %15, 3
%17 = add nuw nsw i64 %16, 1
%18 = and i64 %17, 1
%19 = icmp eq i64 %15, 0
br i1 %19, label %69, label %20
20: ; preds = %9
%21 = and i64 %17, 4611686018427387902
br label %22
22: ; preds = %22, %20
%23 = phi i64 [ 0, %20 ], [ %66, %22 ]
%24 = phi i64 [ 0, %20 ], [ %67, %22 ]
%25 = getelementptr inbounds i16, i16* %1, i64 %23
%26 = bitcast i16* %25 to <4 x i16>*
%27 = load <4 x i16>, <4 x i16>* %26
%28 = getelementptr inbounds i16, i16* %25, i64 4
%29 = bitcast i16* %28 to <4 x i16>*
%30 = load <4 x i16>, <4 x i16>* %29
%31 = sext <4 x i16> %27 to <4 x i32>
%32 = sext <4 x i16> %30 to <4 x i32>
%33 = mul nsw <4 x i32> %12, %31
%34 = mul nsw <4 x i32> %14, %32
%35 = getelementptr inbounds i32, i32* %0, i64 %23
%36 = bitcast i32* %35 to <4 x i32>*
%37 = load <4 x i32>, <4 x i32>* %36
%38 = getelementptr inbounds i32, i32* %35, i64 4
%39 = bitcast i32* %38 to <4 x i32>*
%40 = load <4 x i32>, <4 x i32>* %39
%41 = add nsw <4 x i32> %33, %37
%42 = add nsw <4 x i32> %34, %40
%43 = bitcast i32* %35 to <4 x i32>*
store <4 x i32> %41, <4 x i32>* %43
%44 = bitcast i32* %38 to <4 x i32>*
store <4 x i32> %42, <4 x i32>* %44
%45 = or i64 %23, 8
%46 = getelementptr inbounds i16, i16* %1, i64 %45
%47 = bitcast i16* %46 to <4 x i16>*
%48 = load <4 x i16>, <4 x i16>* %47
%49 = getelementptr inbounds i16, i16* %46, i64 4
%50 = bitcast i16* %49 to <4 x i16>*
%51 = load <4 x i16>, <4 x i16>* %50
%52 = sext <4 x i16> %48 to <4 x i32>
%53 = sext <4 x i16> %51 to <4 x i32>
%54 = mul nsw <4 x i32> %12, %52
%55 = mul nsw <4 x i32> %14, %53
%56 = getelementptr inbounds i32, i32* %0, i64 %45
%57 = bitcast i32* %56 to <4 x i32>*
%58 = load <4 x i32>, <4 x i32>* %57
%59 = getelementptr inbounds i32, i32* %56, i64 4
%60 = bitcast i32* %59 to <4 x i32>*
%61 = load <4 x i32>, <4 x i32>* %60
%62 = add nsw <4 x i32> %54, %58
%63 = add nsw <4 x i32> %55, %61
%64 = bitcast i32* %56 to <4 x i32>*
store <4 x i32> %62, <4 x i32>* %64
%65 = bitcast i32* %59 to <4 x i32>*
store <4 x i32> %63, <4 x i32>* %65
%66 = add nuw i64 %23, 16
%67 = add i64 %24, 2
%68 = icmp eq i64 %67, %21
br i1 %68, label %69, label %22
69: ; preds = %22, %9
%70 = phi i64 [ 0, %9 ], [ %66, %22 ]
%71 = icmp eq i64 %18, 0
br i1 %71, label %93, label %72
72: ; preds = %69
%73 = getelementptr inbounds i16, i16* %1, i64 %70
%74 = bitcast i16* %73 to <4 x i16>*
%75 = load <4 x i16>, <4 x i16>* %74
%76 = getelementptr inbounds i16, i16* %73, i64 4
%77 = bitcast i16* %76 to <4 x i16>*
%78 = load <4 x i16>, <4 x i16>* %77
%79 = sext <4 x i16> %75 to <4 x i32>
%80 = sext <4 x i16> %78 to <4 x i32>
%81 = mul nsw <4 x i32> %12, %79
%82 = mul nsw <4 x i32> %14, %80
%83 = getelementptr inbounds i32, i32* %0, i64 %70
%84 = bitcast i32* %83 to <4 x i32>*
%85 = load <4 x i32>, <4 x i32>* %84
%86 = getelementptr inbounds i32, i32* %83, i64 4
%87 = bitcast i32* %86 to <4 x i32>*
%88 = load <4 x i32>, <4 x i32>* %87
%89 = add nsw <4 x i32> %81, %85
%90 = add nsw <4 x i32> %82, %88
%91 = bitcast i32* %83 to <4 x i32>*
store <4 x i32> %89, <4 x i32>* %91
%92 = bitcast i32* %86 to <4 x i32>*
store <4 x i32> %90, <4 x i32>* %92
br label %93
93: ; preds = %69, %72
%94 = icmp eq i64 %10, %7
br i1 %94, label %108, label %95
95: ; preds = %6, %93
%96 = phi i64 [ 0, %6 ], [ %10, %93 ]
br label %97
97: ; preds = %95, %97
%98 = phi i64 [ %106, %97 ], [ %96, %95 ]
%99 = getelementptr inbounds i16, i16* %1, i64 %98
%100 = load i16, i16* %99
%101 = sext i16 %100 to i32
%102 = mul nsw i32 %101, %2
%103 = getelementptr inbounds i32, i32* %0, i64 %98
%104 = load i32, i32* %103
%105 = add nsw i32 %102, %104
store i32 %105, i32* %103
%106 = add nuw nsw i64 %98, 1
%107 = icmp eq i64 %106, %7
br i1 %107, label %108, label %97
108: ; preds = %97, %93, %4
ret void
}