void UPDATE2(int *taps, short *taps_short, 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;
    taps_short[i] = taps[i] >> 16;
  }
}


x86 assembly (clang, -O0):
<UPDATE2>:
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)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    1179 <UPDATE2+0x79>
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movswl (%rax,%rcx,2),%eax
imul   -0x1c(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%edx
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
add    (%rax,%rcx,4),%edx
mov    %edx,(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
sar    $0x10,%eax
mov    %ax,%dx
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    %dx,(%rax,%rcx,2)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    111e <UPDATE2+0x1e>
pop    %rbp
ret

x86 assembly (clang, -O2):
<UPDATE2>:
push   %rbx
test   %r8d,%r8d
jle    12d4 <UPDATE2+0x1d4>
mov    %r8d,%r9d
cmp    $0x4,%r8d
jb     1125 <UPDATE2+0x25>
lea    (%rdx,%r9,2),%rax
cmp    %rsi,%rax
jbe    119a <UPDATE2+0x9a>
lea    (%rsi,%r9,2),%rax
cmp    %rdx,%rax
jbe    119a <UPDATE2+0x9a>
xor    %ebx,%ebx
mov    %rbx,%r8
not    %r8
test   $0x1,%r9b
je     1151 <UPDATE2+0x51>
movswl (%rdx,%rbx,2),%r10d
imul   %ecx,%r10d
add    (%rdi,%rbx,4),%r10d
mov    %r10d,(%rdi,%rbx,4)
shr    $0x10,%r10d
mov    %r10w,(%rsi,%rbx,2)
or     $0x1,%rbx
add    %r9,%r8
je     12d4 <UPDATE2+0x1d4>
nopw   0x0(%rax,%rax,1)
movswl (%rdx,%rbx,2),%eax
imul   %ecx,%eax
add    (%rdi,%rbx,4),%eax
mov    %eax,(%rdi,%rbx,4)
shr    $0x10,%eax
mov    %ax,(%rsi,%rbx,2)
movswl 0x2(%rdx,%rbx,2),%eax
imul   %ecx,%eax
add    0x4(%rdi,%rbx,4),%eax
mov    %eax,0x4(%rdi,%rbx,4)
shr    $0x10,%eax
mov    %ax,0x2(%rsi,%rbx,2)
add    $0x2,%rbx
cmp    %rbx,%r9
jne    1160 <UPDATE2+0x60>
jmp    12d4 <UPDATE2+0x1d4>
mov    %r9d,%ebx
and    $0xfffffffc,%ebx
movd   %ecx,%xmm0
pshufd $0x0,%xmm0,%xmm0
lea    -0x4(%rbx),%r8
mov    %r8,%r10
shr    $0x2,%r10
add    $0x1,%r10
pshufd $0xf5,%xmm0,%xmm1
test   %r8,%r8
je     12d6 <UPDATE2+0x1d6>
mov    %r10,%r11
and    $0xfffffffffffffffe,%r11
xor    %r8d,%r8d
movq   (%rdx,%r8,2),%xmm2
punpcklwd %xmm2,%xmm2
psrad  $0x10,%xmm2
movdqa %xmm0,%xmm3
pmuludq %xmm2,%xmm3
pshufd $0xe8,%xmm3,%xmm3
pshufd $0xf5,%xmm2,%xmm2
pmuludq %xmm1,%xmm2
pshufd $0xe8,%xmm2,%xmm2
punpckldq %xmm2,%xmm3
movdqu (%rdi,%r8,4),%xmm2
paddd  %xmm3,%xmm2
movdqu %xmm2,(%rdi,%r8,4)
psrad  $0x10,%xmm2
packssdw %xmm2,%xmm2
movq   %xmm2,(%rsi,%r8,2)
movq   0x8(%rdx,%r8,2),%xmm2
punpcklwd %xmm2,%xmm2
psrad  $0x10,%xmm2
movdqa %xmm0,%xmm3
pmuludq %xmm2,%xmm3
pshufd $0xe8,%xmm3,%xmm3
pshufd $0xf5,%xmm2,%xmm2
pmuludq %xmm1,%xmm2
pshufd $0xe8,%xmm2,%xmm2
punpckldq %xmm2,%xmm3
movdqu 0x10(%rdi,%r8,4),%xmm2
paddd  %xmm3,%xmm2
movdqu %xmm2,0x10(%rdi,%r8,4)
psrad  $0x10,%xmm2
packssdw %xmm2,%xmm2
movq   %xmm2,0x8(%rsi,%r8,2)
add    $0x8,%r8
add    $0xfffffffffffffffe,%r11
jne    11d0 <UPDATE2+0xd0>
test   $0x1,%r10b
je     12cb <UPDATE2+0x1cb>
movq   (%rdx,%r8,2),%xmm2
punpcklwd %xmm2,%xmm2
psrad  $0x10,%xmm2
pmuludq %xmm2,%xmm0
pshufd $0xe8,%xmm0,%xmm0
pshufd $0xf5,%xmm2,%xmm2
pmuludq %xmm2,%xmm1
pshufd $0xe8,%xmm1,%xmm1
punpckldq %xmm1,%xmm0
movdqu (%rdi,%r8,4),%xmm1
paddd  %xmm0,%xmm1
movdqu %xmm1,(%rdi,%r8,4)
psrad  $0x10,%xmm1
packssdw %xmm1,%xmm1
movq   %xmm1,(%rsi,%r8,2)
cmp    %r9,%rbx
jne    1127 <UPDATE2+0x27>
pop    %rbx
ret
xor    %r8d,%r8d
test   $0x1,%r10b
jne    1282 <UPDATE2+0x182>
jmp    12cb <UPDATE2+0x1cb>

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

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