int mulin(b, ac3)
unsigned short b[], ac3[];
{
  register unsigned short *p, *r;
  unsigned short act[14 + 1];
  unsigned short carry, y;
  unsigned long lp, a;
  int i;
  qclear(act);
  act[0] = ac3[0];
  act[1] = ac3[1];
  act[14] = 0;
  r = &act[14];
  y = b[3];
  p = &ac3[14 - 1];
  for (i = 14 - 1; i >= 3; i--) {
    if (*p == 0) {
      --p;
      --r;
      continue;
    }
    lp = (unsigned long)(*p--) * y;
    a = (unsigned long)(*r) + (lp & 0xffff);
    if (a & 0x10000)
      carry = 1;
    else
      carry = 0;
    *r-- = a;
    a = (unsigned long)(*r) + ((lp >> 16) & 0xffff) + carry;
    *r-- = a;
    if (a & 0x10000)
      *r += 1;
    ++r;
  }
  mdnorm(act);
  qmov(act, ac3);
  return 0;
}


x86 assembly (clang, -O0):
<mulin>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
lea    -0x40(%rbp),%rdi
mov    $0x0,%al
call   1050 <qclear@plt>
mov    -0x10(%rbp),%rax
mov    (%rax),%ax
mov    %ax,-0x40(%rbp)
mov    -0x10(%rbp),%rax
mov    0x2(%rax),%ax
mov    %ax,-0x3e(%rbp)
movw   $0x0,-0x24(%rbp)
lea    -0x40(%rbp),%rax
add    $0x1c,%rax
mov    %rax,-0x20(%rbp)
mov    -0x8(%rbp),%rax
mov    0x6(%rax),%ax
mov    %ax,-0x44(%rbp)
mov    -0x10(%rbp),%rax
add    $0x1a,%rax
mov    %rax,-0x18(%rbp)
movl   $0xd,-0x5c(%rbp)
cmpl   $0x3,-0x5c(%rbp)
jl     12b7 <mulin+0x187>
mov    -0x18(%rbp),%rax
movzwl (%rax),%eax
cmp    $0x0,%eax
jne    11ca <mulin+0x9a>
mov    -0x18(%rbp),%rax
add    $0xfffffffffffffffe,%rax
mov    %rax,-0x18(%rbp)
mov    -0x20(%rbp),%rax
add    $0xfffffffffffffffe,%rax
mov    %rax,-0x20(%rbp)
jmp    12a9 <mulin+0x179>
mov    -0x18(%rbp),%rax
mov    %rax,%rcx
add    $0xfffffffffffffffe,%rcx
mov    %rcx,-0x18(%rbp)
movzwl (%rax),%eax
movzwl -0x44(%rbp),%ecx
imul   %rcx,%rax
mov    %rax,-0x50(%rbp)
mov    -0x20(%rbp),%rax
movzwl (%rax),%eax
mov    -0x50(%rbp),%rcx
and    $0xffff,%rcx
add    %rcx,%rax
mov    %rax,-0x58(%rbp)
mov    -0x58(%rbp),%rax
and    $0x10000,%rax
cmp    $0x0,%rax
je     1220 <mulin+0xf0>
movw   $0x1,-0x42(%rbp)
jmp    1226 <mulin+0xf6>
movw   $0x0,-0x42(%rbp)
mov    -0x58(%rbp),%rax
mov    %ax,%cx
mov    -0x20(%rbp),%rax
mov    %rax,%rdx
add    $0xfffffffffffffffe,%rdx
mov    %rdx,-0x20(%rbp)
mov    %cx,(%rax)
mov    -0x20(%rbp),%rax
movzwl (%rax),%eax
mov    -0x50(%rbp),%rcx
shr    $0x10,%rcx
and    $0xffff,%rcx
add    %rcx,%rax
movzwl -0x42(%rbp),%ecx
add    %rcx,%rax
mov    %rax,-0x58(%rbp)
mov    -0x58(%rbp),%rax
mov    %ax,%cx
mov    -0x20(%rbp),%rax
mov    %rax,%rdx
add    $0xfffffffffffffffe,%rdx
mov    %rdx,-0x20(%rbp)
mov    %cx,(%rax)
mov    -0x58(%rbp),%rax
and    $0x10000,%rax
cmp    $0x0,%rax
je     129d <mulin+0x16d>
mov    -0x20(%rbp),%rax
movzwl (%rax),%ecx
add    $0x1,%ecx
mov    %cx,(%rax)
mov    -0x20(%rbp),%rax
add    $0x2,%rax
mov    %rax,-0x20(%rbp)
mov    -0x5c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x5c(%rbp)
jmp    1193 <mulin+0x63>
lea    -0x40(%rbp),%rdi
mov    $0x0,%al
call   1030 <mdnorm@plt>
lea    -0x40(%rbp),%rdi
mov    -0x10(%rbp),%rsi
mov    $0x0,%al
call   1040 <qmov@plt>
xor    %eax,%eax
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<mulin>:
push   %r14
push   %rbx
sub    $0x28,%rsp
mov    %rsi,%rbx
mov    %rdi,%r14
mov    %rsp,%rdi
xor    %eax,%eax
call   1050 <qclear@plt>
movzwl (%rbx),%eax
mov    %ax,(%rsp)
movzwl 0x2(%rbx),%eax
mov    %ax,0x2(%rsp)
movw   $0x0,0x1c(%rsp)
movzwl 0x6(%r14),%eax
mov    $0xe,%ecx
jmp    1179 <mulin+0x49>
nopw   0x0(%rax,%rax,1)
add    $0xffffffffffffffff,%rcx
cmp    $0x3,%ecx
jbe    11bd <mulin+0x8d>
movzwl -0x2(%rbx,%rcx,2),%edx
test   %rdx,%rdx
je     1170 <mulin+0x40>
imul   %rax,%rdx
movzwl (%rsp,%rcx,2),%esi
movzwl %dx,%edi
add    %rsi,%rdi
mov    %di,(%rsp,%rcx,2)
shr    $0x10,%rdi
movzwl -0x2(%rsp,%rcx,2),%esi
shr    $0x10,%rdx
add    %rsi,%rdx
add    %rdi,%rdx
mov    %dx,-0x2(%rsp,%rcx,2)
test   $0x10000,%edx
je     1170 <mulin+0x40>
addw   $0x1,-0x4(%rsp,%rcx,2)
jmp    1170 <mulin+0x40>
mov    %rsp,%r14
mov    %r14,%rdi
xor    %eax,%eax
call   1030 <mdnorm@plt>
mov    %r14,%rdi
mov    %rbx,%rsi
xor    %eax,%eax
call   1040 <qmov@plt>
xor    %eax,%eax
add    $0x28,%rsp
pop    %rbx
pop    %r14
ret

