#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 pow(double, int);
double quantile(double *, int, double);
void distestim(double *x, int *n, double *smoothing, double *value) {
  int i, j;
  int count;
  int peakCount = 0;
  double lim[2];
  double range;
  double peakRange = 0;
  double M;
  double D[*n];
  double sumD = 0;
  double xBin[*n];
  if (*smoothing == 0) {
    *smoothing = quantile(x, *n, 0.975);
    *smoothing -= quantile(x, *n, 0.025);
    *smoothing /= pow(*n, 0.5);
  }
  *smoothing *= 0.5;
  for (i = 0; i < *n; i++) {
    lim[0] = x[i] - (*smoothing);
    lim[1] = x[i] + (*smoothing);
    sumD = 0;
    count = 0;
    for (j = 0; j < *n; j++) {
      if ((x[j] >= lim[0]) & (x[j] <= lim[1])) {
        xBin[count] = x[j];
        D[count] = pow(x[j] - x[i], 2);
        if (D[count] < 0)
          D[count] *= -1;
        sumD += D[count];
        count++;
      }
    }
    if (count >= peakCount) {
      range = sumD;
      for (j = 0; j < count; j++)
        D[j] = 1 - D[j] / sumD;
      sumD = 0;
      M = 0;
      for (j = 0; j < count; j++) {
        M += xBin[j] * D[j];
        sumD += D[j];
      }
      M /= sumD;
      if ((count > peakCount) || (range < peakRange)) {
        *value = M;
        peakCount = count;
        peakRange = range;
      }
    }
  }
}


