float compute_percentile_f(float x, float *xv, int n) {
  int i;
  float p = 0.0;
  for (i = 1; i <= n; i++) {
    if (xv[i] >= x)
      p++;
  }
  p = p / (float)n;
  return p;
}


x86 assembly (clang, -O0):
<compute_percentile_f>:
push   %rbp
mov    %rsp,%rbp
movss  %xmm0,-0x4(%rbp)
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x1c(%rbp)
movl   $0x1,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     1167 <compute_percentile_f+0x67>
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
ucomiss -0x4(%rbp),%xmm0
jb     1154 <compute_percentile_f+0x54>
movss  0xeb6(%rip),%xmm0
addss  -0x1c(%rbp),%xmm0
movss  %xmm0,-0x1c(%rbp)
jmp    1159 <compute_percentile_f+0x59>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    111f <compute_percentile_f+0x1f>
movss  -0x1c(%rbp),%xmm0
cvtsi2ssl -0x14(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x1c(%rbp)
movss  -0x1c(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<compute_percentile_f>:
test   %esi,%esi
jle    112b <compute_percentile_f+0x2b>
mov    %esi,%edx
lea    -0x1(%rdx),%rcx
mov    %edx,%eax
and    $0x3,%eax
cmp    $0x3,%rcx
jae    1133 <compute_percentile_f+0x33>
xorps  %xmm1,%xmm1
mov    $0x1,%ecx
test   %rax,%rax
jne    11ca <compute_percentile_f+0xca>
jmp    1202 <compute_percentile_f+0x102>
xorps  %xmm1,%xmm1
jmp    1202 <compute_percentile_f+0x102>
and    $0xfffffffc,%edx
xorps  %xmm1,%xmm1
xor    %ecx,%ecx
movss  0xebd(%rip),%xmm2
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movaps %xmm0,%xmm3
cmpnless 0x4(%rdi,%rcx,4),%xmm3
movaps %xmm3,%xmm4
andps  %xmm1,%xmm4
addss  %xmm2,%xmm1
andnps %xmm1,%xmm3
orps   %xmm4,%xmm3
movaps %xmm0,%xmm1
cmpnless 0x8(%rdi,%rcx,4),%xmm1
movaps %xmm1,%xmm4
andps  %xmm3,%xmm4
addss  %xmm2,%xmm3
andnps %xmm3,%xmm1
orps   %xmm4,%xmm1
movaps %xmm0,%xmm3
cmpnless 0xc(%rdi,%rcx,4),%xmm3
movaps %xmm3,%xmm4
andps  %xmm1,%xmm4
addss  %xmm2,%xmm1
andnps %xmm1,%xmm3
orps   %xmm4,%xmm3
movaps %xmm0,%xmm1
cmpnless 0x10(%rdi,%rcx,4),%xmm1
movaps %xmm1,%xmm4
andps  %xmm3,%xmm4
addss  %xmm2,%xmm3
andnps %xmm3,%xmm1
orps   %xmm4,%xmm1
add    $0x4,%rcx
cmp    %rcx,%rdx
jne    1150 <compute_percentile_f+0x50>
add    $0x1,%rcx
test   %rax,%rax
je     1202 <compute_percentile_f+0x102>
lea    (%rdi,%rcx,4),%rcx
xor    %edx,%edx
movss  0xe28(%rip),%xmm2
nopl   0x0(%rax,%rax,1)
movaps %xmm1,%xmm3
addss  %xmm2,%xmm3
movaps %xmm0,%xmm4
cmpnless (%rcx,%rdx,4),%xmm4
andps  %xmm4,%xmm1
andnps %xmm3,%xmm4
orps   %xmm4,%xmm1
add    $0x1,%rdx
cmp    %rdx,%rax
jne    11e0 <compute_percentile_f+0xe0>
xorps  %xmm0,%xmm0
cvtsi2ss %esi,%xmm0
divss  %xmm0,%xmm1
movaps %xmm1,%xmm0
ret

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

LLVM IR (-O2):
define float @compute_percentile_f(float %0, float* nocapture %1, i32 %2) {
%4 = icmp slt i32 %2, 1
br i1 %4, label %61, label %5
5: ; preds = %3
%6 = add nuw i32 %2, 1
%7 = zext i32 %6 to i64
%8 = add nsw i64 %7, -1
%9 = add nsw i64 %7, -2
%10 = and i64 %8, 3
%11 = icmp ult i64 %9, 3
br i1 %11, label %44, label %12
12: ; preds = %5
%13 = and i64 %8, -4
br label %14
14: ; preds = %14, %12
%15 = phi i64 [ 1, %12 ], [ %41, %14 ]
%16 = phi float [ 0.000000e+00, %12 ], [ %40, %14 ]
%17 = phi i64 [ 0, %12 ], [ %42, %14 ]
%18 = getelementptr inbounds float, float* %1, i64 %15
%19 = load float, float* %18
%20 = fcmp ult float %19, %0
%21 = fadd float %16, 1.000000e+00
%22 = select i1 %20, float %16, float %21
%23 = add nuw nsw i64 %15, 1
%24 = getelementptr inbounds float, float* %1, i64 %23
%25 = load float, float* %24
%26 = fcmp ult float %25, %0
%27 = fadd float %22, 1.000000e+00
%28 = select i1 %26, float %22, float %27
%29 = add nuw nsw i64 %15, 2
%30 = getelementptr inbounds float, float* %1, i64 %29
%31 = load float, float* %30
%32 = fcmp ult float %31, %0
%33 = fadd float %28, 1.000000e+00
%34 = select i1 %32, float %28, float %33
%35 = add nuw nsw i64 %15, 3
%36 = getelementptr inbounds float, float* %1, i64 %35
%37 = load float, float* %36
%38 = fcmp ult float %37, %0
%39 = fadd float %34, 1.000000e+00
%40 = select i1 %38, float %34, float %39
%41 = add nuw nsw i64 %15, 4
%42 = add i64 %17, 4
%43 = icmp eq i64 %42, %13
br i1 %43, label %44, label %14
44: ; preds = %14, %5
%45 = phi float [ undef, %5 ], [ %40, %14 ]
%46 = phi i64 [ 1, %5 ], [ %41, %14 ]
%47 = phi float [ 0.000000e+00, %5 ], [ %40, %14 ]
%48 = icmp eq i64 %10, 0
br i1 %48, label %61, label %49
49: ; preds = %44, %49
%50 = phi i64 [ %58, %49 ], [ %46, %44 ]
%51 = phi float [ %57, %49 ], [ %47, %44 ]
%52 = phi i64 [ %59, %49 ], [ 0, %44 ]
%53 = getelementptr inbounds float, float* %1, i64 %50
%54 = load float, float* %53
%55 = fcmp ult float %54, %0
%56 = fadd float %51, 1.000000e+00
%57 = select i1 %55, float %51, float %56
%58 = add nuw nsw i64 %50, 1
%59 = add i64 %52, 1
%60 = icmp eq i64 %59, %10
br i1 %60, label %61, label %49
61: ; preds = %44, %49, %3
%62 = phi float [ 0.000000e+00, %3 ], [ %45, %44 ], [ %57, %49 ]
%63 = sitofp i32 %2 to float
%64 = fdiv float %62, %63
ret float %64
}