int getnbestdiagsum(row, col, score, ndots, mincol, maxcol, minrow, maxrow,
                    level, region, minscore, bestscore)
int *col;
int *row;
double *score;
int ndots;
int mincol, maxcol, minrow, maxrow;
int level;
int region;
double minscore;
double *bestscore;
{
  int i;
  double *count;
  int *index;
  int d, d1, n;
  int maxdiag, maxindex;
  int bestdiag;
  int dmin, dmax;
  maxdiag = ((maxrow - mincol) > (maxcol - mincol)) ? (maxrow - mincol)
                                                    : (maxcol - minrow);
  if (maxdiag < 0)
    maxdiag = -maxdiag;
  maxindex = maxdiag;
  count = (double *)mymalloc((long)(maxindex + 1) * sizeof(double));
  index = (int *)mymalloc((long)(maxindex + 1) * sizeof(int));
  for (i = 0; i < maxindex; i++)
    count[i] = 0;
  for (i = 0; i < maxindex; i++)
    index[i] = i;
  for (i = 0; i < ndots; i++) {
    if (row[i] > col[i]) {
      continue;
    }
    if (row[i] >= minrow && row[i] <= maxrow && col[i] >= mincol &&
        col[i] <= maxcol) {
      d = col[i] - row[i];
      dmin = ((d - region) < 0) ? 0 : (d - region);
      dmax = ((d + region) >= maxindex) ? maxindex : (d + region);
      for (d1 = dmin; d1 <= dmax; d1++)
        if (score[i] > minscore)
          count[d1] += score[i];
    }
  }
  QuickSortIndicesDouble(0, maxindex - 1, index, count);
  bestdiag = index[maxindex - level - 1];
  d = bestdiag;
  n = 1;
  i = bestdiag + 1;
  while ((i < maxindex) && (count[i] == count[bestdiag])) {
    d += i;
    n++;
    i++;
  }
  i = bestdiag - 1;
  while ((i >= 0) && (count[i] == count[bestdiag])) {
    d += i;
    n++;
    i--;
  }
  bestdiag = rint(d / n);
  *bestscore = count[bestdiag];
  free(count);
  free(index);
  return bestdiag;
}