x86 assembly (clang, -O0):
<distestim>:
push   %rbp
mov    %rsp,%rbp
sub    $0x90,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movl   $0x0,-0x30(%rbp)
movq   $0x0,-0x50(%rbp)
mov    -0x10(%rbp),%rax
mov    (%rax),%eax
mov    %rsp,%rcx
mov    %rcx,-0x60(%rbp)
lea    0xf(,%rax,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rcx
sub    %rdx,%rcx
mov    %rcx,-0x88(%rbp)
mov    %rcx,%rsp
mov    %rax,-0x68(%rbp)
movq   $0x0,-0x70(%rbp)
mov    -0x10(%rbp),%rax
mov    (%rax),%eax
lea    0xf(,%rax,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rcx
sub    %rdx,%rcx
mov    %rcx,-0x80(%rbp)
mov    %rcx,%rsp
mov    %rax,-0x78(%rbp)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1228 <distestim+0x108>
jp     1228 <distestim+0x108>
mov    -0x8(%rbp),%rdi
mov    -0x10(%rbp),%rax
mov    (%rax),%esi
movsd  0xe39(%rip),%xmm0
call   1040 <quantile@plt>
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x8(%rbp),%rdi
mov    -0x10(%rbp),%rax
mov    (%rax),%esi
movsd  0xe12(%rip),%xmm0
call   1040 <quantile@plt>
movaps %xmm0,%xmm1
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
mov    -0x10(%rbp),%rax
cvtsi2sdl (%rax),%xmm0
xor    %edi,%edi
call   1030 <pow@plt>
movaps %xmm0,%xmm1
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
mov    -0x18(%rbp),%rax
movsd  0xddc(%rip),%xmm0
mulsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x10(%rbp),%rcx
cmp    (%rcx),%eax
jge    14de <distestim+0x3be>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
subsd  (%rax),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
addsd  (%rax),%xmm0
movsd  %xmm0,-0x38(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x70(%rbp)
movl   $0x0,-0x2c(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x10(%rbp),%rcx
cmp    (%rcx),%eax
jge    13a2 <distestim+0x282>
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x40(%rbp),%xmm0
setae  %al
and    $0x1,%al
movzbl %al,%eax
mov    -0x8(%rbp),%rcx
movslq -0x28(%rbp),%rdx
movsd  (%rcx,%rdx,8),%xmm1
movsd  -0x38(%rbp),%xmm0
ucomisd %xmm1,%xmm0
setae  %cl
and    $0x1,%cl
movzbl %cl,%ecx
and    %ecx,%eax
cmp    $0x0,%eax
je     138f <distestim+0x26f>
mov    -0x80(%rbp),%rax
mov    -0x8(%rbp),%rcx
movslq -0x28(%rbp),%rdx
movsd  (%rcx,%rdx,8),%xmm0
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    $0x2,%edi
call   1030 <pow@plt>
mov    -0x88(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movslq -0x2c(%rbp),%rcx
xorps  %xmm0,%xmm0
ucomisd (%rax,%rcx,8),%xmm0
jbe    136c <distestim+0x24c>
mov    -0x88(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  0xcbe(%rip),%xmm0
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x88(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1394 <distestim+0x274>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    129c <distestim+0x17c>
mov    -0x2c(%rbp),%eax
cmp    -0x30(%rbp),%eax
jl     14cb <distestim+0x3ab>
movsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    1403 <distestim+0x2e3>
mov    -0x88(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
divsd  -0x70(%rbp),%xmm1
movsd  0xc30(%rip),%xmm0
subsd  %xmm1,%xmm0
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    13bf <distestim+0x29f>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x70(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    1476 <distestim+0x356>
mov    -0x88(%rbp),%rax
mov    -0x80(%rbp),%rcx
movslq -0x28(%rbp),%rdx
movsd  (%rcx,%rdx,8),%xmm0
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x58(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    141a <distestim+0x2fa>
movsd  -0x70(%rbp),%xmm1
movsd  -0x58(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x30(%rbp),%eax
jg     14a9 <distestim+0x389>
movsd  -0x48(%rbp),%xmm1
movsd  -0x50(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    14c6 <distestim+0x3a6>
movsd  -0x58(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x30(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
jmp    14cb <distestim+0x3ab>
jmp    14d0 <distestim+0x3b0>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1243 <distestim+0x123>
mov    -0x60(%rbp),%rax
mov    %rax,%rsp
mov    %rbp,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<distestim>:
push   %rbp
mov    %rsp,%rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
mov    %rcx,-0x50(%rbp)
mov    %rdx,%r14
mov    %rsi,%r12
mov    %rdi,%r13
mov    (%rsi),%esi
mov    %rsp,%r15
lea    0xf(,%rsi,8),%rax
and    $0xfffffffffffffff0,%rax
sub    %rax,%r15
mov    %r15,%rsp
mov    %rsp,%rbx
sub    %rax,%rbx
mov    %rbx,%rsp
movsd  (%rdx),%xmm1
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm1
mov    %r12,-0x58(%rbp)
jne    11c6 <distestim+0xa6>
jp     11c6 <distestim+0xa6>
movsd  0xe86(%rip),%xmm0
mov    %r13,%rdi
call   1040 <quantile@plt>
movsd  %xmm0,(%r14)
mov    (%r12),%esi
movsd  0xe75(%rip),%xmm0
mov    %r13,%rdi
call   1040 <quantile@plt>
movsd  (%r14),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,(%r14)
xorps  %xmm0,%xmm0
cvtsi2sdl (%r12),%xmm0
xor    %edi,%edi
call   1030 <pow@plt>
movsd  (%r14),%xmm1
divsd  %xmm0,%xmm1
mov    (%r12),%esi
mulsd  0xe42(%rip),%xmm1
movsd  %xmm1,(%r14)
test   %esi,%esi
jle    156e <distestim+0x44e>
xor    %eax,%eax
mov    %rax,-0x40(%rbp)
movsd  0xe2f(%rip),%xmm5
movapd 0xe4f(%rip),%xmm6
xorpd  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
xor    %eax,%eax
mov    %r14,-0x48(%rbp)
test   %esi,%esi
jg     1230 <distestim+0x110>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
test   %eax,%eax
jg     1512 <distestim+0x3f2>
xor    %r12d,%r12d
xorpd  %xmm7,%xmm7
xorpd  %xmm0,%xmm0
jmp    1350 <distestim+0x230>
nopl   0x0(%rax,%rax,1)
mov    %eax,-0x2c(%rbp)
mov    -0x40(%rbp),%rax
movsd  0x0(%r13,%rax,8),%xmm0
movapd %xmm0,%xmm2
subsd  %xmm1,%xmm2
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
xor    %r14d,%r14d
xorpd  %xmm7,%xmm7
xor    %r12d,%r12d
movsd  %xmm2,-0x60(%rbp)
jmp    1270 <distestim+0x150>
add    $0x1,%r14
movslq %esi,%rax
cmp    %rax,%r14
jge    1310 <distestim+0x1f0>
movsd  0x0(%r13,%r14,8),%xmm0
ucomisd %xmm2,%xmm0
jb     1260 <distestim+0x140>
movsd  -0x68(%rbp),%xmm1
ucomisd %xmm0,%xmm1
jb     1260 <distestim+0x140>
movslq %r12d,%r12
movsd  %xmm0,(%rbx,%r12,8)
mov    -0x40(%rbp),%rax
subsd  0x0(%r13,%rax,8),%xmm0
mov    $0x2,%edi
movapd %xmm7,-0x80(%rbp)
call   1030 <pow@plt>
movapd -0x80(%rbp),%xmm7
movapd 0xd88(%rip),%xmm6
movsd  0xd58(%rip),%xmm5
movapd %xmm0,%xmm1
xorpd  0xd64(%rip),%xmm1
movapd %xmm0,%xmm2
cmpltsd 0xd47(%rip),%xmm2
andpd  %xmm2,%xmm1
andnpd %xmm0,%xmm2
orpd   %xmm1,%xmm2
movlpd %xmm2,(%r15,%r12,8)
addsd  %xmm2,%xmm7
movsd  -0x60(%rbp),%xmm2
add    $0x1,%r12d
mov    -0x58(%rbp),%rax
mov    (%rax),%esi
jmp    1260 <distestim+0x140>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    -0x2c(%rbp),%eax
cmp    %eax,%r12d
jge    1330 <distestim+0x210>
mov    %eax,%r12d
movsd  -0x38(%rbp),%xmm7
mov    -0x48(%rbp),%r14
jmp    152c <distestim+0x40c>
nopl   0x0(%rax)
xorpd  %xmm0,%xmm0
test   %r12d,%r12d
mov    -0x48(%rbp),%r14
jle    1350 <distestim+0x230>
mov    %r12d,%eax
cmp    $0x1,%r12d
jne    1370 <distestim+0x250>
xor    %ecx,%ecx
jmp    1400 <distestim+0x2e0>
nopl   (%rax)
xorpd  %xmm1,%xmm1
cmp    %eax,%r12d
jle    1507 <distestim+0x3e7>
jmp    1520 <distestim+0x400>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    %eax,%ecx
and    $0xfffffffe,%ecx
movapd %xmm7,%xmm1
unpcklpd %xmm7,%xmm1
lea    -0x2(%rcx),%rdx
mov    %rdx,%r8
shr    %r8
add    $0x1,%r8
test   %rdx,%rdx
je     155d <distestim+0x43d>
mov    %r8,%rdx
and    $0xfffffffffffffffe,%rdx
xor    %edi,%edi
nopl   (%rax)
movapd (%r15,%rdi,8),%xmm2
movapd 0x10(%r15,%rdi,8),%xmm3
divpd  %xmm1,%xmm2
movapd %xmm6,%xmm4
subpd  %xmm2,%xmm4
movapd %xmm4,(%r15,%rdi,8)
divpd  %xmm1,%xmm3
movapd %xmm6,%xmm2
subpd  %xmm3,%xmm2
movapd %xmm2,0x10(%r15,%rdi,8)
add    $0x4,%rdi
add    $0xfffffffffffffffe,%rdx
jne    13a0 <distestim+0x280>
test   $0x1,%r8b
je     13fa <distestim+0x2da>
movapd (%r15,%rdi,8),%xmm2
divpd  %xmm1,%xmm2
movapd %xmm6,%xmm1
subpd  %xmm2,%xmm1
movapd %xmm1,(%r15,%rdi,8)
cmp    %rax,%rcx
je     1421 <distestim+0x301>
nop
movsd  (%r15,%rcx,8),%xmm1
divsd  %xmm7,%xmm1
movapd %xmm5,%xmm2
subsd  %xmm1,%xmm2
movsd  %xmm2,(%r15,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%rax
jne    1400 <distestim+0x2e0>
test   %r12d,%r12d
jle    1444 <distestim+0x324>
lea    -0x1(%rax),%rdx
mov    %eax,%ecx
and    $0x3,%ecx
cmp    $0x3,%rdx
jae    144d <distestim+0x32d>
xorpd  %xmm1,%xmm1
xor    %edx,%edx
xorpd  %xmm0,%xmm0
jmp    14c3 <distestim+0x3a3>
xorpd  %xmm1,%xmm1
jmp    14ff <distestim+0x3df>
and    $0xfffffffc,%eax
xorpd  %xmm1,%xmm1
xor    %edx,%edx
xorpd  %xmm0,%xmm0
nopw   0x0(%rax,%rax,1)
movsd  (%r15,%rdx,8),%xmm2
movsd  0x8(%r15,%rdx,8),%xmm3
addsd  %xmm2,%xmm1
mulsd  (%rbx,%rdx,8),%xmm2
addsd  %xmm0,%xmm2
addsd  %xmm3,%xmm1
mulsd  0x8(%rbx,%rdx,8),%xmm3
addsd  %xmm2,%xmm3
movsd  0x10(%r15,%rdx,8),%xmm0
movsd  0x10(%rbx,%rdx,8),%xmm2
mulsd  %xmm0,%xmm2
addsd  %xmm3,%xmm2
addsd  %xmm0,%xmm1
movsd  0x18(%r15,%rdx,8),%xmm3
movsd  0x18(%rbx,%rdx,8),%xmm0
mulsd  %xmm3,%xmm0
addsd  %xmm2,%xmm0
addsd  %xmm3,%xmm1
add    $0x4,%rdx
cmp    %rdx,%rax
jne    1460 <distestim+0x340>
test   %rcx,%rcx
je     14ff <distestim+0x3df>
lea    (%r15,%rdx,8),%rax
lea    (%rbx,%rdx,8),%rdx
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  (%rax,%rdi,8),%xmm2
movsd  (%rdx,%rdi,8),%xmm3
mulsd  %xmm2,%xmm3
addsd  %xmm3,%xmm0
addsd  %xmm2,%xmm1
add    $0x1,%rdi
cmp    %rdi,%rcx
jne    14e0 <distestim+0x3c0>
mov    -0x2c(%rbp),%eax
cmp    %eax,%r12d
jg     1520 <distestim+0x400>
movsd  -0x38(%rbp),%xmm2
ucomisd %xmm7,%xmm2
ja     1520 <distestim+0x400>
mov    %eax,%r12d
movsd  -0x38(%rbp),%xmm7
jmp    152c <distestim+0x40c>
nopl   0x0(%rax)
divsd  %xmm1,%xmm0
mov    -0x50(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x40(%rbp),%rdx
add    $0x1,%rdx
movslq %esi,%rax
mov    %rdx,%rcx
mov    %rdx,-0x40(%rbp)
cmp    %rax,%rdx
jge    156e <distestim+0x44e>
movsd  (%r14),%xmm1
movsd  %xmm7,-0x38(%rbp)
mov    %r12d,%eax
test   %esi,%esi
jg     1230 <distestim+0x110>
jmp    1210 <distestim+0xf0>
xor    %edi,%edi
test   $0x1,%r8b
jne    13e2 <distestim+0x2c2>
jmp    13fa <distestim+0x2da>
lea    -0x28(%rbp),%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @distestim(double* %0, i32* %1, double* %2, double* %3) {
%5 = alloca double*
%6 = alloca i32*
%7 = alloca double*
%8 = alloca double*
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca [2 x double]
%14 = alloca double
%15 = alloca double
%16 = alloca double
%17 = alloca i8*
%18 = alloca i64
%19 = alloca double
%20 = alloca i64
store double* %0, double** %5
store i32* %1, i32** %6
store double* %2, double** %7
store double* %3, double** %8
store i32 0, i32* %12
store double 0.000000e+00, double* %15
%21 = load i32*, i32** %6
%22 = load i32, i32* %21
%23 = zext i32 %22 to i64
%24 = call i8* @llvm.stacksave()
store i8* %24, i8** %17
%25 = alloca double, i64 %23
store i64 %23, i64* %18
store double 0.000000e+00, double* %19
%26 = load i32*, i32** %6
%27 = load i32, i32* %26
%28 = zext i32 %27 to i64
%29 = alloca double, i64 %28
store i64 %28, i64* %20
%30 = load double*, double** %7
%31 = load double, double* %30
%32 = fcmp oeq double %31, 0.000000e+00
br i1 %32, label %33, label %53
33: ; preds = %4
%34 = load double*, double** %5
%35 = load i32*, i32** %6
%36 = load i32, i32* %35
%37 = call double @quantile(double* %34, i32 %36, double 0x3FEF333333333333)
%38 = load double*, double** %7
store double %37, double* %38
%39 = load double*, double** %5
%40 = load i32*, i32** %6
%41 = load i32, i32* %40
%42 = call double @quantile(double* %39, i32 %41, double 2.500000e-02)
%43 = load double*, double** %7
%44 = load double, double* %43
%45 = fsub double %44, %42
store double %45, double* %43
%46 = load i32*, i32** %6
%47 = load i32, i32* %46
%48 = sitofp i32 %47 to double
%49 = call double @pow(double %48, i32 0)
%50 = load double*, double** %7
%51 = load double, double* %50
%52 = fdiv double %51, %49
store double %52, double* %50
br label %53
53: ; preds = %33, %4
%54 = load double*, double** %7
%55 = load double, double* %54
%56 = fmul double %55, 5.000000e-01
store double %56, double* %54
store i32 0, i32* %9
br label %57
57: ; preds = %222, %53
%58 = load i32, i32* %9
%59 = load i32*, i32** %6
%60 = load i32, i32* %59
%61 = icmp slt i32 %58, %60
br i1 %61, label %62, label %225
62: ; preds = %57
%63 = load double*, double** %5
%64 = load i32, i32* %9
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds double, double* %63, i64 %65
%67 = load double, double* %66
%68 = load double*, double** %7
%69 = load double, double* %68
%70 = fsub double %67, %69
%71 = getelementptr inbounds [2 x double], [2 x double]* %13, i64 0, i64 0
store double %70, double* %71
%72 = load double*, double** %5
%73 = load i32, i32* %9
%74 = sext i32 %73 to i64
%75 = getelementptr inbounds double, double* %72, i64 %74
%76 = load double, double* %75
%77 = load double*, double** %7
%78 = load double, double* %77
%79 = fadd double %76, %78
%80 = getelementptr inbounds [2 x double], [2 x double]* %13, i64 0, i64 1
store double %79, double* %80
store double 0.000000e+00, double* %19
store i32 0, i32* %11
store i32 0, i32* %10
br label %81
81: ; preds = %152, %62
%82 = load i32, i32* %10
%83 = load i32*, i32** %6
%84 = load i32, i32* %83
%85 = icmp slt i32 %82, %84
br i1 %85, label %86, label %155
86: ; preds = %81
%87 = load double*, double** %5
%88 = load i32, i32* %10
%89 = sext i32 %88 to i64
%90 = getelementptr inbounds double, double* %87, i64 %89
%91 = load double, double* %90
%92 = getelementptr inbounds [2 x double], [2 x double]* %13, i64 0, i64 0
%93 = load double, double* %92
%94 = fcmp oge double %91, %93
%95 = zext i1 %94 to i32
%96 = load double*, double** %5
%97 = load i32, i32* %10
%98 = sext i32 %97 to i64
%99 = getelementptr inbounds double, double* %96, i64 %98
%100 = load double, double* %99
%101 = getelementptr inbounds [2 x double], [2 x double]* %13, i64 0, i64 1
%102 = load double, double* %101
%103 = fcmp ole double %100, %102
%104 = zext i1 %103 to i32
%105 = and i32 %95, %104
%106 = icmp ne i32 %105, 0
br i1 %106, label %107, label %151
107: ; preds = %86
%108 = load double*, double** %5
%109 = load i32, i32* %10
%110 = sext i32 %109 to i64
%111 = getelementptr inbounds double, double* %108, i64 %110
%112 = load double, double* %111
%113 = load i32, i32* %11
%114 = sext i32 %113 to i64
%115 = getelementptr inbounds double, double* %29, i64 %114
store double %112, double* %115
%116 = load double*, double** %5
%117 = load i32, i32* %10
%118 = sext i32 %117 to i64
%119 = getelementptr inbounds double, double* %116, i64 %118
%120 = load double, double* %119
%121 = load double*, double** %5
%122 = load i32, i32* %9
%123 = sext i32 %122 to i64
%124 = getelementptr inbounds double, double* %121, i64 %123
%125 = load double, double* %124
%126 = fsub double %120, %125
%127 = call double @pow(double %126, i32 2)
%128 = load i32, i32* %11
%129 = sext i32 %128 to i64
%130 = getelementptr inbounds double, double* %25, i64 %129
store double %127, double* %130
%131 = load i32, i32* %11
%132 = sext i32 %131 to i64
%133 = getelementptr inbounds double, double* %25, i64 %132
%134 = load double, double* %133
%135 = fcmp olt double %134, 0.000000e+00
br i1 %135, label %136, label %142
136: ; preds = %107
%137 = load i32, i32* %11
%138 = sext i32 %137 to i64
%139 = getelementptr inbounds double, double* %25, i64 %138
%140 = load double, double* %139
%141 = fmul double %140, -1.000000e+00
store double %141, double* %139
br label %142
142: ; preds = %136, %107
%143 = load i32, i32* %11
%144 = sext i32 %143 to i64
%145 = getelementptr inbounds double, double* %25, i64 %144
%146 = load double, double* %145
%147 = load double, double* %19
%148 = fadd double %147, %146
store double %148, double* %19
%149 = load i32, i32* %11
%150 = add nsw i32 %149, 1
store i32 %150, i32* %11
br label %151
151: ; preds = %142, %86
br label %152
152: ; preds = %151
%153 = load i32, i32* %10
%154 = add nsw i32 %153, 1
store i32 %154, i32* %10
br label %81
155: ; preds = %81
%156 = load i32, i32* %11
%157 = load i32, i32* %12
%158 = icmp sge i32 %156, %157
br i1 %158, label %159, label %221
159: ; preds = %155
%160 = load double, double* %19
store double %160, double* %14
store i32 0, i32* %10
br label %161
161: ; preds = %176, %159
%162 = load i32, i32* %10
%163 = load i32, i32* %11
%164 = icmp slt i32 %162, %163
br i1 %164, label %165, label %179
165: ; preds = %161
%166 = load i32, i32* %10
%167 = sext i32 %166 to i64
%168 = getelementptr inbounds double, double* %25, i64 %167
%169 = load double, double* %168
%170 = load double, double* %19
%171 = fdiv double %169, %170
%172 = fsub double 1.000000e+00, %171
%173 = load i32, i32* %10
%174 = sext i32 %173 to i64
%175 = getelementptr inbounds double, double* %25, i64 %174
store double %172, double* %175
br label %176
176: ; preds = %165
%177 = load i32, i32* %10
%178 = add nsw i32 %177, 1
store i32 %178, i32* %10
br label %161
179: ; preds = %161
store double 0.000000e+00, double* %19
store double 0.000000e+00, double* %16
store i32 0, i32* %10
br label %180
180: ; preds = %201, %179
%181 = load i32, i32* %10
%182 = load i32, i32* %11
%183 = icmp slt i32 %181, %182
br i1 %183, label %184, label %204
184: ; preds = %180
%185 = load i32, i32* %10
%186 = sext i32 %185 to i64
%187 = getelementptr inbounds double, double* %29, i64 %186
%188 = load double, double* %187
%189 = load i32, i32* %10
%190 = sext i32 %189 to i64
%191 = getelementptr inbounds double, double* %25, i64 %190
%192 = load double, double* %191
%193 = load double, double* %16
%194 = call double @llvm.fmuladd.f64(double %188, double %192, double %193)
store double %194, double* %16
%195 = load i32, i32* %10
%196 = sext i32 %195 to i64
%197 = getelementptr inbounds double, double* %25, i64 %196
%198 = load double, double* %197
%199 = load double, double* %19
%200 = fadd double %199, %198
store double %200, double* %19
br label %201
201: ; preds = %184
%202 = load i32, i32* %10
%203 = add nsw i32 %202, 1
store i32 %203, i32* %10
br label %180
204: ; preds = %180
%205 = load double, double* %19
%206 = load double, double* %16
%207 = fdiv double %206, %205
store double %207, double* %16
%208 = load i32, i32* %11
%209 = load i32, i32* %12
%210 = icmp sgt i32 %208, %209
br i1 %210, label %215, label %211
211: ; preds = %204
%212 = load double, double* %14
%213 = load double, double* %15
%214 = fcmp olt double %212, %213
br i1 %214, label %215, label %220
215: ; preds = %211, %204
%216 = load double, double* %16
%217 = load double*, double** %8
store double %216, double* %217
%218 = load i32, i32* %11
store i32 %218, i32* %12
%219 = load double, double* %14
store double %219, double* %15
br label %220
220: ; preds = %215, %211
br label %221
221: ; preds = %220, %155
br label %222
222: ; preds = %221
%223 = load i32, i32* %9
%224 = add nsw i32 %223, 1
store i32 %224, i32* %9
br label %57
225: ; preds = %57
%226 = load i8*, i8** %17
ret void
}

LLVM IR (-O2):
define void @distestim(double* %0, i32* nocapture %1, double* nocapture %2, double* nocapture writeonly %3) {
%5 = load i32, i32* %1
%6 = zext i32 %5 to i64
%7 = alloca double, i64 %6
%8 = alloca double, i64 %6
%9 = load double, double* %2
%10 = fcmp oeq double %9, 0.000000e+00
br i1 %10, label %11, label %23
11: ; preds = %4
%12 = tail call double @quantile(double* %0, i32 %5, double 0x3FEF333333333333)
store double %12, double* %2
%13 = load i32, i32* %1
%14 = tail call double @quantile(double* %0, i32 %13, double 2.500000e-02)
%15 = load double, double* %2
%16 = fsub double %15, %14
store double %16, double* %2
%17 = load i32, i32* %1
%18 = sitofp i32 %17 to double
%19 = tail call double @pow(double %18, i32 0)
%20 = load double, double* %2
%21 = fdiv double %20, %19
%22 = load i32, i32* %1
br label %23
23: ; preds = %11, %4
%24 = phi i32 [ %22, %11 ], [ %5, %4 ]
%25 = phi double [ %21, %11 ], [ %9, %4 ]
%26 = fmul double %25, 5.000000e-01
store double %26, double* %2
%27 = icmp sgt i32 %24, 0
br i1 %27, label %28, label %215
28: ; preds = %23, %213
%29 = phi i32 [ %207, %213 ], [ %24, %23 ]
%30 = phi double [ %214, %213 ], [ %26, %23 ]
%31 = phi i64 [ %210, %213 ], [ 0, %23 ]
%32 = phi double [ %209, %213 ], [ 0.000000e+00, %23 ]
%33 = phi i32 [ %208, %213 ], [ 0, %23 ]
%34 = getelementptr inbounds double, double* %0, i64 %31
%35 = load double, double* %34
%36 = fsub double %35, %30
%37 = fadd double %35, %30
%38 = icmp sgt i32 %29, 0
br i1 %38, label %39, label %73
39: ; preds = %28, %63
%40 = phi i32 [ %64, %63 ], [ %29, %28 ]
%41 = phi i32 [ %65, %63 ], [ %29, %28 ]
%42 = phi i64 [ %68, %63 ], [ 0, %28 ]
%43 = phi double [ %67, %63 ], [ 0.000000e+00, %28 ]
%44 = phi i32 [ %66, %63 ], [ 0, %28 ]
%45 = getelementptr inbounds double, double* %0, i64 %42
%46 = load double, double* %45
%47 = fcmp oge double %46, %36
%48 = fcmp ole double %46, %37
%49 = and i1 %47, %48
br i1 %49, label %50, label %63
50: ; preds = %39
%51 = sext i32 %44 to i64
%52 = getelementptr inbounds double, double* %8, i64 %51
store double %46, double* %52
%53 = load double, double* %34
%54 = fsub double %46, %53
%55 = tail call double @pow(double %54, i32 2)
%56 = getelementptr inbounds double, double* %7, i64 %51
%57 = fcmp olt double %55, 0.000000e+00
%58 = fneg double %55
%59 = select i1 %57, double %58, double %55
store double %59, double* %56
%60 = fadd double %43, %59
%61 = add nsw i32 %44, 1
%62 = load i32, i32* %1
br label %63
63: ; preds = %39, %50
%64 = phi i32 [ %62, %50 ], [ %40, %39 ]
%65 = phi i32 [ %62, %50 ], [ %41, %39 ]
%66 = phi i32 [ %61, %50 ], [ %44, %39 ]
%67 = phi double [ %60, %50 ], [ %43, %39 ]
%68 = add nuw nsw i64 %42, 1
%69 = sext i32 %65 to i64
%70 = icmp slt i64 %68, %69
br i1 %70, label %39, label %71
71: ; preds = %63
%72 = icmp slt i32 %66, %33
br i1 %72, label %206, label %75
73: ; preds = %28
%74 = icmp sgt i32 %33, 0
br i1 %74, label %206, label %195
75: ; preds = %71
%76 = icmp sgt i32 %66, 0
br i1 %76, label %77, label %195
77: ; preds = %75
%78 = zext i32 %66 to i64
%79 = icmp eq i32 %66, 1
br i1 %79, label %122, label %80
80: ; preds = %77
%81 = and i64 %78, 4294967294
%82 = insertelement <2 x double> poison, double %67, i64 0
%83 = shufflevector <2 x double> %82, <2 x double> poison, <2 x i32> zeroinitializer
%84 = add nsw i64 %81, -2
%85 = lshr exact i64 %84, 1
%86 = add nuw i64 %85, 1
%87 = and i64 %86, 1
%88 = icmp eq i64 %84, 0
br i1 %88, label %110, label %89
89: ; preds = %80
%90 = and i64 %86, -2
br label %91
91: ; preds = %91, %89
%92 = phi i64 [ 0, %89 ], [ %107, %91 ]
%93 = phi i64 [ 0, %89 ], [ %108, %91 ]
%94 = getelementptr inbounds double, double* %7, i64 %92
%95 = bitcast double* %94 to <2 x double>*
%96 = load <2 x double>, <2 x double>* %95
%97 = fdiv <2 x double> %96, %83
%98 = fsub <2 x double> <double 1.000000e+00, double 1.000000e+00>, %97
%99 = bitcast double* %94 to <2 x double>*
store <2 x double> %98, <2 x double>* %99
%100 = or i64 %92, 2
%101 = getelementptr inbounds double, double* %7, i64 %100
%102 = bitcast double* %101 to <2 x double>*
%103 = load <2 x double>, <2 x double>* %102
%104 = fdiv <2 x double> %103, %83
%105 = fsub <2 x double> <double 1.000000e+00, double 1.000000e+00>, %104
%106 = bitcast double* %101 to <2 x double>*
store <2 x double> %105, <2 x double>* %106
%107 = add nuw i64 %92, 4
%108 = add nuw i64 %93, 2
%109 = icmp eq i64 %108, %90
br i1 %109, label %110, label %91
110: ; preds = %91, %80
%111 = phi i64 [ 0, %80 ], [ %107, %91 ]
%112 = icmp eq i64 %87, 0
br i1 %112, label %120, label %113
113: ; preds = %110
%114 = getelementptr inbounds double, double* %7, i64 %111
%115 = bitcast double* %114 to <2 x double>*
%116 = load <2 x double>, <2 x double>* %115
%117 = fdiv <2 x double> %116, %83
%118 = fsub <2 x double> <double 1.000000e+00, double 1.000000e+00>, %117
%119 = bitcast double* %114 to <2 x double>*
store <2 x double> %118, <2 x double>* %119
br label %120
120: ; preds = %110, %113
%121 = icmp eq i64 %81, %78
br i1 %121, label %124, label %122
122: ; preds = %77, %120
%123 = phi i64 [ 0, %77 ], [ %81, %120 ]
br label %131
124: ; preds = %131, %120
br i1 %76, label %125, label %195
125: ; preds = %124
%126 = add nsw i64 %78, -1
%127 = and i64 %78, 3
%128 = icmp ult i64 %126, 3
br i1 %128, label %174, label %129
129: ; preds = %125
%130 = and i64 %78, 4294967292
br label %139
131: ; preds = %122, %131
%132 = phi i64 [ %137, %131 ], [ %123, %122 ]
%133 = getelementptr inbounds double, double* %7, i64 %132
%134 = load double, double* %133
%135 = fdiv double %134, %67
%136 = fsub double 1.000000e+00, %135
store double %136, double* %133
%137 = add nuw nsw i64 %132, 1
%138 = icmp eq i64 %137, %78
br i1 %138, label %124, label %131
139: ; preds = %139, %129
%140 = phi i64 [ 0, %129 ], [ %171, %139 ]
%141 = phi double [ 0.000000e+00, %129 ], [ %170, %139 ]
%142 = phi double [ 0.000000e+00, %129 ], [ %169, %139 ]
%143 = phi i64 [ 0, %129 ], [ %172, %139 ]
%144 = getelementptr inbounds double, double* %8, i64 %140
%145 = load double, double* %144
%146 = getelementptr inbounds double, double* %7, i64 %140
%147 = load double, double* %146
%148 = tail call double @llvm.fmuladd.f64(double %145, double %147, double %142)
%149 = fadd double %141, %147
%150 = or i64 %140, 1
%151 = getelementptr inbounds double, double* %8, i64 %150
%152 = load double, double* %151
%153 = getelementptr inbounds double, double* %7, i64 %150
%154 = load double, double* %153
%155 = tail call double @llvm.fmuladd.f64(double %152, double %154, double %148)
%156 = fadd double %149, %154
%157 = or i64 %140, 2
%158 = getelementptr inbounds double, double* %8, i64 %157
%159 = load double, double* %158
%160 = getelementptr inbounds double, double* %7, i64 %157
%161 = load double, double* %160
%162 = tail call double @llvm.fmuladd.f64(double %159, double %161, double %155)
%163 = fadd double %156, %161
%164 = or i64 %140, 3
%165 = getelementptr inbounds double, double* %8, i64 %164
%166 = load double, double* %165
%167 = getelementptr inbounds double, double* %7, i64 %164
%168 = load double, double* %167
%169 = tail call double @llvm.fmuladd.f64(double %166, double %168, double %162)
%170 = fadd double %163, %168
%171 = add nuw nsw i64 %140, 4
%172 = add i64 %143, 4
%173 = icmp eq i64 %172, %130
br i1 %173, label %174, label %139
174: ; preds = %139, %125
%175 = phi double [ undef, %125 ], [ %169, %139 ]
%176 = phi double [ undef, %125 ], [ %170, %139 ]
%177 = phi i64 [ 0, %125 ], [ %171, %139 ]
%178 = phi double [ 0.000000e+00, %125 ], [ %170, %139 ]
%179 = phi double [ 0.000000e+00, %125 ], [ %169, %139 ]
%180 = icmp eq i64 %127, 0
br i1 %180, label %195, label %181
181: ; preds = %174, %181
%182 = phi i64 [ %192, %181 ], [ %177, %174 ]
%183 = phi double [ %191, %181 ], [ %178, %174 ]
%184 = phi double [ %190, %181 ], [ %179, %174 ]
%185 = phi i64 [ %193, %181 ], [ 0, %174 ]
%186 = getelementptr inbounds double, double* %8, i64 %182
%187 = load double, double* %186
%188 = getelementptr inbounds double, double* %7, i64 %182
%189 = load double, double* %188
%190 = tail call double @llvm.fmuladd.f64(double %187, double %189, double %184)
%191 = fadd double %183, %189
%192 = add nuw nsw i64 %182, 1
%193 = add i64 %185, 1
%194 = icmp eq i64 %193, %127
br i1 %194, label %195, label %181
195: ; preds = %174, %181, %73, %75, %124
%196 = phi double [ %67, %124 ], [ %67, %75 ], [ 0.000000e+00, %73 ], [ %67, %181 ], [ %67, %174 ]
%197 = phi i32 [ %66, %124 ], [ %66, %75 ], [ 0, %73 ], [ %66, %181 ], [ %66, %174 ]
%198 = phi i32 [ %64, %124 ], [ %64, %75 ], [ %29, %73 ], [ %64, %181 ], [ %64, %174 ]
%199 = phi double [ 0.000000e+00, %124 ], [ 0.000000e+00, %75 ], [ 0.000000e+00, %73 ], [ %175, %174 ], [ %190, %181 ]
%200 = phi double [ 0.000000e+00, %124 ], [ 0.000000e+00, %75 ], [ 0.000000e+00, %73 ], [ %176, %174 ], [ %191, %181 ]
%201 = icmp sgt i32 %197, %33
%202 = fcmp olt double %196, %32
%203 = select i1 %201, i1 true, i1 %202
br i1 %203, label %204, label %206
204: ; preds = %195
%205 = fdiv double %199, %200
store double %205, double* %3
br label %206
206: ; preds = %73, %195, %71, %204
%207 = phi i32 [ %198, %204 ], [ %64, %71 ], [ %198, %195 ], [ %29, %73 ]
%208 = phi i32 [ %197, %204 ], [ %33, %71 ], [ %33, %195 ], [ %33, %73 ]
%209 = phi double [ %196, %204 ], [ %32, %71 ], [ %32, %195 ], [ %32, %73 ]
%210 = add nuw nsw i64 %31, 1
%211 = sext i32 %207 to i64
%212 = icmp slt i64 %210, %211
br i1 %212, label %213, label %215
213: ; preds = %206
%214 = load double, double* %2
br label %28
215: ; preds = %206, %23
ret void
}