int inlidxso(const double *err, const double *sgn, int len, double th,
             int *inl_buff, int **inls) {
  int i, po = 0, ne = 0;
  for (i = 0; i < len; i++)
    if (err[i] <= th) {
      if (sgn[i] > 0) {
        inl_buff[po] = i;
        po++;
      } else {
        ne++;
        inl_buff[len - ne] = i;
      }
    }
  if (po >= ne) {
    *inls = inl_buff;
    return po;
  }
  *inls = inl_buff + len - ne;
  return ne;
}


x86 assembly (clang, -O0):
<inlidxso>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
movsd  %xmm0,-0x28(%rbp)
mov    %rcx,-0x30(%rbp)
mov    %r8,-0x38(%rbp)
movl   $0x0,-0x40(%rbp)
movl   $0x0,-0x44(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    11c3 <inlidxso+0xc3>
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x28(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jb     11b0 <inlidxso+0xb0>
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    118f <inlidxso+0x8f>
mov    -0x3c(%rbp),%edx
mov    -0x30(%rbp),%rax
movslq -0x40(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    11ab <inlidxso+0xab>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
mov    -0x3c(%rbp),%edx
mov    -0x30(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    -0x44(%rbp),%ecx
movslq %ecx,%rcx
mov    %edx,(%rax,%rcx,4)
jmp    11b0 <inlidxso+0xb0>
jmp    11b5 <inlidxso+0xb5>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1131 <inlidxso+0x31>
mov    -0x40(%rbp),%eax
cmp    -0x44(%rbp),%eax
jl     11e5 <inlidxso+0xe5>
mov    -0x30(%rbp),%rcx
mov    -0x38(%rbp),%rax
mov    %rcx,(%rax)
mov    -0x40(%rbp),%eax
mov    %eax,-0x4(%rbp)
jmp    1211 <inlidxso+0x111>
mov    -0x30(%rbp),%rcx
movslq -0x1c(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rcx
movslq -0x44(%rbp),%rdx
xor    %eax,%eax
sub    %rdx,%rax
shl    $0x2,%rax
add    %rax,%rcx
mov    -0x38(%rbp),%rax
mov    %rcx,(%rax)
mov    -0x44(%rbp),%eax
mov    %eax,-0x4(%rbp)
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<inlidxso>:
push   %rbx
test   %edx,%edx
jle    1158 <inlidxso+0x58>
mov    %edx,%r10d
xor    %r11d,%r11d
xorpd  %xmm1,%xmm1
xor    %r9d,%r9d
xor    %eax,%eax
jmp    1139 <inlidxso+0x39>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%r9d
mov    %edx,%ebx
sub    %r9d,%ebx
movslq %ebx,%rbx
mov    %r11d,(%rcx,%rbx,4)
add    $0x1,%r11
cmp    %r11,%r10
je     115d <inlidxso+0x5d>
ucomisd (%rdi,%r11,8),%xmm0
jb     1130 <inlidxso+0x30>
movsd  (%rsi,%r11,8),%xmm2
ucomisd %xmm1,%xmm2
jbe    1120 <inlidxso+0x20>
cltq
mov    %r11d,(%rcx,%rax,4)
add    $0x1,%eax
jmp    1130 <inlidxso+0x30>
xor    %eax,%eax
xor    %r9d,%r9d
movslq %edx,%rdx
lea    (%rcx,%rdx,4),%rdx
movslq %r9d,%rsi
lea    0x0(,%rsi,4),%rdi
sub    %rdi,%rdx
cmp    %esi,%eax
cmovge %rcx,%rdx
cmovl  %esi,%eax
mov    %rdx,(%r8)
pop    %rbx
ret

LLVM IR (-O0):
define i32 @inlidxso(double* %0, double* %1, i32 %2, double %3, i32* %4, i32** %5) {
%7 = alloca i32
%8 = alloca double*
%9 = alloca double*
%10 = alloca i32
%11 = alloca double
%12 = alloca i32*
%13 = alloca i32**
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
store double* %0, double** %8
store double* %1, double** %9
store i32 %2, i32* %10
store double %3, double* %11
store i32* %4, i32** %12
store i32** %5, i32*** %13
store i32 0, i32* %15
store i32 0, i32* %16
store i32 0, i32* %14
br label %17
17: ; preds = %56, %6
%18 = load i32, i32* %14
%19 = load i32, i32* %10
%20 = icmp slt i32 %18, %19
br i1 %20, label %21, label %59
21: ; preds = %17
%22 = load double*, double** %8
%23 = load i32, i32* %14
%24 = sext i32 %23 to i64
%25 = getelementptr inbounds double, double* %22, i64 %24
%26 = load double, double* %25
%27 = load double, double* %11
%28 = fcmp ole double %26, %27
br i1 %28, label %29, label %55
29: ; preds = %21
%30 = load double*, double** %9
%31 = load i32, i32* %14
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds double, double* %30, i64 %32
%34 = load double, double* %33
%35 = fcmp ogt double %34, 0.000000e+00
br i1 %35, label %36, label %44
36: ; preds = %29
%37 = load i32, i32* %14
%38 = load i32*, i32** %12
%39 = load i32, i32* %15
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds i32, i32* %38, i64 %40
store i32 %37, i32* %41
%42 = load i32, i32* %15
%43 = add nsw i32 %42, 1
store i32 %43, i32* %15
br label %54
44: ; preds = %29
%45 = load i32, i32* %16
%46 = add nsw i32 %45, 1
store i32 %46, i32* %16
%47 = load i32, i32* %14
%48 = load i32*, i32** %12
%49 = load i32, i32* %10
%50 = load i32, i32* %16
%51 = sub nsw i32 %49, %50
%52 = sext i32 %51 to i64
%53 = getelementptr inbounds i32, i32* %48, i64 %52
store i32 %47, i32* %53
br label %54
54: ; preds = %44, %36
br label %55
55: ; preds = %54, %21
br label %56
56: ; preds = %55
%57 = load i32, i32* %14
%58 = add nsw i32 %57, 1
store i32 %58, i32* %14
br label %17
59: ; preds = %17
%60 = load i32, i32* %15
%61 = load i32, i32* %16
%62 = icmp sge i32 %60, %61
br i1 %62, label %63, label %67
63: ; preds = %59
%64 = load i32*, i32** %12
%65 = load i32**, i32*** %13
store i32* %64, i32** %65
%66 = load i32, i32* %15
store i32 %66, i32* %7
br label %78
67: ; preds = %59
%68 = load i32*, i32** %12
%69 = load i32, i32* %10
%70 = sext i32 %69 to i64
%71 = getelementptr inbounds i32, i32* %68, i64 %70
%72 = load i32, i32* %16
%73 = sext i32 %72 to i64
%74 = sub i64 0, %73
%75 = getelementptr inbounds i32, i32* %71, i64 %74
%76 = load i32**, i32*** %13
store i32* %75, i32** %76
%77 = load i32, i32* %16
store i32 %77, i32* %7
br label %78
78: ; preds = %67, %63
%79 = load i32, i32* %7
ret i32 %79
}

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