x86 assembly (clang, -O0):
<getnbestdiagsum>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    0x30(%rbp),%rax
mov    0x28(%rbp),%eax
mov    0x20(%rbp),%eax
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
mov    %r9d,-0x24(%rbp)
movsd  %xmm0,-0x30(%rbp)
mov    0x18(%rbp),%eax
sub    -0x20(%rbp),%eax
mov    -0x24(%rbp),%ecx
sub    -0x20(%rbp),%ecx
cmp    %ecx,%eax
jle    1199 <getnbestdiagsum+0x59>
mov    0x18(%rbp),%eax
sub    -0x20(%rbp),%eax
mov    %eax,-0x6c(%rbp)
jmp    11a2 <getnbestdiagsum+0x62>
mov    -0x24(%rbp),%eax
sub    0x10(%rbp),%eax
mov    %eax,-0x6c(%rbp)
mov    -0x6c(%rbp),%eax
mov    %eax,-0x58(%rbp)
cmpl   $0x0,-0x58(%rbp)
jge    11ba <getnbestdiagsum+0x7a>
xor    %eax,%eax
sub    -0x58(%rbp),%eax
mov    %eax,-0x58(%rbp)
mov    -0x58(%rbp),%eax
mov    %eax,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
add    $0x1,%eax
movslq %eax,%rdi
shl    $0x3,%rdi
mov    $0x0,%al
call   1060 <mymalloc@plt>
cltq
mov    %rax,-0x40(%rbp)
mov    -0x5c(%rbp),%eax
add    $0x1,%eax
movslq %eax,%rdi
shl    $0x2,%rdi
mov    $0x0,%al
call   1060 <mymalloc@plt>
cltq
mov    %rax,-0x48(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x5c(%rbp),%eax
jge    1225 <getnbestdiagsum+0xe5>
mov    -0x40(%rbp),%rax
movslq -0x34(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    11fb <getnbestdiagsum+0xbb>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x5c(%rbp),%eax
jge    1254 <getnbestdiagsum+0x114>
mov    -0x34(%rbp),%edx
mov    -0x48(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    122c <getnbestdiagsum+0xec>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    13b6 <getnbestdiagsum+0x276>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    -0x10(%rbp),%rcx
movslq -0x34(%rbp),%rdx
cmp    (%rcx,%rdx,4),%eax
jle    1288 <getnbestdiagsum+0x148>
jmp    13a8 <getnbestdiagsum+0x268>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
cmp    0x10(%rbp),%eax
jl     13a3 <getnbestdiagsum+0x263>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
cmp    0x18(%rbp),%eax
jg     13a3 <getnbestdiagsum+0x263>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
cmp    -0x20(%rbp),%eax
jl     13a3 <getnbestdiagsum+0x263>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
cmp    -0x24(%rbp),%eax
jg     13a3 <getnbestdiagsum+0x263>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    -0x8(%rbp),%rcx
movslq -0x34(%rbp),%rdx
sub    (%rcx,%rdx,4),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
sub    0x28(%rbp),%eax
cmp    $0x0,%eax
jge    130a <getnbestdiagsum+0x1ca>
xor    %eax,%eax
mov    %eax,-0x70(%rbp)
jmp    1313 <getnbestdiagsum+0x1d3>
mov    -0x4c(%rbp),%eax
sub    0x28(%rbp),%eax
mov    %eax,-0x70(%rbp)
mov    -0x70(%rbp),%eax
mov    %eax,-0x64(%rbp)
mov    -0x4c(%rbp),%eax
add    0x28(%rbp),%eax
cmp    -0x5c(%rbp),%eax
jl     1333 <getnbestdiagsum+0x1f3>
mov    -0x5c(%rbp),%eax
mov    %eax,-0x74(%rbp)
jmp    133c <getnbestdiagsum+0x1fc>
mov    -0x4c(%rbp),%eax
add    0x28(%rbp),%eax
mov    %eax,-0x74(%rbp)
mov    -0x74(%rbp),%eax
mov    %eax,-0x68(%rbp)
mov    -0x64(%rbp),%eax
mov    %eax,-0x50(%rbp)
mov    -0x50(%rbp),%eax
cmp    -0x68(%rbp),%eax
jg     139e <getnbestdiagsum+0x25e>
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x30(%rbp),%xmm0
jbe    138b <getnbestdiagsum+0x24b>
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x40(%rbp),%rax
movslq -0x50(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1390 <getnbestdiagsum+0x250>
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    1348 <getnbestdiagsum+0x208>
jmp    13a3 <getnbestdiagsum+0x263>
jmp    13a8 <getnbestdiagsum+0x268>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    125b <getnbestdiagsum+0x11b>
mov    -0x5c(%rbp),%esi
sub    $0x1,%esi
mov    -0x48(%rbp),%rdx
mov    -0x40(%rbp),%rcx
xor    %edi,%edi
mov    $0x0,%al
call   1040 <QuickSortIndicesDouble@plt>
mov    -0x48(%rbp),%rax
mov    -0x5c(%rbp),%ecx
sub    0x20(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x60(%rbp)
mov    -0x60(%rbp),%eax
mov    %eax,-0x4c(%rbp)
movl   $0x1,-0x54(%rbp)
mov    -0x60(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%ecx
xor    %eax,%eax
cmp    -0x5c(%rbp),%ecx
mov    %al,-0x75(%rbp)
jge    142f <getnbestdiagsum+0x2ef>
mov    -0x40(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x40(%rbp),%rax
movslq -0x60(%rbp),%rcx
ucomisd (%rax,%rcx,8),%xmm0
sete   %al
setnp  %cl
and    %cl,%al
mov    %al,-0x75(%rbp)
mov    -0x75(%rbp),%al
test   $0x1,%al
jne    143f <getnbestdiagsum+0x2ff>
jmp    145f <getnbestdiagsum+0x31f>
mov    -0x34(%rbp),%eax
add    -0x4c(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    13f9 <getnbestdiagsum+0x2b9>
mov    -0x60(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x34(%rbp)
xor    %eax,%eax
cmpl   $0x0,-0x34(%rbp)
mov    %al,-0x76(%rbp)
jl     149c <getnbestdiagsum+0x35c>
mov    -0x40(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x40(%rbp),%rax
movslq -0x60(%rbp),%rcx
ucomisd (%rax,%rcx,8),%xmm0
sete   %al
setnp  %cl
and    %cl,%al
mov    %al,-0x76(%rbp)
mov    -0x76(%rbp),%al
test   $0x1,%al
jne    14ac <getnbestdiagsum+0x36c>
jmp    14cc <getnbestdiagsum+0x38c>
mov    -0x34(%rbp),%eax
add    -0x4c(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
mov    -0x34(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x34(%rbp)
jmp    1468 <getnbestdiagsum+0x328>
mov    -0x4c(%rbp),%eax
cltd
idivl  -0x54(%rbp)
cvtsi2sd %eax,%xmm0
call   1050 <rint@plt>
cvttsd2si %xmm0,%eax
mov    %eax,-0x60(%rbp)
mov    -0x40(%rbp),%rax
movslq -0x60(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x30(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x40(%rbp),%rdi
call   1030 <free@plt>
mov    -0x48(%rbp),%rdi
call   1030 <free@plt>
mov    -0x60(%rbp),%eax
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<getnbestdiagsum>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
movq   %xmm0,0x18(%rsp)
mov    %ecx,0x4(%rsp)
mov    %rdx,%rbp
mov    %rsi,0x10(%rsp)
mov    %rdi,0x20(%rsp)
mov    0x68(%rsp),%ecx
mov    %ecx,%eax
mov    %r8d,0x8(%rsp)
sub    %r8d,%eax
mov    %r9d,%r15d
sub    0x60(%rsp),%r15d
mov    %r9d,0xc(%rsp)
cmp    %r9d,%ecx
cmovg  %eax,%r15d
mov    %r15d,%r13d
neg    %r13d
cmovs  %r15d,%r13d
lea    0x8(,%r13,8),%rdi
xor    %eax,%eax
call   1060 <mymalloc@plt>
movslq %eax,%rbx
lea    0x4(,%r13,4),%rdi
xor    %eax,%eax
call   1060 <mymalloc@plt>
movslq %eax,%r14
test   %r15d,%r15d
je     135d <getnbestdiagsum+0x21d>
cmp    $0x2,%r13d
mov    $0x1,%r15d
cmovae %r13d,%r15d
lea    0x0(,%r15,8),%rdx
xor    %r12d,%r12d
mov    %rbx,%rdi
xor    %esi,%esi
call   1040 <memset@plt>
cmp    $0x8,%r15d
jb     1350 <getnbestdiagsum+0x210>
mov    %r15d,%r12d
and    $0x7ffffff8,%r12d
lea    -0x8(%r12),%rdx
mov    %rdx,%rcx
shr    $0x3,%rcx
add    $0x1,%rcx
mov    %ecx,%eax
and    $0x3,%eax
cmp    $0x18,%rdx
jae    121f <getnbestdiagsum+0xdf>
movdqa 0xde8(%rip),%xmm0
xor    %edx,%edx
jmp    12f0 <getnbestdiagsum+0x1b0>
and    $0xfffffffffffffffc,%rcx
movdqa 0xdd5(%rip),%xmm0
xor    %edx,%edx
movdqa 0xdda(%rip),%xmm8
movdqa 0xde1(%rip),%xmm9
movdqa 0xde8(%rip),%xmm10
movdqa 0xdf0(%rip),%xmm4
movdqa 0xdf8(%rip),%xmm5
movdqa 0xe00(%rip),%xmm6
movdqa 0xe08(%rip),%xmm7
movdqa 0xe10(%rip),%xmm1
movdqa %xmm0,%xmm2
paddd  %xmm8,%xmm2
movdqu %xmm0,(%r14,%rdx,4)
movdqu %xmm2,0x10(%r14,%rdx,4)
movdqa %xmm0,%xmm2
paddd  %xmm9,%xmm2
movdqa %xmm0,%xmm3
paddd  %xmm10,%xmm3
movdqu %xmm2,0x20(%r14,%rdx,4)
movdqu %xmm3,0x30(%r14,%rdx,4)
movdqa %xmm0,%xmm2
paddd  %xmm4,%xmm2
movdqa %xmm0,%xmm3
paddd  %xmm5,%xmm3
movdqu %xmm2,0x40(%r14,%rdx,4)
movdqu %xmm3,0x50(%r14,%rdx,4)
movdqa %xmm0,%xmm2
paddd  %xmm6,%xmm2
movdqa %xmm0,%xmm3
paddd  %xmm7,%xmm3
movdqu %xmm2,0x60(%r14,%rdx,4)
movdqu %xmm3,0x70(%r14,%rdx,4)
add    $0x20,%rdx
paddd  %xmm1,%xmm0
add    $0xfffffffffffffffc,%rcx
jne    1270 <getnbestdiagsum+0x130>
test   %rax,%rax
je     1340 <getnbestdiagsum+0x200>
lea    (%r14,%rdx,4),%rcx
add    $0x10,%rcx
shl    $0x5,%rax
xor    %edx,%edx
movdqa 0xd05(%rip),%xmm1
movdqa 0xd0d(%rip),%xmm2
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movdqa %xmm0,%xmm3
paddd  %xmm1,%xmm3
movdqu %xmm0,-0x10(%rcx,%rdx,1)
movdqu %xmm3,(%rcx,%rdx,1)
paddd  %xmm2,%xmm0
add    $0x20,%rdx
cmp    %rdx,%rax
jne    1320 <getnbestdiagsum+0x1e0>
cmp    %r15,%r12
je     135d <getnbestdiagsum+0x21d>
cs nopw 0x0(%rax,%rax,1)
nop
mov    %r12d,(%r14,%r12,4)
add    $0x1,%r12
cmp    %r12,%r15
jne    1350 <getnbestdiagsum+0x210>
mov    0x70(%rsp),%r15d
mov    0x4(%rsp),%eax
test   %eax,%eax
movsd  0x18(%rsp),%xmm1
mov    0xc(%rsp),%r11d
mov    0x8(%rsp),%r12d
mov    0x10(%rsp),%rcx
jle    146d <getnbestdiagsum+0x32d>
mov    0x78(%rsp),%r9d
mov    %eax,%r10d
xor    %r8d,%r8d
xor    %esi,%esi
jmp    13ad <getnbestdiagsum+0x26d>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%rsi
cmp    %r10,%rsi
je     146d <getnbestdiagsum+0x32d>
mov    0x20(%rsp),%rax
mov    (%rax,%rsi,4),%eax
cmp    0x68(%rsp),%eax
jg     13a0 <getnbestdiagsum+0x260>
cmp    0x60(%rsp),%eax
jl     13a0 <getnbestdiagsum+0x260>
mov    (%rcx,%rsi,4),%edx
mov    %edx,%edi
sub    %eax,%edi
jl     13a0 <getnbestdiagsum+0x260>
cmp    %r12d,%edx
jl     13a0 <getnbestdiagsum+0x260>
cmp    %r11d,%edx
jg     13a0 <getnbestdiagsum+0x260>
mov    %edi,%edx
sub    %r9d,%edx
cmovs  %r8d,%edx
add    %r9d,%edi
cmp    %r13d,%edi
cmovge %r13d,%edi
cmp    %edi,%edx
jg     13a0 <getnbestdiagsum+0x260>
mov    %edx,%eax
add    $0x1,%edi
mov    %edi,%ecx
sub    %edx,%ecx
mov    %rax,%rdx
test   $0x1,%cl
je     1416 <getnbestdiagsum+0x2d6>
movsd  0x0(%rbp,%rsi,8),%xmm0
ucomisd %xmm1,%xmm0
jbe    1412 <getnbestdiagsum+0x2d2>
addsd  (%rbx,%rax,8),%xmm0
movsd  %xmm0,(%rbx,%rax,8)
lea    0x1(%rax),%rdx
not    %rax
add    %rdi,%rax
mov    0x10(%rsp),%rcx
jne    143d <getnbestdiagsum+0x2fd>
jmp    13a0 <getnbestdiagsum+0x260>
nopl   0x0(%rax,%rax,1)
add    $0x2,%rdx
cmp    %rdx,%rdi
je     13a0 <getnbestdiagsum+0x260>
movsd  0x0(%rbp,%rsi,8),%xmm0
ucomisd %xmm1,%xmm0
jbe    1453 <getnbestdiagsum+0x313>
addsd  (%rbx,%rdx,8),%xmm0
movsd  %xmm0,(%rbx,%rdx,8)
movsd  0x0(%rbp,%rsi,8),%xmm0
ucomisd %xmm1,%xmm0
jbe    1430 <getnbestdiagsum+0x2f0>
addsd  0x8(%rbx,%rdx,8),%xmm0
movsd  %xmm0,0x8(%rbx,%rdx,8)
jmp    1430 <getnbestdiagsum+0x2f0>
lea    -0x1(%r13),%esi
xor    %edi,%edi
mov    %r14,%rdx
mov    %rbx,%rcx
xor    %eax,%eax
call   1050 <QuickSortIndicesDouble@plt>
not    %r15d
add    %r13d,%r15d
movslq %r15d,%rax
mov    (%r14,%rax,4),%r8d
lea    0x1(%r8),%edx
mov    $0x1,%ecx
mov    %r8d,%eax
cmp    %r13d,%edx
jge    14e3 <getnbestdiagsum+0x3a3>
movslq %r8d,%rax
movsd  (%rbx,%rax,8),%xmm0
movslq %edx,%rax
mov    %r13d,%ecx
sub    %r8d,%ecx
lea    (%rbx,%rax,8),%rsi
mov    %edx,%edi
neg    %r13d
add    %rdi,%r13
xor    %ebp,%ebp
mov    %r8d,%eax
movsd  (%rsi,%rbp,8),%xmm1
ucomisd %xmm0,%xmm1
jne    14de <getnbestdiagsum+0x39e>
jp     14de <getnbestdiagsum+0x39e>
add    %edi,%eax
add    %ebp,%eax
add    $0x1,%rbp
mov    %r13d,%edx
add    %ebp,%edx
jne    14c0 <getnbestdiagsum+0x380>
jmp    14e3 <getnbestdiagsum+0x3a3>
add    $0x1,%ebp
mov    %ebp,%ecx
mov    0x80(%rsp),%rsi
test   %r8d,%r8d
jle    152b <getnbestdiagsum+0x3eb>
mov    %r8,%rdx
movsd  (%rbx,%r8,8),%xmm0
add    %ecx,%r8d
nopl   0x0(%rax)
lea    -0x1(%rdx),%rdi
mov    %edi,%ebp
movsd  (%rbx,%rbp,8),%xmm1
ucomisd %xmm0,%xmm1
jne    152b <getnbestdiagsum+0x3eb>
jp     152b <getnbestdiagsum+0x3eb>
add    %edx,%eax
add    $0xffffffff,%eax
add    $0x1,%ecx
lea    0x1(%rdi),%rbp
mov    %rdi,%rdx
cmp    $0x1,%rbp
ja     1500 <getnbestdiagsum+0x3c0>
mov    %r8d,%ecx
cltd
idiv   %ecx
movslq %eax,%rbp
movsd  (%rbx,%rbp,8),%xmm0
movsd  %xmm0,(%rsi)
mov    %rbx,%rdi
call   1030 <free@plt>
mov    %r14,%rdi
call   1030 <free@plt>
mov    %ebp,%eax
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define i32 @getnbestdiagsum(i32* %0, i32* %1, double* %2, i32 %3, i32 %4, i32 %5, i32 %6, i32 %7, i32 %8, i32 %9, double %10, double* %11) {
%13 = alloca i32*
%14 = alloca i32*
%15 = alloca double*
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
%20 = alloca i32
%21 = alloca i32
%22 = alloca i32
%23 = alloca double
%24 = alloca double*
%25 = alloca i32
%26 = alloca double*
%27 = alloca i32*
%28 = alloca i32
%29 = alloca i32
%30 = alloca i32
%31 = alloca i32
%32 = alloca i32
%33 = alloca i32
%34 = alloca i32
%35 = alloca i32
store i32* %0, i32** %13
store i32* %1, i32** %14
store double* %2, double** %15
store i32 %3, i32* %16
store i32 %4, i32* %17
store i32 %5, i32* %18
store i32 %6, i32* %19
store i32 %7, i32* %20
store i32 %8, i32* %21
store i32 %9, i32* %22
store double %10, double* %23
store double* %11, double** %24
%36 = load i32, i32* %20
%37 = load i32, i32* %17
%38 = sub nsw i32 %36, %37
%39 = load i32, i32* %18
%40 = load i32, i32* %17
%41 = sub nsw i32 %39, %40
%42 = icmp sgt i32 %38, %41
br i1 %42, label %43, label %47
43: ; preds = %12
%44 = load i32, i32* %20
%45 = load i32, i32* %17
%46 = sub nsw i32 %44, %45
br label %51
47: ; preds = %12
%48 = load i32, i32* %18
%49 = load i32, i32* %19
%50 = sub nsw i32 %48, %49
br label %51
51: ; preds = %47, %43
%52 = phi i32 [ %46, %43 ], [ %50, %47 ]
store i32 %52, i32* %31
%53 = load i32, i32* %31
%54 = icmp slt i32 %53, 0
br i1 %54, label %55, label %58
55: ; preds = %51
%56 = load i32, i32* %31
%57 = sub nsw i32 0, %56
store i32 %57, i32* %31
br label %58
58: ; preds = %55, %51
%59 = load i32, i32* %31
store i32 %59, i32* %32
%60 = load i32, i32* %32
%61 = add nsw i32 %60, 1
%62 = sext i32 %61 to i64
%63 = mul i64 %62, 8
%64 = call i32 (i64, ...) bitcast (i32 (...)* @mymalloc to i32 (i64, ...)*)(i64 %63)
%65 = sext i32 %64 to i64
%66 = inttoptr i64 %65 to double*
store double* %66, double** %26
%67 = load i32, i32* %32
%68 = add nsw i32 %67, 1
%69 = sext i32 %68 to i64
%70 = mul i64 %69, 4
%71 = call i32 (i64, ...) bitcast (i32 (...)* @mymalloc to i32 (i64, ...)*)(i64 %70)
%72 = sext i32 %71 to i64
%73 = inttoptr i64 %72 to i32*
store i32* %73, i32** %27
store i32 0, i32* %25
br label %74
74: ; preds = %83, %58
%75 = load i32, i32* %25
%76 = load i32, i32* %32
%77 = icmp slt i32 %75, %76
br i1 %77, label %78, label %86
78: ; preds = %74
%79 = load double*, double** %26
%80 = load i32, i32* %25
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds double, double* %79, i64 %81
store double 0.000000e+00, double* %82
br label %83
83: ; preds = %78
%84 = load i32, i32* %25
%85 = add nsw i32 %84, 1
store i32 %85, i32* %25
br label %74
86: ; preds = %74
store i32 0, i32* %25
br label %87
87: ; preds = %97, %86
%88 = load i32, i32* %25
%89 = load i32, i32* %32
%90 = icmp slt i32 %88, %89
br i1 %90, label %91, label %100
91: ; preds = %87
%92 = load i32, i32* %25
%93 = load i32*, i32** %27
%94 = load i32, i32* %25
%95 = sext i32 %94 to i64
%96 = getelementptr inbounds i32, i32* %93, i64 %95
store i32 %92, i32* %96
br label %97
97: ; preds = %91
%98 = load i32, i32* %25
%99 = add nsw i32 %98, 1
store i32 %99, i32* %25
br label %87
100: ; preds = %87
store i32 0, i32* %25
br label %101
101: ; preds = %217, %100
%102 = load i32, i32* %25
%103 = load i32, i32* %16
%104 = icmp slt i32 %102, %103
br i1 %104, label %105, label %220
105: ; preds = %101
%106 = load i32*, i32** %13
%107 = load i32, i32* %25
%108 = sext i32 %107 to i64
%109 = getelementptr inbounds i32, i32* %106, i64 %108
%110 = load i32, i32* %109
%111 = load i32*, i32** %14
%112 = load i32, i32* %25
%113 = sext i32 %112 to i64
%114 = getelementptr inbounds i32, i32* %111, i64 %113
%115 = load i32, i32* %114
%116 = icmp sgt i32 %110, %115
br i1 %116, label %117, label %118
117: ; preds = %105
br label %217
118: ; preds = %105
%119 = load i32*, i32** %13
%120 = load i32, i32* %25
%121 = sext i32 %120 to i64
%122 = getelementptr inbounds i32, i32* %119, i64 %121
%123 = load i32, i32* %122
%124 = load i32, i32* %19
%125 = icmp sge i32 %123, %124
br i1 %125, label %126, label %216
126: ; preds = %118
%127 = load i32*, i32** %13
%128 = load i32, i32* %25
%129 = sext i32 %128 to i64
%130 = getelementptr inbounds i32, i32* %127, i64 %129
%131 = load i32, i32* %130
%132 = load i32, i32* %20
%133 = icmp sle i32 %131, %132
br i1 %133, label %134, label %216
134: ; preds = %126
%135 = load i32*, i32** %14
%136 = load i32, i32* %25
%137 = sext i32 %136 to i64
%138 = getelementptr inbounds i32, i32* %135, i64 %137
%139 = load i32, i32* %138
%140 = load i32, i32* %17
%141 = icmp sge i32 %139, %140
br i1 %141, label %142, label %216
142: ; preds = %134
%143 = load i32*, i32** %14
%144 = load i32, i32* %25
%145 = sext i32 %144 to i64
%146 = getelementptr inbounds i32, i32* %143, i64 %145
%147 = load i32, i32* %146
%148 = load i32, i32* %18
%149 = icmp sle i32 %147, %148
br i1 %149, label %150, label %216
150: ; preds = %142
%151 = load i32*, i32** %14
%152 = load i32, i32* %25
%153 = sext i32 %152 to i64
%154 = getelementptr inbounds i32, i32* %151, i64 %153
%155 = load i32, i32* %154
%156 = load i32*, i32** %13
%157 = load i32, i32* %25
%158 = sext i32 %157 to i64
%159 = getelementptr inbounds i32, i32* %156, i64 %158
%160 = load i32, i32* %159
%161 = sub nsw i32 %155, %160
store i32 %161, i32* %28
%162 = load i32, i32* %28
%163 = load i32, i32* %22
%164 = sub nsw i32 %162, %163
%165 = icmp slt i32 %164, 0
br i1 %165, label %166, label %167
166: ; preds = %150
br label %171
167: ; preds = %150
%168 = load i32, i32* %28
%169 = load i32, i32* %22
%170 = sub nsw i32 %168, %169
br label %171
171: ; preds = %167, %166
%172 = phi i32 [ 0, %166 ], [ %170, %167 ]
store i32 %172, i32* %34
%173 = load i32, i32* %28
%174 = load i32, i32* %22
%175 = add nsw i32 %173, %174
%176 = load i32, i32* %32
%177 = icmp sge i32 %175, %176
br i1 %177, label %178, label %180
178: ; preds = %171
%179 = load i32, i32* %32
br label %184
180: ; preds = %171
%181 = load i32, i32* %28
%182 = load i32, i32* %22
%183 = add nsw i32 %181, %182
br label %184
184: ; preds = %180, %178
%185 = phi i32 [ %179, %178 ], [ %183, %180 ]
store i32 %185, i32* %35
%186 = load i32, i32* %34
store i32 %186, i32* %29
br label %187
187: ; preds = %212, %184
%188 = load i32, i32* %29
%189 = load i32, i32* %35
%190 = icmp sle i32 %188, %189
br i1 %190, label %191, label %215
191: ; preds = %187
%192 = load double*, double** %15
%193 = load i32, i32* %25
%194 = sext i32 %193 to i64
%195 = getelementptr inbounds double, double* %192, i64 %194
%196 = load double, double* %195
%197 = load double, double* %23
%198 = fcmp ogt double %196, %197
br i1 %198, label %199, label %211
199: ; preds = %191
%200 = load double*, double** %15
%201 = load i32, i32* %25
%202 = sext i32 %201 to i64
%203 = getelementptr inbounds double, double* %200, i64 %202
%204 = load double, double* %203
%205 = load double*, double** %26
%206 = load i32, i32* %29
%207 = sext i32 %206 to i64
%208 = getelementptr inbounds double, double* %205, i64 %207
%209 = load double, double* %208
%210 = fadd double %209, %204
store double %210, double* %208
br label %211
211: ; preds = %199, %191
br label %212
212: ; preds = %211
%213 = load i32, i32* %29
%214 = add nsw i32 %213, 1
store i32 %214, i32* %29
br label %187
215: ; preds = %187
br label %216
216: ; preds = %215, %142, %134, %126, %118
br label %217
217: ; preds = %216, %117
%218 = load i32, i32* %25
%219 = add nsw i32 %218, 1
store i32 %219, i32* %25
br label %101
220: ; preds = %101
%221 = load i32, i32* %32
%222 = sub nsw i32 %221, 1
%223 = load i32*, i32** %27
%224 = load double*, double** %26
%225 = call i32 (i32, i32, i32*, double*, ...) bitcast (i32 (...)* @QuickSortIndicesDouble to i32 (i32, i32, i32*, double*, ...)*)(i32 0, i32 %222, i32* %223, double* %224)
%226 = load i32*, i32** %27
%227 = load i32, i32* %32
%228 = load i32, i32* %21
%229 = sub nsw i32 %227, %228
%230 = sub nsw i32 %229, 1
%231 = sext i32 %230 to i64
%232 = getelementptr inbounds i32, i32* %226, i64 %231
%233 = load i32, i32* %232
store i32 %233, i32* %33
%234 = load i32, i32* %33
store i32 %234, i32* %28
store i32 1, i32* %30
%235 = load i32, i32* %33
%236 = add nsw i32 %235, 1
store i32 %236, i32* %25
br label %237
237: ; preds = %255, %220
%238 = load i32, i32* %25
%239 = load i32, i32* %32
%240 = icmp slt i32 %238, %239
br i1 %240, label %241, label %253
241: ; preds = %237
%242 = load double*, double** %26
%243 = load i32, i32* %25
%244 = sext i32 %243 to i64
%245 = getelementptr inbounds double, double* %242, i64 %244
%246 = load double, double* %245
%247 = load double*, double** %26
%248 = load i32, i32* %33
%249 = sext i32 %248 to i64
%250 = getelementptr inbounds double, double* %247, i64 %249
%251 = load double, double* %250
%252 = fcmp oeq double %246, %251
br label %253
253: ; preds = %241, %237
%254 = phi i1 [ false, %237 ], [ %252, %241 ]
br i1 %254, label %255, label %263
255: ; preds = %253
%256 = load i32, i32* %25
%257 = load i32, i32* %28
%258 = add nsw i32 %257, %256
store i32 %258, i32* %28
%259 = load i32, i32* %30
%260 = add nsw i32 %259, 1
store i32 %260, i32* %30
%261 = load i32, i32* %25
%262 = add nsw i32 %261, 1
store i32 %262, i32* %25
br label %237
263: ; preds = %253
%264 = load i32, i32* %33
%265 = sub nsw i32 %264, 1
store i32 %265, i32* %25
br label %266
266: ; preds = %283, %263
%267 = load i32, i32* %25
%268 = icmp sge i32 %267, 0
br i1 %268, label %269, label %281
269: ; preds = %266
%270 = load double*, double** %26
%271 = load i32, i32* %25
%272 = sext i32 %271 to i64
%273 = getelementptr inbounds double, double* %270, i64 %272
%274 = load double, double* %273
%275 = load double*, double** %26
%276 = load i32, i32* %33
%277 = sext i32 %276 to i64
%278 = getelementptr inbounds double, double* %275, i64 %277
%279 = load double, double* %278
%280 = fcmp oeq double %274, %279
br label %281
281: ; preds = %269, %266
%282 = phi i1 [ false, %266 ], [ %280, %269 ]
br i1 %282, label %283, label %291
283: ; preds = %281
%284 = load i32, i32* %25
%285 = load i32, i32* %28
%286 = add nsw i32 %285, %284
store i32 %286, i32* %28
%287 = load i32, i32* %30
%288 = add nsw i32 %287, 1
store i32 %288, i32* %30
%289 = load i32, i32* %25
%290 = add nsw i32 %289, -1
store i32 %290, i32* %25
br label %266
291: ; preds = %281
%292 = load i32, i32* %28
%293 = load i32, i32* %30
%294 = sdiv i32 %292, %293
%295 = sitofp i32 %294 to double
%296 = call double @llvm.rint.f64(double %295)
%297 = fptosi double %296 to i32
store i32 %297, i32* %33
%298 = load double*, double** %26
%299 = load i32, i32* %33
%300 = sext i32 %299 to i64
%301 = getelementptr inbounds double, double* %298, i64 %300
%302 = load double, double* %301
%303 = load double*, double** %24
store double %302, double* %303
%304 = load double*, double** %26
%305 = bitcast double* %304 to i8*
%306 = load i32*, i32** %27
%307 = bitcast i32* %306 to i8*
%308 = load i32, i32* %33
ret i32 %308
}

LLVM IR (-O2):
define i32 @getnbestdiagsum(i32* nocapture %0, i32* nocapture %1, double* nocapture %2, i32 %3, i32 %4, i32 %5, i32 %6, i32 %7, i32 %8, i32 %9, double %10, double* nocapture writeonly %11) {
%13 = icmp sgt i32 %7, %5
%14 = sub nsw i32 %7, %4
%15 = sub nsw i32 %5, %6
%16 = select i1 %13, i32 %14, i32 %15
%17 = tail call i32 @llvm.abs.i32(i32 %16, i1 true)
%18 = add nuw nsw i32 %17, 1
%19 = zext i32 %18 to i64
%20 = shl nuw nsw i64 %19, 3
%21 = tail call i32 (i64, ...) bitcast (i32 (...)* @mymalloc to i32 (i64, ...)*)(i64 %20)
%22 = sext i32 %21 to i64
%23 = inttoptr i64 %22 to double*
%24 = shl nuw nsw i64 %19, 2
%25 = tail call i32 (i64, ...) bitcast (i32 (...)* @mymalloc to i32 (i64, ...)*)(i64 %24)
%26 = sext i32 %25 to i64
%27 = inttoptr i64 %26 to i32*
%28 = icmp eq i32 %16, 0
br i1 %28, label %101, label %29
29: ; preds = %12
%30 = inttoptr i64 %22 to i8*
%31 = call i32 @llvm.umax.i32(i32 %17, i32 1)
%32 = zext i32 %31 to i64
%33 = shl nuw nsw i64 %32, 3
%34 = call i32 @llvm.umax.i32(i32 %17, i32 1)
%35 = zext i32 %34 to i64
%36 = icmp ult i32 %34, 8
br i1 %36, label %99, label %37
37: ; preds = %29
%38 = and i64 %35, 2147483640
%39 = add nsw i64 %38, -8
%40 = lshr exact i64 %39, 3
%41 = add nuw nsw i64 %40, 1
%42 = and i64 %41, 3
%43 = icmp ult i64 %39, 24
br i1 %43, label %80, label %44
44: ; preds = %37
%45 = and i64 %41, 4611686018427387900
br label %46
46: ; preds = %46, %44
%47 = phi i64 [ 0, %44 ], [ %76, %46 ]
%48 = phi <4 x i32> [ <i32 0, i32 1, i32 2, i32 3>, %44 ], [ %77, %46 ]
%49 = phi i64 [ 0, %44 ], [ %78, %46 ]
%50 = add <4 x i32> %48, <i32 4, i32 4, i32 4, i32 4>
%51 = getelementptr inbounds i32, i32* %27, i64 %47
%52 = bitcast i32* %51 to <4 x i32>*
store <4 x i32> %48, <4 x i32>* %52
%53 = getelementptr inbounds i32, i32* %51, i64 4
%54 = bitcast i32* %53 to <4 x i32>*
store <4 x i32> %50, <4 x i32>* %54
%55 = or i64 %47, 8
%56 = add <4 x i32> %48, <i32 8, i32 8, i32 8, i32 8>
%57 = add <4 x i32> %48, <i32 12, i32 12, i32 12, i32 12>
%58 = getelementptr inbounds i32, i32* %27, i64 %55
%59 = bitcast i32* %58 to <4 x i32>*
store <4 x i32> %56, <4 x i32>* %59
%60 = getelementptr inbounds i32, i32* %58, i64 4
%61 = bitcast i32* %60 to <4 x i32>*
store <4 x i32> %57, <4 x i32>* %61
%62 = or i64 %47, 16
%63 = add <4 x i32> %48, <i32 16, i32 16, i32 16, i32 16>
%64 = add <4 x i32> %48, <i32 20, i32 20, i32 20, i32 20>
%65 = getelementptr inbounds i32, i32* %27, i64 %62
%66 = bitcast i32* %65 to <4 x i32>*
store <4 x i32> %63, <4 x i32>* %66
%67 = getelementptr inbounds i32, i32* %65, i64 4
%68 = bitcast i32* %67 to <4 x i32>*
store <4 x i32> %64, <4 x i32>* %68
%69 = or i64 %47, 24
%70 = add <4 x i32> %48, <i32 24, i32 24, i32 24, i32 24>
%71 = add <4 x i32> %48, <i32 28, i32 28, i32 28, i32 28>
%72 = getelementptr inbounds i32, i32* %27, i64 %69
%73 = bitcast i32* %72 to <4 x i32>*
store <4 x i32> %70, <4 x i32>* %73
%74 = getelementptr inbounds i32, i32* %72, i64 4
%75 = bitcast i32* %74 to <4 x i32>*
store <4 x i32> %71, <4 x i32>* %75
%76 = add nuw i64 %47, 32
%77 = add <4 x i32> %48, <i32 32, i32 32, i32 32, i32 32>
%78 = add i64 %49, 4
%79 = icmp eq i64 %78, %45
br i1 %79, label %80, label %46
80: ; preds = %46, %37
%81 = phi i64 [ 0, %37 ], [ %76, %46 ]
%82 = phi <4 x i32> [ <i32 0, i32 1, i32 2, i32 3>, %37 ], [ %77, %46 ]
%83 = icmp eq i64 %42, 0
br i1 %83, label %97, label %84
84: ; preds = %80, %84
%85 = phi i64 [ %93, %84 ], [ %81, %80 ]
%86 = phi <4 x i32> [ %94, %84 ], [ %82, %80 ]
%87 = phi i64 [ %95, %84 ], [ 0, %80 ]
%88 = add <4 x i32> %86, <i32 4, i32 4, i32 4, i32 4>
%89 = getelementptr inbounds i32, i32* %27, i64 %85
%90 = bitcast i32* %89 to <4 x i32>*
store <4 x i32> %86, <4 x i32>* %90
%91 = getelementptr inbounds i32, i32* %89, i64 4
%92 = bitcast i32* %91 to <4 x i32>*
store <4 x i32> %88, <4 x i32>* %92
%93 = add nuw i64 %85, 8
%94 = add <4 x i32> %86, <i32 8, i32 8, i32 8, i32 8>
%95 = add i64 %87, 1
%96 = icmp eq i64 %95, %42
br i1 %96, label %97, label %84
97: ; preds = %84, %80
%98 = icmp eq i64 %38, %35
br i1 %98, label %101, label %99
99: ; preds = %29, %97
%100 = phi i64 [ 0, %29 ], [ %38, %97 ]
br label %105
101: ; preds = %105, %97, %12
%102 = icmp sgt i32 %3, 0
br i1 %102, label %103, label %179
103: ; preds = %101
%104 = zext i32 %3 to i64
br label %111
105: ; preds = %99, %105
%106 = phi i64 [ %109, %105 ], [ %100, %99 ]
%107 = getelementptr inbounds i32, i32* %27, i64 %106
%108 = trunc i64 %106 to i32
store i32 %108, i32* %107
%109 = add nuw nsw i64 %106, 1
%110 = icmp eq i64 %109, %35
br i1 %110, label %101, label %105
111: ; preds = %103, %176
%112 = phi i64 [ 0, %103 ], [ %177, %176 ]
%113 = getelementptr inbounds i32, i32* %0, i64 %112
%114 = load i32, i32* %113
%115 = getelementptr inbounds i32, i32* %1, i64 %112
%116 = load i32, i32* %115
%117 = icmp slt i32 %116, %114
%118 = icmp slt i32 %114, %6
%119 = or i1 %118, %117
%120 = icmp sgt i32 %114, %7
%121 = or i1 %120, %119
%122 = icmp slt i32 %116, %4
%123 = select i1 %121, i1 true, i1 %122
%124 = icmp sgt i32 %116, %5
%125 = select i1 %123, i1 true, i1 %124
br i1 %125, label %176, label %126
126: ; preds = %111
%127 = sub nsw i32 %116, %114
%128 = sub nsw i32 %127, %9
%129 = icmp sgt i32 %128, 0
%130 = select i1 %129, i32 %128, i32 0
%131 = add nsw i32 %127, %9
%132 = icmp slt i32 %131, %17
%133 = select i1 %132, i32 %131, i32 %17
%134 = icmp sgt i32 %130, %133
br i1 %134, label %176, label %135
135: ; preds = %126
%136 = getelementptr inbounds double, double* %2, i64 %112
%137 = zext i32 %130 to i64
%138 = add i32 %133, 1
%139 = zext i32 %138 to i64
%140 = sub nsw i64 %139, %137
%141 = xor i64 %137, -1
%142 = and i64 %140, 1
%143 = icmp eq i64 %142, 0
br i1 %143, label %153, label %144
144: ; preds = %135
%145 = load double, double* %136
%146 = fcmp ogt double %145, %10
br i1 %146, label %147, label %151
147: ; preds = %144
%148 = getelementptr inbounds double, double* %23, i64 %137
%149 = load double, double* %148
%150 = fadd double %145, %149
store double %150, double* %148
br label %151
151: ; preds = %147, %144
%152 = add nuw nsw i64 %137, 1
br label %153
153: ; preds = %151, %135
%154 = phi i64 [ %137, %135 ], [ %152, %151 ]
%155 = sub nsw i64 0, %139
%156 = icmp eq i64 %141, %155
br i1 %156, label %176, label %157
157: ; preds = %153, %173
%158 = phi i64 [ %174, %173 ], [ %154, %153 ]
%159 = load double, double* %136
%160 = fcmp ogt double %159, %10
br i1 %160, label %161, label %165
161: ; preds = %157
%162 = getelementptr inbounds double, double* %23, i64 %158
%163 = load double, double* %162
%164 = fadd double %159, %163
store double %164, double* %162
br label %165
165: ; preds = %157, %161
%166 = load double, double* %136
%167 = fcmp ogt double %166, %10
br i1 %167, label %168, label %173
168: ; preds = %165
%169 = add nuw nsw i64 %158, 1
%170 = getelementptr inbounds double, double* %23, i64 %169
%171 = load double, double* %170
%172 = fadd double %166, %171
store double %172, double* %170
br label %173
173: ; preds = %168, %165
%174 = add nuw nsw i64 %158, 2
%175 = icmp eq i64 %174, %139
br i1 %175, label %176, label %157
176: ; preds = %153, %173, %126, %111
%177 = add nuw nsw i64 %112, 1
%178 = icmp eq i64 %177, %104
br i1 %178, label %179, label %111
179: ; preds = %176, %101
%180 = add nsw i32 %17, -1
%181 = tail call i32 (i32, i32, i32*, double*, ...) bitcast (i32 (...)* @QuickSortIndicesDouble to i32 (i32, i32, i32*, double*, ...)*)(i32 0, i32 %180, i32* %27, double* %23)
%182 = xor i32 %8, -1
%183 = add i32 %17, %182
%184 = sext i32 %183 to i64
%185 = getelementptr inbounds i32, i32* %27, i64 %184
%186 = load i32, i32* %185
%187 = add i32 %186, 1
%188 = icmp slt i32 %187, %17
br i1 %188, label %189, label %209
189: ; preds = %179
%190 = sext i32 %186 to i64
%191 = getelementptr inbounds double, double* %23, i64 %190
%192 = load double, double* %191
%193 = sext i32 %187 to i64
%194 = sub i32 %17, %186
br label %195
195: ; preds = %189, %202
%196 = phi i64 [ %193, %189 ], [ %206, %202 ]
%197 = phi i32 [ %186, %189 ], [ %204, %202 ]
%198 = phi i32 [ 1, %189 ], [ %205, %202 ]
%199 = getelementptr inbounds double, double* %23, i64 %196
%200 = load double, double* %199
%201 = fcmp oeq double %200, %192
br i1 %201, label %202, label %209
202: ; preds = %195
%203 = trunc i64 %196 to i32
%204 = add nsw i32 %197, %203
%205 = add nuw nsw i32 %198, 1
%206 = add nsw i64 %196, 1
%207 = trunc i64 %206 to i32
%208 = icmp eq i32 %17, %207
br i1 %208, label %209, label %195
209: ; preds = %195, %202, %179
%210 = phi i32 [ 1, %179 ], [ %194, %202 ], [ %198, %195 ]
%211 = phi i32 [ %186, %179 ], [ %204, %202 ], [ %197, %195 ]
%212 = icmp sgt i32 %186, 0
br i1 %212, label %213, label %233
213: ; preds = %209
%214 = zext i32 %186 to i64
%215 = getelementptr inbounds double, double* %23, i64 %214
%216 = load double, double* %215
%217 = zext i32 %186 to i64
%218 = add i32 %210, %186
br label %219
219: ; preds = %213, %228
%220 = phi i64 [ %217, %213 ], [ %223, %228 ]
%221 = phi i32 [ %211, %213 ], [ %230, %228 ]
%222 = phi i32 [ %210, %213 ], [ %231, %228 ]
%223 = add nsw i64 %220, -1
%224 = and i64 %223, 4294967295
%225 = getelementptr inbounds double, double* %23, i64 %224
%226 = load double, double* %225
%227 = fcmp oeq double %226, %216
br i1 %227, label %228, label %233
228: ; preds = %219
%229 = trunc i64 %223 to i32
%230 = add nsw i32 %221, %229
%231 = add nuw nsw i32 %222, 1
%232 = icmp ugt i64 %220, 1
br i1 %232, label %219, label %233
233: ; preds = %219, %228, %209
%234 = phi i32 [ %210, %209 ], [ %218, %228 ], [ %222, %219 ]
%235 = phi i32 [ %211, %209 ], [ %230, %228 ], [ %221, %219 ]
%236 = sdiv i32 %235, %234
%237 = sext i32 %236 to i64
%238 = getelementptr inbounds double, double* %23, i64 %237
%239 = load double, double* %238
store double %239, double* %11
%240 = inttoptr i64 %22 to i8*
tail call void @free(i8* %240)
%241 = inttoptr i64 %26 to i8*
tail call void @free(i8* %241)
ret i32 %236
}