#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 */
int /*<<< orphan*/ dirty_sort(int, int *, double *, int *);
void weighted_KS(int m, int *g, double *r, int n_gene, double *w,
                 double *result, int *iwrk) {
  double pos_score = 0.0;
  double neg_score = 0.0;
  result[0] = 0.0;
  result[1] = 0.0;
  if (m == 0)
    return;
  dirty_sort(m, g, r, iwrk);
  double sum_weight = 0;
  for (int i = 0; i < m; i++)
    sum_weight += w[g[i]];
  double P = 0.0, Q, Qprev;
  double d, dprev;
  for (int i = 1; i <= m; i++) {
    if (g[i - 1] > 0) {
      Qprev = (r[g[i - 1] - 1] - (i - 1)) / (n_gene - m);
    } else {
      Qprev = 0.0;
    }
    dprev = P - Qprev;
    P += w[g[i - 1]] / sum_weight;
    Q = (r[g[i - 1]] - i) / (n_gene - m);
    d = P - Q;
    if (d > pos_score)
      pos_score = d;
    if (dprev < neg_score)
      neg_score = dprev;
  }
  result[0] = pos_score;
  result[1] = neg_score;
  return;
}


x86 assembly (clang, -O0):
<weighted_KS>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x30(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
mov    -0x30(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x8(%rax)
cmpl   $0x0,-0x4(%rbp)
jne    116b <weighted_KS+0x5b>
jmp    1329 <weighted_KS+0x219>
mov    -0x4(%rbp),%edi
mov    -0x10(%rbp),%rsi
mov    -0x18(%rbp),%rdx
mov    0x10(%rbp),%rcx
call   1030 <dirty_sort@plt>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11c7 <weighted_KS+0xb7>
mov    -0x28(%rbp),%rax
mov    -0x10(%rbp),%rcx
movslq -0x4c(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    118e <weighted_KS+0x7e>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
movl   $0x1,-0x7c(%rbp)
mov    -0x7c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     130e <weighted_KS+0x1fe>
mov    -0x10(%rbp),%rax
mov    -0x7c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
cmpl   $0x0,(%rax,%rcx,4)
jle    123e <weighted_KS+0x12e>
mov    -0x18(%rbp),%rax
mov    -0x10(%rbp),%rcx
mov    -0x7c(%rbp),%edx
sub    $0x1,%edx
movslq %edx,%rdx
mov    (%rcx,%rdx,4),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x7c(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
subsd  %xmm1,%xmm0
mov    -0x1c(%rbp),%eax
sub    -0x4(%rbp),%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
jmp    1246 <weighted_KS+0x136>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x58(%rbp),%xmm0
subsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x28(%rbp),%rax
mov    -0x10(%rbp),%rcx
mov    -0x7c(%rbp),%edx
sub    $0x1,%edx
movslq %edx,%rdx
movslq (%rcx,%rdx,4),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  -0x48(%rbp),%xmm0
addsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x10(%rbp),%rcx
mov    -0x7c(%rbp),%edx
sub    $0x1,%edx
movslq %edx,%rdx
movslq (%rcx,%rdx,4),%rcx
movsd  (%rax,%rcx,8),%xmm0
cvtsi2sdl -0x7c(%rbp),%xmm1
subsd  %xmm1,%xmm0
mov    -0x1c(%rbp),%eax
sub    -0x4(%rbp),%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x58(%rbp),%xmm0
subsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x70(%rbp),%xmm0
ucomisd -0x38(%rbp),%xmm0
jbe    12dd <weighted_KS+0x1cd>
movsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x78(%rbp),%xmm1
movsd  -0x40(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    12fb <weighted_KS+0x1eb>
movsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
jmp    1300 <weighted_KS+0x1f0>
mov    -0x7c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x7c(%rbp)
jmp    11d6 <weighted_KS+0xc6>
movsd  -0x38(%rbp),%xmm0
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x40(%rbp),%xmm0
mov    -0x30(%rbp),%rax
movsd  %xmm0,0x8(%rax)
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<weighted_KS>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
push   %rax
xorpd  %xmm0,%xmm0
movupd %xmm0,(%r9)
test   %edi,%edi
je     1285 <weighted_KS+0x175>
mov    %r9,%r14
mov    %r8,%rbx
mov    %ecx,%r13d
mov    %rdx,%r15
mov    %rsi,%r12
mov    %edi,%ebp
mov    0x40(%rsp),%rcx
call   1030 <dirty_sort@plt>
test   %ebp,%ebp
jle    1272 <weighted_KS+0x162>
mov    %ebp,%eax
lea    -0x1(%rax),%rdx
mov    %eax,%ecx
and    $0x3,%ecx
cmp    $0x3,%rdx
jae    1168 <weighted_KS+0x58>
xorpd  %xmm0,%xmm0
xor    %edx,%edx
jmp    11b0 <weighted_KS+0xa0>
mov    %eax,%esi
and    $0xfffffffc,%esi
xorpd  %xmm0,%xmm0
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movslq (%r12,%rdx,4),%rdi
addsd  (%rbx,%rdi,8),%xmm0
movslq 0x4(%r12,%rdx,4),%rdi
addsd  (%rbx,%rdi,8),%xmm0
movslq 0x8(%r12,%rdx,4),%rdi
addsd  (%rbx,%rdi,8),%xmm0
movslq 0xc(%r12,%rdx,4),%rdi
addsd  (%rbx,%rdi,8),%xmm0
add    $0x4,%rdx
cmp    %rdx,%rsi
jne    1180 <weighted_KS+0x70>
test   %rcx,%rcx
je     11d2 <weighted_KS+0xc2>
lea    (%r12,%rdx,4),%rdx
xor    %esi,%esi
nopl   0x0(%rax,%rax,1)
movslq (%rdx,%rsi,4),%rdi
addsd  (%rbx,%rdi,8),%xmm0
add    $0x1,%rsi
cmp    %rsi,%rcx
jne    11c0 <weighted_KS+0xb0>
test   %ebp,%ebp
jle    1272 <weighted_KS+0x162>
sub    %ebp,%r13d
cvtsi2sd %r13d,%xmm1
xorpd  %xmm2,%xmm2
xor    %ecx,%ecx
xorpd  %xmm5,%xmm5
xorpd  %xmm3,%xmm3
jmp    124b <weighted_KS+0x13b>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movapd %xmm3,%xmm4
subsd  %xmm6,%xmm4
minsd  %xmm5,%xmm4
movsd  (%rbx,%rdx,8),%xmm5
divsd  %xmm0,%xmm5
addsd  %xmm5,%xmm3
movsd  (%r15,%rdx,8),%xmm5
add    $0x1,%rcx
xorps  %xmm6,%xmm6
cvtsi2sd %ecx,%xmm6
subsd  %xmm6,%xmm5
divsd  %xmm1,%xmm5
movapd %xmm3,%xmm6
subsd  %xmm5,%xmm6
maxsd  %xmm2,%xmm6
movapd %xmm6,%xmm2
movapd %xmm4,%xmm5
cmp    %rcx,%rax
je     127a <weighted_KS+0x16a>
movslq (%r12,%rcx,4),%rdx
xorpd  %xmm6,%xmm6
test   %rdx,%rdx
jle    1200 <weighted_KS+0xf0>
lea    -0x1(%rdx),%esi
movsd  (%r15,%rsi,8),%xmm6
xorps  %xmm4,%xmm4
cvtsi2sd %ecx,%xmm4
subsd  %xmm4,%xmm6
divsd  %xmm1,%xmm6
jmp    1200 <weighted_KS+0xf0>
xorpd  %xmm4,%xmm4
xorpd  %xmm6,%xmm6
movsd  %xmm6,(%r14)
movsd  %xmm4,0x8(%r14)
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @weighted_KS(i32 %0, i32* %1, double* %2, i32 %3, double* %4, double* %5, i32* %6) {
%8 = alloca i32
%9 = alloca i32*
%10 = alloca double*
%11 = alloca i32
%12 = alloca double*
%13 = alloca double*
%14 = alloca i32*
%15 = alloca double
%16 = alloca double
%17 = alloca double
%18 = alloca i32
%19 = alloca double
%20 = alloca double
%21 = alloca double
%22 = alloca double
%23 = alloca double
%24 = alloca i32
store i32 %0, i32* %8
store i32* %1, i32** %9
store double* %2, double** %10
store i32 %3, i32* %11
store double* %4, double** %12
store double* %5, double** %13
store i32* %6, i32** %14
store double 0.000000e+00, double* %15
store double 0.000000e+00, double* %16
%25 = load double*, double** %13
%26 = getelementptr inbounds double, double* %25, i64 0
store double 0.000000e+00, double* %26
%27 = load double*, double** %13
%28 = getelementptr inbounds double, double* %27, i64 1
store double 0.000000e+00, double* %28
%29 = load i32, i32* %8
%30 = icmp eq i32 %29, 0
br i1 %30, label %31, label %32
31: ; preds = %7
br label %153
32: ; preds = %7
%33 = load i32, i32* %8
%34 = load i32*, i32** %9
%35 = load double*, double** %10
%36 = load i32*, i32** %14
%37 = call i32 @dirty_sort(i32 %33, i32* %34, double* %35, i32* %36)
store double 0.000000e+00, double* %17
store i32 0, i32* %18
br label %38
38: ; preds = %54, %32
%39 = load i32, i32* %18
%40 = load i32, i32* %8
%41 = icmp slt i32 %39, %40
br i1 %41, label %42, label %57
42: ; preds = %38
%43 = load double*, double** %12
%44 = load i32*, i32** %9
%45 = load i32, i32* %18
%46 = sext i32 %45 to i64
%47 = getelementptr inbounds i32, i32* %44, i64 %46
%48 = load i32, i32* %47
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds double, double* %43, i64 %49
%51 = load double, double* %50
%52 = load double, double* %17
%53 = fadd double %52, %51
store double %53, double* %17
br label %54
54: ; preds = %42
%55 = load i32, i32* %18
%56 = add nsw i32 %55, 1
store i32 %56, i32* %18
br label %38
57: ; preds = %38
store double 0.000000e+00, double* %19
store i32 1, i32* %24
br label %58
58: ; preds = %143, %57
%59 = load i32, i32* %24
%60 = load i32, i32* %8
%61 = icmp sle i32 %59, %60
br i1 %61, label %62, label %146
62: ; preds = %58
%63 = load i32*, i32** %9
%64 = load i32, i32* %24
%65 = sub nsw i32 %64, 1
%66 = sext i32 %65 to i64
%67 = getelementptr inbounds i32, i32* %63, i64 %66
%68 = load i32, i32* %67
%69 = icmp sgt i32 %68, 0
br i1 %69, label %70, label %91
70: ; preds = %62
%71 = load double*, double** %10
%72 = load i32*, i32** %9
%73 = load i32, i32* %24
%74 = sub nsw i32 %73, 1
%75 = sext i32 %74 to i64
%76 = getelementptr inbounds i32, i32* %72, i64 %75
%77 = load i32, i32* %76
%78 = sub nsw i32 %77, 1
%79 = sext i32 %78 to i64
%80 = getelementptr inbounds double, double* %71, i64 %79
%81 = load double, double* %80
%82 = load i32, i32* %24
%83 = sub nsw i32 %82, 1
%84 = sitofp i32 %83 to double
%85 = fsub double %81, %84
%86 = load i32, i32* %11
%87 = load i32, i32* %8
%88 = sub nsw i32 %86, %87
%89 = sitofp i32 %88 to double
%90 = fdiv double %85, %89
store double %90, double* %21
br label %92
91: ; preds = %62
store double 0.000000e+00, double* %21
br label %92
92: ; preds = %91, %70
%93 = load double, double* %19
%94 = load double, double* %21
%95 = fsub double %93, %94
store double %95, double* %23
%96 = load double*, double** %12
%97 = load i32*, i32** %9
%98 = load i32, i32* %24
%99 = sub nsw i32 %98, 1
%100 = sext i32 %99 to i64
%101 = getelementptr inbounds i32, i32* %97, i64 %100
%102 = load i32, i32* %101
%103 = sext i32 %102 to i64
%104 = getelementptr inbounds double, double* %96, i64 %103
%105 = load double, double* %104
%106 = load double, double* %17
%107 = fdiv double %105, %106
%108 = load double, double* %19
%109 = fadd double %108, %107
store double %109, double* %19
%110 = load double*, double** %10
%111 = load i32*, i32** %9
%112 = load i32, i32* %24
%113 = sub nsw i32 %112, 1
%114 = sext i32 %113 to i64
%115 = getelementptr inbounds i32, i32* %111, i64 %114
%116 = load i32, i32* %115
%117 = sext i32 %116 to i64
%118 = getelementptr inbounds double, double* %110, i64 %117
%119 = load double, double* %118
%120 = load i32, i32* %24
%121 = sitofp i32 %120 to double
%122 = fsub double %119, %121
%123 = load i32, i32* %11
%124 = load i32, i32* %8
%125 = sub nsw i32 %123, %124
%126 = sitofp i32 %125 to double
%127 = fdiv double %122, %126
store double %127, double* %20
%128 = load double, double* %19
%129 = load double, double* %20
%130 = fsub double %128, %129
store double %130, double* %22
%131 = load double, double* %22
%132 = load double, double* %15
%133 = fcmp ogt double %131, %132
br i1 %133, label %134, label %136
134: ; preds = %92
%135 = load double, double* %22
store double %135, double* %15
br label %136
136: ; preds = %134, %92
%137 = load double, double* %23
%138 = load double, double* %16
%139 = fcmp olt double %137, %138
br i1 %139, label %140, label %142
140: ; preds = %136
%141 = load double, double* %23
store double %141, double* %16
br label %142
142: ; preds = %140, %136
br label %143
143: ; preds = %142
%144 = load i32, i32* %24
%145 = add nsw i32 %144, 1
store i32 %145, i32* %24
br label %58
146: ; preds = %58
%147 = load double, double* %15
%148 = load double*, double** %13
%149 = getelementptr inbounds double, double* %148, i64 0
store double %147, double* %149
%150 = load double, double* %16
%151 = load double*, double** %13
%152 = getelementptr inbounds double, double* %151, i64 1
store double %150, double* %152
br label %153
153: ; preds = %146, %31
ret void
}

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