void mpsad(unsigned char w[], unsigned char u[], int n, int iv) {
  int j;
  unsigned short ireg;
  ireg = 256 * iv;
  for (j = n; j >= 1; j--) {
    ireg = u[j] + ((unsigned char)((ireg) >> 8 & 0xff));
    w[j + 1] = ((unsigned char)((ireg)&0xff));
  }
  w[1] = ((unsigned char)((ireg) >> 8 & 0xff));
}


x86 assembly (clang, -O0):
<mpsad>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    -0x18(%rbp),%eax
shl    $0x8,%eax
mov    %ax,-0x1e(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0x1c(%rbp)
cmpl   $0x1,-0x1c(%rbp)
jl     1177 <mpsad+0x77>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movzbl (%rax,%rcx,1),%eax
movzwl -0x1e(%rbp),%ecx
sar    $0x8,%ecx
and    $0xff,%ecx
movzbl %cl,%ecx
add    %ecx,%eax
mov    %ax,-0x1e(%rbp)
movzwl -0x1e(%rbp),%eax
and    $0xff,%eax
mov    %al,%dl
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    %dl,(%rax,%rcx,1)
mov    -0x1c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x1c(%rbp)
jmp    1122 <mpsad+0x22>
movzwl -0x1e(%rbp),%eax
sar    $0x8,%eax
and    $0xff,%eax
mov    %al,%cl
mov    -0x8(%rbp),%rax
mov    %cl,0x1(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<mpsad>:
push   %rbx
mov    %ecx,%ebx
shl    $0x8,%ebx
test   %edx,%edx
jle    1158 <mpsad+0x58>
mov    %edx,%eax
test   $0x1,%al
je     1123 <mpsad+0x23>
movzbl (%rsi,%rax,1),%r8d
movzbl %cl,%ebx
add    %r8d,%ebx
mov    %bl,0x1(%rax,%rdi,1)
add    $0xffffffffffffffff,%rax
cmp    $0x1,%edx
je     1158 <mpsad+0x58>
add    $0x1,%rax
nopl   0x0(%rax)
movzbl -0x1(%rsi,%rax,1),%edx
movzwl %bx,%ecx
shr    $0x8,%ecx
add    %edx,%ecx
mov    %cl,(%rdi,%rax,1)
movzbl -0x2(%rsi,%rax,1),%edx
movzbl %ch,%ebx
add    %edx,%ebx
mov    %bl,-0x1(%rdi,%rax,1)
add    $0xfffffffffffffffe,%rax
cmp    $0x1,%rax
ja     1130 <mpsad+0x30>
mov    %bh,0x1(%rdi)
pop    %rbx
ret

LLVM IR (-O0):
define void @mpsad(i8* %0, i8* %1, i32 %2, i32 %3) {
%5 = alloca i8*
%6 = alloca i8*
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i16
store i8* %0, i8** %5
store i8* %1, i8** %6
store i32 %2, i32* %7
store i32 %3, i32* %8
%11 = load i32, i32* %8
%12 = mul nsw i32 256, %11
%13 = trunc i32 %12 to i16
store i16 %13, i16* %10
%14 = load i32, i32* %7
store i32 %14, i32* %9
br label %15
15: ; preds = %42, %4
%16 = load i32, i32* %9
%17 = icmp sge i32 %16, 1
br i1 %17, label %18, label %45
18: ; preds = %15
%19 = load i8*, i8** %6
%20 = load i32, i32* %9
%21 = sext i32 %20 to i64
%22 = getelementptr inbounds i8, i8* %19, i64 %21
%23 = load i8, i8* %22
%24 = zext i8 %23 to i32
%25 = load i16, i16* %10
%26 = zext i16 %25 to i32
%27 = ashr i32 %26, 8
%28 = and i32 %27, 255
%29 = trunc i32 %28 to i8
%30 = zext i8 %29 to i32
%31 = add nsw i32 %24, %30
%32 = trunc i32 %31 to i16
store i16 %32, i16* %10
%33 = load i16, i16* %10
%34 = zext i16 %33 to i32
%35 = and i32 %34, 255
%36 = trunc i32 %35 to i8
%37 = load i8*, i8** %5
%38 = load i32, i32* %9
%39 = add nsw i32 %38, 1
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds i8, i8* %37, i64 %40
store i8 %36, i8* %41
br label %42
42: ; preds = %18
%43 = load i32, i32* %9
%44 = add nsw i32 %43, -1
store i32 %44, i32* %9
br label %15
45: ; preds = %15
%46 = load i16, i16* %10
%47 = zext i16 %46 to i32
%48 = ashr i32 %47, 8
%49 = and i32 %48, 255
%50 = trunc i32 %49 to i8
%51 = load i8*, i8** %5
%52 = getelementptr inbounds i8, i8* %51, i64 1
store i8 %50, i8* %52
ret void
}

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