void Subsample_Horizontal(short *s, unsigned char *d, int x0, int lx, int lxs,
                          int lxd, int ly, int m, int n) {
  int i;
  int i1;
  int j;
  int id;
  int c1;
  int c2;
  int v;
  short *s1;
  short *s2;
  unsigned char *d1;
  int tmp;
  {
    i1 = 0;
    while (i1 < lx) {
      d1 = d + i1;
      i = x0 + i1;
      id = (i * m) / n;
      s1 = s + id;
      if (id < lxs - 1) {
        s2 = s1 + 1;
      } else {
        s2 = s1;
      }
      c2 = (16 * ((i * m) % n) + (n >> 1)) / n;
      c1 = 16 - c2;
      j = 0;
      while (j < ly) {
        v = c1 * (int)*s1 + c2 * (int)*s2;
        if (v >= 0) {
          tmp = 128;
        } else {
          tmp = 127;
        }
        *d1 = (unsigned char)((v + tmp) >> 8);
        d1 += lxd;
        s1 += lxs;
        s2 += lxs;
        j++;
      }
      i1++;
    }
    return;
  }
}


x86 assembly (clang, -O0):
<Subsample_Horizontal>:
push   %rbp
mov    %rsp,%rbp
mov    0x20(%rbp),%eax
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8d,-0x1c(%rbp)
mov    %r9d,-0x20(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    127a <Subsample_Horizontal+0x17a>
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
add    %rcx,%rax
mov    %rax,-0x58(%rbp)
mov    -0x14(%rbp),%eax
add    -0x28(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
imul   0x18(%rbp),%eax
cltd
idivl  0x20(%rbp)
mov    %eax,-0x30(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
shl    $0x1,%rcx
add    %rcx,%rax
mov    %rax,-0x48(%rbp)
mov    -0x30(%rbp),%eax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1191 <Subsample_Horizontal+0x91>
mov    -0x48(%rbp),%rax
add    $0x2,%rax
mov    %rax,-0x50(%rbp)
jmp    1199 <Subsample_Horizontal+0x99>
mov    -0x48(%rbp),%rax
mov    %rax,-0x50(%rbp)
mov    -0x24(%rbp),%eax
imul   0x18(%rbp),%eax
cltd
idivl  0x20(%rbp)
mov    %edx,%eax
shl    $0x4,%eax
mov    0x20(%rbp),%ecx
sar    $0x1,%ecx
add    %ecx,%eax
cltd
idivl  0x20(%rbp)
mov    %eax,-0x38(%rbp)
mov    $0x10,%eax
sub    -0x38(%rbp),%eax
mov    %eax,-0x34(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    0x10(%rbp),%eax
jge    126c <Subsample_Horizontal+0x16c>
mov    -0x34(%rbp),%eax
mov    -0x48(%rbp),%rcx
movswl (%rcx),%ecx
imul   %ecx,%eax
mov    -0x38(%rbp),%ecx
mov    -0x50(%rbp),%rdx
movswl (%rdx),%edx
imul   %edx,%ecx
add    %ecx,%eax
mov    %eax,-0x3c(%rbp)
cmpl   $0x0,-0x3c(%rbp)
jl     120b <Subsample_Horizontal+0x10b>
movl   $0x80,-0x5c(%rbp)
jmp    1212 <Subsample_Horizontal+0x112>
movl   $0x7f,-0x5c(%rbp)
mov    -0x3c(%rbp),%eax
add    -0x5c(%rbp),%eax
sar    $0x8,%eax
mov    %al,%cl
mov    -0x58(%rbp),%rax
mov    %cl,(%rax)
mov    -0x20(%rbp),%ecx
mov    -0x58(%rbp),%rax
movslq %ecx,%rcx
add    %rcx,%rax
mov    %rax,-0x58(%rbp)
mov    -0x1c(%rbp),%ecx
mov    -0x48(%rbp),%rax
movslq %ecx,%rcx
shl    $0x1,%rcx
add    %rcx,%rax
mov    %rax,-0x48(%rbp)
mov    -0x1c(%rbp),%ecx
mov    -0x50(%rbp),%rax
movslq %ecx,%rcx
shl    $0x1,%rcx
add    %rcx,%rax
mov    %rax,-0x50(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11ca <Subsample_Horizontal+0xca>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    112a <Subsample_Horizontal+0x2a>
pop    %rbp
ret

x86 assembly (clang, -O2):
<Subsample_Horizontal>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %ecx,%ecx
jle    129a <Subsample_Horizontal+0x19a>
mov    %edx,%r11d
mov    0x48(%rsp),%r12d
mov    0x38(%rsp),%eax
mov    %r12d,%edx
sar    %edx
mov    %rdx,-0x8(%rsp)
movslq %r9d,%rbp
movslq %r8d,%r15
add    $0xffffffff,%r8d
mov    %ecx,%ecx
mov    %rcx,-0x10(%rsp)
and    $0xfffffffe,%eax
mov    %eax,-0x3c(%rsp)
lea    0x0(,%r15,4),%r13
lea    (%r15,%r15,1),%rax
mov    %rax,-0x38(%rsp)
xor    %r14d,%r14d
mov    %rbp,-0x18(%rsp)
mov    %r11,-0x20(%rsp)
mov    %rsi,-0x28(%rsp)
mov    %rdi,-0x30(%rsp)
jmp    117f <Subsample_Horizontal+0x7f>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r14
cmp    -0x10(%rsp),%r14
je     129a <Subsample_Horizontal+0x19a>
lea    (%r14,%r11,1),%eax
imul   0x40(%rsp),%eax
cltd
idiv   %r12d
mov    %eax,%ecx
shl    $0x4,%edx
mov    -0x8(%rsp),%rax
add    %edx,%eax
cltd
idiv   %r12d
cmpl   $0x0,0x38(%rsp)
jle    1170 <Subsample_Horizontal+0x70>
mov    $0x10,%edx
sub    %eax,%edx
movslq %ecx,%rcx
lea    (%rdi,%rcx,2),%rbx
xor    %r10d,%r10d
cmp    %r8d,%ecx
setl   %r10b
lea    (%rsi,%r14,1),%rcx
cmpl   $0x1,0x38(%rsp)
jne    11d0 <Subsample_Horizontal+0xd0>
lea    (%rbx,%r10,2),%r9
jmp    1265 <Subsample_Horizontal+0x165>
nop
mov    %r8d,%r12d
mov    -0x38(%rsp),%rsi
lea    (%rsi,%r10,2),%r8
add    %r10,%r10
mov    -0x3c(%rsp),%r11d
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    %rbx,%rsi
movswl (%rbx),%edi
imul   %edx,%edi
movswl (%rbx,%r10,1),%r9d
imul   %eax,%r9d
add    %edi,%r9d
mov    %r9d,%edi
not    %edi
shr    $0x1f,%edi
lea    (%rdi,%r9,1),%ebx
add    $0x7f,%ebx
mov    %bh,(%rcx)
movswl (%rsi,%r15,2),%edi
imul   %edx,%edi
movswl (%rsi,%r8,1),%ebx
imul   %eax,%ebx
add    %edi,%ebx
mov    %ebx,%edi
not    %edi
shr    $0x1f,%edi
add    %edi,%ebx
add    $0x7f,%ebx
mov    %bh,(%rcx,%rbp,1)
add    %rbp,%rcx
add    %rbp,%rcx
lea    (%rsi,%r13,1),%rbx
add    $0xfffffffe,%r11d
jne    11f0 <Subsample_Horizontal+0xf0>
add    %r8,%rsi
lea    (%rsi,%r15,2),%r9
mov    %r12d,%r8d
mov    -0x20(%rsp),%r11
mov    -0x28(%rsp),%rsi
mov    -0x30(%rsp),%rdi
mov    0x48(%rsp),%r12d
mov    0x38(%rsp),%ebp
test   $0x1,%bpl
mov    -0x18(%rsp),%rbp
je     1170 <Subsample_Horizontal+0x70>
movswl (%rbx),%ebx
imul   %ebx,%edx
movswl (%r9),%ebx
imul   %ebx,%eax
add    %edx,%eax
mov    %eax,%edx
not    %edx
shr    $0x1f,%edx
add    %edx,%eax
add    $0x7f,%eax
mov    %ah,(%rcx)
jmp    1170 <Subsample_Horizontal+0x70>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @Subsample_Horizontal(i16* %0, i8* %1, i32 %2, i32 %3, i32 %4, i32 %5, i32 %6, i32 %7, i32 %8) {
%10 = alloca i16*
%11 = alloca i8*
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
%20 = alloca i32
%21 = alloca i32
%22 = alloca i32
%23 = alloca i32
%24 = alloca i32
%25 = alloca i32
%26 = alloca i16*
%27 = alloca i16*
%28 = alloca i8*
%29 = alloca i32
store i16* %0, i16** %10
store i8* %1, i8** %11
store i32 %2, i32* %12
store i32 %3, i32* %13
store i32 %4, i32* %14
store i32 %5, i32* %15
store i32 %6, i32* %16
store i32 %7, i32* %17
store i32 %8, i32* %18
store i32 0, i32* %20
br label %30
30: ; preds = %115, %9
%31 = load i32, i32* %20
%32 = load i32, i32* %13
%33 = icmp slt i32 %31, %32
br i1 %33, label %34, label %118
34: ; preds = %30
%35 = load i8*, i8** %11
%36 = load i32, i32* %20
%37 = sext i32 %36 to i64
%38 = getelementptr inbounds i8, i8* %35, i64 %37
store i8* %38, i8** %28
%39 = load i32, i32* %12
%40 = load i32, i32* %20
%41 = add nsw i32 %39, %40
store i32 %41, i32* %19
%42 = load i32, i32* %19
%43 = load i32, i32* %17
%44 = mul nsw i32 %42, %43
%45 = load i32, i32* %18
%46 = sdiv i32 %44, %45
store i32 %46, i32* %22
%47 = load i16*, i16** %10
%48 = load i32, i32* %22
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds i16, i16* %47, i64 %49
store i16* %50, i16** %26
%51 = load i32, i32* %22
%52 = load i32, i32* %14
%53 = sub nsw i32 %52, 1
%54 = icmp slt i32 %51, %53
br i1 %54, label %55, label %58
55: ; preds = %34
%56 = load i16*, i16** %26
%57 = getelementptr inbounds i16, i16* %56, i64 1
store i16* %57, i16** %27
br label %60
58: ; preds = %34
%59 = load i16*, i16** %26
store i16* %59, i16** %27
br label %60
60: ; preds = %58, %55
%61 = load i32, i32* %19
%62 = load i32, i32* %17
%63 = mul nsw i32 %61, %62
%64 = load i32, i32* %18
%65 = srem i32 %63, %64
%66 = mul nsw i32 16, %65
%67 = load i32, i32* %18
%68 = ashr i32 %67, 1
%69 = add nsw i32 %66, %68
%70 = load i32, i32* %18
%71 = sdiv i32 %69, %70
store i32 %71, i32* %24
%72 = load i32, i32* %24
%73 = sub nsw i32 16, %72
store i32 %73, i32* %23
store i32 0, i32* %21
br label %74
74: ; preds = %94, %60
%75 = load i32, i32* %21
%76 = load i32, i32* %16
%77 = icmp slt i32 %75, %76
br i1 %77, label %78, label %115
78: ; preds = %74
%79 = load i32, i32* %23
%80 = load i16*, i16** %26
%81 = load i16, i16* %80
%82 = sext i16 %81 to i32
%83 = mul nsw i32 %79, %82
%84 = load i32, i32* %24
%85 = load i16*, i16** %27
%86 = load i16, i16* %85
%87 = sext i16 %86 to i32
%88 = mul nsw i32 %84, %87
%89 = add nsw i32 %83, %88
store i32 %89, i32* %25
%90 = load i32, i32* %25
%91 = icmp sge i32 %90, 0
br i1 %91, label %92, label %93
92: ; preds = %78
store i32 128, i32* %29
br label %94
93: ; preds = %78
store i32 127, i32* %29
br label %94
94: ; preds = %93, %92
%95 = load i32, i32* %25
%96 = load i32, i32* %29
%97 = add nsw i32 %95, %96
%98 = ashr i32 %97, 8
%99 = trunc i32 %98 to i8
%100 = load i8*, i8** %28
store i8 %99, i8* %100
%101 = load i32, i32* %15
%102 = load i8*, i8** %28
%103 = sext i32 %101 to i64
%104 = getelementptr inbounds i8, i8* %102, i64 %103
store i8* %104, i8** %28
%105 = load i32, i32* %14
%106 = load i16*, i16** %26
%107 = sext i32 %105 to i64
%108 = getelementptr inbounds i16, i16* %106, i64 %107
store i16* %108, i16** %26
%109 = load i32, i32* %14
%110 = load i16*, i16** %27
%111 = sext i32 %109 to i64
%112 = getelementptr inbounds i16, i16* %110, i64 %111
store i16* %112, i16** %27
%113 = load i32, i32* %21
%114 = add nsw i32 %113, 1
store i32 %114, i32* %21
br label %74
115: ; preds = %74
%116 = load i32, i32* %20
%117 = add nsw i32 %116, 1
store i32 %117, i32* %20
br label %30
118: ; preds = %30
ret void
}

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