int mulv(a, b, c, prec)
unsigned short a[], b[], c[];
int prec;
{
  register unsigned short *p, *q, *r;
  register unsigned long u, lp;
  unsigned short carry;
  int k, i;
  k = prec + 2;
  p = &c[2];
  do
    *p++ = 0;
  while (--k);
  r = &c[prec + 3];
  for (k = prec + 2; k >= 3; k--) {
    q = &b[3];
    p = &a[k];
    for (i = k; i >= 3; i--) {
      if ((*p == 0) || (*q == 0)) {
        --p;
        ++q;
        continue;
      }
      lp = (unsigned long)(*p--) * (unsigned long)(*q++);
      u = (unsigned long)(*r) + (lp & 0xffff);
      if (u & 0x10000)
        carry = 1;
      else
        carry = 0;
      *r-- = u;
      u = (unsigned long)(*r) + ((lp >> 16) & 0xffff) + carry;
      *r-- = u;
      if (u & 0x10000)
        *r += 1;
      r += 2;
    }
    --r;
  }
  return 0;
}


x86 assembly (clang, -O0):
<mulv>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x50(%rbp)
mov    -0x18(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x28(%rbp)
mov    -0x28(%rbp),%rax
mov    %rax,%rcx
add    $0x2,%rcx
mov    %rcx,-0x28(%rbp)
movw   $0x0,(%rax)
mov    -0x50(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x50(%rbp)
cmp    $0x0,%eax
jne    1128 <mulv+0x28>
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
shl    $0x1,%rcx
add    %rcx,%rax
mov    %rax,-0x38(%rbp)
mov    -0x1c(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x50(%rbp)
cmpl   $0x3,-0x50(%rbp)
jl     12fa <mulv+0x1fa>
mov    -0x10(%rbp),%rax
add    $0x6,%rax
mov    %rax,-0x30(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x50(%rbp),%rcx
shl    $0x1,%rcx
add    %rcx,%rax
mov    %rax,-0x28(%rbp)
mov    -0x50(%rbp),%eax
mov    %eax,-0x54(%rbp)
cmpl   $0x3,-0x54(%rbp)
jl     12e0 <mulv+0x1e0>
mov    -0x28(%rbp),%rax
movzwl (%rax),%eax
cmp    $0x0,%eax
je     11c8 <mulv+0xc8>
mov    -0x30(%rbp),%rax
movzwl (%rax),%eax
cmp    $0x0,%eax
jne    11e5 <mulv+0xe5>
mov    -0x28(%rbp),%rax
add    $0xfffffffffffffffe,%rax
mov    %rax,-0x28(%rbp)
mov    -0x30(%rbp),%rax
add    $0x2,%rax
mov    %rax,-0x30(%rbp)
jmp    12d2 <mulv+0x1d2>
mov    -0x28(%rbp),%rax
mov    %rax,%rcx
add    $0xfffffffffffffffe,%rcx
mov    %rcx,-0x28(%rbp)
movzwl (%rax),%eax
mov    -0x30(%rbp),%rcx
mov    %rcx,%rdx
add    $0x2,%rdx
mov    %rdx,-0x30(%rbp)
movzwl (%rcx),%ecx
imul   %rcx,%rax
mov    %rax,-0x48(%rbp)
mov    -0x38(%rbp),%rax
movzwl (%rax),%eax
mov    -0x48(%rbp),%rcx
and    $0xffff,%rcx
add    %rcx,%rax
mov    %rax,-0x40(%rbp)
mov    -0x40(%rbp),%rax
and    $0x10000,%rax
cmp    $0x0,%rax
je     1249 <mulv+0x149>
movw   $0x1,-0x4a(%rbp)
jmp    124f <mulv+0x14f>
movw   $0x0,-0x4a(%rbp)
mov    -0x40(%rbp),%rax
mov    %ax,%cx
mov    -0x38(%rbp),%rax
mov    %rax,%rdx
add    $0xfffffffffffffffe,%rdx
mov    %rdx,-0x38(%rbp)
mov    %cx,(%rax)
mov    -0x38(%rbp),%rax
movzwl (%rax),%eax
mov    -0x48(%rbp),%rcx
shr    $0x10,%rcx
and    $0xffff,%rcx
add    %rcx,%rax
movzwl -0x4a(%rbp),%ecx
add    %rcx,%rax
mov    %rax,-0x40(%rbp)
mov    -0x40(%rbp),%rax
mov    %ax,%cx
mov    -0x38(%rbp),%rax
mov    %rax,%rdx
add    $0xfffffffffffffffe,%rdx
mov    %rdx,-0x38(%rbp)
mov    %cx,(%rax)
mov    -0x40(%rbp),%rax
and    $0x10000,%rax
cmp    $0x0,%rax
je     12c6 <mulv+0x1c6>
mov    -0x38(%rbp),%rax
movzwl (%rax),%ecx
add    $0x1,%ecx
mov    %cx,(%rax)
mov    -0x38(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x38(%rbp)
mov    -0x54(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x54(%rbp)
jmp    119e <mulv+0x9e>
mov    -0x38(%rbp),%rax
add    $0xfffffffffffffffe,%rax
mov    %rax,-0x38(%rbp)
mov    -0x50(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x50(%rbp)
jmp    116f <mulv+0x6f>
xor    %eax,%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<mulv>:
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %ecx,%r13d
mov    %rdx,%r12
mov    %rsi,%r14
mov    %rdi,%r15
lea    0x4(%rdx),%rdi
lea    0x1(%r13),%eax
lea    (%rax,%rax,1),%rdx
add    $0x2,%rdx
xor    %esi,%esi
call   1030 <memset@plt>
test   %r13d,%r13d
jle    11ea <mulv+0xda>
lea    0x2(%r13),%r8d
add    $0x3,%r13d
lea    (%r12,%r13,2),%r9
add    $0x6,%r14
jmp    1178 <mulv+0x68>
nopw   0x0(%rax,%rax,1)
add    $0xfffffffffffffffe,%r9
mov    %r8,%rax
add    $0xffffffffffffffff,%rax
add    $0xffffffff,%r13d
cmp    $0x3,%r8d
mov    %rax,%r8
jle    11ea <mulv+0xda>
cmp    $0x3,%r8d
jl     11ea <mulv+0xda>
lea    (%r15,%r8,2),%rdx
mov    %r13d,%esi
mov    %r14,%rdi
jmp    11a0 <mulv+0x90>
nopw   0x0(%rax,%rax,1)
add    $0x2,%rdi
add    $0xfffffffffffffffe,%rdx
add    $0xffffffff,%esi
cmp    $0x3,%esi
jle    1160 <mulv+0x50>
movzwl (%rdx),%eax
test   %rax,%rax
je     1190 <mulv+0x80>
movzwl (%rdi),%ebx
test   %rbx,%rbx
je     1190 <mulv+0x80>
imul   %rax,%rbx
movzwl (%r9),%eax
movzwl %bx,%ecx
add    %rax,%rcx
mov    %cx,(%r9)
shr    $0x10,%rcx
movzwl -0x2(%r9),%eax
shr    $0x10,%rbx
add    %rax,%rbx
add    %rcx,%rbx
mov    %bx,-0x2(%r9)
test   $0x10000,%ebx
je     1190 <mulv+0x80>
addw   $0x1,-0x4(%r9)
jmp    1190 <mulv+0x80>
xor    %eax,%eax
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret

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

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