LLVM IR (-O0):
define i32 @mulin(i16* %0, i16* %1) {
%3 = alloca i16*
%4 = alloca i16*
%5 = alloca i16*
%6 = alloca i16*
%7 = alloca [15 x i16]
%8 = alloca i16
%9 = alloca i16
%10 = alloca i64
%11 = alloca i64
%12 = alloca i32
store i16* %0, i16** %3
store i16* %1, i16** %4
%13 = getelementptr inbounds [15 x i16], [15 x i16]* %7, i64 0, i64 0
%14 = call i32 (i16*, ...) bitcast (i32 (...)* @qclear to i32 (i16*, ...)*)(i16* %13)
%15 = load i16*, i16** %4
%16 = getelementptr inbounds i16, i16* %15, i64 0
%17 = load i16, i16* %16
%18 = getelementptr inbounds [15 x i16], [15 x i16]* %7, i64 0, i64 0
store i16 %17, i16* %18
%19 = load i16*, i16** %4
%20 = getelementptr inbounds i16, i16* %19, i64 1
%21 = load i16, i16* %20
%22 = getelementptr inbounds [15 x i16], [15 x i16]* %7, i64 0, i64 1
store i16 %21, i16* %22
%23 = getelementptr inbounds [15 x i16], [15 x i16]* %7, i64 0, i64 14
store i16 0, i16* %23
%24 = getelementptr inbounds [15 x i16], [15 x i16]* %7, i64 0, i64 14
store i16* %24, i16** %6
%25 = load i16*, i16** %3
%26 = getelementptr inbounds i16, i16* %25, i64 3
%27 = load i16, i16* %26
store i16 %27, i16* %9
%28 = load i16*, i16** %4
%29 = getelementptr inbounds i16, i16* %28, i64 13
store i16* %29, i16** %5
store i32 13, i32* %12
br label %30
30: ; preds = %93, %2
%31 = load i32, i32* %12
%32 = icmp sge i32 %31, 3
br i1 %32, label %33, label %96
33: ; preds = %30
%34 = load i16*, i16** %5
%35 = load i16, i16* %34
%36 = zext i16 %35 to i32
%37 = icmp eq i32 %36, 0
br i1 %37, label %38, label %43
38: ; preds = %33
%39 = load i16*, i16** %5
%40 = getelementptr inbounds i16, i16* %39, i32 -1
store i16* %40, i16** %5
%41 = load i16*, i16** %6
%42 = getelementptr inbounds i16, i16* %41, i32 -1
store i16* %42, i16** %6
br label %93
43: ; preds = %33
%44 = load i16*, i16** %5
%45 = getelementptr inbounds i16, i16* %44, i32 -1
store i16* %45, i16** %5
%46 = load i16, i16* %44
%47 = zext i16 %46 to i64
%48 = load i16, i16* %9
%49 = zext i16 %48 to i64
%50 = mul i64 %47, %49
store i64 %50, i64* %10
%51 = load i16*, i16** %6
%52 = load i16, i16* %51
%53 = zext i16 %52 to i64
%54 = load i64, i64* %10
%55 = and i64 %54, 65535
%56 = add i64 %53, %55
store i64 %56, i64* %11
%57 = load i64, i64* %11
%58 = and i64 %57, 65536
%59 = icmp ne i64 %58, 0
br i1 %59, label %60, label %61
60: ; preds = %43
store i16 1, i16* %8
br label %62
61: ; preds = %43
store i16 0, i16* %8
br label %62
62: ; preds = %61, %60
%63 = load i64, i64* %11
%64 = trunc i64 %63 to i16
%65 = load i16*, i16** %6
%66 = getelementptr inbounds i16, i16* %65, i32 -1
store i16* %66, i16** %6
store i16 %64, i16* %65
%67 = load i16*, i16** %6
%68 = load i16, i16* %67
%69 = zext i16 %68 to i64
%70 = load i64, i64* %10
%71 = lshr i64 %70, 16
%72 = and i64 %71, 65535
%73 = add i64 %69, %72
%74 = load i16, i16* %8
%75 = zext i16 %74 to i64
%76 = add i64 %73, %75
store i64 %76, i64* %11
%77 = load i64, i64* %11
%78 = trunc i64 %77 to i16
%79 = load i16*, i16** %6
%80 = getelementptr inbounds i16, i16* %79, i32 -1
store i16* %80, i16** %6
store i16 %78, i16* %79
%81 = load i64, i64* %11
%82 = and i64 %81, 65536
%83 = icmp ne i64 %82, 0
br i1 %83, label %84, label %90
84: ; preds = %62
%85 = load i16*, i16** %6
%86 = load i16, i16* %85
%87 = zext i16 %86 to i32
%88 = add nsw i32 %87, 1
%89 = trunc i32 %88 to i16
store i16 %89, i16* %85
br label %90
90: ; preds = %84, %62
%91 = load i16*, i16** %6
%92 = getelementptr inbounds i16, i16* %91, i32 1
store i16* %92, i16** %6
br label %93
93: ; preds = %90, %38
%94 = load i32, i32* %12
%95 = add nsw i32 %94, -1
store i32 %95, i32* %12
br label %30
96: ; preds = %30
%97 = getelementptr inbounds [15 x i16], [15 x i16]* %7, i64 0, i64 0
%98 = call i32 (i16*, ...) bitcast (i32 (...)* @mdnorm to i32 (i16*, ...)*)(i16* %97)
%99 = getelementptr inbounds [15 x i16], [15 x i16]* %7, i64 0, i64 0
%100 = load i16*, i16** %4
%101 = call i32 (i16*, i16*, ...) bitcast (i32 (...)* @qmov to i32 (i16*, i16*, ...)*)(i16* %99, i16* %100)
ret i32 0
}

