#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);
double c_canberra_location(long *x, long *y, long n, long k) {
  double d = 0.0;
  double xx, yy;
  long i;
  for (i = 0; i < n; i++) {
    xx = ((x[i] + 1) < k + 1) ? (x[i] + 1) : k + 1;
    yy = ((y[i] + 1) < k + 1) ? (y[i] + 1) : k + 1;
    d += fabs(xx - yy) / (xx + yy);
  }
  return d;
}


x86 assembly (clang, -O0):
<c_canberra_location>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
movq   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%rax
cmp    -0x18(%rbp),%rax
jge    1218 <c_canberra_location+0x118>
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
add    $0x1,%rax
mov    -0x20(%rbp),%rcx
add    $0x1,%rcx
cmp    %rcx,%rax
jge    116c <c_canberra_location+0x6c>
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
add    $0x1,%rax
mov    %rax,-0x48(%rbp)
jmp    1178 <c_canberra_location+0x78>
mov    -0x20(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x48(%rbp)
mov    -0x48(%rbp),%rax
cvtsi2sd %rax,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x40(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
add    $0x1,%rax
mov    -0x20(%rbp),%rcx
add    $0x1,%rcx
cmp    %rcx,%rax
jge    11c0 <c_canberra_location+0xc0>
mov    -0x10(%rbp),%rax
mov    -0x40(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
add    $0x1,%rax
mov    %rax,-0x50(%rbp)
jmp    11cc <c_canberra_location+0xcc>
mov    -0x20(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x50(%rbp)
mov    -0x50(%rbp),%rax
cvtsi2sd %rax,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x30(%rbp),%xmm0
subsd  -0x38(%rbp),%xmm0
movaps 0xe15(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  -0x30(%rbp),%xmm1
addsd  -0x38(%rbp),%xmm1
divsd  %xmm1,%xmm0
addsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    1124 <c_canberra_location+0x24>
movsd  -0x28(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<c_canberra_location>:
test   %rdx,%rdx
jle    1170 <c_canberra_location+0x70>
lea    0x1(%rcx),%r8
xorpd  %xmm0,%xmm0
xor    %eax,%eax
movapd 0xee9(%rip),%xmm1
nopw   0x0(%rax,%rax,1)
mov    (%rdi,%rax,8),%r9
lea    0x1(%r9),%r10
cmp    %rcx,%r9
cmovge %r8,%r10
xorps  %xmm2,%xmm2
cvtsi2sd %r10,%xmm2
mov    (%rsi,%rax,8),%r9
lea    0x1(%r9),%r10
cmp    %rcx,%r9
cmovge %r8,%r10
xorps  %xmm3,%xmm3
cvtsi2sd %r10,%xmm3
movapd %xmm2,%xmm4
subsd  %xmm3,%xmm4
andpd  %xmm1,%xmm4
addsd  %xmm3,%xmm2
divsd  %xmm2,%xmm4
addsd  %xmm4,%xmm0
add    $0x1,%rax
cmp    %rax,%rdx
jne    1120 <c_canberra_location+0x20>
ret
xorps  %xmm0,%xmm0
ret

LLVM IR (-O0):
define double @c_canberra_location(i64* %0, i64* %1, i64 %2, i64 %3) {
%5 = alloca i64*
%6 = alloca i64*
%7 = alloca i64
%8 = alloca i64
%9 = alloca double
%10 = alloca double
%11 = alloca double
%12 = alloca i64
store i64* %0, i64** %5
store i64* %1, i64** %6
store i64 %2, i64* %7
store i64 %3, i64* %8
store double 0.000000e+00, double* %9
store i64 0, i64* %12
br label %13
13: ; preds = %68, %4
%14 = load i64, i64* %12
%15 = load i64, i64* %7
%16 = icmp slt i64 %14, %15
br i1 %16, label %17, label %71
17: ; preds = %13
%18 = load i64*, i64** %5
%19 = load i64, i64* %12
%20 = getelementptr inbounds i64, i64* %18, i64 %19
%21 = load i64, i64* %20
%22 = add nsw i64 %21, 1
%23 = load i64, i64* %8
%24 = add nsw i64 %23, 1
%25 = icmp slt i64 %22, %24
br i1 %25, label %26, label %32
26: ; preds = %17
%27 = load i64*, i64** %5
%28 = load i64, i64* %12
%29 = getelementptr inbounds i64, i64* %27, i64 %28
%30 = load i64, i64* %29
%31 = add nsw i64 %30, 1
br label %35
32: ; preds = %17
%33 = load i64, i64* %8
%34 = add nsw i64 %33, 1
br label %35
35: ; preds = %32, %26
%36 = phi i64 [ %31, %26 ], [ %34, %32 ]
%37 = sitofp i64 %36 to double
store double %37, double* %10
%38 = load i64*, i64** %6
%39 = load i64, i64* %12
%40 = getelementptr inbounds i64, i64* %38, i64 %39
%41 = load i64, i64* %40
%42 = add nsw i64 %41, 1
%43 = load i64, i64* %8
%44 = add nsw i64 %43, 1
%45 = icmp slt i64 %42, %44
br i1 %45, label %46, label %52
46: ; preds = %35
%47 = load i64*, i64** %6
%48 = load i64, i64* %12
%49 = getelementptr inbounds i64, i64* %47, i64 %48
%50 = load i64, i64* %49
%51 = add nsw i64 %50, 1
br label %55
52: ; preds = %35
%53 = load i64, i64* %8
%54 = add nsw i64 %53, 1
br label %55
55: ; preds = %52, %46
%56 = phi i64 [ %51, %46 ], [ %54, %52 ]
%57 = sitofp i64 %56 to double
store double %57, double* %11
%58 = load double, double* %10
%59 = load double, double* %11
%60 = fsub double %58, %59
%61 = call double @llvm.fabs.f64(double %60)
%62 = load double, double* %10
%63 = load double, double* %11
%64 = fadd double %62, %63
%65 = fdiv double %61, %64
%66 = load double, double* %9
%67 = fadd double %66, %65
store double %67, double* %9
br label %68
68: ; preds = %55
%69 = load i64, i64* %12
%70 = add nsw i64 %69, 1
store i64 %70, i64* %12
br label %13
71: ; preds = %13
%72 = load double, double* %9
ret double %72
}

LLVM IR (-O2):
define double @c_canberra_location(i64* nocapture %0, i64* nocapture %1, i64 %2, i64 %3) {
%5 = icmp sgt i64 %2, 0
br i1 %5, label %6, label %30
6: ; preds = %4
%7 = add nsw i64 %3, 1
br label %8
8: ; preds = %6, %8
%9 = phi i64 [ 0, %6 ], [ %28, %8 ]
%10 = phi double [ 0.000000e+00, %6 ], [ %27, %8 ]
%11 = getelementptr inbounds i64, i64* %0, i64 %9
%12 = load i64, i64* %11
%13 = icmp slt i64 %12, %3
%14 = add nsw i64 %12, 1
%15 = select i1 %13, i64 %14, i64 %7
%16 = sitofp i64 %15 to double
%17 = getelementptr inbounds i64, i64* %1, i64 %9
%18 = load i64, i64* %17
%19 = icmp slt i64 %18, %3
%20 = add nsw i64 %18, 1
%21 = select i1 %19, i64 %20, i64 %7
%22 = sitofp i64 %21 to double
%23 = fsub double %16, %22
%24 = tail call double @llvm.fabs.f64(double %23)
%25 = fadd double %16, %22
%26 = fdiv double %24, %25
%27 = fadd double %10, %26
%28 = add nuw nsw i64 %9, 1
%29 = icmp eq i64 %28, %2
br i1 %29, label %30, label %8
30: ; preds = %8, %4
%31 = phi double [ 0.000000e+00, %4 ], [ %27, %8 ]
ret double %31
}