int subSample(a, b, npix, step)
float *a;
float *b;
int npix, step;
{
  register int ip, i, got;
  if (step < 1) {
    step = 1;
  }
  got = 0;
  ip = 0;
  for (i = 0; i < npix; i++) {
    if (isfinite(a[ip])) {
      b[got++] = a[ip];
    }
    ip += step;
  }
  return got;
}


x86 assembly (clang, -O0):
<subSample>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
cmpl   $0x1,-0x18(%rbp)
jge    1137 <subSample+0x27>
movl   $0x1,-0x18(%rbp)
movl   $0x0,-0x24(%rbp)
movl   $0x0,-0x1c(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11b4 <subSample+0xa4>
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    $0x1,%al
call   1030 <isfinite@plt>
cmp    $0x0,%eax
je     119d <subSample+0x8d>
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x24(%rbp)
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x18(%rbp),%eax
add    -0x1c(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    114c <subSample+0x3c>
mov    -0x24(%rbp),%eax
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<subSample>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
test   %edx,%edx
jle    116e <subSample+0x5e>
mov    %edx,%ebp
mov    %rsi,%r14
mov    %rdi,%rbx
test   %ecx,%ecx
mov    $0x1,%r12d
cmovg  %ecx,%r12d
shl    $0x2,%r12
xor    %r15d,%r15d
jmp    1148 <subSample+0x38>
nopl   0x0(%rax)
add    %r12,%rbx
add    $0xffffffff,%ebp
je     1171 <subSample+0x61>
movss  (%rbx),%xmm0
cvtss2sd %xmm0,%xmm0
mov    $0x1,%al
call   1030 <isfinite@plt>
test   %eax,%eax
je     1140 <subSample+0x30>
movss  (%rbx),%xmm0
movslq %r15d,%rax
add    $0x1,%r15d
movss  %xmm0,(%r14,%rax,4)
jmp    1140 <subSample+0x30>
xor    %r15d,%r15d
mov    %r15d,%eax
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define i32 @subSample(float* nocapture %0, float* nocapture writeonly %1, i32 %2, i32 %3) {
%5 = icmp sgt i32 %2, 0
br i1 %5, label %6, label %29
6: ; preds = %4
%7 = icmp sgt i32 %3, 1
%8 = select i1 %7, i32 %3, i32 1
%9 = zext i32 %8 to i64
br label %10
10: ; preds = %6, %24
%11 = phi i64 [ 0, %6 ], [ %26, %24 ]
%12 = phi i32 [ 0, %6 ], [ %25, %24 ]
%13 = phi i32 [ 0, %6 ], [ %27, %24 ]
%14 = getelementptr inbounds float, float* %0, i64 %11
%15 = load float, float* %14
%16 = fpext float %15 to double
%17 = tail call i32 (double, ...) bitcast (i32 (...)* @isfinite to i32 (double, ...)*)(double %16)
%18 = icmp eq i32 %17, 0
br i1 %18, label %24, label %19
19: ; preds = %10
%20 = load float, float* %14
%21 = add nsw i32 %12, 1
%22 = sext i32 %12 to i64
%23 = getelementptr inbounds float, float* %1, i64 %22
store float %20, float* %23
br label %24
24: ; preds = %19, %10
%25 = phi i32 [ %21, %19 ], [ %12, %10 ]
%26 = add nuw nsw i64 %11, %9
%27 = add nuw nsw i32 %13, 1
%28 = icmp eq i32 %27, %2
br i1 %28, label %29, label %10
29: ; preds = %24, %4
%30 = phi i32 [ 0, %4 ], [ %25, %24 ]
ret i32 %30
}