LLVM IR (-O2):
define i32 @mulin(i16* nocapture %0, i16* %1) {
%3 = alloca [15 x i16]
%4 = bitcast [15 x i16]* %3 to i8*
%5 = getelementptr inbounds [15 x i16], [15 x i16]* %3, i64 0, i64 0
%6 = call i32 (i16*, ...) bitcast (i32 (...)* @qclear to i32 (i16*, ...)*)(i16* %5)
%7 = load i16, i16* %1
store i16 %7, i16* %5
%8 = getelementptr inbounds i16, i16* %1, i64 1
%9 = load i16, i16* %8
%10 = getelementptr inbounds [15 x i16], [15 x i16]* %3, i64 0, i64 1
store i16 %9, i16* %10
%11 = getelementptr inbounds [15 x i16], [15 x i16]* %3, i64 0, i64 14
store i16 0, i16* %11
%12 = getelementptr inbounds i16, i16* %0, i64 3
%13 = load i16, i16* %12
%14 = getelementptr inbounds i16, i16* %1, i64 13
%15 = zext i16 %13 to i64
br label %16
16: ; preds = %2, %44
%17 = phi i32 [ 13, %2 ], [ %47, %44 ]
%18 = phi i16* [ %14, %2 ], [ %46, %44 ]
%19 = phi i16* [ %11, %2 ], [ %45, %44 ]
%20 = load i16, i16* %18
%21 = icmp eq i16 %20, 0
br i1 %21, label %44, label %22
22: ; preds = %16
%23 = zext i16 %20 to i64
%24 = mul nuw nsw i64 %23, %15
%25 = load i16, i16* %19
%26 = zext i16 %25 to i64
%27 = and i64 %24, 65535
%28 = add nuw nsw i64 %27, %26
%29 = lshr i64 %28, 16
%30 = trunc i64 %28 to i16
%31 = getelementptr inbounds i16, i16* %19, i64 -1
store i16 %30, i16* %19
%32 = load i16, i16* %31
%33 = zext i16 %32 to i64
%34 = lshr i64 %24, 16
%35 = add nuw nsw i64 %34, %33
%36 = add nuw nsw i64 %35, %29
%37 = trunc i64 %36 to i16
%38 = getelementptr inbounds i16, i16* %19, i64 -2
store i16 %37, i16* %31
%39 = and i64 %36, 65536
%40 = icmp eq i64 %39, 0
br i1 %40, label %44, label %41
41: ; preds = %22
%42 = load i16, i16* %38
%43 = add i16 %42, 1
store i16 %43, i16* %38
br label %44
44: ; preds = %22, %41, %16
%45 = getelementptr inbounds i16, i16* %19, i64 -1
%46 = getelementptr inbounds i16, i16* %18, i64 -1
%47 = add nsw i32 %17, -1
%48 = icmp ugt i32 %17, 3
br i1 %48, label %16, label %49
49: ; preds = %44
%50 = call i32 (i16*, ...) bitcast (i32 (...)* @mdnorm to i32 (i16*, ...)*)(i16* %5)
%51 = call i32 (i16*, i16*, ...) bitcast (i32 (...)* @qmov to i32 (i16*, i16*, ...)*)(i16* %5, i16* %1)
ret i32 0
}