#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
double fabs(double);
long neigh(double *y, long n, long d, long m, long kmax, long npmax, double eps,
           long *nlist) {
  long nn, i, nfound;
  for (nfound = 0, nn = n - 1; nn >= (m - 1) * d && nfound < kmax; nn--) {
    for (i = 0; i < m; i++)
      if (fabs((double)(y[n - i * d] - y[nn - i * d])) >= eps)
        break;
    if (i == m)
      nlist[nfound++] = nn;
    if (npmax && nn < n - npmax)
      break;
  }
  return nfound;
}


x86 assembly (clang, -O0):
<neigh>:
push   %rbp
mov    %rsp,%rbp
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movsd  %xmm0,-0x38(%rbp)
movq   $0x0,-0x50(%rbp)
mov    -0x10(%rbp),%rax
sub    $0x1,%rax
mov    %rax,-0x40(%rbp)
mov    -0x40(%rbp),%rcx
mov    -0x20(%rbp),%rdx
sub    $0x1,%rdx
imul   -0x18(%rbp),%rdx
xor    %eax,%eax
cmp    %rdx,%rcx
mov    %al,-0x51(%rbp)
jl     1166 <neigh+0x66>
mov    -0x50(%rbp),%rax
cmp    -0x28(%rbp),%rax
setl   %al
mov    %al,-0x51(%rbp)
mov    -0x51(%rbp),%al
test   $0x1,%al
jne    1176 <neigh+0x76>
jmp    1253 <neigh+0x153>
movq   $0x0,-0x48(%rbp)
mov    -0x48(%rbp),%rax
cmp    -0x20(%rbp),%rax
jge    11ef <neigh+0xef>
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%rcx
mov    -0x48(%rbp),%rdx
imul   -0x18(%rbp),%rdx
sub    %rdx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%rcx
mov    -0x48(%rbp),%rdx
imul   -0x18(%rbp),%rdx
sub    %rdx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movaps 0xe3b(%rip),%xmm1
pand   %xmm1,%xmm0
ucomisd -0x38(%rbp),%xmm0
jb     11d9 <neigh+0xd9>
jmp    11ef <neigh+0xef>
jmp    11de <neigh+0xde>
mov    -0x48(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x48(%rbp)
jmp    117e <neigh+0x7e>
mov    -0x48(%rbp),%rax
cmp    -0x20(%rbp),%rax
jne    1218 <neigh+0x118>
mov    -0x40(%rbp),%rdx
mov    0x10(%rbp),%rax
mov    -0x50(%rbp),%rcx
mov    %rcx,%rsi
add    $0x1,%rsi
mov    %rsi,-0x50(%rbp)
mov    %rdx,(%rax,%rcx,8)
cmpq   $0x0,-0x30(%rbp)
je     123d <neigh+0x13d>
mov    -0x40(%rbp),%rax
mov    -0x10(%rbp),%rcx
sub    -0x30(%rbp),%rcx
cmp    %rcx,%rax
jge    123d <neigh+0x13d>
jmp    1253 <neigh+0x153>
jmp    1242 <neigh+0x142>
mov    -0x40(%rbp),%rax
add    $0xffffffffffffffff,%rax
mov    %rax,-0x40(%rbp)
jmp    1139 <neigh+0x39>
mov    -0x50(%rbp),%rax
pop    %rbp
ret

x86 assembly (clang, -O2):
<neigh>:
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    0x30(%rsp),%r10
lea    -0x1(%rcx),%r11
imul   %rdx,%r11
mov    %rsi,%r15
sub    %r9,%r15
lea    (%rdi,%rsi,8),%r14
shl    $0x3,%rdx
neg    %rdx
mov    $0xfffffffffffffff8,%r13
xor    %eax,%eax
movapd 0xec8(%rip),%xmm1
nopl   0x0(%rax,%rax,1)
cmp    %r11,%rsi
jle    11ad <neigh+0xad>
cmp    %r8,%rax
jge    11ad <neigh+0xad>
mov    %rsi,%r12
add    $0xffffffffffffffff,%rsi
test   %rcx,%rcx
jle    1190 <neigh+0x90>
mov    %r14,%rdi
xor    %ebx,%ebx
nopl   0x0(%rax,%rax,1)
movsd  (%rdi),%xmm2
subsd  (%rdi,%r13,1),%xmm2
andpd  %xmm1,%xmm2
ucomisd %xmm0,%xmm2
jae    1192 <neigh+0x92>
add    $0x1,%rbx
add    %rdx,%rdi
cmp    %rbx,%rcx
jne    1160 <neigh+0x60>
jmp    1197 <neigh+0x97>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
xor    %ebx,%ebx
cmp    %rcx,%rbx
jne    119f <neigh+0x9f>
mov    %rsi,(%r10,%rax,8)
add    $0x1,%rax
add    $0xfffffffffffffff8,%r13
test   %r9,%r9
je     1140 <neigh+0x40>
cmp    %r15,%r12
jg     1140 <neigh+0x40>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define i64 @neigh(double* %0, i64 %1, i64 %2, i64 %3, i64 %4, i64 %5, double %6, i64* %7) {
%9 = alloca double*
%10 = alloca i64
%11 = alloca i64
%12 = alloca i64
%13 = alloca i64
%14 = alloca i64
%15 = alloca double
%16 = alloca i64*
%17 = alloca i64
%18 = alloca i64
%19 = alloca i64
store double* %0, double** %9
store i64 %1, i64* %10
store i64 %2, i64* %11
store i64 %3, i64* %12
store i64 %4, i64* %13
store i64 %5, i64* %14
store double %6, double* %15
store i64* %7, i64** %16
store i64 0, i64* %19
%20 = load i64, i64* %10
%21 = sub nsw i64 %20, 1
store i64 %21, i64* %17
br label %22
22: ; preds = %87, %8
%23 = load i64, i64* %17
%24 = load i64, i64* %12
%25 = sub nsw i64 %24, 1
%26 = load i64, i64* %11
%27 = mul nsw i64 %25, %26
%28 = icmp sge i64 %23, %27
br i1 %28, label %29, label %33
29: ; preds = %22
%30 = load i64, i64* %19
%31 = load i64, i64* %13
%32 = icmp slt i64 %30, %31
br label %33
33: ; preds = %29, %22
%34 = phi i1 [ false, %22 ], [ %32, %29 ]
br i1 %34, label %35, label %90
35: ; preds = %33
store i64 0, i64* %18
br label %36
36: ; preds = %63, %35
%37 = load i64, i64* %18
%38 = load i64, i64* %12
%39 = icmp slt i64 %37, %38
br i1 %39, label %40, label %66
40: ; preds = %36
%41 = load double*, double** %9
%42 = load i64, i64* %10
%43 = load i64, i64* %18
%44 = load i64, i64* %11
%45 = mul nsw i64 %43, %44
%46 = sub nsw i64 %42, %45
%47 = getelementptr inbounds double, double* %41, i64 %46
%48 = load double, double* %47
%49 = load double*, double** %9
%50 = load i64, i64* %17
%51 = load i64, i64* %18
%52 = load i64, i64* %11
%53 = mul nsw i64 %51, %52
%54 = sub nsw i64 %50, %53
%55 = getelementptr inbounds double, double* %49, i64 %54
%56 = load double, double* %55
%57 = fsub double %48, %56
%58 = call double @llvm.fabs.f64(double %57)
%59 = load double, double* %15
%60 = fcmp oge double %58, %59
br i1 %60, label %61, label %62
61: ; preds = %40
br label %66
62: ; preds = %40
br label %63
63: ; preds = %62
%64 = load i64, i64* %18
%65 = add nsw i64 %64, 1
store i64 %65, i64* %18
br label %36
66: ; preds = %61, %36
%67 = load i64, i64* %18
%68 = load i64, i64* %12
%69 = icmp eq i64 %67, %68
br i1 %69, label %70, label %76
70: ; preds = %66
%71 = load i64, i64* %17
%72 = load i64*, i64** %16
%73 = load i64, i64* %19
%74 = add nsw i64 %73, 1
store i64 %74, i64* %19
%75 = getelementptr inbounds i64, i64* %72, i64 %73
store i64 %71, i64* %75
br label %76
76: ; preds = %70, %66
%77 = load i64, i64* %14
%78 = icmp ne i64 %77, 0
br i1 %78, label %79, label %86
79: ; preds = %76
%80 = load i64, i64* %17
%81 = load i64, i64* %10
%82 = load i64, i64* %14
%83 = sub nsw i64 %81, %82
%84 = icmp slt i64 %80, %83
br i1 %84, label %85, label %86
85: ; preds = %79
br label %90
86: ; preds = %79, %76
br label %87
87: ; preds = %86
%88 = load i64, i64* %17
%89 = add nsw i64 %88, -1
store i64 %89, i64* %17
br label %22
90: ; preds = %85, %33
%91 = load i64, i64* %19
ret i64 %91
}

LLVM IR (-O2):
define i64 @neigh(double* nocapture %0, i64 %1, i64 %2, i64 %3, i64 %4, i64 %5, double %6, i64* nocapture writeonly %7) {
%9 = add nsw i64 %3, -1
%10 = mul nsw i64 %9, %2
%11 = icmp sgt i64 %3, 0
%12 = icmp eq i64 %5, 0
%13 = sub nsw i64 %1, %5
br label %14
14: ; preds = %43, %8
%15 = phi i64 [ %1, %8 ], [ %17, %43 ]
%16 = phi i64 [ 0, %8 ], [ %44, %43 ]
%17 = add nsw i64 %15, -1
%18 = icmp sgt i64 %15, %10
%19 = icmp slt i64 %16, %4
%20 = select i1 %18, i1 %19, i1 false
br i1 %20, label %21, label %47
21: ; preds = %14
br i1 %11, label %22, label %37
22: ; preds = %21, %34
%23 = phi i64 [ %35, %34 ], [ 0, %21 ]
%24 = mul nsw i64 %23, %2
%25 = sub nsw i64 %1, %24
%26 = getelementptr inbounds double, double* %0, i64 %25
%27 = load double, double* %26
%28 = sub nsw i64 %17, %24
%29 = getelementptr inbounds double, double* %0, i64 %28
%30 = load double, double* %29
%31 = fsub double %27, %30
%32 = tail call double @llvm.fabs.f64(double %31)
%33 = fcmp ult double %32, %6
br i1 %33, label %34, label %37
34: ; preds = %22
%35 = add nuw nsw i64 %23, 1
%36 = icmp eq i64 %35, %3
br i1 %36, label %40, label %22
37: ; preds = %22, %21
%38 = phi i64 [ 0, %21 ], [ %23, %22 ]
%39 = icmp eq i64 %38, %3
br i1 %39, label %40, label %43
40: ; preds = %34, %37
%41 = add nsw i64 %16, 1
%42 = getelementptr inbounds i64, i64* %7, i64 %16
store i64 %17, i64* %42
br label %43
43: ; preds = %40, %37
%44 = phi i64 [ %41, %40 ], [ %16, %37 ]
%45 = icmp sgt i64 %15, %13
%46 = select i1 %12, i1 true, i1 %45
br i1 %46, label %14, label %47
47: ; preds = %43, %14
%48 = phi i64 [ %16, %14 ], [ %44, %43 ]
ret i64 %48
}