#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 */
typedef size_t int64_t;

/* Variables and functions */
float chi2_scale(float, float *, float *, size_t);
double fabs(float);
float calc_scale_from_bins(float rs, float *bin_r, float *weights,
                           int64_t num_bins) {
  if (!(rs > bin_r[0]) || !(rs < bin_r[num_bins]))
    rs = bin_r[(num_bins + 1) / 2];
  float initial_rs = rs;
  float chi2 = chi2_scale(rs, bin_r, weights, num_bins);
  float drs = rs / 10.0;
  int64_t iter = 0;
  float last_chi2 = 1e30;
  if (!(rs > 0))
    return initial_rs;
  while (fabs(chi2 - last_chi2) > 0.005 * chi2 && iter < 15) {
    last_chi2 = chi2;
    if (rs + drs > bin_r[num_bins])
      drs = 0.1 * (bin_r[num_bins] - rs);
    if (rs - drs < bin_r[0])
      drs = 0.1 * (rs - bin_r[0]);
    float chi2_right = chi2_scale(rs + drs, bin_r, weights, num_bins);
    float chi2_left = chi2_scale(rs - drs, bin_r, weights, num_bins);
    float dx = 0.5 * (chi2_right - chi2_left) / drs;
    float dx2 = (chi2_right + chi2_left - 2.0 * chi2) / (drs * drs);
    float move = 0;
    if (dx2 != 0)
      move = -dx / dx2;
    if (!move)
      return rs;
    if (rs + move > 4 * rs)
      move = 3.0 * rs;
    if (rs + move < 0.25 * rs)
      move = -0.75 * rs;
    float new_rs = rs + move;
    if (new_rs > bin_r[num_bins])
      new_rs = rs + 0.8 * (bin_r[num_bins] - rs);
    else if (new_rs < bin_r[0])
      new_rs = rs + 0.8 * (bin_r[0] - rs);
    float new_chi2 = chi2_scale(new_rs, bin_r, weights, num_bins);
    if (chi2_right < new_chi2) {
      new_chi2 = chi2_right;
      new_rs = rs + drs;
    }
    if (chi2_left < new_chi2) {
      new_chi2 = chi2_left;
      new_rs = rs - drs;
    }
    if (last_chi2 < new_chi2) {
      drs *= 0.5;
      continue;
    }
    drs = fabs(rs - new_rs) / 10.0;
    rs = new_rs;
    chi2 = new_chi2;
    iter++;
  }
  return rs;
}


