#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(float);
float sqrt(float);
void kendl1(float data1[], float data2[], unsigned long n, float *tau, float *z,
            float *prob) {
  float erfcc(float x);
  unsigned long n2 = 0, n1 = 0, k, j;
  long is = 0;
  float svar, aa, a2, a1;
  for (j = 1; j < n; j++) {
    for (k = (j + 1); k <= n; k++) {
      a1 = data1[j] - data1[k];
      a2 = data2[j] - data2[k];
      aa = a1 * a2;
      if (aa) {
        ++n1;
        ++n2;
        aa > 0.0 ? ++is : --is;
      } else {
        if (a1)
          ++n1;
        if (a2)
          ++n2;
      }
    }
  }
  *tau = is / (sqrt((double)n1) * sqrt((double)n2));
  svar = (4.0 * n + 10.0) / (9.0 * n * (n - 1.0));
  *z = (*tau) / sqrt(svar);
  *prob = erfcc(fabs(*z) / 1.4142136);
}


x86 assembly (clang, -O0):
<kendl1>:
push   %rbp
mov    %rsp,%rbp
sub    $0xa0,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movq   $0x0,-0x38(%rbp)
movq   $0x0,-0x40(%rbp)
movq   $0x0,-0x58(%rbp)
movq   $0x1,-0x50(%rbp)
mov    -0x50(%rbp),%rax
cmp    -0x18(%rbp),%rax
jae    12d5 <kendl1+0x1a5>
mov    -0x50(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x48(%rbp)
mov    -0x48(%rbp),%rax
cmp    -0x18(%rbp),%rax
ja     12bf <kendl1+0x18f>
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x68(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x50(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x48(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x64(%rbp)
movss  -0x68(%rbp),%xmm0
mulss  -0x64(%rbp),%xmm0
movss  %xmm0,-0x60(%rbp)
movss  -0x60(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jne    1204 <kendl1+0xd4>
jp     1204 <kendl1+0xd4>
jmp    1254 <kendl1+0x124>
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
movss  -0x60(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    1243 <kendl1+0x113>
mov    -0x58(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x58(%rbp)
jmp    124f <kendl1+0x11f>
mov    -0x58(%rbp),%rax
add    $0xffffffffffffffff,%rax
mov    %rax,-0x58(%rbp)
jmp    12a9 <kendl1+0x179>
movss  -0x68(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jne    1270 <kendl1+0x140>
jp     1270 <kendl1+0x140>
jmp    127c <kendl1+0x14c>
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
movss  -0x64(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jne    1298 <kendl1+0x168>
jp     1298 <kendl1+0x168>
jmp    12a4 <kendl1+0x174>
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    12a9 <kendl1+0x179>
jmp    12ae <kendl1+0x17e>
mov    -0x48(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x48(%rbp)
jmp    118d <kendl1+0x5d>
jmp    12c4 <kendl1+0x194>
mov    -0x50(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x50(%rbp)
jmp    1173 <kendl1+0x43>
cvtsi2ssq -0x58(%rbp),%xmm0
movss  %xmm0,-0x94(%rbp)
ff
movq   -0x40(%rbp),%xmm0
movaps 0xd41(%rip),%xmm1
movaps %xmm1,-0x90(%rbp)
punpckldq %xmm1,%xmm0
movapd 0xd3e(%rip),%xmm1
movaps %xmm1,-0x80(%rbp)
subpd  %xmm1,%xmm0
movaps %xmm0,%xmm1
unpckhpd %xmm0,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
call   1050 <sqrt@plt>
movaps -0x90(%rbp),%xmm2
movaps -0x80(%rbp),%xmm1
movss  %xmm0,-0x98(%rbp)
ff
movq   -0x38(%rbp),%xmm0
punpckldq %xmm2,%xmm0
subpd  %xmm1,%xmm0
movaps %xmm0,%xmm1
unpckhpd %xmm0,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
call   1050 <sqrt@plt>
movss  -0x98(%rbp),%xmm3
ff
movss  -0x94(%rbp),%xmm2
ff
movaps -0x90(%rbp),%xmm1
movaps %xmm0,%xmm4
movaps -0x80(%rbp),%xmm0
mulss  %xmm4,%xmm3
divss  %xmm3,%xmm2
mov    -0x20(%rbp),%rax
movss  %xmm2,(%rax)
movq   -0x18(%rbp),%xmm2
punpckldq %xmm1,%xmm2
subpd  %xmm0,%xmm2
movaps %xmm2,%xmm0
unpckhpd %xmm2,%xmm2
addsd  %xmm0,%xmm2
movsd  0xc70(%rip),%xmm1
movaps %xmm2,%xmm0
mulsd  %xmm1,%xmm0
movsd  0xc69(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  0xc65(%rip),%xmm3
movaps %xmm2,%xmm1
mulsd  %xmm3,%xmm1
movsd  0xc3e(%rip),%xmm3
subsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x5c(%rbp)
mov    -0x20(%rbp),%rax
movss  (%rax),%xmm0
movss  %xmm0,-0x6c(%rbp)
movss  -0x5c(%rbp),%xmm0
call   1050 <sqrt@plt>
movaps %xmm0,%xmm1
movss  -0x6c(%rbp),%xmm0
divss  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
movss  (%rax),%xmm0
call   1030 <fabs@plt>
movsd  0xbe1(%rip),%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
call   1040 <erfcc@plt>
mov    -0x30(%rbp),%rax
movss  %xmm0,(%rax)
add    $0xa0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<kendl1>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %r9,0x10(%rsp)
mov    %r8,%r15
mov    %rcx,%r12
mov    %rdx,%r13
cmp    $0x2,%rdx
jb     11ec <kendl1+0xbc>
mov    $0x1,%ecx
xor    %ebx,%ebx
xorps  %xmm0,%xmm0
xor    %eax,%eax
xor    %r14d,%r14d
jmp    1175 <kendl1+0x45>
nopw   0x0(%rax,%rax,1)
cmp    %r13,%rcx
je     11f3 <kendl1+0xc3>
movss  (%rdi,%rcx,4),%xmm1
movss  (%rsi,%rcx,4),%xmm2
add    $0x1,%rcx
mov    %rcx,%rdx
jmp    119d <kendl1+0x6d>
nopl   0x0(%rax,%rax,1)
add    $0xffffffffffffffff,%r14
add    $0x1,%rdx
cmp    %r13,%rdx
ja     1170 <kendl1+0x40>
movaps %xmm1,%xmm4
subss  (%rdi,%rdx,4),%xmm4
movaps %xmm2,%xmm3
subss  (%rsi,%rdx,4),%xmm3
movaps %xmm4,%xmm5
mulss  %xmm3,%xmm5
ucomiss %xmm0,%xmm5
jne    11bb <kendl1+0x8b>
jnp    11d0 <kendl1+0xa0>
add    $0x1,%rax
add    $0x1,%rbx
ucomiss %xmm0,%xmm5
jbe    1190 <kendl1+0x60>
add    $0x1,%r14
jmp    1194 <kendl1+0x64>
xchg   %ax,%ax
cmpneqss %xmm0,%xmm4
movd   %xmm4,%ebp
and    $0x1,%ebp
add    %rbp,%rax
ucomiss %xmm0,%xmm3
jne    11e6 <kendl1+0xb6>
jnp    1194 <kendl1+0x64>
add    $0x1,%rbx
jmp    1194 <kendl1+0x64>
xor    %r14d,%r14d
xor    %eax,%eax
xor    %ebx,%ebx
movq   %rax,%xmm0
punpckldq 0xe00(%rip),%xmm0
subpd  0xe08(%rip),%xmm0
movapd %xmm0,%xmm1
unpckhpd %xmm0,%xmm1
addsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
xorps  %xmm4,%xmm4
ucomiss %xmm4,%xmm0
jb     122c <kendl1+0xfc>
xorps  %xmm2,%xmm2
sqrtss %xmm0,%xmm2
jmp    1237 <kendl1+0x107>
call   1050 <sqrt@plt>
xorps  %xmm4,%xmm4
movaps %xmm0,%xmm2
xorps  %xmm3,%xmm3
cvtsi2ss %r14,%xmm3
movq   %rbx,%xmm0
punpckldq 0xdb4(%rip),%xmm0
subpd  0xdbc(%rip),%xmm0
movapd %xmm0,%xmm1
unpckhpd %xmm0,%xmm1
addsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
ucomiss %xmm4,%xmm0
jb     1272 <kendl1+0x142>
sqrtss %xmm0,%xmm0
jmp    128f <kendl1+0x15f>
movss  %xmm3,0x8(%rsp)
movss  %xmm2,0xc(%rsp)
call   1050 <sqrt@plt>
movss  0xc(%rsp),%xmm2
movss  0x8(%rsp),%xmm3
mulss  %xmm0,%xmm2
divss  %xmm2,%xmm3
movq   %r13,%xmm0
punpckldq 0xd5c(%rip),%xmm0
subpd  0xd64(%rip),%xmm0
movapd %xmm0,%xmm1
unpckhpd %xmm0,%xmm1
addsd  %xmm0,%xmm1
movsd  0xd60(%rip),%xmm0
mulsd  %xmm1,%xmm0
addsd  0xd5c(%rip),%xmm0
movsd  0xd5c(%rip),%xmm2
mulsd  %xmm1,%xmm2
addsd  0xd58(%rip),%xmm1
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm3,(%r12)
xorpd  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jb     1301 <kendl1+0x1d1>
sqrtss %xmm0,%xmm0
jmp    1312 <kendl1+0x1e2>
movss  %xmm3,0x8(%rsp)
call   1050 <sqrt@plt>
movss  0x8(%rsp),%xmm3
divss  %xmm0,%xmm3
movss  %xmm3,(%r15)
movaps %xmm3,%xmm0
call   1030 <fabs@plt>
divsd  0xd15(%rip),%xmm0
cvtsd2ss %xmm0,%xmm0
call   1040 <erfcc@plt>
mov    0x10(%rsp),%rax
movss  %xmm0,(%rax)
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @kendl1(float* nocapture %0, float* nocapture %1, i64 %2, float* nocapture writeonly %3, float* nocapture writeonly %4, float* nocapture writeonly %5) {
%7 = icmp ugt i64 %2, 1
br i1 %7, label %10, label %54
8: ; preds = %48
%9 = icmp eq i64 %15, %2
br i1 %9, label %54, label %10
10: ; preds = %6, %8
%11 = phi i64 [ %51, %8 ], [ 0, %6 ]
%12 = phi i64 [ %50, %8 ], [ 0, %6 ]
%13 = phi i64 [ %15, %8 ], [ 1, %6 ]
%14 = phi i64 [ %49, %8 ], [ 0, %6 ]
%15 = add nuw i64 %13, 1
%16 = getelementptr inbounds float, float* %0, i64 %13
%17 = load float, float* %16
%18 = getelementptr inbounds float, float* %1, i64 %13
%19 = load float, float* %18
br label %20
20: ; preds = %10, %48
%21 = phi i64 [ %11, %10 ], [ %51, %48 ]
%22 = phi i64 [ %12, %10 ], [ %50, %48 ]
%23 = phi i64 [ %15, %10 ], [ %52, %48 ]
%24 = phi i64 [ %14, %10 ], [ %49, %48 ]
%25 = getelementptr inbounds float, float* %0, i64 %23
%26 = load float, float* %25
%27 = fsub float %17, %26
%28 = getelementptr inbounds float, float* %1, i64 %23
%29 = load float, float* %28
%30 = fsub float %19, %29
%31 = fmul float %27, %30
%32 = fcmp une float %31, 0.000000e+00
br i1 %32, label %33, label %41
33: ; preds = %20
%34 = add i64 %22, 1
%35 = add i64 %21, 1
%36 = fcmp ogt float %31, 0.000000e+00
br i1 %36, label %37, label %39
37: ; preds = %33
%38 = add nsw i64 %24, 1
br label %48
39: ; preds = %33
%40 = add nsw i64 %24, -1
br label %48
41: ; preds = %20
%42 = fcmp une float %27, 0.000000e+00
%43 = zext i1 %42 to i64
%44 = add i64 %22, %43
%45 = fcmp une float %30, 0.000000e+00
br i1 %45, label %46, label %48
46: ; preds = %41
%47 = add i64 %21, 1
br label %48
48: ; preds = %39, %37, %46, %41
%49 = phi i64 [ %38, %37 ], [ %40, %39 ], [ %24, %46 ], [ %24, %41 ]
%50 = phi i64 [ %34, %37 ], [ %34, %39 ], [ %44, %46 ], [ %44, %41 ]
%51 = phi i64 [ %35, %37 ], [ %35, %39 ], [ %47, %46 ], [ %21, %41 ]
%52 = add i64 %23, 1
%53 = icmp ugt i64 %52, %2
br i1 %53, label %8, label %20
54: ; preds = %8, %6
%55 = phi i64 [ 0, %6 ], [ %49, %8 ]
%56 = phi i64 [ 0, %6 ], [ %50, %8 ]
%57 = phi i64 [ 0, %6 ], [ %51, %8 ]
%58 = sitofp i64 %55 to float
%59 = uitofp i64 %56 to double
%60 = fptrunc double %59 to float
%61 = tail call float @sqrt(float %60)
%62 = uitofp i64 %57 to double
%63 = fptrunc double %62 to float
%64 = tail call float @sqrt(float %63)
%65 = fmul float %61, %64
%66 = fdiv float %58, %65
store float %66, float* %3
%67 = uitofp i64 %2 to double
%68 = tail call double @llvm.fmuladd.f64(double %67, double 4.000000e+00, double 1.000000e+01)
%69 = fmul double %67, 9.000000e+00
%70 = fadd double %67, -1.000000e+00
%71 = fmul double %69, %70
%72 = fdiv double %68, %71
%73 = fptrunc double %72 to float
%74 = tail call float @sqrt(float %73)
%75 = fdiv float %66, %74
store float %75, float* %4
%76 = tail call double @fabs(float %75)
%77 = fdiv double %76, 0x3FF6A09E7098EF50
%78 = fptrunc double %77 to float
%79 = tail call float @erfcc(float %78)
store float %79, float* %5
ret void
}