int num_src_points(const int *restrict mask, const int src_add[4],
                   double src_lats[4]) {
  int icount = 0;
  for (int n = 0; n < 4; ++n) {
    if (mask[src_add[n]])
      icount++;
    else
      src_lats[n] = 0.;
  }
  return (icount);
}


x86 assembly (clang, -O0):
<num_src_points>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
movl   $0x0,-0x20(%rbp)
cmpl   $0x4,-0x20(%rbp)
jge    1173 <num_src_points+0x73>
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%rcx
movslq -0x20(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
cmpl   $0x0,(%rax,%rcx,4)
je     1150 <num_src_points+0x50>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1160 <num_src_points+0x60>
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1165 <num_src_points+0x65>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    111e <num_src_points+0x1e>
mov    -0x1c(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<num_src_points>:
movslq (%rsi),%rcx
mov    $0x1,%eax
cmpl   $0x0,(%rdi,%rcx,4)
jne    1117 <num_src_points+0x17>
movq   $0x0,(%rdx)
xor    %eax,%eax
movslq 0x4(%rsi),%rcx
cmpl   $0x0,(%rdi,%rcx,4)
je     1126 <num_src_points+0x26>
add    $0x1,%eax
jmp    112e <num_src_points+0x2e>
movq   $0x0,0x8(%rdx)
movslq 0x8(%rsi),%rcx
cmpl   $0x0,(%rdi,%rcx,4)
je     113d <num_src_points+0x3d>
add    $0x1,%eax
jmp    1145 <num_src_points+0x45>
movq   $0x0,0x10(%rdx)
movslq 0xc(%rsi),%rcx
cmpl   $0x0,(%rdi,%rcx,4)
je     1153 <num_src_points+0x53>
add    $0x1,%eax
ret
movq   $0x0,0x18(%rdx)
ret

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

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