void cleandata(double *p) {
  int i, j = 0;
  for (i = 0; i < 5; i++) {
    if (*(p + i) >= 0.001) {
      *(p + j) = *(p + i);
      j++;
    }
  }
  for (; j < 5; j++)
    *(p + j) = 0;
}


x86 assembly (clang, -O0):
<cleandata>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
movl   $0x0,-0x10(%rbp)
movl   $0x0,-0xc(%rbp)
cmpl   $0x5,-0xc(%rbp)
jge    1175 <cleandata+0x75>
mov    -0x8(%rbp),%rax
movslq -0xc(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xecb(%rip),%xmm1
ucomisd %xmm1,%xmm0
jb     1162 <cleandata+0x62>
mov    -0x8(%rbp),%rax
movslq -0xc(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1167 <cleandata+0x67>
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xc(%rbp)
jmp    1116 <cleandata+0x16>
jmp    117a <cleandata+0x7a>
cmpl   $0x5,-0x10(%rbp)
jge    11a2 <cleandata+0xa2>
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    117a <cleandata+0x7a>
pop    %rbp
ret

x86 assembly (clang, -O2):
<cleandata>:
push   %rax
movsd  (%rdi),%xmm1
movsd  0xee3(%rip),%xmm2
xor    %eax,%eax
ucomisd %xmm2,%xmm1
movsd  0x8(%rdi),%xmm0
setae  %cl
ucomisd %xmm2,%xmm0
jae    1135 <cleandata+0x25>
mov    %cl,%al
jmp    114a <cleandata+0x3a>
xor    %eax,%eax
ucomisd 0xec1(%rip),%xmm1
setae  %al
movsd  %xmm0,(%rdi,%rax,8)
add    $0x1,%eax
movsd  0x10(%rdi),%xmm0
ucomisd 0xea9(%rip),%xmm0
jb     1163 <cleandata+0x53>
mov    %eax,%ecx
movsd  %xmm0,(%rdi,%rcx,8)
add    $0x1,%eax
movsd  0x18(%rdi),%xmm0
ucomisd 0xe90(%rip),%xmm0
jb     117c <cleandata+0x6c>
mov    %eax,%ecx
movsd  %xmm0,(%rdi,%rcx,8)
add    $0x1,%eax
movsd  0x20(%rdi),%xmm0
ucomisd 0xe77(%rip),%xmm0
jb     119a <cleandata+0x8a>
mov    %eax,%ecx
movsd  %xmm0,(%rdi,%rcx,8)
cmp    $0x4,%eax
jae    11b6 <cleandata+0xa6>
add    $0x1,%eax
mov    %eax,%ecx
lea    (%rdi,%rcx,8),%rdi
mov    $0x4,%ecx
sub    %eax,%ecx
lea    0x8(,%rcx,8),%rdx
xor    %esi,%esi
call   1030 <memset@plt>
pop    %rax
ret

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

LLVM IR (-O2):
define void @cleandata(double* nocapture %0) {
%2 = load double, double* %0
%3 = fcmp oge double %2, 1.000000e-03
%4 = zext i1 %3 to i32
%5 = getelementptr inbounds double, double* %0, i64 1
%6 = load double, double* %5
%7 = fcmp ult double %6, 1.000000e-03
br i1 %7, label %21, label %17
8: ; preds = %39, %44
%9 = phi i32 [ %47, %44 ], [ %40, %39 ]
%10 = zext i32 %9 to i64
%11 = getelementptr double, double* %0, i64 %10
%12 = bitcast double* %11 to i8*
%13 = sub nuw nsw i32 4, %9
%14 = zext i32 %13 to i64
%15 = shl nuw nsw i64 %14, 3
%16 = add nuw nsw i64 %15, 8
br label %49
17: ; preds = %1
%18 = zext i1 %3 to i64
%19 = getelementptr inbounds double, double* %0, i64 %18
store double %6, double* %19
%20 = select i1 %3, i32 2, i32 1
br label %21
21: ; preds = %17, %1
%22 = phi i32 [ %20, %17 ], [ %4, %1 ]
%23 = getelementptr inbounds double, double* %0, i64 2
%24 = load double, double* %23
%25 = fcmp ult double %24, 1.000000e-03
br i1 %25, label %30, label %26
26: ; preds = %21
%27 = zext i32 %22 to i64
%28 = getelementptr inbounds double, double* %0, i64 %27
store double %24, double* %28
%29 = add nuw nsw i32 %22, 1
br label %30
30: ; preds = %26, %21
%31 = phi i32 [ %29, %26 ], [ %22, %21 ]
%32 = getelementptr inbounds double, double* %0, i64 3
%33 = load double, double* %32
%34 = fcmp ult double %33, 1.000000e-03
br i1 %34, label %39, label %35
35: ; preds = %30
%36 = zext i32 %31 to i64
%37 = getelementptr inbounds double, double* %0, i64 %36
store double %33, double* %37
%38 = add nuw nsw i32 %31, 1
br label %39
39: ; preds = %35, %30
%40 = phi i32 [ %38, %35 ], [ %31, %30 ]
%41 = getelementptr inbounds double, double* %0, i64 4
%42 = load double, double* %41
%43 = fcmp ult double %42, 1.000000e-03
br i1 %43, label %8, label %44
44: ; preds = %39
%45 = zext i32 %40 to i64
%46 = getelementptr inbounds double, double* %0, i64 %45
store double %42, double* %46
%47 = add nuw nsw i32 %40, 1
%48 = icmp ult i32 %40, 4
br i1 %48, label %8, label %49
49: ; preds = %8, %44
ret void
}