x86 assembly (clang, -O0):
<calc_scale_from_bins>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
movss  %xmm0,-0x8(%rbp)
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
movss  -0x8(%rbp),%xmm0
mov    -0x10(%rbp),%rax
ucomiss (%rax),%xmm0
jbe    1166 <calc_scale_from_bins+0x46>
movss  -0x8(%rbp),%xmm1
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
ucomiss %xmm1,%xmm0
ja     1180 <calc_scale_from_bins+0x60>
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%rcx
add    $0x1,%rcx
shr    $0x1,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x8(%rbp)
movss  -0x8(%rbp),%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x8(%rbp),%xmm0
mov    -0x10(%rbp),%rdi
mov    -0x18(%rbp),%rsi
mov    -0x20(%rbp),%rdx
call   1040 <chi2_scale@plt>
movss  %xmm0,-0x28(%rbp)
movss  -0x8(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe52(%rip),%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x2c(%rbp)
movq   $0x0,-0x38(%rbp)
movss  0xe2d(%rip),%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x8(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
ja     11f8 <calc_scale_from_bins+0xd8>
movss  -0x24(%rbp),%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    1623 <calc_scale_from_bins+0x503>
jmp    11fd <calc_scale_from_bins+0xdd>
movss  -0x28(%rbp),%xmm0
subss  -0x3c(%rbp),%xmm0
call   1030 <fabs@plt>
movss  -0x28(%rbp),%xmm1
cvtss2sd %xmm1,%xmm2
movsd  0xdf3(%rip),%xmm1
mulsd  %xmm2,%xmm1
xor    %eax,%eax
ucomisd %xmm1,%xmm0
mov    %al,-0x59(%rbp)
jbe    123b <calc_scale_from_bins+0x11b>
cmpq   $0xf,-0x38(%rbp)
setb   %al
mov    %al,-0x59(%rbp)
mov    -0x59(%rbp),%al
test   $0x1,%al
jne    124b <calc_scale_from_bins+0x12b>
jmp    1619 <calc_scale_from_bins+0x4f9>
movss  -0x28(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x8(%rbp),%xmm0
addss  -0x2c(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%rcx
ucomiss (%rax,%rcx,4),%xmm0
jbe    129c <calc_scale_from_bins+0x17c>
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x8(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xd89(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  -0x8(%rbp),%xmm1
subss  -0x2c(%rbp),%xmm1
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm0
ucomiss %xmm1,%xmm0
jbe    12dd <calc_scale_from_bins+0x1bd>
movss  -0x8(%rbp),%xmm0
mov    -0x10(%rbp),%rax
subss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xd48(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  -0x8(%rbp),%xmm0
addss  -0x2c(%rbp),%xmm0
mov    -0x10(%rbp),%rdi
mov    -0x18(%rbp),%rsi
mov    -0x20(%rbp),%rdx
call   1040 <chi2_scale@plt>
movss  %xmm0,-0x40(%rbp)
movss  -0x8(%rbp),%xmm0
subss  -0x2c(%rbp),%xmm0
mov    -0x10(%rbp),%rdi
mov    -0x18(%rbp),%rsi
mov    -0x20(%rbp),%rdx
call   1040 <chi2_scale@plt>
movss  %xmm0,-0x44(%rbp)
movss  -0x40(%rbp),%xmm0
subss  -0x44(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xced(%rip),%xmm0
mulsd  %xmm1,%xmm0
movss  -0x2c(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x48(%rbp)
movss  -0x40(%rbp),%xmm0
addss  -0x44(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movss  -0x28(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
addsd  %xmm1,%xmm1
subsd  %xmm1,%xmm0
movss  -0x2c(%rbp),%xmm1
mulss  -0x2c(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x4c(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x50(%rbp)
movss  -0x4c(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jne    13ab <calc_scale_from_bins+0x28b>
jp     13ab <calc_scale_from_bins+0x28b>
jmp    13c7 <calc_scale_from_bins+0x2a7>
movss  -0x48(%rbp),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
divss  -0x4c(%rbp),%xmm0
movss  %xmm0,-0x50(%rbp)
movss  -0x50(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jne    13ed <calc_scale_from_bins+0x2cd>
jp     13ed <calc_scale_from_bins+0x2cd>
movss  -0x8(%rbp),%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    1623 <calc_scale_from_bins+0x503>
movss  -0x8(%rbp),%xmm0
addss  -0x50(%rbp),%xmm0
movss  0xc05(%rip),%xmm1
mulss  -0x8(%rbp),%xmm1
ucomiss %xmm1,%xmm0
jbe    142b <calc_scale_from_bins+0x30b>
movss  -0x8(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xc0a(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x50(%rbp)
movss  -0x8(%rbp),%xmm0
addss  -0x50(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movss  -0x8(%rbp),%xmm0
cvtss2sd %xmm0,%xmm2
movsd  0xbe6(%rip),%xmm0
mulsd  %xmm2,%xmm0
ucomisd %xmm1,%xmm0
jbe    1476 <calc_scale_from_bins+0x356>
movss  -0x8(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xbcf(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x50(%rbp)
movss  -0x8(%rbp),%xmm0
addss  -0x50(%rbp),%xmm0
movss  %xmm0,-0x54(%rbp)
movss  -0x54(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%rcx
ucomiss (%rax,%rcx,4),%xmm0
jbe    14d9 <calc_scale_from_bins+0x3b9>
movss  -0x8(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x8(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xb7d(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x54(%rbp)
jmp    1527 <calc_scale_from_bins+0x407>
movss  -0x54(%rbp),%xmm1
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm0
ucomiss %xmm1,%xmm0
jbe    1522 <calc_scale_from_bins+0x402>
movss  -0x8(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm0
subss  -0x8(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xb2f(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x54(%rbp)
jmp    1527 <calc_scale_from_bins+0x407>
movss  -0x54(%rbp),%xmm0
mov    -0x10(%rbp),%rdi
mov    -0x18(%rbp),%rsi
mov    -0x20(%rbp),%rdx
call   1040 <chi2_scale@plt>
movss  %xmm0,-0x58(%rbp)
movss  -0x40(%rbp),%xmm1
movss  -0x58(%rbp),%xmm0
ucomiss %xmm1,%xmm0
jbe    156e <calc_scale_from_bins+0x44e>
movss  -0x40(%rbp),%xmm0
movss  %xmm0,-0x58(%rbp)
movss  -0x8(%rbp),%xmm0
addss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x54(%rbp)
movss  -0x44(%rbp),%xmm1
movss  -0x58(%rbp),%xmm0
ucomiss %xmm1,%xmm0
jbe    159a <calc_scale_from_bins+0x47a>
movss  -0x44(%rbp),%xmm0
movss  %xmm0,-0x58(%rbp)
movss  -0x8(%rbp),%xmm0
subss  -0x2c(%rbp),%xmm0
movss  %xmm0,-0x54(%rbp)
movss  -0x3c(%rbp),%xmm1
movss  -0x58(%rbp),%xmm0
ucomiss %xmm1,%xmm0
jbe    15d0 <calc_scale_from_bins+0x4b0>
movss  -0x2c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xa62(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x2c(%rbp)
jmp    11fd <calc_scale_from_bins+0xdd>
movss  -0x8(%rbp),%xmm0
subss  -0x54(%rbp),%xmm0
call   1030 <fabs@plt>
movsd  0xa21(%rip),%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  -0x54(%rbp),%xmm0
movss  %xmm0,-0x8(%rbp)
movss  -0x58(%rbp),%xmm0
movss  %xmm0,-0x28(%rbp)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    11fd <calc_scale_from_bins+0xdd>
movss  -0x8(%rbp),%xmm0
movss  %xmm0,-0x4(%rbp)
movss  -0x4(%rbp),%xmm0
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<calc_scale_from_bins>:
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x78,%rsp
mov    %rdx,%r15
mov    %rsi,%r14
mov    %rdi,%rbx
ucomiss (%rdi),%xmm0
jbe    1144 <calc_scale_from_bins+0x24>
movss  (%rbx,%r15,4),%xmm1
ucomiss %xmm0,%xmm1
ja     1155 <calc_scale_from_bins+0x35>
lea    (%r15,%r15,1),%rax
add    $0x2,%rax
and    $0xfffffffffffffffc,%rax
movss  (%rbx,%rax,1),%xmm0
movss  %xmm0,0x8(%rsp)
mov    %rbx,%rdi
mov    %r14,%rsi
mov    %r15,%rdx
call   1040 <chi2_scale@plt>
movss  0x8(%rsp),%xmm4
movaps %xmm0,%xmm1
xorps  %xmm0,%xmm0
ucomiss %xmm0,%xmm4
jbe    14be <calc_scale_from_bins+0x39e>
movss  0xe7e(%rip),%xmm0
addss  %xmm1,%xmm0
movss  %xmm1,0xc(%rsp)
call   1030 <fabs@plt>
movss  0xc(%rsp),%xmm2
cvtss2sd %xmm2,%xmm3
movsd  0xe71(%rip),%xmm1
mulsd  %xmm3,%xmm1
ucomisd %xmm1,%xmm0
jbe    14b8 <calc_scale_from_bins+0x398>
movss  0x8(%rsp),%xmm4
movaps %xmm4,%xmm0
divss  0xe3a(%rip),%xmm0
xor    %r12d,%r12d
nopl   0x0(%rax)
movss  %xmm4,0x8(%rsp)
movsd  %xmm3,0x30(%rsp)
movss  %xmm2,0xc(%rsp)
movaps %xmm0,%xmm1
addss  %xmm4,%xmm1
movss  (%rbx,%r15,4),%xmm2
movaps %xmm2,%xmm3
subss  %xmm4,%xmm3
cvtss2sd %xmm3,%xmm3
movsd  0xe1e(%rip),%xmm5
mulsd  %xmm5,%xmm3
cvtsd2ss %xmm3,%xmm3
cmpltss %xmm1,%xmm2
movaps %xmm2,%xmm1
andnps %xmm0,%xmm1
andps  %xmm3,%xmm2
orps   %xmm1,%xmm2
movaps %xmm4,%xmm3
subss  %xmm2,%xmm3
movss  (%rbx),%xmm0
movaps %xmm4,%xmm1
subss  %xmm0,%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  %xmm5,%xmm1
cvtsd2ss %xmm1,%xmm1
cmpltss %xmm0,%xmm3
movaps %xmm3,%xmm0
andnps %xmm2,%xmm0
andps  %xmm1,%xmm3
orps   %xmm0,%xmm3
movaps %xmm3,0x10(%rsp)
addss  %xmm3,%xmm4
movaps %xmm4,0x60(%rsp)
movaps %xmm4,%xmm0
mov    %rbx,%rdi
mov    %r14,%rsi
mov    %r15,%rdx
call   1040 <chi2_scale@plt>
movaps %xmm0,0x20(%rsp)
movss  0x8(%rsp),%xmm0
subss  0x10(%rsp),%xmm0
movaps %xmm0,0x50(%rsp)
mov    %rbx,%rdi
mov    %r14,%rsi
mov    %r15,%rdx
call   1040 <chi2_scale@plt>
movaps 0x10(%rsp),%xmm3
movaps 0x20(%rsp),%xmm4
movaps %xmm4,%xmm1
addss  %xmm0,%xmm1
xorps  %xmm2,%xmm2
cvtss2sd %xmm1,%xmm2
movsd  0x30(%rsp),%xmm1
addsd  %xmm1,%xmm1
subsd  %xmm1,%xmm2
movaps %xmm3,%xmm1
mulss  %xmm3,%xmm1
cvtss2sd %xmm1,%xmm1
divsd  %xmm1,%xmm2
xorps  %xmm1,%xmm1
cvtsd2ss %xmm2,%xmm1
xorps  %xmm2,%xmm2
ucomiss %xmm2,%xmm1
jne    12d7 <calc_scale_from_bins+0x1b7>
jnp    14b8 <calc_scale_from_bins+0x398>
cvtss2sd %xmm3,%xmm7
movaps %xmm4,%xmm3
subss  %xmm0,%xmm3
cvtss2sd %xmm3,%xmm3
mulsd  0xd3a(%rip),%xmm3
divsd  %xmm7,%xmm3
xorps  %xmm5,%xmm5
cvtsd2ss %xmm3,%xmm5
divss  %xmm1,%xmm5
ucomiss %xmm2,%xmm5
jne    1308 <calc_scale_from_bins+0x1e8>
jnp    14b8 <calc_scale_from_bins+0x398>
movss  0x8(%rsp),%xmm4
movaps %xmm4,%xmm2
addss  %xmm5,%xmm2
movaps %xmm4,%xmm1
mulss  0xce8(%rip),%xmm1
movaps %xmm4,%xmm3
mulss  0xce1(%rip),%xmm3
cmpltss %xmm2,%xmm1
andps  %xmm1,%xmm3
andnps %xmm5,%xmm1
orps   %xmm3,%xmm1
movaps %xmm4,%xmm2
addss  %xmm1,%xmm2
cvtss2sd %xmm2,%xmm2
xorps  %xmm5,%xmm5
cvtss2sd %xmm4,%xmm5
movaps %xmm5,%xmm3
mulsd  0xcda(%rip),%xmm3
ucomisd %xmm2,%xmm3
movaps %xmm4,%xmm2
mulss  0xcab(%rip),%xmm2
ja     136a <calc_scale_from_bins+0x24a>
movaps %xmm1,%xmm2
addss  %xmm4,%xmm2
movss  (%rbx,%r15,4),%xmm1
ucomiss %xmm1,%xmm2
movaps %xmm0,0x40(%rsp)
movsd  %xmm7,0x10(%rsp)
ja     138d <calc_scale_from_bins+0x26d>
movss  (%rbx),%xmm1
ucomiss %xmm2,%xmm1
jbe    13a8 <calc_scale_from_bins+0x288>
subss  %xmm4,%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  0xc9b(%rip),%xmm1
addsd  %xmm5,%xmm1
xorps  %xmm2,%xmm2
cvtsd2ss %xmm1,%xmm2
movaps %xmm2,0x30(%rsp)
movaps %xmm2,%xmm0
mov    %rbx,%rdi
mov    %r14,%rsi
mov    %r15,%rdx
call   1040 <chi2_scale@plt>
movaps 0x20(%rsp),%xmm4
movaps %xmm4,%xmm1
minss  %xmm0,%xmm1
movaps 0x40(%rsp),%xmm3
movaps %xmm3,%xmm5
minss  %xmm1,%xmm5
movss  0xc(%rsp),%xmm2
ucomiss %xmm2,%xmm5
jbe    1400 <calc_scale_from_bins+0x2e0>
movsd  0x10(%rsp),%xmm0
mulsd  0xc59(%rip),%xmm0
movsd  %xmm0,0x10(%rsp)
movaps %xmm2,%xmm1
jmp    146c <calc_scale_from_bins+0x34c>
nopw   0x0(%rax,%rax,1)
cmpltss %xmm0,%xmm4
movaps 0x60(%rsp),%xmm0
andps  %xmm4,%xmm0
andnps 0x30(%rsp),%xmm4
orps   %xmm0,%xmm4
cmpltss %xmm1,%xmm3
movaps %xmm3,%xmm0
andnps %xmm4,%xmm0
andps  0x50(%rsp),%xmm3
orps   %xmm0,%xmm3
movaps %xmm3,0x40(%rsp)
movss  0x8(%rsp),%xmm0
subss  %xmm3,%xmm0
movss  %xmm5,0x20(%rsp)
call   1030 <fabs@plt>
movss  0x20(%rsp),%xmm1
movss  0xc(%rsp),%xmm2
divsd  0xbea(%rip),%xmm0
movsd  %xmm0,0x10(%rsp)
add    $0x1,%r12
movapd 0x40(%rsp),%xmm0
movss  %xmm0,0x8(%rsp)
movss  %xmm1,0x20(%rsp)
movaps %xmm1,%xmm0
subss  %xmm2,%xmm0
call   1030 <fabs@plt>
movss  0x20(%rsp),%xmm2
xorps  %xmm3,%xmm3
cvtss2sd %xmm2,%xmm3
movaps %xmm3,%xmm1
mulsd  0xb82(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    14b8 <calc_scale_from_bins+0x398>
movsd  0x10(%rsp),%xmm0
cvtsd2ss %xmm0,%xmm0
cmp    $0xf,%r12
movss  0x8(%rsp),%xmm4
jb     11d0 <calc_scale_from_bins+0xb0>
jmp    14be <calc_scale_from_bins+0x39e>
movss  0x8(%rsp),%xmm4
movaps %xmm4,%xmm0
add    $0x78,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define float @calc_scale_from_bins(float %0, float* %1, float* %2, i64 %3) {
%5 = alloca float
%6 = alloca float
%7 = alloca float*
%8 = alloca float*
%9 = alloca i64
%10 = alloca float
%11 = alloca float
%12 = alloca float
%13 = alloca i64
%14 = alloca float
%15 = alloca float
%16 = alloca float
%17 = alloca float
%18 = alloca float
%19 = alloca float
%20 = alloca float
%21 = alloca float
store float %0, float* %6
store float* %1, float** %7
store float* %2, float** %8
store i64 %3, i64* %9
%22 = load float, float* %6
%23 = load float*, float** %7
%24 = getelementptr inbounds float, float* %23, i64 0
%25 = load float, float* %24
%26 = fcmp ogt float %22, %25
br i1 %26, label %27, label %34
27: ; preds = %4
%28 = load float, float* %6
%29 = load float*, float** %7
%30 = load i64, i64* %9
%31 = getelementptr inbounds float, float* %29, i64 %30
%32 = load float, float* %31
%33 = fcmp olt float %28, %32
br i1 %33, label %41, label %34
34: ; preds = %27, %4
%35 = load float*, float** %7
%36 = load i64, i64* %9
%37 = add i64 %36, 1
%38 = udiv i64 %37, 2
%39 = getelementptr inbounds float, float* %35, i64 %38
%40 = load float, float* %39
store float %40, float* %6
br label %41
41: ; preds = %34, %27
%42 = load float, float* %6
store float %42, float* %10
%43 = load float, float* %6
%44 = load float*, float** %7
%45 = load float*, float** %8
%46 = load i64, i64* %9
%47 = call float @chi2_scale(float %43, float* %44, float* %45, i64 %46)
store float %47, float* %11
%48 = load float, float* %6
%49 = fpext float %48 to double
%50 = fdiv double %49, 1.000000e+01
%51 = fptrunc double %50 to float
store float %51, float* %12
store i64 0, i64* %13
store float 0x46293E5940000000, float* %14
%52 = load float, float* %6
%53 = fcmp ogt float %52, 0.000000e+00
br i1 %53, label %56, label %54
54: ; preds = %41
%55 = load float, float* %10
store float %55, float* %5
br label %268
56: ; preds = %41
br label %57
57: ; preds = %255, %250, %56
%58 = load float, float* %11
%59 = load float, float* %14
%60 = fsub float %58, %59
%61 = call double @fabs(float %60)
%62 = load float, float* %11
%63 = fpext float %62 to double
%64 = fmul double 5.000000e-03, %63
%65 = fcmp ogt double %61, %64
br i1 %65, label %66, label %69
66: ; preds = %57
%67 = load i64, i64* %13
%68 = icmp ult i64 %67, 15
br label %69
69: ; preds = %66, %57
%70 = phi i1 [ false, %57 ], [ %68, %66 ]
br i1 %70, label %71, label %266
71: ; preds = %69
%72 = load float, float* %11
store float %72, float* %14
%73 = load float, float* %6
%74 = load float, float* %12
%75 = fadd float %73, %74
%76 = load float*, float** %7
%77 = load i64, i64* %9
%78 = getelementptr inbounds float, float* %76, i64 %77
%79 = load float, float* %78
%80 = fcmp ogt float %75, %79
br i1 %80, label %81, label %91
81: ; preds = %71
%82 = load float*, float** %7
%83 = load i64, i64* %9
%84 = getelementptr inbounds float, float* %82, i64 %83
%85 = load float, float* %84
%86 = load float, float* %6
%87 = fsub float %85, %86
%88 = fpext float %87 to double
%89 = fmul double 1.000000e-01, %88
%90 = fptrunc double %89 to float
store float %90, float* %12
br label %91
91: ; preds = %81, %71
%92 = load float, float* %6
%93 = load float, float* %12
%94 = fsub float %92, %93
%95 = load float*, float** %7
%96 = getelementptr inbounds float, float* %95, i64 0
%97 = load float, float* %96
%98 = fcmp olt float %94, %97
br i1 %98, label %99, label %108
99: ; preds = %91
%100 = load float, float* %6
%101 = load float*, float** %7
%102 = getelementptr inbounds float, float* %101, i64 0
%103 = load float, float* %102
%104 = fsub float %100, %103
%105 = fpext float %104 to double
%106 = fmul double 1.000000e-01, %105
%107 = fptrunc double %106 to float
store float %107, float* %12
br label %108
108: ; preds = %99, %91
%109 = load float, float* %6
%110 = load float, float* %12
%111 = fadd float %109, %110
%112 = load float*, float** %7
%113 = load float*, float** %8
%114 = load i64, i64* %9
%115 = call float @chi2_scale(float %111, float* %112, float* %113, i64 %114)
store float %115, float* %15
%116 = load float, float* %6
%117 = load float, float* %12
%118 = fsub float %116, %117
%119 = load float*, float** %7
%120 = load float*, float** %8
%121 = load i64, i64* %9
%122 = call float @chi2_scale(float %118, float* %119, float* %120, i64 %121)
store float %122, float* %16
%123 = load float, float* %15
%124 = load float, float* %16
%125 = fsub float %123, %124
%126 = fpext float %125 to double
%127 = fmul double 5.000000e-01, %126
%128 = load float, float* %12
%129 = fpext float %128 to double
%130 = fdiv double %127, %129
%131 = fptrunc double %130 to float
store float %131, float* %17
%132 = load float, float* %15
%133 = load float, float* %16
%134 = fadd float %132, %133
%135 = fpext float %134 to double
%136 = load float, float* %11
%137 = fpext float %136 to double
%138 = call double @llvm.fmuladd.f64(double -2.000000e+00, double %137, double %135)
%139 = load float, float* %12
%140 = load float, float* %12
%141 = fmul float %139, %140
%142 = fpext float %141 to double
%143 = fdiv double %138, %142
%144 = fptrunc double %143 to float
store float %144, float* %18
store float 0.000000e+00, float* %19
%145 = load float, float* %18
%146 = fcmp une float %145, 0.000000e+00
br i1 %146, label %147, label %152
147: ; preds = %108
%148 = load float, float* %17
%149 = fneg float %148
%150 = load float, float* %18
%151 = fdiv float %149, %150
store float %151, float* %19
br label %152
152: ; preds = %147, %108
%153 = load float, float* %19
%154 = fcmp une float %153, 0.000000e+00
br i1 %154, label %157, label %155
155: ; preds = %152
%156 = load float, float* %6
store float %156, float* %5
br label %268
157: ; preds = %152
%158 = load float, float* %6
%159 = load float, float* %19
%160 = fadd float %158, %159
%161 = load float, float* %6
%162 = fmul float 4.000000e+00, %161
%163 = fcmp ogt float %160, %162
br i1 %163, label %164, label %169
164: ; preds = %157
%165 = load float, float* %6
%166 = fpext float %165 to double
%167 = fmul double 3.000000e+00, %166
%168 = fptrunc double %167 to float
store float %168, float* %19
br label %169
169: ; preds = %164, %157
%170 = load float, float* %6
%171 = load float, float* %19
%172 = fadd float %170, %171
%173 = fpext float %172 to double
%174 = load float, float* %6
%175 = fpext float %174 to double
%176 = fmul double 2.500000e-01, %175
%177 = fcmp olt double %173, %176
br i1 %177, label %178, label %183
178: ; preds = %169
%179 = load float, float* %6
%180 = fpext float %179 to double
%181 = fmul double -7.500000e-01, %180
%182 = fptrunc double %181 to float
store float %182, float* %19
br label %183
183: ; preds = %178, %169
%184 = load float, float* %6
%185 = load float, float* %19
%186 = fadd float %184, %185
store float %186, float* %20
%187 = load float, float* %20
%188 = load float*, float** %7
%189 = load i64, i64* %9
%190 = getelementptr inbounds float, float* %188, i64 %189
%191 = load float, float* %190
%192 = fcmp ogt float %187, %191
br i1 %192, label %193, label %205
193: ; preds = %183
%194 = load float, float* %6
%195 = fpext float %194 to double
%196 = load float*, float** %7
%197 = load i64, i64* %9
%198 = getelementptr inbounds float, float* %196, i64 %197
%199 = load float, float* %198
%200 = load float, float* %6
%201 = fsub float %199, %200
%202 = fpext float %201 to double
%203 = call double @llvm.fmuladd.f64(double 8.000000e-01, double %202, double %195)
%204 = fptrunc double %203 to float
store float %204, float* %20
br label %223
205: ; preds = %183
%206 = load float, float* %20
%207 = load float*, float** %7
%208 = getelementptr inbounds float, float* %207, i64 0
%209 = load float, float* %208
%210 = fcmp olt float %206, %209
br i1 %210, label %211, label %222
211: ; preds = %205
%212 = load float, float* %6
%213 = fpext float %212 to double
%214 = load float*, float** %7
%215 = getelementptr inbounds float, float* %214, i64 0
%216 = load float, float* %215
%217 = load float, float* %6
%218 = fsub float %216, %217
%219 = fpext float %218 to double
%220 = call double @llvm.fmuladd.f64(double 8.000000e-01, double %219, double %213)
%221 = fptrunc double %220 to float
store float %221, float* %20
br label %222
222: ; preds = %211, %205
br label %223
223: ; preds = %222, %193
%224 = load float, float* %20
%225 = load float*, float** %7
%226 = load float*, float** %8
%227 = load i64, i64* %9
%228 = call float @chi2_scale(float %224, float* %225, float* %226, i64 %227)
store float %228, float* %21
%229 = load float, float* %15
%230 = load float, float* %21
%231 = fcmp olt float %229, %230
br i1 %231, label %232, label %237
232: ; preds = %223
%233 = load float, float* %15
store float %233, float* %21
%234 = load float, float* %6
%235 = load float, float* %12
%236 = fadd float %234, %235
store float %236, float* %20
br label %237
237: ; preds = %232, %223
%238 = load float, float* %16
%239 = load float, float* %21
%240 = fcmp olt float %238, %239
br i1 %240, label %241, label %246
241: ; preds = %237
%242 = load float, float* %16
store float %242, float* %21
%243 = load float, float* %6
%244 = load float, float* %12
%245 = fsub float %243, %244
store float %245, float* %20
br label %246
246: ; preds = %241, %237
%247 = load float, float* %14
%248 = load float, float* %21
%249 = fcmp olt float %247, %248
br i1 %249, label %250, label %255
250: ; preds = %246
%251 = load float, float* %12
%252 = fpext float %251 to double
%253 = fmul double %252, 5.000000e-01
%254 = fptrunc double %253 to float
store float %254, float* %12
br label %57
255: ; preds = %246
%256 = load float, float* %6
%257 = load float, float* %20
%258 = fsub float %256, %257
%259 = call double @fabs(float %258)
%260 = fdiv double %259, 1.000000e+01
%261 = fptrunc double %260 to float
store float %261, float* %12
%262 = load float, float* %20
store float %262, float* %6
%263 = load float, float* %21
store float %263, float* %11
%264 = load i64, i64* %13
%265 = add i64 %264, 1
store i64 %265, i64* %13
br label %57
266: ; preds = %69
%267 = load float, float* %6
store float %267, float* %5
br label %268
268: ; preds = %266, %155, %54
%269 = load float, float* %5
ret float %269
}

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