void train(short *t, short *w, int n, int err) {
  n = (n + 7) & -8;
  for (int i = 0; i < n; ++i) {
    int wt = w[i] + ((t[i] * err * 2 >> 16) + 1 >> 1);
    if (wt < -32768)
      wt = -32768;
    if (wt > 32767)
      wt = 32767;
    w[i] = wt;
  }
}


x86 assembly (clang, -O0):
<train>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    -0x14(%rbp),%eax
add    $0x7,%eax
and    $0xfffffff8,%eax
mov    %eax,-0x14(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11a6 <train+0xa6>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movswl (%rax,%rcx,2),%eax
mov    -0x8(%rbp),%rcx
movslq -0x1c(%rbp),%rdx
movswl (%rcx,%rdx,2),%ecx
imul   -0x18(%rbp),%ecx
shl    $0x1,%ecx
sar    $0x10,%ecx
add    $0x1,%ecx
sar    $0x1,%ecx
add    %ecx,%eax
mov    %eax,-0x20(%rbp)
cmpl   $0xffff8000,-0x20(%rbp)
jge    1172 <train+0x72>
movl   $0xffff8000,-0x20(%rbp)
cmpl   $0x7fff,-0x20(%rbp)
jle    1186 <train+0x86>
movl   $0x7fff,-0x20(%rbp)
mov    -0x20(%rbp),%eax
mov    %ax,%dx
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mov    %dx,(%rax,%rcx,2)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1125 <train+0x25>
pop    %rbp
ret

x86 assembly (clang, -O2):
<train>:
add    $0x7,%edx
and    $0xfffffff8,%edx
jle    1236 <train+0x136>
mov    %edx,%r8d
lea    (%rdi,%r8,2),%rax
cmp    %rsi,%rax
jbe    1172 <train+0x72>
lea    (%rsi,%r8,2),%rax
cmp    %rdi,%rax
jbe    1172 <train+0x72>
xor    %edx,%edx
jmp    1141 <train+0x41>
cs nopw 0x0(%rax,%rax,1)
nop
mov    %ax,(%rsi,%rdx,2)
add    $0x1,%rdx
cmp    %rdx,%r8
je     1236 <train+0x136>
movswl (%rsi,%rdx,2),%r9d
movswl (%rdi,%rdx,2),%eax
imul   %ecx,%eax
sar    $0xf,%eax
add    $0x1,%eax
sar    %eax
add    %r9d,%eax
cmp    $0xffff8001,%eax
jge    1164 <train+0x64>
mov    $0xffff8000,%eax
cmp    $0x7fff,%eax
jl     1130 <train+0x30>
mov    $0x7fff,%eax
jmp    1130 <train+0x30>
movd   %ecx,%xmm0
pshufd $0x0,%xmm0,%xmm0
xor    %ecx,%ecx
pshufd $0xf5,%xmm0,%xmm1
pcmpeqd %xmm2,%xmm2
cs nopw 0x0(%rax,%rax,1)
movdqu (%rsi,%rcx,2),%xmm3
punpcklwd %xmm3,%xmm4
psrad  $0x10,%xmm4
punpckhwd %xmm3,%xmm3
psrad  $0x10,%xmm3
movdqu (%rdi,%rcx,2),%xmm5
punpcklwd %xmm5,%xmm6
psrad  $0x10,%xmm6
punpckhwd %xmm5,%xmm7
psrad  $0x10,%xmm7
movdqa %xmm0,%xmm5
pmuludq %xmm7,%xmm5
pshufd $0xe8,%xmm5,%xmm5
pshufd $0xf5,%xmm7,%xmm7
pmuludq %xmm1,%xmm7
pshufd $0xe8,%xmm7,%xmm7
punpckldq %xmm7,%xmm5
movdqa %xmm0,%xmm7
pmuludq %xmm6,%xmm7
pshufd $0xe8,%xmm7,%xmm7
pshufd $0xf5,%xmm6,%xmm6
pmuludq %xmm1,%xmm6
pshufd $0xe8,%xmm6,%xmm6
punpckldq %xmm6,%xmm7
psrad  $0xf,%xmm7
psrad  $0xf,%xmm5
psubd  %xmm2,%xmm5
psubd  %xmm2,%xmm7
psrad  $0x1,%xmm7
paddd  %xmm4,%xmm7
psrad  $0x1,%xmm5
paddd  %xmm3,%xmm5
packssdw %xmm5,%xmm7
movdqu %xmm7,(%rsi,%rcx,2)
add    $0x8,%rcx
cmp    %rcx,%r8
jne    1190 <train+0x90>
ret

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

LLVM IR (-O2):
define void @train(i16* nocapture %0, i16* nocapture %1, i32 %2, i32 %3) {
%5 = add nsw i32 %2, 7
%6 = and i32 %5, -8
%7 = icmp sgt i32 %6, 0
br i1 %7, label %8, label %41
8: ; preds = %4
%9 = zext i32 %6 to i64
%10 = getelementptr i16, i16* %1, i64 %9
%11 = getelementptr i16, i16* %0, i64 %9
%12 = icmp ugt i16* %11, %1
%13 = icmp ugt i16* %10, %0
%14 = and i1 %12, %13
br i1 %14, label %42, label %15
15: ; preds = %8
%16 = insertelement <8 x i32> poison, i32 %3, i64 0
%17 = shufflevector <8 x i32> %16, <8 x i32> poison, <8 x i32> zeroinitializer
br label %18
18: ; preds = %18, %15
%19 = phi i64 [ 0, %15 ], [ %39, %18 ]
%20 = getelementptr inbounds i16, i16* %1, i64 %19
%21 = bitcast i16* %20 to <8 x i16>*
%22 = load <8 x i16>, <8 x i16>* %21
%23 = sext <8 x i16> %22 to <8 x i32>
%24 = getelementptr inbounds i16, i16* %0, i64 %19
%25 = bitcast i16* %24 to <8 x i16>*
%26 = load <8 x i16>, <8 x i16>* %25
%27 = sext <8 x i16> %26 to <8 x i32>
%28 = mul nsw <8 x i32> %17, %27
%29 = ashr <8 x i32> %28, <i32 15, i32 15, i32 15, i32 15, i32 15, i32 15, i32 15, i32 15>
%30 = add nsw <8 x i32> %29, <i32 1, i32 1, i32 1, i32 1, i32 1, i32 1, i32 1, i32 1>
%31 = ashr <8 x i32> %30, <i32 1, i32 1, i32 1, i32 1, i32 1, i32 1, i32 1, i32 1>
%32 = add nsw <8 x i32> %31, %23
%33 = icmp sgt <8 x i32> %32, <i32 -32768, i32 -32768, i32 -32768, i32 -32768, i32 -32768, i32 -32768, i32 -32768, i32 -32768>
%34 = select <8 x i1> %33, <8 x i32> %32, <8 x i32> <i32 -32768, i32 -32768, i32 -32768, i32 -32768, i32 -32768, i32 -32768, i32 -32768, i32 -32768>
%35 = icmp slt <8 x i32> %34, <i32 32767, i32 32767, i32 32767, i32 32767, i32 32767, i32 32767, i32 32767, i32 32767>
%36 = select <8 x i1> %35, <8 x i32> %34, <8 x i32> <i32 32767, i32 32767, i32 32767, i32 32767, i32 32767, i32 32767, i32 32767, i32 32767>
%37 = trunc <8 x i32> %36 to <8 x i16>
%38 = bitcast i16* %20 to <8 x i16>*
store <8 x i16> %37, <8 x i16>* %38
%39 = add nuw i64 %19, 8
%40 = icmp eq i64 %39, %9
br i1 %40, label %41, label %18
41: ; preds = %18, %42, %4
ret void
42: ; preds = %8, %42
%43 = phi i64 [ %60, %42 ], [ 0, %8 ]
%44 = getelementptr inbounds i16, i16* %1, i64 %43
%45 = load i16, i16* %44
%46 = sext i16 %45 to i32
%47 = getelementptr inbounds i16, i16* %0, i64 %43
%48 = load i16, i16* %47
%49 = sext i16 %48 to i32
%50 = mul nsw i32 %49, %3
%51 = ashr i32 %50, 15
%52 = add nsw i32 %51, 1
%53 = ashr i32 %52, 1
%54 = add nsw i32 %53, %46
%55 = icmp sgt i32 %54, -32768
%56 = select i1 %55, i32 %54, i32 -32768
%57 = icmp slt i32 %56, 32767
%58 = select i1 %57, i32 %56, i32 32767
%59 = trunc i32 %58 to i16
store i16 %59, i16* %44
%60 = add nuw nsw i64 %43, 1
%61 = icmp eq i64 %60, %9
br i1 %61, label %41, label %42
}