#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 km_ks_dist(int na, const double a[], int nb, const double b[]) {
  int ia = 0, ib = 0;
  double fa = 0, fb = 0, sup = 0, na1 = 1. / na, nb1 = 1. / nb;
  while (ia < na || ib < nb) {
    if (ia == na)
      fb += nb1, ++ib;
    else if (ib == nb)
      fa += na1, ++ia;
    else if (a[ia] < b[ib])
      fa += na1, ++ia;
    else if (a[ia] > b[ib])
      fb += nb1, ++ib;
    else
      fa += na1, fb += nb1, ++ia, ++ib;
    if (sup < fabs(fa - fb))
      sup = fabs(fa - fb);
  }
  return sup;
}


x86 assembly (clang, -O0):
<km_ks_dist>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
movl   $0x0,-0x28(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
cvtsi2sdl -0x4(%rbp),%xmm1
movsd  0xebb(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
cvtsi2sdl -0x14(%rbp),%xmm1
movsd  0xea5(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x24(%rbp),%ecx
mov    $0x1,%al
cmp    -0x4(%rbp),%ecx
mov    %al,-0x51(%rbp)
jl     1181 <km_ks_dist+0x81>
mov    -0x28(%rbp),%eax
cmp    -0x14(%rbp),%eax
setl   %al
mov    %al,-0x51(%rbp)
mov    -0x51(%rbp),%al
test   $0x1,%al
jne    1191 <km_ks_dist+0x91>
jmp    12e3 <km_ks_dist+0x1e3>
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jne    11ba <km_ks_dist+0xba>
movsd  -0x50(%rbp),%xmm0
addsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    12a0 <km_ks_dist+0x1a0>
mov    -0x28(%rbp),%eax
cmp    -0x14(%rbp),%eax
jne    11e3 <km_ks_dist+0xe3>
movsd  -0x48(%rbp),%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    129b <km_ks_dist+0x19b>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd %xmm1,%xmm0
jbe    1224 <km_ks_dist+0x124>
movsd  -0x48(%rbp),%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1296 <km_ks_dist+0x196>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x28(%rbp),%rcx
ucomisd (%rax,%rcx,8),%xmm0
jbe    1261 <km_ks_dist+0x161>
movsd  -0x50(%rbp),%xmm0
addsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1291 <km_ks_dist+0x191>
movsd  -0x48(%rbp),%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x50(%rbp),%xmm0
addsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1296 <km_ks_dist+0x196>
jmp    129b <km_ks_dist+0x19b>
jmp    12a0 <km_ks_dist+0x1a0>
movsd  -0x40(%rbp),%xmm1
movsd  -0x30(%rbp),%xmm0
subsd  -0x38(%rbp),%xmm0
movaps 0xd5a(%rip),%xmm2
pand   %xmm2,%xmm0
ucomisd %xmm1,%xmm0
jbe    12de <km_ks_dist+0x1de>
movsd  -0x30(%rbp),%xmm0
subsd  -0x38(%rbp),%xmm0
movaps 0xd3b(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
jmp    1164 <km_ks_dist+0x64>
movsd  -0x40(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<km_ks_dist>:
test   %edi,%edi
jg     1110 <km_ks_dist+0x10>
xorpd  %xmm0,%xmm0
test   %edx,%edx
jle    120a <km_ks_dist+0x10a>
movd   %edx,%xmm0
movd   %edi,%xmm1
punpckldq %xmm0,%xmm1
cvtdq2pd %xmm1,%xmm0
movapd 0xed8(%rip),%xmm1
divpd  %xmm0,%xmm1
xorpd  %xmm2,%xmm2
xorpd  %xmm0,%xmm0
xor    %r9d,%r9d
movapd 0xed1(%rip),%xmm3
xor    %r8d,%r8d
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movapd %xmm0,%xmm4
cmp    %edi,%r9d
jne    1170 <km_ks_dist+0x70>
movapd %xmm1,%xmm0
addpd  %xmm2,%xmm0
add    $0x1,%r8d
shufpd $0x2,%xmm0,%xmm2
mov    %edi,%r9d
jmp    11e0 <km_ks_dist+0xe0>
nop
cmp    %edx,%r8d
jne    1190 <km_ks_dist+0x90>
add    $0x1,%r9d
addsd  %xmm1,%xmm2
mov    %edx,%r8d
jmp    11e0 <km_ks_dist+0xe0>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movslq %r9d,%rax
movsd  (%rsi,%rax,8),%xmm0
movslq %r8d,%rax
movsd  (%rcx,%rax,8),%xmm5
ucomisd %xmm0,%xmm5
jbe    11b0 <km_ks_dist+0xb0>
add    $0x1,%r9d
addsd  %xmm1,%xmm2
jmp    11e0 <km_ks_dist+0xe0>
ucomisd %xmm5,%xmm0
jbe    11c9 <km_ks_dist+0xc9>
movapd %xmm1,%xmm0
addpd  %xmm2,%xmm0
add    $0x1,%r8d
shufpd $0x2,%xmm0,%xmm2
jmp    11e0 <km_ks_dist+0xe0>
addpd  %xmm1,%xmm2
add    $0x1,%r9d
add    $0x1,%r8d
cs nopw 0x0(%rax,%rax,1)
nop
movapd %xmm2,%xmm5
unpckhpd %xmm2,%xmm5
movapd %xmm2,%xmm0
subsd  %xmm5,%xmm0
andpd  %xmm3,%xmm0
maxsd  %xmm4,%xmm0
cmp    %edi,%r9d
jl     1150 <km_ks_dist+0x50>
cmp    %edx,%r8d
jl     1150 <km_ks_dist+0x50>
ret

LLVM IR (-O0):
define double @km_ks_dist(i32 %0, double* %1, i32 %2, double* %3) {
%5 = alloca i32
%6 = alloca double*
%7 = alloca i32
%8 = alloca double*
%9 = alloca i32
%10 = alloca i32
%11 = alloca double
%12 = alloca double
%13 = alloca double
%14 = alloca double
%15 = alloca double
store i32 %0, i32* %5
store double* %1, double** %6
store i32 %2, i32* %7
store double* %3, double** %8
store i32 0, i32* %9
store i32 0, i32* %10
store double 0.000000e+00, double* %11
store double 0.000000e+00, double* %12
store double 0.000000e+00, double* %13
%16 = load i32, i32* %5
%17 = sitofp i32 %16 to double
%18 = fdiv double 1.000000e+00, %17
store double %18, double* %14
%19 = load i32, i32* %7
%20 = sitofp i32 %19 to double
%21 = fdiv double 1.000000e+00, %20
store double %21, double* %15
br label %22
22: ; preds = %114, %4
%23 = load i32, i32* %9
%24 = load i32, i32* %5
%25 = icmp slt i32 %23, %24
br i1 %25, label %30, label %26
26: ; preds = %22
%27 = load i32, i32* %10
%28 = load i32, i32* %7
%29 = icmp slt i32 %27, %28
br label %30
30: ; preds = %26, %22
%31 = phi i1 [ true, %22 ], [ %29, %26 ]
br i1 %31, label %32, label %115
32: ; preds = %30
%33 = load i32, i32* %9
%34 = load i32, i32* %5
%35 = icmp eq i32 %33, %34
br i1 %35, label %36, label %42
36: ; preds = %32
%37 = load double, double* %15
%38 = load double, double* %12
%39 = fadd double %38, %37
store double %39, double* %12
%40 = load i32, i32* %10
%41 = add nsw i32 %40, 1
store i32 %41, i32* %10
br label %102
42: ; preds = %32
%43 = load i32, i32* %10
%44 = load i32, i32* %7
%45 = icmp eq i32 %43, %44
br i1 %45, label %46, label %52
46: ; preds = %42
%47 = load double, double* %14
%48 = load double, double* %11
%49 = fadd double %48, %47
store double %49, double* %11
%50 = load i32, i32* %9
%51 = add nsw i32 %50, 1
store i32 %51, i32* %9
br label %101
52: ; preds = %42
%53 = load double*, double** %6
%54 = load i32, i32* %9
%55 = sext i32 %54 to i64
%56 = getelementptr inbounds double, double* %53, i64 %55
%57 = load double, double* %56
%58 = load double*, double** %8
%59 = load i32, i32* %10
%60 = sext i32 %59 to i64
%61 = getelementptr inbounds double, double* %58, i64 %60
%62 = load double, double* %61
%63 = fcmp olt double %57, %62
br i1 %63, label %64, label %70
64: ; preds = %52
%65 = load double, double* %14
%66 = load double, double* %11
%67 = fadd double %66, %65
store double %67, double* %11
%68 = load i32, i32* %9
%69 = add nsw i32 %68, 1
store i32 %69, i32* %9
br label %100
70: ; preds = %52
%71 = load double*, double** %6
%72 = load i32, i32* %9
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds double, double* %71, i64 %73
%75 = load double, double* %74
%76 = load double*, double** %8
%77 = load i32, i32* %10
%78 = sext i32 %77 to i64
%79 = getelementptr inbounds double, double* %76, i64 %78
%80 = load double, double* %79
%81 = fcmp ogt double %75, %80
br i1 %81, label %82, label %88
82: ; preds = %70
%83 = load double, double* %15
%84 = load double, double* %12
%85 = fadd double %84, %83
store double %85, double* %12
%86 = load i32, i32* %10
%87 = add nsw i32 %86, 1
store i32 %87, i32* %10
br label %99
88: ; preds = %70
%89 = load double, double* %14
%90 = load double, double* %11
%91 = fadd double %90, %89
store double %91, double* %11
%92 = load double, double* %15
%93 = load double, double* %12
%94 = fadd double %93, %92
store double %94, double* %12
%95 = load i32, i32* %9
%96 = add nsw i32 %95, 1
store i32 %96, i32* %9
%97 = load i32, i32* %10
%98 = add nsw i32 %97, 1
store i32 %98, i32* %10
br label %99
99: ; preds = %88, %82
br label %100
100: ; preds = %99, %64
br label %101
101: ; preds = %100, %46
br label %102
102: ; preds = %101, %36
%103 = load double, double* %13
%104 = load double, double* %11
%105 = load double, double* %12
%106 = fsub double %104, %105
%107 = call double @llvm.fabs.f64(double %106)
%108 = fcmp olt double %103, %107
br i1 %108, label %109, label %114
109: ; preds = %102
%110 = load double, double* %11
%111 = load double, double* %12
%112 = fsub double %110, %111
%113 = call double @llvm.fabs.f64(double %112)
store double %113, double* %13
br label %114
114: ; preds = %109, %102
br label %22
115: ; preds = %30
%116 = load double, double* %13
ret double %116
}

LLVM IR (-O2):
define double @km_ks_dist(i32 %0, double* nocapture %1, i32 %2, double* nocapture %3) {
%5 = insertelement <2 x i32> poison, i32 %0, i64 0
%6 = insertelement <2 x i32> %5, i32 %2, i64 1
%7 = sitofp <2 x i32> %6 to <2 x double>
%8 = fdiv <2 x double> <double 1.000000e+00, double 1.000000e+00>, %7
%9 = icmp sgt i32 %0, 0
%10 = icmp sgt i32 %2, 0
%11 = or i1 %9, %10
br i1 %11, label %12, label %63
12: ; preds = %4, %50
%13 = phi i32 [ %52, %50 ], [ 0, %4 ]
%14 = phi i32 [ %51, %50 ], [ 0, %4 ]
%15 = phi double [ %59, %50 ], [ 0.000000e+00, %4 ]
%16 = phi <2 x double> [ %53, %50 ], [ zeroinitializer, %4 ]
%17 = icmp eq i32 %13, %0
br i1 %17, label %18, label %22
18: ; preds = %12
%19 = fadd <2 x double> %8, %16
%20 = add nsw i32 %14, 1
%21 = shufflevector <2 x double> %16, <2 x double> %19, <2 x i32> <i32 0, i32 3>
br label %50
22: ; preds = %12
%23 = icmp eq i32 %14, %2
br i1 %23, label %24, label %28
24: ; preds = %22
%25 = fadd <2 x double> %8, %16
%26 = add nsw i32 %13, 1
%27 = shufflevector <2 x double> %25, <2 x double> %16, <2 x i32> <i32 0, i32 3>
br label %50
28: ; preds = %22
%29 = sext i32 %13 to i64
%30 = getelementptr inbounds double, double* %1, i64 %29
%31 = load double, double* %30
%32 = sext i32 %14 to i64
%33 = getelementptr inbounds double, double* %3, i64 %32
%34 = load double, double* %33
%35 = fcmp olt double %31, %34
br i1 %35, label %36, label %40
36: ; preds = %28
%37 = fadd <2 x double> %8, %16
%38 = add nsw i32 %13, 1
%39 = shufflevector <2 x double> %37, <2 x double> %16, <2 x i32> <i32 0, i32 3>
br label %50
40: ; preds = %28
%41 = fcmp ogt double %31, %34
br i1 %41, label %42, label %46
42: ; preds = %40
%43 = fadd <2 x double> %8, %16
%44 = add nsw i32 %14, 1
%45 = shufflevector <2 x double> %16, <2 x double> %43, <2 x i32> <i32 0, i32 3>
br label %50
46: ; preds = %40
%47 = fadd <2 x double> %8, %16
%48 = add nsw i32 %13, 1
%49 = add nsw i32 %14, 1
br label %50
50: ; preds = %24, %42, %46, %36, %18
%51 = phi i32 [ %20, %18 ], [ %2, %24 ], [ %14, %36 ], [ %44, %42 ], [ %49, %46 ]
%52 = phi i32 [ %0, %18 ], [ %26, %24 ], [ %38, %36 ], [ %13, %42 ], [ %48, %46 ]
%53 = phi <2 x double> [ %21, %18 ], [ %27, %24 ], [ %39, %36 ], [ %45, %42 ], [ %47, %46 ]
%54 = shufflevector <2 x double> %53, <2 x double> poison, <2 x i32> <i32 1, i32 undef>
%55 = fsub <2 x double> %53, %54
%56 = extractelement <2 x double> %55, i64 0
%57 = tail call double @llvm.fabs.f64(double %56)
%58 = fcmp olt double %15, %57
%59 = select i1 %58, double %57, double %15
%60 = icmp slt i32 %52, %0
%61 = icmp slt i32 %51, %2
%62 = select i1 %60, i1 true, i1 %61
br i1 %62, label %12, label %63
63: ; preds = %50, %4
%64 = phi double [ 0.000000e+00, %4 ], [ %59, %50 ]
ret double %64
}