int FindNextMax(double *pFMagSpectrum, int iHighBinBound, int *pICurrentLoc,
                double *pFMaxVal, double fMinPeakMag) {
  int iCurrentBin = *pICurrentLoc;
  double fPrevVal = pFMagSpectrum[iCurrentBin - 1];
  double fCurrentVal = pFMagSpectrum[iCurrentBin];
  double fNextVal =
      (iCurrentBin >= iHighBinBound) ? 0 : pFMagSpectrum[iCurrentBin + 1];
  while (iCurrentBin <= iHighBinBound) {
    if (fCurrentVal > fMinPeakMag && fCurrentVal >= fPrevVal &&
        fCurrentVal >= fNextVal)
      break;
    iCurrentBin++;
    fPrevVal = fCurrentVal;
    fCurrentVal = fNextVal;
    fNextVal = pFMagSpectrum[1 + iCurrentBin];
  }
  *pICurrentLoc = iCurrentBin + 1;
  *pFMaxVal = fCurrentVal;
  return iCurrentBin;
}


x86 assembly (clang, -O0):
<FindNextMax>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movsd  %xmm0,-0x28(%rbp)
mov    -0x18(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jl     1163 <FindNextMax+0x63>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
jmp    117a <FindNextMax+0x7a>
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     11fe <FindNextMax+0xfe>
movsd  -0x40(%rbp),%xmm0
ucomisd -0x28(%rbp),%xmm0
jbe    11c5 <FindNextMax+0xc5>
movsd  -0x40(%rbp),%xmm0
ucomisd -0x38(%rbp),%xmm0
jb     11c5 <FindNextMax+0xc5>
movsd  -0x40(%rbp),%xmm0
ucomisd -0x48(%rbp),%xmm0
jb     11c5 <FindNextMax+0xc5>
jmp    11fe <FindNextMax+0xfe>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
jmp    1184 <FindNextMax+0x84>
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
mov    -0x18(%rbp),%rax
mov    %ecx,(%rax)
movsd  -0x40(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x2c(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<FindNextMax>:
mov    (%rdx),%eax
movslq %eax,%r8
xorpd  %xmm3,%xmm3
cmp    %esi,%r8d
jge    1115 <FindNextMax+0x15>
movsd  0x8(%rdi,%r8,8),%xmm3
movsd  (%rdi,%r8,8),%xmm1
cmp    %esi,%r8d
jg     1175 <FindNextMax+0x75>
movsd  -0x8(%rdi,%r8,8),%xmm2
add    $0x1,%esi
lea    (%rdi,%r8,8),%rdi
add    $0x10,%rdi
jmp    114f <FindNextMax+0x4f>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rdi),%xmm3
add    $0x1,%eax
add    $0x8,%rdi
cmp    %eax,%esi
je     1173 <FindNextMax+0x73>
movapd %xmm2,%xmm4
movapd %xmm1,%xmm2
movapd %xmm3,%xmm1
ucomisd %xmm0,%xmm2
jbe    1140 <FindNextMax+0x40>
ucomisd %xmm4,%xmm2
jb     1140 <FindNextMax+0x40>
ucomisd %xmm1,%xmm2
jb     1140 <FindNextMax+0x40>
movapd %xmm2,%xmm1
jmp    1175 <FindNextMax+0x75>
mov    %esi,%eax
lea    0x1(%rax),%esi
mov    %esi,(%rdx)
movsd  %xmm1,(%rcx)
ret

LLVM IR (-O0):
define i32 @FindNextMax(double* %0, i32 %1, i32* %2, double* %3, double %4) {
%6 = alloca double*
%7 = alloca i32
%8 = alloca i32*
%9 = alloca double*
%10 = alloca double
%11 = alloca i32
%12 = alloca double
%13 = alloca double
%14 = alloca double
store double* %0, double** %6
store i32 %1, i32* %7
store i32* %2, i32** %8
store double* %3, double** %9
store double %4, double* %10
%15 = load i32*, i32** %8
%16 = load i32, i32* %15
store i32 %16, i32* %11
%17 = load double*, double** %6
%18 = load i32, i32* %11
%19 = sub nsw i32 %18, 1
%20 = sext i32 %19 to i64
%21 = getelementptr inbounds double, double* %17, i64 %20
%22 = load double, double* %21
store double %22, double* %12
%23 = load double*, double** %6
%24 = load i32, i32* %11
%25 = sext i32 %24 to i64
%26 = getelementptr inbounds double, double* %23, i64 %25
%27 = load double, double* %26
store double %27, double* %13
%28 = load i32, i32* %11
%29 = load i32, i32* %7
%30 = icmp sge i32 %28, %29
br i1 %30, label %31, label %32
31: ; preds = %5
br label %39
32: ; preds = %5
%33 = load double*, double** %6
%34 = load i32, i32* %11
%35 = add nsw i32 %34, 1
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds double, double* %33, i64 %36
%38 = load double, double* %37
br label %39
39: ; preds = %32, %31
%40 = phi double [ 0.000000e+00, %31 ], [ %38, %32 ]
store double %40, double* %14
br label %41
41: ; preds = %58, %39
%42 = load i32, i32* %11
%43 = load i32, i32* %7
%44 = icmp sle i32 %42, %43
br i1 %44, label %45, label %69
45: ; preds = %41
%46 = load double, double* %13
%47 = load double, double* %10
%48 = fcmp ogt double %46, %47
br i1 %48, label %49, label %58
49: ; preds = %45
%50 = load double, double* %13
%51 = load double, double* %12
%52 = fcmp oge double %50, %51
br i1 %52, label %53, label %58
53: ; preds = %49
%54 = load double, double* %13
%55 = load double, double* %14
%56 = fcmp oge double %54, %55
br i1 %56, label %57, label %58
57: ; preds = %53
br label %69
58: ; preds = %53, %49, %45
%59 = load i32, i32* %11
%60 = add nsw i32 %59, 1
store i32 %60, i32* %11
%61 = load double, double* %13
store double %61, double* %12
%62 = load double, double* %14
store double %62, double* %13
%63 = load double*, double** %6
%64 = load i32, i32* %11
%65 = add nsw i32 1, %64
%66 = sext i32 %65 to i64
%67 = getelementptr inbounds double, double* %63, i64 %66
%68 = load double, double* %67
store double %68, double* %14
br label %41
69: ; preds = %57, %41
%70 = load i32, i32* %11
%71 = add nsw i32 %70, 1
%72 = load i32*, i32** %8
store i32 %71, i32* %72
%73 = load double, double* %13
%74 = load double*, double** %9
store double %73, double* %74
%75 = load i32, i32* %11
ret i32 %75
}

LLVM IR (-O2):
define i32 @FindNextMax(double* nocapture %0, i32 %1, i32* nocapture %2, double* nocapture writeonly %3, double %4) {
%6 = load i32, i32* %2
%7 = add nsw i32 %6, -1
%8 = sext i32 %7 to i64
%9 = getelementptr inbounds double, double* %0, i64 %8
%10 = load double, double* %9
%11 = sext i32 %6 to i64
%12 = getelementptr inbounds double, double* %0, i64 %11
%13 = load double, double* %12
%14 = icmp slt i32 %6, %1
br i1 %14, label %15, label %20
15: ; preds = %5
%16 = add nsw i32 %6, 1
%17 = sext i32 %16 to i64
%18 = getelementptr inbounds double, double* %0, i64 %17
%19 = load double, double* %18
br label %20
20: ; preds = %5, %15
%21 = phi double [ %19, %15 ], [ 0.000000e+00, %5 ]
%22 = icmp sgt i32 %6, %1
br i1 %22, label %44, label %23
23: ; preds = %20
%24 = add i32 %1, 1
br label %25
25: ; preds = %23, %35
%26 = phi i64 [ %11, %23 ], [ %36, %35 ]
%27 = phi double [ %21, %23 ], [ %39, %35 ]
%28 = phi double [ %13, %23 ], [ %27, %35 ]
%29 = phi double [ %10, %23 ], [ %28, %35 ]
%30 = fcmp ule double %28, %4
%31 = fcmp ult double %28, %29
%32 = select i1 %30, i1 true, i1 %31
%33 = fcmp ult double %28, %27
%34 = select i1 %32, i1 true, i1 %33
br i1 %34, label %35, label %42
35: ; preds = %25
%36 = add nsw i64 %26, 1
%37 = add nsw i64 %26, 2
%38 = getelementptr inbounds double, double* %0, i64 %37
%39 = load double, double* %38
%40 = trunc i64 %36 to i32
%41 = icmp eq i32 %24, %40
br i1 %41, label %44, label %25
42: ; preds = %25
%43 = trunc i64 %26 to i32
br label %44
44: ; preds = %35, %42, %20
%45 = phi i32 [ %6, %20 ], [ %43, %42 ], [ %24, %35 ]
%46 = phi double [ %13, %20 ], [ %28, %42 ], [ %27, %35 ]
%47 = add nsw i32 %45, 1
store i32 %47, i32* %2
store double %46, double* %3
ret i32 %45
}