void test_pmulhrsw_c(short *pDst, short *pSrc, int xmm) {
  int i;
  for (i = 0; i < (4 << xmm); i++) {
    int a = pSrc[i] * pDst[i];
    pDst[i] = (short)(((a >> 14) + 1) >> 1);
  }
}


x86 assembly (clang, -O0):
<test_pmulhrsw_c>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
movl   $0x0,-0x18(%rbp)
mov    -0x18(%rbp),%eax
mov    -0x14(%rbp),%ecx
mov    $0x4,%edx
shl    %cl,%edx
mov    %edx,%ecx
cmp    %ecx,%eax
jge    1174 <test_pmulhrsw_c+0x74>
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
movswl (%rax,%rcx,2),%eax
mov    -0x8(%rbp),%rcx
movslq -0x18(%rbp),%rdx
movswl (%rcx,%rdx,2),%ecx
imul   %ecx,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
sar    $0xe,%eax
add    $0x1,%eax
sar    $0x1,%eax
mov    %ax,%dx
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
mov    %dx,(%rax,%rcx,2)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1116 <test_pmulhrsw_c+0x16>
pop    %rbp
ret

x86 assembly (clang, -O2):
<test_pmulhrsw_c>:
mov    %edx,%ecx
mov    $0x4,%eax
shl    %cl,%eax
test   %eax,%eax
jle    12ca <test_pmulhrsw_c+0x1ca>
mov    %eax,%r9d
cmp    $0x8,%eax
jae    1120 <test_pmulhrsw_c+0x20>
xor    %ecx,%ecx
jmp    1290 <test_pmulhrsw_c+0x190>
lea    (%rsi,%r9,2),%rax
cmp    %rdi,%rax
jbe    1139 <test_pmulhrsw_c+0x39>
lea    (%rdi,%r9,2),%rax
cmp    %rsi,%rax
jbe    1139 <test_pmulhrsw_c+0x39>
xor    %ecx,%ecx
jmp    1290 <test_pmulhrsw_c+0x190>
mov    %r9d,%ecx
and    $0xfffffff8,%ecx
lea    -0x8(%rcx),%rax
mov    %rax,%r8
shr    $0x3,%r8
add    $0x1,%r8
test   %rax,%rax
je     12cb <test_pmulhrsw_c+0x1cb>
mov    %r8,%rax
and    $0xfffffffffffffffe,%rax
xor    %edx,%edx
pcmpeqd %xmm0,%xmm0
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movdqu (%rsi,%rdx,2),%xmm2
movdqu (%rdi,%rdx,2),%xmm3
movdqu 0x10(%rdi,%rdx,2),%xmm1
movdqa %xmm3,%xmm4
pmulhw %xmm2,%xmm4
pmullw %xmm2,%xmm3
movdqa %xmm3,%xmm2
punpcklwd %xmm4,%xmm2
punpckhwd %xmm4,%xmm3
psrld  $0xe,%xmm3
psrld  $0xe,%xmm2
psubd  %xmm0,%xmm2
psubd  %xmm0,%xmm3
pslld  $0xf,%xmm3
psrad  $0x10,%xmm3
pslld  $0xf,%xmm2
psrad  $0x10,%xmm2
packssdw %xmm3,%xmm2
movdqu %xmm2,(%rdi,%rdx,2)
movdqu 0x10(%rsi,%rdx,2),%xmm2
movdqa %xmm1,%xmm3
pmulhw %xmm2,%xmm3
pmullw %xmm2,%xmm1
movdqa %xmm1,%xmm2
punpcklwd %xmm3,%xmm2
punpckhwd %xmm3,%xmm1
psrld  $0xe,%xmm1
psrld  $0xe,%xmm2
psubd  %xmm0,%xmm2
psubd  %xmm0,%xmm1
pslld  $0xf,%xmm1
psrad  $0x10,%xmm1
pslld  $0xf,%xmm2
psrad  $0x10,%xmm2
packssdw %xmm1,%xmm2
movdqu %xmm2,0x10(%rdi,%rdx,2)
add    $0x10,%rdx
add    $0xfffffffffffffffe,%rax
jne    1170 <test_pmulhrsw_c+0x70>
test   $0x1,%r8b
je     127e <test_pmulhrsw_c+0x17e>
movdqu (%rsi,%rdx,2),%xmm0
movdqu (%rdi,%rdx,2),%xmm1
movdqa %xmm1,%xmm2
pmulhw %xmm0,%xmm2
pmullw %xmm0,%xmm1
movdqa %xmm1,%xmm0
punpcklwd %xmm2,%xmm0
punpckhwd %xmm2,%xmm1
psrld  $0xe,%xmm1
psrld  $0xe,%xmm0
pcmpeqd %xmm2,%xmm2
psubd  %xmm2,%xmm0
psubd  %xmm2,%xmm1
pslld  $0xf,%xmm1
psrad  $0x10,%xmm1
pslld  $0xf,%xmm0
psrad  $0x10,%xmm0
packssdw %xmm1,%xmm0
movdqu %xmm0,(%rdi,%rdx,2)
cmp    %r9,%rcx
je     12ca <test_pmulhrsw_c+0x1ca>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movswl (%rsi,%rcx,2),%eax
movswl (%rdi,%rcx,2),%edx
imul   %eax,%edx
shr    $0xe,%edx
add    $0x1,%edx
shr    %edx
mov    %dx,(%rdi,%rcx,2)
movswl 0x2(%rsi,%rcx,2),%eax
movswl 0x2(%rdi,%rcx,2),%edx
imul   %eax,%edx
shr    $0xe,%edx
add    $0x1,%edx
shr    %edx
mov    %dx,0x2(%rdi,%rcx,2)
add    $0x2,%rcx
cmp    %rcx,%r9
jne    1290 <test_pmulhrsw_c+0x190>
ret
xor    %edx,%edx
test   $0x1,%r8b
jne    1229 <test_pmulhrsw_c+0x129>
jmp    127e <test_pmulhrsw_c+0x17e>

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

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