short fix_dot(short *hpa, short *pb, int n) {
  short *pa;
  long sum;
  register short a, b;
  sum = 0L;
  while (n--) {
    a = *pa++;
    b = *pb++;
    a = ((long)(a) * (long)(b)) >> 15;
    sum += a;
  }
  if (sum > 0x7FFF)
    sum = 0x7FFF;
  else if (sum < -0x7FFF)
    sum = -0x7FFF;
  return (short)sum;
}


x86 assembly (clang, -O0):
<fix_dot>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
movq   $0x0,-0x28(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,%ecx
add    $0xffffffff,%ecx
mov    %ecx,-0x14(%rbp)
cmp    $0x0,%eax
je     117f <fix_dot+0x7f>
mov    -0x20(%rbp),%rax
mov    %rax,%rcx
add    $0x2,%rcx
mov    %rcx,-0x20(%rbp)
mov    (%rax),%ax
mov    %ax,-0x2a(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,%rcx
add    $0x2,%rcx
mov    %rcx,-0x10(%rbp)
mov    (%rax),%ax
mov    %ax,-0x2c(%rbp)
movswq -0x2a(%rbp),%rax
movswq -0x2c(%rbp),%rcx
imul   %rcx,%rax
sar    $0xf,%rax
mov    %ax,-0x2a(%rbp)
movswq -0x2a(%rbp),%rax
add    -0x28(%rbp),%rax
mov    %rax,-0x28(%rbp)
jmp    1117 <fix_dot+0x17>
cmpq   $0x7fff,-0x28(%rbp)
jle    119a <fix_dot+0x9a>
movq   $0x7fff,-0x28(%rbp)
jmp    11b5 <fix_dot+0xb5>
cmpq   $0xffffffffffff8001,-0x28(%rbp)
ff
jge    11b0 <fix_dot+0xb0>
movq   $0xffffffffffff8001,-0x28(%rbp)
ff
jmp    11b5 <fix_dot+0xb5>
mov    -0x28(%rbp),%rax
cwtl
pop    %rbp
ret

x86 assembly (clang, -O2):
<fix_dot>:
test   %edx,%edx
je     1114 <fix_dot+0x14>
push   %rbx
cmp    $0x1,%edx
jne    111b <fix_dot+0x1b>
xor    %r11d,%r11d
test   $0x1,%dl
jne    1177 <fix_dot+0x77>
jmp    118a <fix_dot+0x8a>
xor    %edx,%edx
jmp    11ac <fix_dot+0xac>
mov    %edx,%r8d
and    $0xfffffffe,%r8d
xor    %edi,%edi
xor    %r9d,%r9d
xor    %r11d,%r11d
nopw   0x0(%rax,%rax,1)
movswl (%rax),%r10d
movswl (%rsi,%rdi,2),%ecx
imul   %r10d,%ecx
shr    $0xf,%ecx
movswq %cx,%rbx
add    %r11,%rbx
movswl 0x2(%rax),%r10d
add    $0x4,%rax
movswl 0x2(%rsi,%rdi,2),%ecx
imul   %r10d,%ecx
shr    $0xf,%ecx
movswq %cx,%r11
add    %rbx,%r11
add    $0xfffffffffffffffc,%r9
add    $0x2,%rdi
cmp    %edi,%r8d
jne    1130 <fix_dot+0x30>
sub    %r9,%rsi
test   $0x1,%dl
je     118a <fix_dot+0x8a>
movswl (%rax),%eax
movswl (%rsi),%ecx
imul   %eax,%ecx
shr    $0xf,%ecx
movswq %cx,%rax
add    %rax,%r11
mov    $0xffffffffffff8001,%rdx
cmp    $0xffffffffffff8002,%r11
pop    %rbx
jl     11ac <fix_dot+0xac>
mov    $0x7fff,%eax
mov    %r11,%rdx
cmp    $0x7ffe,%r11
jg     11af <fix_dot+0xaf>
mov    %rdx,%rax
ret

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

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