int numgtz(double *a, int n) {
  int num, k;
  num = 0;
  for (k = 0; k < n; ++k) {
    if (a[k] > .000001)
      ++num;
  }
  return num;
}


x86 assembly (clang, -O0):
<numgtz>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
movl   $0x0,-0x10(%rbp)
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1160 <numgtz+0x60>
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xec6(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    114d <numgtz+0x4d>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1152 <numgtz+0x52>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1119 <numgtz+0x19>
mov    -0x10(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<numgtz>:
test   %esi,%esi
jle    1114 <numgtz+0x14>
mov    %esi,%ecx
cmp    $0x4,%esi
jae    1117 <numgtz+0x17>
xor    %edx,%edx
xor    %eax,%eax
jmp    1202 <numgtz+0x102>
xor    %eax,%eax
ret
mov    %ecx,%edx
and    $0xfffffffc,%edx
lea    -0x4(%rdx),%rax
mov    %rax,%r8
shr    $0x2,%r8
add    $0x1,%r8
test   %rax,%rax
je     122b <numgtz+0x12b>
mov    %r8,%rax
and    $0xfffffffffffffffe,%rax
pxor   %xmm0,%xmm0
xor    %esi,%esi
movapd 0xeb7(%rip),%xmm2
pxor   %xmm1,%xmm1
nopl   (%rax)
movupd (%rdi,%rsi,8),%xmm3
movupd 0x10(%rdi,%rsi,8),%xmm4
movupd 0x20(%rdi,%rsi,8),%xmm5
movupd 0x30(%rdi,%rsi,8),%xmm6
movapd %xmm2,%xmm7
cmpltpd %xmm3,%xmm7
shufps $0xe8,%xmm7,%xmm7
psubd  %xmm7,%xmm0
movapd %xmm2,%xmm3
cmpltpd %xmm4,%xmm3
shufps $0xe8,%xmm3,%xmm3
psubd  %xmm3,%xmm1
movapd %xmm2,%xmm3
cmpltpd %xmm5,%xmm3
shufps $0xe8,%xmm3,%xmm3
psubd  %xmm3,%xmm0
movapd %xmm2,%xmm3
cmpltpd %xmm6,%xmm3
shufps $0xe8,%xmm3,%xmm3
psubd  %xmm3,%xmm1
add    $0x8,%rsi
add    $0xfffffffffffffffe,%rax
jne    1150 <numgtz+0x50>
test   $0x1,%r8b
je     11ec <numgtz+0xec>
movupd (%rdi,%rsi,8),%xmm2
movupd 0x10(%rdi,%rsi,8),%xmm3
movapd 0xe32(%rip),%xmm4
movapd %xmm4,%xmm5
cmpltpd %xmm2,%xmm5
shufps $0xe8,%xmm5,%xmm5
psubd  %xmm5,%xmm0
cmpltpd %xmm3,%xmm4
shufps $0xe8,%xmm4,%xmm4
psubd  %xmm4,%xmm1
paddd  %xmm1,%xmm0
pshufd $0x55,%xmm0,%xmm1
paddd  %xmm0,%xmm1
movd   %xmm1,%eax
cmp    %rcx,%rdx
je     122a <numgtz+0x12a>
movsd  0xe06(%rip),%xmm0
nopw   0x0(%rax,%rax,1)
movsd  (%rdi,%rdx,8),%xmm1
xor    %esi,%esi
ucomisd %xmm0,%xmm1
seta   %sil
add    %esi,%eax
add    $0x1,%rdx
cmp    %rdx,%rcx
jne    1210 <numgtz+0x110>
ret
pxor   %xmm0,%xmm0
xor    %esi,%esi
pxor   %xmm1,%xmm1
test   $0x1,%r8b
jne    11bb <numgtz+0xbb>
jmp    11ec <numgtz+0xec>

LLVM IR (-O0):
define i32 @numgtz(double* %0, i32 %1) {
%3 = alloca double*
%4 = alloca i32
%5 = alloca i32
%6 = alloca i32
store double* %0, double** %3
store i32 %1, i32* %4
store i32 0, i32* %5
store i32 0, i32* %6
br label %7
7: ; preds = %22, %2
%8 = load i32, i32* %6
%9 = load i32, i32* %4
%10 = icmp slt i32 %8, %9
br i1 %10, label %11, label %25
11: ; preds = %7
%12 = load double*, double** %3
%13 = load i32, i32* %6
%14 = sext i32 %13 to i64
%15 = getelementptr inbounds double, double* %12, i64 %14
%16 = load double, double* %15
%17 = fcmp ogt double %16, 0x3EB0C6F7A0B5ED8D
br i1 %17, label %18, label %21
18: ; preds = %11
%19 = load i32, i32* %5
%20 = add nsw i32 %19, 1
store i32 %20, i32* %5
br label %21
21: ; preds = %18, %11
br label %22
22: ; preds = %21
%23 = load i32, i32* %6
%24 = add nsw i32 %23, 1
store i32 %24, i32* %6
br label %7
25: ; preds = %7
%26 = load i32, i32* %5
ret i32 %26
}

LLVM IR (-O2):
define i32 @numgtz(double* nocapture %0, i32 %1) {
%3 = icmp sgt i32 %1, 0
br i1 %3, label %4, label %88
4: ; preds = %2
%5 = zext i32 %1 to i64
%6 = icmp ult i32 %1, 4
br i1 %6, label %75, label %7
7: ; preds = %4
%8 = and i64 %5, 4294967292
%9 = add nsw i64 %8, -4
%10 = lshr exact i64 %9, 2
%11 = add nuw nsw i64 %10, 1
%12 = and i64 %11, 1
%13 = icmp eq i64 %9, 0
br i1 %13, label %49, label %14
14: ; preds = %7
%15 = and i64 %11, 9223372036854775806
br label %16
16: ; preds = %16, %14
%17 = phi i64 [ 0, %14 ], [ %46, %16 ]
%18 = phi <2 x i32> [ zeroinitializer, %14 ], [ %44, %16 ]
%19 = phi <2 x i32> [ zeroinitializer, %14 ], [ %45, %16 ]
%20 = phi i64 [ 0, %14 ], [ %47, %16 ]
%21 = getelementptr inbounds double, double* %0, i64 %17
%22 = bitcast double* %21 to <2 x double>*
%23 = load <2 x double>, <2 x double>* %22
%24 = getelementptr inbounds double, double* %21, i64 2
%25 = bitcast double* %24 to <2 x double>*
%26 = load <2 x double>, <2 x double>* %25
%27 = fcmp ogt <2 x double> %23, <double 0x3EB0C6F7A0B5ED8D, double 0x3EB0C6F7A0B5ED8D>
%28 = fcmp ogt <2 x double> %26, <double 0x3EB0C6F7A0B5ED8D, double 0x3EB0C6F7A0B5ED8D>
%29 = zext <2 x i1> %27 to <2 x i32>
%30 = zext <2 x i1> %28 to <2 x i32>
%31 = add <2 x i32> %18, %29
%32 = add <2 x i32> %19, %30
%33 = or i64 %17, 4
%34 = getelementptr inbounds double, double* %0, i64 %33
%35 = bitcast double* %34 to <2 x double>*
%36 = load <2 x double>, <2 x double>* %35
%37 = getelementptr inbounds double, double* %34, i64 2
%38 = bitcast double* %37 to <2 x double>*
%39 = load <2 x double>, <2 x double>* %38
%40 = fcmp ogt <2 x double> %36, <double 0x3EB0C6F7A0B5ED8D, double 0x3EB0C6F7A0B5ED8D>
%41 = fcmp ogt <2 x double> %39, <double 0x3EB0C6F7A0B5ED8D, double 0x3EB0C6F7A0B5ED8D>
%42 = zext <2 x i1> %40 to <2 x i32>
%43 = zext <2 x i1> %41 to <2 x i32>
%44 = add <2 x i32> %31, %42
%45 = add <2 x i32> %32, %43
%46 = add nuw i64 %17, 8
%47 = add i64 %20, 2
%48 = icmp eq i64 %47, %15
br i1 %48, label %49, label %16
49: ; preds = %16, %7
%50 = phi <2 x i32> [ undef, %7 ], [ %44, %16 ]
%51 = phi <2 x i32> [ undef, %7 ], [ %45, %16 ]
%52 = phi i64 [ 0, %7 ], [ %46, %16 ]
%53 = phi <2 x i32> [ zeroinitializer, %7 ], [ %44, %16 ]
%54 = phi <2 x i32> [ zeroinitializer, %7 ], [ %45, %16 ]
%55 = icmp eq i64 %12, 0
br i1 %55, label %69, label %56
56: ; preds = %49
%57 = getelementptr inbounds double, double* %0, i64 %52
%58 = bitcast double* %57 to <2 x double>*
%59 = load <2 x double>, <2 x double>* %58
%60 = getelementptr inbounds double, double* %57, i64 2
%61 = bitcast double* %60 to <2 x double>*
%62 = load <2 x double>, <2 x double>* %61
%63 = fcmp ogt <2 x double> %59, <double 0x3EB0C6F7A0B5ED8D, double 0x3EB0C6F7A0B5ED8D>
%64 = fcmp ogt <2 x double> %62, <double 0x3EB0C6F7A0B5ED8D, double 0x3EB0C6F7A0B5ED8D>
%65 = zext <2 x i1> %63 to <2 x i32>
%66 = zext <2 x i1> %64 to <2 x i32>
%67 = add <2 x i32> %53, %65
%68 = add <2 x i32> %54, %66
br label %69
69: ; preds = %49, %56
%70 = phi <2 x i32> [ %50, %49 ], [ %67, %56 ]
%71 = phi <2 x i32> [ %51, %49 ], [ %68, %56 ]
%72 = add <2 x i32> %71, %70
%73 = call i32 @llvm.vector.reduce.add.v2i32(<2 x i32> %72)
%74 = icmp eq i64 %8, %5
br i1 %74, label %88, label %75
75: ; preds = %4, %69
%76 = phi i64 [ 0, %4 ], [ %8, %69 ]
%77 = phi i32 [ 0, %4 ], [ %73, %69 ]
br label %78
78: ; preds = %75, %78
%79 = phi i64 [ %86, %78 ], [ %76, %75 ]
%80 = phi i32 [ %85, %78 ], [ %77, %75 ]
%81 = getelementptr inbounds double, double* %0, i64 %79
%82 = load double, double* %81
%83 = fcmp ogt double %82, 0x3EB0C6F7A0B5ED8D
%84 = zext i1 %83 to i32
%85 = add nuw nsw i32 %80, %84
%86 = add nuw nsw i64 %79, 1
%87 = icmp eq i64 %86, %5
br i1 %87, label %88, label %78
88: ; preds = %78, %69, %2
%89 = phi i32 [ 0, %2 ], [ %73, %69 ], [ %85, %78 ]
ret i32 %89
}