void quadct(float x, float y, float xx[], float yy[], unsigned long nn,
            float *fa, float *fb, float *fc, float *fd) {
  unsigned long k, na, nb, nc, nd;
  float ff;
  na = nb = nc = nd = 0;
  for (k = 1; k <= nn; k++) {
    if (yy[k] > y) {
      xx[k] > x ? ++na : ++nb;
    } else {
      xx[k] > x ? ++nd : ++nc;
    }
  }
  ff = 1.0 / nn;
  *fa = ff * na;
  *fb = ff * nb;
  *fc = ff * nc;
  *fd = ff * nd;
}


x86 assembly (clang, -O0):
<quadct>:
push   %rbp
mov    %rsp,%rbp
sub    $0x18,%rsp
mov    0x10(%rbp),%rax
movss  %xmm0,-0x4(%rbp)
movss  %xmm1,-0x8(%rbp)
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
mov    %r8,-0x30(%rbp)
mov    %r9,-0x38(%rbp)
movq   $0x0,-0x60(%rbp)
movq   $0x0,-0x58(%rbp)
movq   $0x0,-0x50(%rbp)
movq   $0x0,-0x48(%rbp)
movq   $0x1,-0x40(%rbp)
mov    -0x40(%rbp),%rax
cmp    -0x20(%rbp),%rax
ja     11fe <quadct+0xfe>
mov    -0x18(%rbp),%rax
mov    -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
ucomiss -0x8(%rbp),%xmm0
jbe    11b4 <quadct+0xb4>
mov    -0x10(%rbp),%rax
mov    -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
ucomiss -0x4(%rbp),%xmm0
jbe    11a3 <quadct+0xa3>
mov    -0x48(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x48(%rbp)
jmp    11af <quadct+0xaf>
mov    -0x50(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x50(%rbp)
jmp    11e8 <quadct+0xe8>
mov    -0x10(%rbp),%rax
mov    -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
ucomiss -0x4(%rbp),%xmm0
jbe    11dc <quadct+0xdc>
mov    -0x60(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x60(%rbp)
jmp    11e8 <quadct+0xe8>
mov    -0x58(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x58(%rbp)
jmp    11ed <quadct+0xed>
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    1156 <quadct+0x56>
movaps 0xdfb(%rip),%xmm0
movq   -0x20(%rbp),%xmm1
punpckldq %xmm0,%xmm1
movapd 0xdfa(%rip),%xmm0
subpd  %xmm0,%xmm1
movaps %xmm1,%xmm0
unpckhpd %xmm1,%xmm1
addsd  %xmm0,%xmm1
movsd  0xdf3(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x64(%rbp)
movss  -0x64(%rbp),%xmm0
movss  %xmm0,-0x70(%rbp)
mov    -0x48(%rbp),%rax
mov    %rax,%rdx
shr    %rdx
mov    %eax,%ecx
and    $0x1,%ecx
or     %rdx,%rcx
cvtsi2ss %rcx,%xmm0
addss  %xmm0,%xmm0
cvtsi2ss %rax,%xmm1
movss  %xmm1,-0x6c(%rbp)
test   %rax,%rax
movss  %xmm0,-0x68(%rbp)
js     1281 <quadct+0x181>
movss  -0x6c(%rbp),%xmm0
movss  %xmm0,-0x68(%rbp)
movss  -0x70(%rbp),%xmm0
movss  -0x68(%rbp),%xmm1
mulss  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movss  %xmm0,(%rax)
movss  -0x64(%rbp),%xmm0
movss  %xmm0,-0x7c(%rbp)
mov    -0x50(%rbp),%rax
mov    %rax,%rdx
shr    %rdx
mov    %eax,%ecx
and    $0x1,%ecx
or     %rdx,%rcx
cvtsi2ss %rcx,%xmm0
addss  %xmm0,%xmm0
cvtsi2ss %rax,%xmm1
movss  %xmm1,-0x78(%rbp)
test   %rax,%rax
movss  %xmm0,-0x74(%rbp)
js     12de <quadct+0x1de>
movss  -0x78(%rbp),%xmm0
movss  %xmm0,-0x74(%rbp)
movss  -0x7c(%rbp),%xmm0
movss  -0x74(%rbp),%xmm1
mulss  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
movss  %xmm0,(%rax)
movss  -0x64(%rbp),%xmm0
movss  %xmm0,-0x88(%rbp)
mov    -0x58(%rbp),%rax
mov    %rax,%rdx
shr    %rdx
mov    %eax,%ecx
and    $0x1,%ecx
or     %rdx,%rcx
cvtsi2ss %rcx,%xmm0
addss  %xmm0,%xmm0
cvtsi2ss %rax,%xmm1
movss  %xmm1,-0x84(%rbp)
ff
test   %rax,%rax
movss  %xmm0,-0x80(%rbp)
js     1344 <quadct+0x244>
movss  -0x84(%rbp),%xmm0
ff
movss  %xmm0,-0x80(%rbp)
movss  -0x88(%rbp),%xmm0
ff
movss  -0x80(%rbp),%xmm1
mulss  %xmm1,%xmm0
mov    -0x38(%rbp),%rax
movss  %xmm0,(%rax)
movss  -0x64(%rbp),%xmm0
movss  %xmm0,-0x94(%rbp)
ff
mov    -0x60(%rbp),%rax
mov    %rax,%rdx
shr    %rdx
mov    %eax,%ecx
and    $0x1,%ecx
or     %rdx,%rcx
cvtsi2ss %rcx,%xmm0
addss  %xmm0,%xmm0
cvtsi2ss %rax,%xmm1
movss  %xmm1,-0x90(%rbp)
ff
test   %rax,%rax
movss  %xmm0,-0x8c(%rbp)
ff
js     13b3 <quadct+0x2b3>
movss  -0x90(%rbp),%xmm0
ff
movss  %xmm0,-0x8c(%rbp)
ff
movss  -0x94(%rbp),%xmm0
ff
movss  -0x8c(%rbp),%xmm1
ff
mulss  %xmm1,%xmm0
mov    0x10(%rbp),%rax
movss  %xmm0,(%rax)
add    $0x18,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<quadct>:
push   %r14
push   %rbx
test   %rdx,%rdx
je     1166 <quadct+0x66>
xor    %eax,%eax
xor    %r10d,%r10d
xor    %r11d,%r11d
xor    %r14d,%r14d
xor    %ebx,%ebx
jmp    112d <quadct+0x2d>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r11
add    $0x1,%rax
cmp    %rax,%rdx
je     1171 <quadct+0x71>
movss  0x4(%rsi,%rax,4),%xmm2
ucomiss %xmm1,%xmm2
movss  0x4(%rdi,%rax,4),%xmm2
jbe    1150 <quadct+0x50>
ucomiss %xmm0,%xmm2
jbe    1160 <quadct+0x60>
add    $0x1,%rbx
jmp    1124 <quadct+0x24>
nopl   0x0(%rax)
ucomiss %xmm0,%xmm2
jbe    1120 <quadct+0x20>
add    $0x1,%r10
jmp    1124 <quadct+0x24>
nopl   0x0(%rax,%rax,1)
add    $0x1,%r14
jmp    1124 <quadct+0x24>
xor    %ebx,%ebx
xor    %r14d,%r14d
xor    %r11d,%r11d
xor    %r10d,%r10d
movq   %rdx,%xmm0
punpckldq 0xe82(%rip),%xmm0
subpd  0xe8a(%rip),%xmm0
movapd %xmm0,%xmm1
unpckhpd %xmm0,%xmm1
addsd  %xmm0,%xmm1
movsd  0xe86(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
test   %rbx,%rbx
js     11b1 <quadct+0xb1>
xorps  %xmm1,%xmm1
cvtsi2ss %rbx,%xmm1
jmp    11c9 <quadct+0xc9>
mov    %rbx,%rax
shr    %rax
and    $0x1,%ebx
or     %rax,%rbx
xorps  %xmm1,%xmm1
cvtsi2ss %rbx,%xmm1
addss  %xmm1,%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,(%rcx)
test   %r14,%r14
js     11e0 <quadct+0xe0>
xorps  %xmm1,%xmm1
cvtsi2ss %r14,%xmm1
jmp    11f9 <quadct+0xf9>
mov    %r14,%rax
shr    %rax
and    $0x1,%r14d
or     %rax,%r14
xorps  %xmm1,%xmm1
cvtsi2ss %r14,%xmm1
addss  %xmm1,%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,(%r8)
test   %r11,%r11
js     1211 <quadct+0x111>
xorps  %xmm1,%xmm1
cvtsi2ss %r11,%xmm1
jmp    122a <quadct+0x12a>
mov    %r11,%rax
shr    %rax
and    $0x1,%r11d
or     %rax,%r11
xorps  %xmm1,%xmm1
cvtsi2ss %r11,%xmm1
addss  %xmm1,%xmm1
mov    0x18(%rsp),%rax
mulss  %xmm0,%xmm1
movss  %xmm1,(%r9)
test   %r10,%r10
js     1247 <quadct+0x147>
xorps  %xmm1,%xmm1
cvtsi2ss %r10,%xmm1
jmp    1260 <quadct+0x160>
mov    %r10,%rcx
shr    %rcx
and    $0x1,%r10d
or     %rcx,%r10
xorps  %xmm1,%xmm1
cvtsi2ss %r10,%xmm1
addss  %xmm1,%xmm1
mulss  %xmm1,%xmm0
movss  %xmm0,(%rax)
pop    %rbx
pop    %r14
ret

LLVM IR (-O0):
define void @quadct(float %0, float %1, float* %2, float* %3, i64 %4, float* %5, float* %6, float* %7, float* %8) {
%10 = alloca float
%11 = alloca float
%12 = alloca float*
%13 = alloca float*
%14 = alloca i64
%15 = alloca float*
%16 = alloca float*
%17 = alloca float*
%18 = alloca float*
%19 = alloca i64
%20 = alloca i64
%21 = alloca i64
%22 = alloca i64
%23 = alloca i64
%24 = alloca float
store float %0, float* %10
store float %1, float* %11
store float* %2, float** %12
store float* %3, float** %13
store i64 %4, i64* %14
store float* %5, float** %15
store float* %6, float** %16
store float* %7, float** %17
store float* %8, float** %18
store i64 0, i64* %23
store i64 0, i64* %22
store i64 0, i64* %21
store i64 0, i64* %20
store i64 1, i64* %19
br label %25
25: ; preds = %67, %9
%26 = load i64, i64* %19
%27 = load i64, i64* %14
%28 = icmp ule i64 %26, %27
br i1 %28, label %29, label %70
29: ; preds = %25
%30 = load float*, float** %13
%31 = load i64, i64* %19
%32 = getelementptr inbounds float, float* %30, i64 %31
%33 = load float, float* %32
%34 = load float, float* %11
%35 = fcmp ogt float %33, %34
br i1 %35, label %36, label %51
36: ; preds = %29
%37 = load float*, float** %12
%38 = load i64, i64* %19
%39 = getelementptr inbounds float, float* %37, i64 %38
%40 = load float, float* %39
%41 = load float, float* %10
%42 = fcmp ogt float %40, %41
br i1 %42, label %43, label %46
43: ; preds = %36
%44 = load i64, i64* %20
%45 = add i64 %44, 1
store i64 %45, i64* %20
br label %49
46: ; preds = %36
%47 = load i64, i64* %21
%48 = add i64 %47, 1
store i64 %48, i64* %21
br label %49
49: ; preds = %46, %43
%50 = phi i64 [ %45, %43 ], [ %48, %46 ]
br label %66
51: ; preds = %29
%52 = load float*, float** %12
%53 = load i64, i64* %19
%54 = getelementptr inbounds float, float* %52, i64 %53
%55 = load float, float* %54
%56 = load float, float* %10
%57 = fcmp ogt float %55, %56
br i1 %57, label %58, label %61
58: ; preds = %51
%59 = load i64, i64* %23
%60 = add i64 %59, 1
store i64 %60, i64* %23
br label %64
61: ; preds = %51
%62 = load i64, i64* %22
%63 = add i64 %62, 1
store i64 %63, i64* %22
br label %64
64: ; preds = %61, %58
%65 = phi i64 [ %60, %58 ], [ %63, %61 ]
br label %66
66: ; preds = %64, %49
br label %67
67: ; preds = %66
%68 = load i64, i64* %19
%69 = add i64 %68, 1
store i64 %69, i64* %19
br label %25
70: ; preds = %25
%71 = load i64, i64* %14
%72 = uitofp i64 %71 to double
%73 = fdiv double 1.000000e+00, %72
%74 = fptrunc double %73 to float
store float %74, float* %24
%75 = load float, float* %24
%76 = load i64, i64* %20
%77 = uitofp i64 %76 to float
%78 = fmul float %75, %77
%79 = load float*, float** %15
store float %78, float* %79
%80 = load float, float* %24
%81 = load i64, i64* %21
%82 = uitofp i64 %81 to float
%83 = fmul float %80, %82
%84 = load float*, float** %16
store float %83, float* %84
%85 = load float, float* %24
%86 = load i64, i64* %22
%87 = uitofp i64 %86 to float
%88 = fmul float %85, %87
%89 = load float*, float** %17
store float %88, float* %89
%90 = load float, float* %24
%91 = load i64, i64* %23
%92 = uitofp i64 %91 to float
%93 = fmul float %90, %92
%94 = load float*, float** %18
store float %93, float* %94
ret void
}

LLVM IR (-O2):
define void @quadct(float %0, float %1, float* nocapture %2, float* nocapture %3, i64 %4, float* nocapture writeonly %5, float* nocapture writeonly %6, float* nocapture writeonly %7, float* nocapture writeonly %8) {
%10 = icmp eq i64 %4, 0
br i1 %10, label %40, label %11
11: ; preds = %9, %33
%12 = phi i64 [ %38, %33 ], [ 1, %9 ]
%13 = phi i64 [ %37, %33 ], [ 0, %9 ]
%14 = phi i64 [ %36, %33 ], [ 0, %9 ]
%15 = phi i64 [ %35, %33 ], [ 0, %9 ]
%16 = phi i64 [ %34, %33 ], [ 0, %9 ]
%17 = getelementptr inbounds float, float* %3, i64 %12
%18 = load float, float* %17
%19 = fcmp ogt float %18, %1
%20 = getelementptr inbounds float, float* %2, i64 %12
%21 = load float, float* %20
%22 = fcmp ogt float %21, %0
br i1 %19, label %23, label %28
23: ; preds = %11
br i1 %22, label %24, label %26
24: ; preds = %23
%25 = add i64 %16, 1
br label %33
26: ; preds = %23
%27 = add i64 %15, 1
br label %33
28: ; preds = %11
br i1 %22, label %29, label %31
29: ; preds = %28
%30 = add i64 %13, 1
br label %33
31: ; preds = %28
%32 = add i64 %14, 1
br label %33
33: ; preds = %26, %24, %31, %29
%34 = phi i64 [ %25, %24 ], [ %16, %26 ], [ %16, %29 ], [ %16, %31 ]
%35 = phi i64 [ %15, %24 ], [ %27, %26 ], [ %15, %29 ], [ %15, %31 ]
%36 = phi i64 [ %14, %24 ], [ %14, %26 ], [ %14, %29 ], [ %32, %31 ]
%37 = phi i64 [ %13, %24 ], [ %13, %26 ], [ %30, %29 ], [ %13, %31 ]
%38 = add nuw i64 %12, 1
%39 = icmp eq i64 %12, %4
br i1 %39, label %40, label %11
40: ; preds = %33, %9
%41 = phi i64 [ 0, %9 ], [ %34, %33 ]
%42 = phi i64 [ 0, %9 ], [ %35, %33 ]
%43 = phi i64 [ 0, %9 ], [ %36, %33 ]
%44 = phi i64 [ 0, %9 ], [ %37, %33 ]
%45 = uitofp i64 %4 to double
%46 = fdiv double 1.000000e+00, %45
%47 = fptrunc double %46 to float
%48 = uitofp i64 %41 to float
%49 = fmul float %47, %48
store float %49, float* %5
%50 = uitofp i64 %42 to float
%51 = fmul float %47, %50
store float %51, float* %6
%52 = uitofp i64 %43 to float
%53 = fmul float %47, %52
store float %53, float* %7
%54 = uitofp i64 %44 to float
%55 = fmul float %47, %54
store float %55, float* %8
ret void
}