void decimation(ndata, data) int *ndata;
float data[];
{
  int i;
  for (i = 0; i < *ndata; i += 2) {
    data[i / 2] = data[i];
  }
  *ndata = *ndata / 2 + *ndata % 2;
  return;
}


x86 assembly (clang, -O0):
<decimation>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
mov    -0x8(%rbp),%rcx
cmp    (%rcx),%eax
jge    115e <decimation+0x5e>
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    %rax,-0x20(%rbp)
mov    -0x14(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x20(%rbp),%rax
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x14(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x14(%rbp)
jmp    1113 <decimation+0x13>
mov    -0x8(%rbp),%rax
mov    (%rax),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x8(%rbp),%rax
mov    (%rax),%eax
mov    $0x2,%esi
cltd
idiv   %esi
add    %edx,%ecx
mov    -0x8(%rbp),%rax
mov    %ecx,(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<decimation>:
mov    (%rdi),%r8d
test   %r8d,%r8d
jle    1193 <decimation+0x93>
lea    -0x1(%r8),%rax
mov    %rax,%rdx
shr    %rdx
add    $0x1,%rdx
mov    %edx,%ecx
and    $0x3,%ecx
cmp    $0x6,%rax
jae    1129 <decimation+0x29>
xor    %eax,%eax
jmp    1168 <decimation+0x68>
and    $0xfffffffffffffffc,%rdx
xor    %eax,%eax
nop
movss  (%rsi,%rax,4),%xmm0
movss  %xmm0,(%rsi,%rax,2)
movss  0x8(%rsi,%rax,4),%xmm0
movss  %xmm0,0x4(%rsi,%rax,2)
movss  0x10(%rsi,%rax,4),%xmm0
movss  %xmm0,0x8(%rsi,%rax,2)
movss  0x18(%rsi,%rax,4),%xmm0
movss  %xmm0,0xc(%rsi,%rax,2)
add    $0x8,%rax
add    $0xfffffffffffffffc,%rdx
jne    1130 <decimation+0x30>
test   %rcx,%rcx
je     1193 <decimation+0x93>
mov    %rax,%rdx
and    $0xfffffffffffffffe,%rdx
lea    (%rsi,%rdx,2),%rdx
lea    (%rsi,%rax,4),%rax
xor    %esi,%esi
xchg   %ax,%ax
movss  (%rax,%rsi,8),%xmm0
movss  %xmm0,(%rdx,%rsi,4)
add    $0x1,%rsi
cmp    %rsi,%rcx
jne    1180 <decimation+0x80>
mov    %r8d,%eax
shr    $0x1f,%eax
add    %r8d,%eax
mov    %eax,%ecx
and    $0xfffffffe,%ecx
sub    %ecx,%r8d
sar    %eax
add    %r8d,%eax
mov    %eax,(%rdi)
ret

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

LLVM IR (-O2):
define void @decimation(i32* nocapture %0, float* nocapture %1) {
%3 = load i32, i32* %0
%4 = icmp sgt i32 %3, 0
br i1 %4, label %5, label %52
5: ; preds = %2
%6 = zext i32 %3 to i64
%7 = add nsw i64 %6, -1
%8 = lshr i64 %7, 1
%9 = add nuw i64 %8, 1
%10 = and i64 %9, 3
%11 = icmp ult i64 %7, 6
br i1 %11, label %39, label %12
12: ; preds = %5
%13 = and i64 %9, -4
br label %14
14: ; preds = %14, %12
%15 = phi i64 [ 0, %12 ], [ %36, %14 ]
%16 = phi i64 [ 0, %12 ], [ %37, %14 ]
%17 = getelementptr inbounds float, float* %1, i64 %15
%18 = load float, float* %17
%19 = lshr exact i64 %15, 1
%20 = getelementptr inbounds float, float* %1, i64 %19
store float %18, float* %20
%21 = or i64 %15, 2
%22 = getelementptr inbounds float, float* %1, i64 %21
%23 = load float, float* %22
%24 = lshr exact i64 %21, 1
%25 = getelementptr inbounds float, float* %1, i64 %24
store float %23, float* %25
%26 = or i64 %15, 4
%27 = getelementptr inbounds float, float* %1, i64 %26
%28 = load float, float* %27
%29 = lshr exact i64 %26, 1
%30 = getelementptr inbounds float, float* %1, i64 %29
store float %28, float* %30
%31 = or i64 %15, 6
%32 = getelementptr inbounds float, float* %1, i64 %31
%33 = load float, float* %32
%34 = lshr exact i64 %31, 1
%35 = getelementptr inbounds float, float* %1, i64 %34
store float %33, float* %35
%36 = add nuw nsw i64 %15, 8
%37 = add i64 %16, 4
%38 = icmp eq i64 %37, %13
br i1 %38, label %39, label %14
39: ; preds = %14, %5
%40 = phi i64 [ 0, %5 ], [ %36, %14 ]
%41 = icmp eq i64 %10, 0
br i1 %41, label %52, label %42
42: ; preds = %39, %42
%43 = phi i64 [ %49, %42 ], [ %40, %39 ]
%44 = phi i64 [ %50, %42 ], [ 0, %39 ]
%45 = getelementptr inbounds float, float* %1, i64 %43
%46 = load float, float* %45
%47 = lshr exact i64 %43, 1
%48 = getelementptr inbounds float, float* %1, i64 %47
store float %46, float* %48
%49 = add nuw nsw i64 %43, 2
%50 = add i64 %44, 1
%51 = icmp eq i64 %50, %10
br i1 %51, label %52, label %42
52: ; preds = %39, %42, %2
%53 = sdiv i32 %3, 2
%54 = srem i32 %3, 2
%55 = add nsw i32 %53, %54
store i32 %55, i32* %0
ret void
}