int mulm(b, ac3)
unsigned short b[], ac3[];
{
  register unsigned short *p, *q;
  unsigned short act[14 + 2];
  unsigned short carry, *r;
  unsigned long lp, a;
  int i, k, m, o;
  qclear(act);
  act[0] = ac3[0];
  act[1] = ac3[1];
  act[14] = 0;
  act[14 + 1] = 0;
  r = &act[14 + 1];
  for (k = 14; k >= 3; k--) {
    if (k == 14) {
      m = 14 - 1;
      o = 4;
    } else {
      m = k;
      o = 3;
    }
    q = &ac3[o];
    p = &b[m];
    for (i = m; i >= o; i--) {
      if ((*p == 0) || (*q == 0)) {
        --p;
        ++q;
        continue;
      }
      lp = (unsigned long)(*p--) * (unsigned long)(*q++);
      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 += 2;
    }
    --r;
  }
  mdnorm(act);
  qmov(act, ac3);
  return 0;
}


x86 assembly (clang, -O0):
<mulm>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%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)
movw   $0x0,-0x22(%rbp)
lea    -0x40(%rbp),%rax
add    $0x1e,%rax
mov    %rax,-0x50(%rbp)
movl   $0xe,-0x68(%rbp)
cmpl   $0x3,-0x68(%rbp)
jl     133f <mulm+0x20f>
cmpl   $0xe,-0x68(%rbp)
jne    11a8 <mulm+0x78>
movl   $0xd,-0x6c(%rbp)
movl   $0x4,-0x70(%rbp)
jmp    11b5 <mulm+0x85>
mov    -0x68(%rbp),%eax
mov    %eax,-0x6c(%rbp)
movl   $0x3,-0x70(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x70(%rbp),%rcx
shl    $0x1,%rcx
add    %rcx,%rax
mov    %rax,-0x20(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x6c(%rbp),%rcx
shl    $0x1,%rcx
add    %rcx,%rax
mov    %rax,-0x18(%rbp)
mov    -0x6c(%rbp),%eax
mov    %eax,-0x64(%rbp)
mov    -0x64(%rbp),%eax
cmp    -0x70(%rbp),%eax
jl     1325 <mulm+0x1f5>
mov    -0x18(%rbp),%rax
movzwl (%rax),%eax
cmp    $0x0,%eax
je     120d <mulm+0xdd>
mov    -0x20(%rbp),%rax
movzwl (%rax),%eax
cmp    $0x0,%eax
jne    122a <mulm+0xfa>
mov    -0x18(%rbp),%rax
add    $0xfffffffffffffffe,%rax
mov    %rax,-0x18(%rbp)
mov    -0x20(%rbp),%rax
add    $0x2,%rax
mov    %rax,-0x20(%rbp)
jmp    1317 <mulm+0x1e7>
mov    -0x18(%rbp),%rax
mov    %rax,%rcx
add    $0xfffffffffffffffe,%rcx
mov    %rcx,-0x18(%rbp)
movzwl (%rax),%eax
mov    -0x20(%rbp),%rcx
mov    %rcx,%rdx
add    $0x2,%rdx
mov    %rdx,-0x20(%rbp)
movzwl (%rcx),%ecx
imul   %rcx,%rax
mov    %rax,-0x58(%rbp)
mov    -0x50(%rbp),%rax
movzwl (%rax),%eax
mov    -0x58(%rbp),%rcx
and    $0xffff,%rcx
add    %rcx,%rax
mov    %rax,-0x60(%rbp)
mov    -0x60(%rbp),%rax
and    $0x10000,%rax
cmp    $0x0,%rax
je     128e <mulm+0x15e>
movw   $0x1,-0x42(%rbp)
jmp    1294 <mulm+0x164>
movw   $0x0,-0x42(%rbp)
mov    -0x60(%rbp),%rax
mov    %ax,%cx
mov    -0x50(%rbp),%rax
mov    %rax,%rdx
add    $0xfffffffffffffffe,%rdx
mov    %rdx,-0x50(%rbp)
mov    %cx,(%rax)
mov    -0x50(%rbp),%rax
movzwl (%rax),%eax
mov    -0x58(%rbp),%rcx
shr    $0x10,%rcx
and    $0xffff,%rcx
add    %rcx,%rax
movzwl -0x42(%rbp),%ecx
add    %rcx,%rax
mov    %rax,-0x60(%rbp)
mov    -0x60(%rbp),%rax
mov    %ax,%cx
mov    -0x50(%rbp),%rax
mov    %rax,%rdx
add    $0xfffffffffffffffe,%rdx
mov    %rdx,-0x50(%rbp)
mov    %cx,(%rax)
mov    -0x60(%rbp),%rax
and    $0x10000,%rax
cmp    $0x0,%rax
je     130b <mulm+0x1db>
mov    -0x50(%rbp),%rax
movzwl (%rax),%ecx
add    $0x1,%ecx
mov    %cx,(%rax)
mov    -0x50(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x50(%rbp)
mov    -0x64(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x64(%rbp)
jmp    11e1 <mulm+0xb1>
mov    -0x50(%rbp),%rax
add    $0xfffffffffffffffe,%rax
mov    %rax,-0x50(%rbp)
mov    -0x68(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x68(%rbp)
jmp    1181 <mulm+0x51>
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    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<mulm>:
push   %r15
push   %r14
push   %rbx
sub    $0x20,%rsp
mov    %rsi,%r14
mov    %rdi,%r15
mov    %rsp,%rdi
xor    %eax,%eax
call   1050 <qclear@plt>
movzwl (%r14),%eax
mov    %ax,(%rsp)
movzwl 0x2(%r14),%eax
mov    %ax,0x2(%rsp)
movl   $0x0,0x1c(%rsp)
lea    0x1e(%rsp),%rax
mov    $0xe,%r9d
mov    $0xd,%r8d
jmp    1195 <mulm+0x65>
cs nopw 0x0(%rax,%rax,1)
add    $0xfffffffffffffffe,%rax
lea    -0x1(%r9),%ecx
cmp    $0x3,%r9d
mov    %ecx,%r9d
jbe    1217 <mulm+0xe7>
xor    %r10d,%r10d
cmp    $0xe,%r9d
sete   %r10b
mov    %r9d,%ecx
cmove  %r8d,%ecx
add    $0x3,%r10d
cmp    %r10d,%ecx
jl     1180 <mulm+0x50>
mov    %ecx,%edi
mov    %r10d,%ecx
lea    (%r14,%rcx,2),%rcx
jmp    11d0 <mulm+0xa0>
nopl   0x0(%rax,%rax,1)
add    $0x2,%rcx
add    $0xffffffffffffffff,%rdi
lea    0x1(%rdi),%edx
cmp    %r10d,%edx
jle    1180 <mulm+0x50>
movzwl (%r15,%rdi,2),%ebx
test   %rbx,%rbx
je     11c0 <mulm+0x90>
movzwl (%rcx),%edx
test   %rdx,%rdx
je     11c0 <mulm+0x90>
imul   %rbx,%rdx
movzwl (%rax),%ebx
movzwl %dx,%esi
add    %rbx,%rsi
mov    %si,(%rax)
shr    $0x10,%rsi
movzwl -0x2(%rax),%ebx
shr    $0x10,%rdx
add    %rbx,%rdx
add    %rsi,%rdx
mov    %dx,-0x2(%rax)
test   $0x10000,%edx
je     11c0 <mulm+0x90>
addw   $0x1,-0x4(%rax)
jmp    11c0 <mulm+0x90>
mov    %rsp,%rbx
mov    %rbx,%rdi
xor    %eax,%eax
call   1030 <mdnorm@plt>
mov    %rbx,%rdi
mov    %r14,%rsi
xor    %eax,%eax
call   1040 <qmov@plt>
xor    %eax,%eax
add    $0x20,%rsp
pop    %rbx
pop    %r14
pop    %r15
ret

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

LLVM IR (-O2):
define i32 @mulm(i16* nocapture %0, i16* %1) {
%3 = alloca [16 x i16]
%4 = bitcast [16 x i16]* %3 to i8*
%5 = getelementptr inbounds [16 x i16], [16 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 [16 x i16], [16 x i16]* %3, i64 0, i64 1
store i16 %9, i16* %10
%11 = getelementptr inbounds [16 x i16], [16 x i16]* %3, i64 0, i64 14
store i16 0, i16* %11
%12 = getelementptr inbounds [16 x i16], [16 x i16]* %3, i64 0, i64 15
store i16 0, i16* %12
br label %13
13: ; preds = %2, %62
%14 = phi i32 [ 14, %2 ], [ %64, %62 ]
%15 = phi i16* [ %12, %2 ], [ %63, %62 ]
%16 = icmp eq i32 %14, 14
%17 = select i1 %16, i32 13, i32 %14
%18 = select i1 %16, i32 4, i32 3
%19 = icmp slt i32 %17, %18
br i1 %19, label %62, label %20
20: ; preds = %13
%21 = zext i32 %17 to i64
%22 = getelementptr inbounds i16, i16* %0, i64 %21
%23 = zext i32 %18 to i64
%24 = getelementptr inbounds i16, i16* %1, i64 %23
%25 = getelementptr inbounds i16, i16* %15, i64 -1
%26 = getelementptr inbounds i16, i16* %15, i64 -2
br label %27
27: ; preds = %20, %57
%28 = phi i32 [ %17, %20 ], [ %60, %57 ]
%29 = phi i16* [ %22, %20 ], [ %59, %57 ]
%30 = phi i16* [ %24, %20 ], [ %58, %57 ]
%31 = load i16, i16* %29
%32 = icmp eq i16 %31, 0
br i1 %32, label %57, label %33
33: ; preds = %27
%34 = load i16, i16* %30
%35 = icmp eq i16 %34, 0
br i1 %35, label %57, label %36
36: ; preds = %33
%37 = zext i16 %31 to i64
%38 = zext i16 %34 to i64
%39 = mul nuw nsw i64 %38, %37
%40 = load i16, i16* %15
%41 = zext i16 %40 to i64
%42 = and i64 %39, 65535
%43 = add nuw nsw i64 %42, %41
%44 = lshr i64 %43, 16
%45 = trunc i64 %43 to i16
store i16 %45, i16* %15
%46 = load i16, i16* %25
%47 = zext i16 %46 to i64
%48 = lshr i64 %39, 16
%49 = add nuw nsw i64 %48, %47
%50 = add nuw nsw i64 %49, %44
%51 = trunc i64 %50 to i16
store i16 %51, i16* %25
%52 = and i64 %50, 65536
%53 = icmp eq i64 %52, 0
br i1 %53, label %57, label %54
54: ; preds = %36
%55 = load i16, i16* %26
%56 = add i16 %55, 1
store i16 %56, i16* %26
br label %57
57: ; preds = %36, %54, %27, %33
%58 = getelementptr inbounds i16, i16* %30, i64 1
%59 = getelementptr inbounds i16, i16* %29, i64 -1
%60 = add nsw i32 %28, -1
%61 = icmp sgt i32 %28, %18
br i1 %61, label %27, label %62
62: ; preds = %57, %13
%63 = getelementptr inbounds i16, i16* %15, i64 -1
%64 = add nsw i32 %14, -1
%65 = icmp ugt i32 %14, 3
br i1 %65, label %13, label %66
66: ; preds = %62
%67 = call i32 (i16*, ...) bitcast (i32 (...)* @mdnorm to i32 (i16*, ...)*)(i16* %5)
%68 = call i32 (i16*, i16*, ...) bitcast (i32 (...)* @qmov to i32 (i16*, i16*, ...)*)(i16* %5, i16* %1)
ret i32 0
}