int findPitchPeriodInRange(short *samples, int minPeriod, int maxPeriod,
                           int *retMinDiff, int *retMaxDiff) {
  int period, bestPeriod = 0, worstPeriod = 255;
  short *s, *p, sVal, pVal;
  unsigned long diff, minDiff = 1, maxDiff = 0;
  int i;
  for (period = minPeriod; period <= maxPeriod; period++) {
    diff = 0;
    s = samples;
    p = samples + period;
    for (i = 0; i < period; i++) {
      sVal = *s++;
      pVal = *p++;
      diff += sVal >= pVal ? (unsigned short)(sVal - pVal)
                           : (unsigned short)(pVal - sVal);
    }
    if (bestPeriod == 0 || diff * bestPeriod < minDiff * period) {
      minDiff = diff;
      bestPeriod = period;
    }
    if (diff * worstPeriod > maxDiff * period) {
      maxDiff = diff;
      worstPeriod = period;
    }
  }
  *retMinDiff = minDiff / bestPeriod;
  *retMaxDiff = maxDiff / worstPeriod;
  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    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
movl   $0x0,-0x28(%rbp)
movl   $0xff,-0x2c(%rbp)
movq   $0x1,-0x58(%rbp)
movq   $0x0,-0x60(%rbp)
mov    -0xc(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x10(%rbp),%eax
jg     1273 <findPitchPeriodInRange+0x173>
movq   $0x0,-0x50(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
shl    $0x1,%rcx
add    %rcx,%rax
mov    %rax,-0x40(%rbp)
movl   $0x0,-0x64(%rbp)
mov    -0x64(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    11f8 <findPitchPeriodInRange+0xf8>
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x2,%rcx
mov    %rcx,-0x38(%rbp)
mov    (%rax),%ax
mov    %ax,-0x42(%rbp)
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x2,%rcx
mov    %rcx,-0x40(%rbp)
mov    (%rax),%ax
mov    %ax,-0x44(%rbp)
movswl -0x42(%rbp),%eax
movswl -0x44(%rbp),%ecx
cmp    %ecx,%eax
jl     11cd <findPitchPeriodInRange+0xcd>
movswl -0x42(%rbp),%eax
movswl -0x44(%rbp),%ecx
sub    %ecx,%eax
movzwl %ax,%eax
mov    %eax,-0x68(%rbp)
jmp    11dd <findPitchPeriodInRange+0xdd>
movswl -0x44(%rbp),%eax
movswl -0x42(%rbp),%ecx
sub    %ecx,%eax
movzwl %ax,%eax
mov    %eax,-0x68(%rbp)
mov    -0x68(%rbp),%eax
cltq
add    -0x50(%rbp),%rax
mov    %rax,-0x50(%rbp)
mov    -0x64(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
jmp    1170 <findPitchPeriodInRange+0x70>
cmpl   $0x0,-0x28(%rbp)
je     1223 <findPitchPeriodInRange+0x123>
mov    -0x50(%rbp),%rax
movslq -0x28(%rbp),%rcx
imul   %rcx,%rax
mov    -0x58(%rbp),%rcx
movslq -0x24(%rbp),%rdx
imul   %rdx,%rcx
cmp    %rcx,%rax
jae    1231 <findPitchPeriodInRange+0x131>
mov    -0x50(%rbp),%rax
mov    %rax,-0x58(%rbp)
mov    -0x24(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x50(%rbp),%rax
movslq -0x2c(%rbp),%rcx
imul   %rcx,%rax
mov    -0x60(%rbp),%rcx
movslq -0x24(%rbp),%rdx
imul   %rdx,%rcx
cmp    %rcx,%rax
jbe    1260 <findPitchPeriodInRange+0x160>
mov    -0x50(%rbp),%rax
mov    %rax,-0x60(%rbp)
mov    -0x24(%rbp),%eax
mov    %eax,-0x2c(%rbp)
jmp    1265 <findPitchPeriodInRange+0x165>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    113a <findPitchPeriodInRange+0x3a>
mov    -0x58(%rbp),%rax
movslq -0x28(%rbp),%rcx
xor    %edx,%edx
div    %rcx
mov    %eax,%ecx
mov    -0x18(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x60(%rbp),%rax
movslq -0x2c(%rbp),%rcx
xor    %edx,%edx
div    %rcx
mov    %eax,%ecx
mov    -0x20(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x28(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<findPitchPeriodInRange>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
cmp    %edx,%esi
jle    111c <findPitchPeriodInRange+0x1c>
mov    $0x1,%eax
xor    %ebp,%ebp
xor    %edi,%edi
jmp    1340 <findPitchPeriodInRange+0x240>
mov    %rcx,-0x28(%rsp)
mov    %r8,-0x20(%rsp)
movslq %esi,%r15
add    $0x1,%edx
mov    %edx,-0x30(%rsp)
add    $0xffffffff,%esi
mov    %rdi,-0x10(%rsp)
lea    0x4(%rdi),%rax
mov    %rax,-0x18(%rsp)
lea    (%r15,%r15,1),%r14
mov    $0xff,%r11d
mov    $0x1,%eax
mov    %rax,-0x38(%rsp)
xor    %r12d,%r12d
pxor   %xmm8,%xmm8
mov    %rsi,-0x8(%rsp)
mov    %esi,%r13d
xor    %eax,%eax
xor    %ebp,%ebp
jmp    11af <findPitchPeriodInRange+0xaf>
nopl   0x0(%rax,%rax,1)
mov    %r15d,%ebp
mov    %r9,-0x38(%rsp)
movslq %r11d,%r11
mov    %r9,%rcx
imul   %r11,%rcx
mov    %rax,%rdx
imul   %r15,%rdx
cmp    %rdx,%rcx
cmova  %r15d,%r11d
cmova  %r9,%rax
add    $0x1,%r15
add    $0x1,%r12d
add    $0x2,%r14
add    $0x1,%r13d
cmp    %r15d,-0x30(%rsp)
je     1311 <findPitchPeriodInRange+0x211>
test   %r15,%r15
mov    %ebp,-0x2c(%rsp)
jle    11e0 <findPitchPeriodInRange+0xe0>
mov    -0x8(%rsp),%rcx
lea    (%rcx,%r12,1),%edi
mov    -0x10(%rsp),%r8
lea    (%r8,%r15,2),%r10
cmp    $0x3,%edi
jae    11f0 <findPitchPeriodInRange+0xf0>
xor    %edx,%edx
xor    %r9d,%r9d
jmp    12b8 <findPitchPeriodInRange+0x1b8>
nopl   0x0(%rax)
xor    %r9d,%r9d
jmp    12e7 <findPitchPeriodInRange+0x1e7>
nopl   0x0(%rax,%rax,1)
mov    %r13d,%esi
add    $0x1,%rsi
and    $0xfffffffffffffffc,%rsi
add    $0x1,%rdi
mov    %rdi,%rdx
and    $0xfffffffffffffffc,%rdx
lea    (%r10,%rdx,2),%r10
lea    (%r8,%rdx,2),%r8
pxor   %xmm9,%xmm9
mov    -0x18(%rsp),%rcx
pxor   %xmm10,%xmm10
nopl   (%rax)
movd   -0x4(%rcx),%xmm1
movd   (%rcx),%xmm5
movd   -0x4(%rcx,%r14,1),%xmm0
movd   (%rcx,%r14,1),%xmm7
movdqa %xmm0,%xmm6
pcmpgtw %xmm1,%xmm6
movdqa %xmm7,%xmm3
pcmpgtw %xmm5,%xmm3
movdqa %xmm1,%xmm4
psubw  %xmm0,%xmm4
movdqa %xmm5,%xmm2
psubw  %xmm7,%xmm2
psubw  %xmm1,%xmm0
psubw  %xmm5,%xmm7
pand   %xmm6,%xmm0
pandn  %xmm4,%xmm6
por    %xmm0,%xmm6
pand   %xmm3,%xmm7
pandn  %xmm2,%xmm3
por    %xmm7,%xmm3
punpcklwd %xmm8,%xmm6
punpckldq %xmm8,%xmm6
paddq  %xmm6,%xmm9
punpcklwd %xmm8,%xmm3
punpckldq %xmm8,%xmm3
paddq  %xmm3,%xmm10
add    $0x8,%rcx
add    $0xfffffffffffffffc,%rsi
jne    1220 <findPitchPeriodInRange+0x120>
paddq  %xmm9,%xmm10
pshufd $0xee,%xmm10,%xmm0
paddq  %xmm10,%xmm0
movq   %xmm0,%r9
cmp    %rdx,%rdi
je     12e7 <findPitchPeriodInRange+0x1e7>
mov    %r15d,%ecx
sub    %edx,%ecx
xor    %esi,%esi
nop
movzwl (%r8,%rsi,2),%edx
movzwl (%r10,%rsi,2),%ebp
mov    %edx,%ebx
sub    %bp,%bx
mov    %ebx,%edi
neg    %edi
cmp    %bp,%dx
cmovge %ebx,%edi
movzwl %di,%edx
add    %rdx,%r9
add    $0x1,%rsi
cmp    %esi,%ecx
jne    12c0 <findPitchPeriodInRange+0x1c0>
mov    -0x2c(%rsp),%ebp
test   %ebp,%ebp
je     1170 <findPitchPeriodInRange+0x70>
movslq %ebp,%rcx
imul   %r9,%rcx
mov    -0x38(%rsp),%rdx
imul   %r15,%rdx
cmp    %rdx,%rcx
jb     1170 <findPitchPeriodInRange+0x70>
jmp    1178 <findPitchPeriodInRange+0x78>
movslq %r11d,%rcx
mov    %rax,%rdx
or     %rcx,%rdx
shr    $0x20,%rdx
je     132a <findPitchPeriodInRange+0x22a>
xor    %edx,%edx
div    %rcx
mov    %rax,%rdi
jmp    1331 <findPitchPeriodInRange+0x231>
xor    %edx,%edx
div    %r11d
mov    %eax,%edi
mov    -0x20(%rsp),%r8
mov    -0x28(%rsp),%rcx
mov    -0x38(%rsp),%rax
movslq %ebp,%rsi
mov    %rax,%rdx
or     %rsi,%rdx
shr    $0x20,%rdx
je     1356 <findPitchPeriodInRange+0x256>
xor    %edx,%edx
div    %rsi
jmp    135a <findPitchPeriodInRange+0x25a>
xor    %edx,%edx
div    %ebp
mov    %eax,(%rcx)
mov    %edi,(%r8)
mov    %ebp,%eax
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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