int findPitchPeriodInRange(short *samples, int numChannels, int minPeriod,
                           int maxPeriod, int skip) {
  int period, bestPeriod = 0;
  short *s, *p, sVal, pVal;
  unsigned long diff, minDiff = 0;
  unsigned int numSamples = (minPeriod + skip - 1) / skip;
  unsigned int bestNumSamples = 0;
  int sampleSkip = skip * numChannels;
  int i;
  for (period = minPeriod; period <= maxPeriod; period += skip) {
    diff = 0;
    s = samples;
    p = samples + period * numChannels;
    for (i = 0; i < numSamples; i++) {
      sVal = *s;
      pVal = *p;
      s += sampleSkip;
      p += sampleSkip;
      diff += sVal >= pVal ? (unsigned short)(sVal - pVal)
                           : (unsigned short)(pVal - sVal);
    }
    if (bestPeriod == 0 || diff * bestNumSamples < minDiff * numSamples) {
      minDiff = diff;
      bestPeriod = period;
      bestNumSamples = numSamples;
    }
    numSamples++;
  }
  return bestPeriod;
}


x86 assembly (clang, -O0):
<findPitchPeriodInRange>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8d,-0x18(%rbp)
movl   $0x0,-0x20(%rbp)
movq   $0x0,-0x48(%rbp)
mov    -0x10(%rbp),%eax
add    -0x18(%rbp),%eax
sub    $0x1,%eax
cltd
idivl  -0x18(%rbp)
mov    %eax,-0x4c(%rbp)
movl   $0x0,-0x50(%rbp)
mov    -0x18(%rbp),%eax
imul   -0xc(%rbp),%eax
mov    %eax,-0x54(%rbp)
mov    -0x10(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     1277 <findPitchPeriodInRange+0x177>
movq   $0x0,-0x40(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-0x28(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x1,%rcx
add    %rcx,%rax
mov    %rax,-0x30(%rbp)
movl   $0x0,-0x58(%rbp)
mov    -0x58(%rbp),%eax
cmp    -0x4c(%rbp),%eax
jae    1223 <findPitchPeriodInRange+0x123>
mov    -0x28(%rbp),%rax
mov    (%rax),%ax
mov    %ax,-0x32(%rbp)
mov    -0x30(%rbp),%rax
mov    (%rax),%ax
mov    %ax,-0x34(%rbp)
mov    -0x54(%rbp),%ecx
mov    -0x28(%rbp),%rax
movslq %ecx,%rcx
shl    $0x1,%rcx
add    %rcx,%rax
mov    %rax,-0x28(%rbp)
mov    -0x54(%rbp),%ecx
mov    -0x30(%rbp),%rax
movslq %ecx,%rcx
shl    $0x1,%rcx
add    %rcx,%rax
mov    %rax,-0x30(%rbp)
movswl -0x32(%rbp),%eax
movswl -0x34(%rbp),%ecx
cmp    %ecx,%eax
jl     11f8 <findPitchPeriodInRange+0xf8>
movswl -0x32(%rbp),%eax
movswl -0x34(%rbp),%ecx
sub    %ecx,%eax
movzwl %ax,%eax
mov    %eax,-0x5c(%rbp)
jmp    1208 <findPitchPeriodInRange+0x108>
movswl -0x34(%rbp),%eax
movswl -0x32(%rbp),%ecx
sub    %ecx,%eax
movzwl %ax,%eax
mov    %eax,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
cltq
add    -0x40(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    1187 <findPitchPeriodInRange+0x87>
cmpl   $0x0,-0x20(%rbp)
je     124c <findPitchPeriodInRange+0x14c>
mov    -0x40(%rbp),%rax
mov    -0x50(%rbp),%ecx
imul   %rcx,%rax
mov    -0x48(%rbp),%rcx
mov    -0x4c(%rbp),%edx
imul   %rdx,%rcx
cmp    %rcx,%rax
jae    1260 <findPitchPeriodInRange+0x160>
mov    -0x40(%rbp),%rax
mov    %rax,-0x48(%rbp)
mov    -0x1c(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x4c(%rbp),%eax
mov    %eax,-0x50(%rbp)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
mov    -0x18(%rbp),%eax
add    -0x1c(%rbp),%eax
mov    %eax,-0x1c(%rbp)
jmp    114b <findPitchPeriodInRange+0x4b>
mov    -0x20(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<findPitchPeriodInRange>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
xor    %r11d,%r11d
cmp    %ecx,%edx
jle    111f <findPitchPeriodInRange+0x1f>
mov    %r11d,%eax
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    %edx,%ebp
mov    %r8d,%ebx
imul   %esi,%ebx
lea    (%r8,%rbp,1),%eax
add    $0xffffffff,%eax
cltd
idiv   %r8d
movslq %ebx,%r9
movslq %ebp,%rbx
movslq %r8d,%rdx
mov    %rdx,-0x20(%rsp)
movslq %esi,%rdx
mov    %rdx,-0x38(%rsp)
movslq %ecx,%rcx
mov    %rcx,-0x28(%rsp)
mov    %eax,%ecx
add    $0xffffffff,%ecx
mov    %ecx,-0x48(%rsp)
xor    %r10d,%r10d
xor    %r14d,%r14d
xor    %r15d,%r15d
xor    %r11d,%r11d
mov    %rdi,-0x30(%rsp)
jmp    118c <findPitchPeriodInRange+0x8c>
nopl   0x0(%rax)
mov    %ebx,%r11d
mov    %rsi,%r15
mov    %eax,%r14d
add    $0x1,%eax
add    -0x20(%rsp),%rbx
add    $0x1,%r10d
cmp    -0x28(%rsp),%rbx
jg     1111 <findPitchPeriodInRange+0x11>
test   %eax,%eax
je     11d0 <findPitchPeriodInRange+0xd0>
mov    %rbx,%rcx
imul   -0x38(%rsp),%rcx
lea    (%rdi,%rcx,2),%rbp
mov    -0x48(%rsp),%ecx
add    %r10d,%ecx
mov    %eax,-0x3c(%rsp)
mov    %rbx,-0x8(%rsp)
mov    %r10,-0x10(%rsp)
mov    %r11d,-0x40(%rsp)
mov    %r14d,-0x44(%rsp)
mov    %r15,-0x18(%rsp)
jne    11e0 <findPitchPeriodInRange+0xe0>
xor    %esi,%esi
mov    %rdi,%r13
jmp    1252 <findPitchPeriodInRange+0x152>
nopl   (%rax)
xor    %esi,%esi
jmp    1298 <findPitchPeriodInRange+0x198>
nopw   0x0(%rax,%rax,1)
mov    %eax,%r10d
and    $0xfffffffe,%r10d
xor    %esi,%esi
mov    %rdi,%r13
xor    %r14d,%r14d
nop
movzwl 0x0(%r13,%r9,2),%r8d
movzwl 0x0(%rbp,%r9,2),%r15d
mov    %r8d,%r11d
sub    %r15w,%r11w
mov    %r11d,%r12d
neg    %r12d
cmp    %r15w,%r8w
movzwl 0x0(%r13),%ebx
cmovge %r11d,%r12d
movzwl 0x0(%rbp),%ecx
mov    %ebx,%eax
sub    %cx,%ax
mov    %eax,%edx
neg    %edx
cmp    %cx,%bx
lea    0x0(%r13,%r9,2),%rcx
lea    0x0(%rbp,%r9,2),%rbp
cmovge %eax,%edx
movzwl %dx,%eax
add    %rsi,%rax
lea    (%rcx,%r9,2),%r13
lea    0x0(%rbp,%r9,2),%rbp
movzwl %r12w,%esi
add    %rax,%rsi
add    $0x2,%r14d
cmp    %r10d,%r14d
jne    11f0 <findPitchPeriodInRange+0xf0>
mov    -0x3c(%rsp),%eax
test   $0x1,%al
mov    -0x30(%rsp),%rdi
mov    -0x8(%rsp),%rbx
mov    -0x10(%rsp),%r10
mov    -0x40(%rsp),%r11d
mov    -0x44(%rsp),%r14d
mov    -0x18(%rsp),%r15
je     1298 <findPitchPeriodInRange+0x198>
movzwl 0x0(%r13),%r8d
movzwl 0x0(%rbp),%ecx
mov    %r8d,%edx
sub    %cx,%dx
mov    %edx,%ebp
neg    %ebp
cmp    %cx,%r8w
cmovge %edx,%ebp
movzwl %bp,%ecx
add    %rcx,%rsi
test   %r11d,%r11d
je     1170 <findPitchPeriodInRange+0x70>
mov    %r14d,%edx
imul   %rsi,%rdx
mov    %eax,%ecx
imul   %r15,%rcx
cmp    %rcx,%rdx
jb     1170 <findPitchPeriodInRange+0x70>
jmp    1179 <findPitchPeriodInRange+0x79>

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

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