double dMeanp(double *pta, double l, double h, int n) {
  int i, cnt;
  double a, sum, mean;
  sum = 0.0;
  cnt = 0;
  for (i = 0; i < n; i++) {
    a = *(pta + i);
    if ((a > l) & (a < h)) {
      sum += a;
      cnt++;
    }
  }
  mean = sum / cnt;
  return mean;
}


x86 assembly (clang, -O0):
<dMeanp>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
movsd  %xmm0,-0x10(%rbp)
movsd  %xmm1,-0x18(%rbp)
mov    %esi,-0x1c(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x0,-0x24(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    11a8 <dMeanp+0xa8>
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
ucomisd -0x10(%rbp),%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
movsd  -0x30(%rbp),%xmm1
movsd  -0x18(%rbp),%xmm0
ucomisd %xmm1,%xmm0
seta   %cl
and    $0x1,%cl
movzbl %cl,%ecx
and    %ecx,%eax
cmp    $0x0,%eax
je     1195 <dMeanp+0x95>
movsd  -0x30(%rbp),%xmm0
addsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    119a <dMeanp+0x9a>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    112b <dMeanp+0x2b>
movsd  -0x38(%rbp),%xmm0
cvtsi2sdl -0x24(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<dMeanp>:
test   %esi,%esi
jle    1154 <dMeanp+0x54>
movapd %xmm0,%xmm2
mov    %esi,%r8d
cmp    $0x1,%esi
jne    115d <dMeanp+0x5d>
xor    %edx,%edx
xorpd  %xmm0,%xmm0
xor    %esi,%esi
test   $0x1,%r8b
je     1148 <dMeanp+0x48>
movsd  (%rdi,%rdx,8),%xmm3
ucomisd %xmm2,%xmm3
seta   %al
ucomisd %xmm3,%xmm1
seta   %cl
and    %al,%cl
movzbl %cl,%ecx
test   %cl,%cl
jne    1142 <dMeanp+0x42>
movsd  0xec6(%rip),%xmm3
addsd  %xmm3,%xmm0
add    %ecx,%esi
xorps  %xmm1,%xmm1
cvtsi2sd %esi,%xmm1
divsd  %xmm1,%xmm0
ret
movsd  0xea4(%rip),%xmm0
ret
mov    %r8d,%r9d
and    $0xfffffffe,%r9d
xor    %edx,%edx
xorpd  %xmm0,%xmm0
movsd  0xe96(%rip),%xmm3
xor    %esi,%esi
jmp    1196 <dMeanp+0x96>
cs nopw 0x0(%rax,%rax,1)
addsd  %xmm4,%xmm0
addsd  %xmm5,%xmm0
movzbl %cl,%eax
add    %eax,%esi
add    $0x2,%rdx
cmp    %rdx,%r9
je     1118 <dMeanp+0x18>
movsd  (%rdi,%rdx,8),%xmm4
ucomisd %xmm2,%xmm4
seta   %al
ucomisd %xmm4,%xmm1
seta   %cl
and    %al,%cl
jne    11b1 <dMeanp+0xb1>
movapd %xmm3,%xmm4
movzbl %cl,%eax
add    %eax,%esi
movsd  0x8(%rdi,%rdx,8),%xmm5
ucomisd %xmm2,%xmm5
seta   %al
ucomisd %xmm5,%xmm1
seta   %cl
and    %al,%cl
jne    1180 <dMeanp+0x80>
movapd %xmm3,%xmm5
jmp    1180 <dMeanp+0x80>

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

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