#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 */
float RBR_ipowdist(int *, int *, int, int);
float pow(double, int);
float sqrt(double);
int RBR_dist_stats(int **coords, int ncoords, float *mean, float *sd,
                   float *longest, float *shortest, int PRECISION,
                   int DIST_POWER) {
  int i, j, k;
  int total;
  float sum, sumpow;
  float powdist, dist;
  sum = 0.0;
  sumpow = 0.0;
  (*longest) = 0.0;
  (*shortest) = 100.0;
  total = 0.0;
  for (j = 0; j < ncoords; ++j) {
    for (k = j + 1; k < ncoords; ++k) {
      total++;
      powdist = RBR_ipowdist(coords[j], coords[k], PRECISION, DIST_POWER);
      if (DIST_POWER == 2)
        dist = sqrt((double)powdist);
      else
        dist = pow((double)powdist, (1 / (float)DIST_POWER));
      sum += dist;
      sumpow += powdist;
      if (dist > (*longest))
        (*longest) = dist;
      if (dist < (*shortest))
        (*shortest) = dist;
    }
  }
  (*mean) = (sum / (float)total);
  (*sd) = sqrt((double)((sumpow - (sum * sum) / (float)total) / (float)total));
  return 0;
}


x86 assembly (clang, -O0):
<RBR_dist_stats>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x44(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x48(%rbp)
mov    -0x28(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax)
mov    -0x30(%rbp),%rax
movss  0xe84(%rip),%xmm0
movss  %xmm0,(%rax)
movl   $0x0,-0x40(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12b0 <RBR_dist_stats+0x180>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    129d <RBR_dist_stats+0x16d>
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rdi
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mov    (%rax,%rcx,8),%rsi
mov    0x10(%rbp),%edx
mov    0x18(%rbp),%ecx
call   1030 <RBR_ipowdist@plt>
movss  %xmm0,-0x4c(%rbp)
cmpl   $0x2,0x18(%rbp)
jne    1202 <RBR_dist_stats+0xd2>
movss  -0x4c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
call   1050 <sqrt@plt>
movss  %xmm0,-0x50(%rbp)
jmp    122a <RBR_dist_stats+0xfa>
movss  -0x4c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
cvtsi2ssl 0x18(%rbp),%xmm2
movss  0xdec(%rip),%xmm1
divss  %xmm2,%xmm1
cvttss2si %xmm1,%edi
call   1040 <pow@plt>
movss  %xmm0,-0x50(%rbp)
movss  -0x50(%rbp),%xmm0
addss  -0x44(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
movss  -0x4c(%rbp),%xmm0
addss  -0x48(%rbp),%xmm0
movss  %xmm0,-0x48(%rbp)
movss  -0x50(%rbp),%xmm0
mov    -0x28(%rbp),%rax
ucomiss (%rax),%xmm0
jbe    1267 <RBR_dist_stats+0x137>
movss  -0x50(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movss  %xmm0,(%rax)
movss  -0x50(%rbp),%xmm1
mov    -0x30(%rbp),%rax
movss  (%rax),%xmm0
ucomiss %xmm1,%xmm0
jbe    128a <RBR_dist_stats+0x15a>
movss  -0x50(%rbp),%xmm0
mov    -0x30(%rbp),%rax
movss  %xmm0,(%rax)
jmp    128f <RBR_dist_stats+0x15f>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    11a3 <RBR_dist_stats+0x73>
jmp    12a2 <RBR_dist_stats+0x172>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    118e <RBR_dist_stats+0x5e>
movss  -0x44(%rbp),%xmm0
cvtsi2ssl -0x40(%rbp),%xmm1
divss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movss  %xmm0,(%rax)
movss  -0x48(%rbp),%xmm0
movss  -0x44(%rbp),%xmm1
mulss  -0x44(%rbp),%xmm1
cvtsi2ssl -0x40(%rbp),%xmm2
divss  %xmm2,%xmm1
subss  %xmm1,%xmm0
cvtsi2ssl -0x40(%rbp),%xmm1
divss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
call   1050 <sqrt@plt>
mov    -0x20(%rbp),%rax
movss  %xmm0,(%rax)
xor    %eax,%eax
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<RBR_dist_stats>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
mov    %rcx,0x30(%rsp)
mov    %rdx,0x28(%rsp)
mov    %r8,0x20(%rsp)
movl   $0x0,(%r8)
movl   $0x42c80000,(%r9)
test   %esi,%esi
jle    1291 <RBR_dist_stats+0x161>
mov    %r9,%rbx
cvtsi2ssl 0x88(%rsp),%xmm0
movss  0xe89(%rip),%xmm1
divss  %xmm0,%xmm1
cvttss2si %xmm1,%eax
mov    %eax,0x8(%rsp)
mov    %esi,%eax
lea    0x8(%rdi),%rcx
mov    %rcx,0x18(%rsp)
lea    -0x1(%rax),%r15
xor    %edx,%edx
xorps  %xmm2,%xmm2
xorps  %xmm3,%xmm3
xor    %ecx,%ecx
mov    %rax,0x38(%rsp)
jmp    11d3 <RBR_dist_stats+0xa3>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    0x4(%rsp),%ecx
add    %ebp,%ecx
mov    0x38(%rsp),%rax
mov    0x40(%rsp),%rdx
addq   $0x8,0x18(%rsp)
add    $0xffffffffffffffff,%r15
cmp    %rax,%rdx
je     1299 <RBR_dist_stats+0x169>
mov    %rdx,%r13
add    $0x1,%rdx
cmp    %rax,%rdx
jae    11c0 <RBR_dist_stats+0x90>
mov    %rdx,0x40(%rsp)
mov    %ecx,0x4(%rsp)
xor    %ebp,%ebp
jmp    1203 <RBR_dist_stats+0xd3>
nopl   0x0(%rax)
addss  %xmm0,%xmm3
addss  0xc(%rsp),%xmm2
add    $0x1,%rbp
cmp    %rbp,%r15
je     11b0 <RBR_dist_stats+0x80>
movss  %xmm3,0x10(%rsp)
movss  %xmm2,0x14(%rsp)
mov    %rdi,%r12
mov    (%rdi,%r13,8),%rdi
mov    0x18(%rsp),%rax
mov    (%rax,%rbp,8),%rsi
mov    0x80(%rsp),%edx
mov    0x88(%rsp),%r14d
mov    %r14d,%ecx
call   1030 <RBR_ipowdist@plt>
movss  %xmm0,0xc(%rsp)
cvtss2sd %xmm0,%xmm0
cmp    $0x2,%r14d
jne    1250 <RBR_dist_stats+0x120>
call   1050 <sqrt@plt>
jmp    1259 <RBR_dist_stats+0x129>
nopl   (%rax)
mov    0x8(%rsp),%edi
call   1040 <pow@plt>
mov    0x20(%rsp),%rax
ucomiss (%rax),%xmm0
movss  0x14(%rsp),%xmm2
mov    %r12,%rdi
movss  0x10(%rsp),%xmm3
jbe    127b <RBR_dist_stats+0x14b>
mov    0x20(%rsp),%rax
movss  %xmm0,(%rax)
movss  (%rbx),%xmm1
ucomiss %xmm0,%xmm1
jbe    11f0 <RBR_dist_stats+0xc0>
movss  %xmm0,(%rbx)
jmp    11f0 <RBR_dist_stats+0xc0>
xorps  %xmm3,%xmm3
xor    %ecx,%ecx
xorps  %xmm2,%xmm2
xorps  %xmm0,%xmm0
cvtsi2ss %ecx,%xmm0
movaps %xmm3,%xmm1
divss  %xmm0,%xmm1
mov    0x28(%rsp),%rax
movss  %xmm1,(%rax)
mulss  %xmm3,%xmm3
divss  %xmm0,%xmm3
subss  %xmm3,%xmm2
divss  %xmm0,%xmm2
xorps  %xmm0,%xmm0
cvtss2sd %xmm2,%xmm0
call   1050 <sqrt@plt>
mov    0x30(%rsp),%rax
movss  %xmm0,(%rax)
xor    %eax,%eax
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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