int getNextKMerDiscrepancy(double *kmerDiscrepancies, int total) {
  double score = -1;
  int kmerLocation = -1;
  for (int i = 0; i < total; i++) {
    if (kmerDiscrepancies[i] > 0 && kmerDiscrepancies[i] > score) {
      score = kmerDiscrepancies[i];
      kmerLocation = i;
    }
  }
  return kmerLocation;
}


x86 assembly (clang, -O0):
<getNextKMerDiscrepancy>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
movsd  0xeed(%rip),%xmm0
movsd  %xmm0,-0x18(%rbp)
movl   $0xffffffff,-0x1c(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    118f <getNextKMerDiscrepancy+0x8f>
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    117c <getNextKMerDiscrepancy+0x7c>
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x18(%rbp),%xmm0
jbe    117c <getNextKMerDiscrepancy+0x7c>
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x20(%rbp),%eax
mov    %eax,-0x1c(%rbp)
jmp    1181 <getNextKMerDiscrepancy+0x81>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1126 <getNextKMerDiscrepancy+0x26>
mov    -0x1c(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<getNextKMerDiscrepancy>:
test   %esi,%esi
jle    113b <getNextKMerDiscrepancy+0x3b>
mov    %esi,%r8d
cmp    $0x1,%esi
jne    1141 <getNextKMerDiscrepancy+0x41>
movsd  0xeec(%rip),%xmm0
mov    $0xffffffff,%eax
xor    %edx,%edx
test   $0x1,%r8b
je     113a <getNextKMerDiscrepancy+0x3a>
movsd  (%rdi,%rdx,8),%xmm1
ucomisd %xmm0,%xmm1
cmovbe %eax,%edx
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm1
cmovbe %eax,%edx
mov    %edx,%eax
ret
mov    $0xffffffff,%eax
ret
mov    %r8d,%r10d
and    $0xfffffffe,%r10d
movsd  0xeb0(%rip),%xmm0
mov    $0xffffffff,%eax
xor    %edx,%edx
xorpd  %xmm1,%xmm1
jmp    116c <getNextKMerDiscrepancy+0x6c>
nopl   (%rax)
cmovne %esi,%eax
add    $0x2,%rdx
cmp    %r10,%rdx
je     111b <getNextKMerDiscrepancy+0x1b>
movsd  (%rdi,%rdx,8),%xmm2
ucomisd %xmm1,%xmm2
seta   %r9b
ucomisd %xmm0,%xmm2
seta   %cl
test   %cl,%r9b
jne    1189 <getNextKMerDiscrepancy+0x89>
movapd %xmm0,%xmm2
movsd  0x8(%rdi,%rdx,8),%xmm0
cmovne %edx,%eax
ucomisd %xmm1,%xmm0
seta   %r9b
ucomisd %xmm2,%xmm0
seta   %cl
lea    0x1(%rdx),%esi
test   %cl,%r9b
jne    1160 <getNextKMerDiscrepancy+0x60>
movapd %xmm2,%xmm0
jmp    1160 <getNextKMerDiscrepancy+0x60>

LLVM IR (-O0):
define i32 @getNextKMerDiscrepancy(double* %0, i32 %1) {
%3 = alloca double*
%4 = alloca i32
%5 = alloca double
%6 = alloca i32
%7 = alloca i32
store double* %0, double** %3
store i32 %1, i32* %4
store double -1.000000e+00, double* %5
store i32 -1, i32* %6
store i32 0, i32* %7
br label %8
8: ; preds = %35, %2
%9 = load i32, i32* %7
%10 = load i32, i32* %4
%11 = icmp slt i32 %9, %10
br i1 %11, label %12, label %38
12: ; preds = %8
%13 = load double*, double** %3
%14 = load i32, i32* %7
%15 = sext i32 %14 to i64
%16 = getelementptr inbounds double, double* %13, i64 %15
%17 = load double, double* %16
%18 = fcmp ogt double %17, 0.000000e+00
br i1 %18, label %19, label %34
19: ; preds = %12
%20 = load double*, double** %3
%21 = load i32, i32* %7
%22 = sext i32 %21 to i64
%23 = getelementptr inbounds double, double* %20, i64 %22
%24 = load double, double* %23
%25 = load double, double* %5
%26 = fcmp ogt double %24, %25
br i1 %26, label %27, label %34
27: ; preds = %19
%28 = load double*, double** %3
%29 = load i32, i32* %7
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds double, double* %28, i64 %30
%32 = load double, double* %31
store double %32, double* %5
%33 = load i32, i32* %7
store i32 %33, i32* %6
br label %34
34: ; preds = %27, %19, %12
br label %35
35: ; preds = %34
%36 = load i32, i32* %7
%37 = add nsw i32 %36, 1
store i32 %37, i32* %7
br label %8
38: ; preds = %8
%39 = load i32, i32* %6
ret i32 %39
}

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