#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 */
int /*<<< orphan*/ free(size_t *);
scalar_t__ malloc(int);
int calculate_ngtdm(int *image, char *mask, int *size, int *bb, int *strides,
                    int *angles, int Na, int Nd, double *ngtdm, int Ng) {
  size_t Ni;
  size_t i, j;
  size_t *cur_idx = (size_t *)malloc(sizeof *cur_idx * Nd);
  size_t a, d;
  int gl;
  double count, sum, diff;
  size_t ngtdm_idx, ngtdm_idx_max = Ng * 3;
  for (gl = 0; gl < Ng; gl++) {
    ngtdm[gl * 3 + 2] = gl + 1;
  }
  Ni = size[0];
  i = bb[0] * strides[0];
  for (d = 1; d < Nd; d++) {
    i += bb[d] * strides[d];
    Ni *= size[d];
  }
  for (; i < Ni; i++) {
    for (d = Nd - 1; d > 0; d--) {
      cur_idx[d] = (i % strides[d - 1]) / strides[d];
      if (cur_idx[d] > bb[Nd + d]) {
        i += (size[d] - cur_idx[d] + bb[d]) * strides[d];
        cur_idx[d] = bb[d];
      } else if (cur_idx[d] < bb[d]) {
        i += (bb[d] - cur_idx[d]) * strides[d];
        cur_idx[d] = bb[d];
      }
    }
    cur_idx[0] = i / strides[0];
    if (cur_idx[0] > bb[Nd])
      break;
    if (mask[i]) {
      count = 0;
      sum = 0;
      for (a = 0; a < Na; a++) {
        j = i;
        for (d = 0; d < Nd; d++) {
          if (cur_idx[d] + angles[a * Nd + d] < bb[d] ||
              cur_idx[d] + angles[a * Nd + d] > bb[Nd + d]) {
            j = i;
            break;
          }
          j += angles[a * Nd + d] * strides[d];
        }
        if (j != i && mask[j]) {
          count++;
          sum += image[j];
        }
      }
      if (count == 0)
        diff = 0;
      else
        diff = (double)image[i] - (sum / count);
      if (diff < 0)
        diff *= -1;
      ngtdm_idx = (image[i] - 1) * 3;
      if (image[i] <= 0 || ngtdm_idx >= ngtdm_idx_max) {
        free(cur_idx);
        return 0;
      }
      ngtdm[ngtdm_idx]++;
      ngtdm[ngtdm_idx + 1] += diff;
    }
  }
  free(cur_idx);
  return 1;
}


x86 assembly (clang, -O0):
<calculate_ngtdm>:
push   %rbp
mov    %rsp,%rbp
sub    $0xa0,%rsp
mov    0x28(%rbp),%eax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%eax
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
mov    %r8,-0x30(%rbp)
mov    %r9,-0x38(%rbp)
movslq 0x18(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   1040 <malloc@plt>
mov    %rax,-0x58(%rbp)
imul   $0x3,0x28(%rbp),%eax
cltq
mov    %rax,-0x98(%rbp)
movl   $0x0,-0x6c(%rbp)
mov    -0x6c(%rbp),%eax
cmp    0x28(%rbp),%eax
jge    11ae <calculate_ngtdm+0x8e>
mov    -0x6c(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm0
mov    0x20(%rbp),%rax
imul   $0x3,-0x6c(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x6c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x6c(%rbp)
jmp    1177 <calculate_ngtdm+0x57>
mov    -0x20(%rbp),%rax
movslq (%rax),%rax
mov    %rax,-0x40(%rbp)
mov    -0x28(%rbp),%rax
mov    (%rax),%eax
mov    -0x30(%rbp),%rcx
imul   (%rcx),%eax
cltq
mov    %rax,-0x48(%rbp)
movq   $0x1,-0x68(%rbp)
mov    -0x68(%rbp),%rax
movslq 0x18(%rbp),%rcx
cmp    %rcx,%rax
jae    122c <calculate_ngtdm+0x10c>
mov    -0x28(%rbp),%rax
mov    -0x68(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    -0x30(%rbp),%rcx
mov    -0x68(%rbp),%rdx
imul   (%rcx,%rdx,4),%eax
cltq
add    -0x48(%rbp),%rax
mov    %rax,-0x48(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x68(%rbp),%rcx
movslq (%rax,%rcx,4),%rax
imul   -0x40(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x68(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x68(%rbp)
jmp    11d4 <calculate_ngtdm+0xb4>
jmp    1231 <calculate_ngtdm+0x111>
mov    -0x48(%rbp),%rax
cmp    -0x40(%rbp),%rax
jae    1672 <calculate_ngtdm+0x552>
mov    0x18(%rbp),%eax
sub    $0x1,%eax
cltq
mov    %rax,-0x68(%rbp)
cmpq   $0x0,-0x68(%rbp)
jbe    1397 <calculate_ngtdm+0x277>
mov    -0x48(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    -0x68(%rbp),%rdx
sub    $0x1,%rdx
movslq (%rcx,%rdx,4),%rcx
xor    %edx,%edx
div    %rcx
mov    %rdx,%rax
mov    -0x30(%rbp),%rcx
mov    -0x68(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
xor    %edx,%edx
div    %rcx
mov    %rax,%rdx
mov    -0x58(%rbp),%rax
mov    -0x68(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
mov    -0x58(%rbp),%rax
mov    -0x68(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x28(%rbp),%rcx
movslq 0x18(%rbp),%rdx
add    -0x68(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
cmp    %rcx,%rax
jbe    1313 <calculate_ngtdm+0x1f3>
mov    -0x20(%rbp),%rax
mov    -0x68(%rbp),%rcx
movslq (%rax,%rcx,4),%rax
mov    -0x58(%rbp),%rcx
mov    -0x68(%rbp),%rdx
sub    (%rcx,%rdx,8),%rax
mov    -0x28(%rbp),%rcx
mov    -0x68(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
add    %rcx,%rax
mov    -0x30(%rbp),%rcx
mov    -0x68(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
imul   %rcx,%rax
add    -0x48(%rbp),%rax
mov    %rax,-0x48(%rbp)
mov    -0x28(%rbp),%rax
mov    -0x68(%rbp),%rcx
movslq (%rax,%rcx,4),%rdx
mov    -0x58(%rbp),%rax
mov    -0x68(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
jmp    1381 <calculate_ngtdm+0x261>
mov    -0x58(%rbp),%rax
mov    -0x68(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x28(%rbp),%rcx
mov    -0x68(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
cmp    %rcx,%rax
jae    137c <calculate_ngtdm+0x25c>
mov    -0x28(%rbp),%rax
mov    -0x68(%rbp),%rcx
movslq (%rax,%rcx,4),%rax
mov    -0x58(%rbp),%rcx
mov    -0x68(%rbp),%rdx
sub    (%rcx,%rdx,8),%rax
mov    -0x30(%rbp),%rcx
mov    -0x68(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
imul   %rcx,%rax
add    -0x48(%rbp),%rax
mov    %rax,-0x48(%rbp)
mov    -0x28(%rbp),%rax
mov    -0x68(%rbp),%rcx
movslq (%rax,%rcx,4),%rdx
mov    -0x58(%rbp),%rax
mov    -0x68(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
jmp    1381 <calculate_ngtdm+0x261>
jmp    1386 <calculate_ngtdm+0x266>
mov    -0x68(%rbp),%rax
add    $0xffffffffffffffff,%rax
mov    %rax,-0x68(%rbp)
jmp    124b <calculate_ngtdm+0x12b>
mov    -0x48(%rbp),%rax
mov    -0x30(%rbp),%rcx
movslq (%rcx),%rcx
xor    %edx,%edx
div    %rcx
mov    %rax,%rcx
mov    -0x58(%rbp),%rax
mov    %rcx,(%rax)
mov    -0x58(%rbp),%rax
mov    (%rax),%rax
mov    -0x28(%rbp),%rcx
movslq 0x18(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
cmp    %rcx,%rax
jbe    13d2 <calculate_ngtdm+0x2b2>
jmp    1672 <calculate_ngtdm+0x552>
mov    -0x18(%rbp),%rax
mov    -0x48(%rbp),%rcx
cmpb   $0x0,(%rax,%rcx,1)
je     165c <calculate_ngtdm+0x53c>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x78(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x80(%rbp)
movq   $0x0,-0x60(%rbp)
mov    -0x60(%rbp),%rax
movslq 0x10(%rbp),%rcx
cmp    %rcx,%rax
jae    1557 <calculate_ngtdm+0x437>
mov    -0x48(%rbp),%rax
mov    %rax,-0x50(%rbp)
movq   $0x0,-0x68(%rbp)
mov    -0x68(%rbp),%rax
movslq 0x18(%rbp),%rcx
cmp    %rcx,%rax
jae    14f5 <calculate_ngtdm+0x3d5>
mov    -0x58(%rbp),%rax
mov    -0x68(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x38(%rbp),%rcx
mov    -0x60(%rbp),%rdx
movslq 0x18(%rbp),%rsi
imul   %rsi,%rdx
add    -0x68(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
add    %rcx,%rax
mov    -0x28(%rbp),%rcx
mov    -0x68(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
cmp    %rcx,%rax
jb     14aa <calculate_ngtdm+0x38a>
mov    -0x58(%rbp),%rax
mov    -0x68(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x38(%rbp),%rcx
mov    -0x60(%rbp),%rdx
movslq 0x18(%rbp),%rsi
imul   %rsi,%rdx
add    -0x68(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
add    %rcx,%rax
mov    -0x28(%rbp),%rcx
movslq 0x18(%rbp),%rdx
add    -0x68(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
cmp    %rcx,%rax
jbe    14b7 <calculate_ngtdm+0x397>
mov    -0x48(%rbp),%rax
mov    %rax,-0x50(%rbp)
jmp    14f5 <calculate_ngtdm+0x3d5>
mov    -0x38(%rbp),%rax
mov    -0x60(%rbp),%rcx
movslq 0x18(%rbp),%rdx
imul   %rdx,%rcx
add    -0x68(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    -0x30(%rbp),%rcx
mov    -0x68(%rbp),%rdx
imul   (%rcx,%rdx,4),%eax
cltq
add    -0x50(%rbp),%rax
mov    %rax,-0x50(%rbp)
mov    -0x68(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x68(%rbp)
jmp    141d <calculate_ngtdm+0x2fd>
mov    -0x50(%rbp),%rax
cmp    -0x48(%rbp),%rax
je     1541 <calculate_ngtdm+0x421>
mov    -0x18(%rbp),%rax
mov    -0x50(%rbp),%rcx
movsbl (%rax,%rcx,1),%eax
cmp    $0x0,%eax
je     1541 <calculate_ngtdm+0x421>
movsd  0xae8(%rip),%xmm0
addsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x50(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm0
addsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
jmp    1546 <calculate_ngtdm+0x426>
mov    -0x60(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x60(%rbp)
jmp    13fc <calculate_ngtdm+0x2dc>
movsd  -0x78(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    157f <calculate_ngtdm+0x45f>
jp     157f <calculate_ngtdm+0x45f>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
jmp    15a2 <calculate_ngtdm+0x482>
mov    -0x10(%rbp),%rax
mov    -0x48(%rbp),%rcx
cvtsi2sdl (%rax,%rcx,4),%xmm0
movsd  -0x80(%rbp),%xmm1
divsd  -0x78(%rbp),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
xorps  %xmm0,%xmm0
ucomisd -0x88(%rbp),%xmm0
ff
jbe    15cb <calculate_ngtdm+0x4ab>
movsd  0xa45(%rip),%xmm0
mulsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0x88(%rbp)
ff
mov    -0x10(%rbp),%rax
mov    -0x48(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
sub    $0x1,%eax
imul   $0x3,%eax,%eax
cltq
mov    %rax,-0x90(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x48(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
jle    160b <calculate_ngtdm+0x4eb>
mov    -0x90(%rbp),%rax
cmp    -0x98(%rbp),%rax
jb     1620 <calculate_ngtdm+0x500>
mov    -0x58(%rbp),%rdi
call   1030 <free@plt>
movl   $0x0,-0x4(%rbp)
jmp    1682 <calculate_ngtdm+0x562>
mov    0x20(%rbp),%rax
mov    -0x90(%rbp),%rcx
movsd  0x9d5(%rip),%xmm0
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x88(%rbp),%xmm0
ff
mov    0x20(%rbp),%rax
mov    -0x90(%rbp),%rcx
addsd  0x8(%rax,%rcx,8),%xmm0
movsd  %xmm0,0x8(%rax,%rcx,8)
jmp    1661 <calculate_ngtdm+0x541>
mov    -0x48(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x48(%rbp)
jmp    1231 <calculate_ngtdm+0x111>
mov    -0x58(%rbp),%rdi
call   1030 <free@plt>
movl   $0x1,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0xa0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<calculate_ngtdm>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x78,%rsp
mov    %r9,0x30(%rsp)
mov    %r8,0x18(%rsp)
mov    %rcx,%r13
mov    %rdx,%r12
mov    %rsi,0x20(%rsp)
mov    %rdi,0x8(%rsp)
mov    0xc8(%rsp),%ebx
mov    0xc0(%rsp),%rbp
mov    0xb8(%rsp),%eax
lea    0x0(,%rax,8),%edi
call   1040 <malloc@plt>
mov    %rbx,%r10
mov    %rax,%r14
test   %r10d,%r10d
jle    1238 <calculate_ngtdm+0x118>
mov    %r10d,%edx
lea    -0x1(%rdx),%rax
mov    %edx,%ecx
and    $0x3,%ecx
cmp    $0x3,%rax
jae    118f <calculate_ngtdm+0x6f>
xor    %eax,%eax
jmp    11f6 <calculate_ngtdm+0xd6>
and    $0xfffffffc,%edx
mov    $0x9,%esi
xor    %eax,%eax
nopl   0x0(%rax)
lea    0x1(%rax),%edi
xorps  %xmm0,%xmm0
cvtsi2sd %edi,%xmm0
lea    0x0(,%rsi,8),%rdi
movsd  %xmm0,-0x38(%rbp,%rsi,8)
lea    0x2(%rax),%ebx
xorps  %xmm0,%xmm0
cvtsi2sd %ebx,%xmm0
movsd  %xmm0,-0x20(%rbp,%rsi,8)
lea    0x3(%rax),%ebx
xorps  %xmm0,%xmm0
cvtsi2sd %ebx,%xmm0
movsd  %xmm0,-0x8(%rbp,%rsi,8)
add    $0x4,%rax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
or     $0x10,%rdi
movsd  %xmm0,0x0(%rbp,%rdi,1)
add    $0xc,%rsi
cmp    %rax,%rdx
jne    11a0 <calculate_ngtdm+0x80>
test   %rcx,%rcx
je     1238 <calculate_ngtdm+0x118>
lea    (%rax,%rax,2),%rdx
lea    0x10(,%rdx,8),%rdx
add    %rbp,%rdx
add    $0x1,%eax
shl    $0x3,%rcx
lea    (%rcx,%rcx,2),%rcx
xor    %esi,%esi
nopw   0x0(%rax,%rax,1)
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,(%rdx,%rsi,1)
add    $0x18,%rsi
add    $0x1,%eax
cmp    %rsi,%rcx
jne    1220 <calculate_ngtdm+0x100>
mov    0xb8(%rsp),%esi
movslq %esi,%rdi
movslq (%r12),%rbp
movslq 0x0(%r13),%r15
mov    0x18(%rsp),%r8
movslq (%r8),%rax
mov    %rax,0x50(%rsp)
imul   %rax,%r15
cmp    $0x2,%esi
jb     12f4 <calculate_ngtdm+0x1d4>
lea    -0x1(%rdi),%rdx
mov    %rdi,%r9
cmp    $0x2,%esi
jne    1277 <calculate_ngtdm+0x157>
mov    $0x1,%eax
jmp    12d4 <calculate_ngtdm+0x1b4>
mov    %rdx,%rcx
and    $0xfffffffffffffffe,%rcx
neg    %rcx
mov    $0x1,%eax
cs nopw 0x0(%rax,%rax,1)
movslq 0x0(%r13,%rax,4),%rdi
movslq (%r8,%rax,4),%rbx
imul   %rdi,%rbx
add    %r15,%rbx
movslq (%r12,%rax,4),%rdi
imul   %rbp,%rdi
movslq 0x4(%r13,%rax,4),%rbp
movslq 0x4(%r8,%rax,4),%r15
imul   %rbp,%r15
add    %rbx,%r15
movslq 0x4(%r12,%rax,4),%rbp
imul   %rdi,%rbp
lea    (%rcx,%rax,1),%rdi
add    $0x2,%rdi
add    $0x2,%rax
cmp    $0x1,%rdi
jne    1290 <calculate_ngtdm+0x170>
test   $0x1,%dl
mov    %r9,%rdi
je     12f4 <calculate_ngtdm+0x1d4>
movslq 0x0(%r13,%rax,4),%rcx
movslq (%r8,%rax,4),%rdx
imul   %rcx,%rdx
add    %rdx,%r15
movslq (%r12,%rax,4),%rax
imul   %rax,%rbp
mov    %rbp,0x48(%rsp)
cmp    %rbp,%r15
jae    164e <calculate_ngtdm+0x52e>
mov    0xb0(%rsp),%ecx
lea    (%r10,%r10,2),%eax
mov    %eax,0x14(%rsp)
lea    -0x1(%rdi),%rax
mov    %rax,0x40(%rsp)
lea    0x0(,%rdi,4),%r9
add    %r13,%r9
movslq 0x0(%r13,%rdi,4),%rax
mov    %rax,0x38(%rsp)
cmp    $0x2,%rdi
mov    $0x1,%ebx
mov    $0x1,%r11d
cmovae %rdi,%r11
cmp    $0x2,%ecx
mov    $0x1,%eax
cmovae %ecx,%eax
cltq
mov    %rax,0x60(%rsp)
lea    -0x4(%r12,%rdi,4),%rax
mov    %rax,0x58(%rsp)
lea    -0x1(%rdi,%rdi,1),%rax
lea    0x0(%r13,%rax,4),%rax
mov    %rax,0x70(%rsp)
lea    (%r14,%rdi,8),%rax
add    $0xfffffffffffffff8,%rax
mov    %rax,0x28(%rsp)
sub    %rdi,%rbx
lea    (%r8,%rdi,4),%rax
add    $0xfffffffffffffff8,%rax
mov    %rax,0x68(%rsp)
shl    $0x2,%rdi
movapd 0xc7a(%rip),%xmm0
xorpd  %xmm1,%xmm1
movapd 0xc7e(%rip),%xmm2
movsd  0xc56(%rip),%xmm3
jmp    1431 <calculate_ngtdm+0x311>
nop
ucomisd %xmm1,%xmm5
jne    1622 <calculate_ngtdm+0x502>
jp     1622 <calculate_ngtdm+0x502>
mov    0x8(%rsp),%rax
mov    (%rax,%r15,4),%eax
xorpd  %xmm6,%xmm6
test   %eax,%eax
jle    1642 <calculate_ngtdm+0x522>
lea    (%rax,%rax,2),%eax
add    $0xfffffffd,%eax
cmp    0x14(%rsp),%eax
jae    1642 <calculate_ngtdm+0x522>
cltq
movapd %xmm6,%xmm4
xorpd  %xmm0,%xmm4
movapd %xmm6,%xmm5
cmpltsd %xmm1,%xmm5
andpd  %xmm5,%xmm4
andnpd %xmm6,%xmm5
orpd   %xmm4,%xmm5
mov    0xc0(%rsp),%rcx
movupd (%rcx,%rax,8),%xmm4
movapd %xmm2,%xmm6
unpcklpd %xmm5,%xmm6
addpd  %xmm4,%xmm6
movupd %xmm6,(%rcx,%rax,8)
add    $0x1,%r15
cmp    0x48(%rsp),%r15
jae    164e <calculate_ngtdm+0x52e>
mov    0x40(%rsp),%rax
test   %eax,%eax
je     1510 <calculate_ngtdm+0x3f0>
mov    (%r8,%rax,4),%ebp
xor    %r12d,%r12d
jmp    1484 <calculate_ngtdm+0x364>
nopw   0x0(%rax,%rax,1)
mov    0x58(%rsp),%rcx
movslq (%rcx,%r12,4),%rdx
sub    %rax,%rdx
movslq -0x4(%r9,%r12,4),%rcx
add    %rcx,%rdx
imul   %r8,%rdx
add    %rdx,%r15
mov    0x28(%rsp),%rax
mov    %rcx,(%rax,%r12,8)
add    $0xffffffffffffffff,%r12
mov    %r10d,%ebp
cmp    %r12,%rbx
je     1510 <calculate_ngtdm+0x3f0>
mov    0x68(%rsp),%rax
mov    (%rax,%r12,4),%r10d
movslq %r10d,%rcx
mov    %r15,%rax
or     %rcx,%rax
shr    $0x20,%rax
je     14b0 <calculate_ngtdm+0x390>
mov    %r15,%rax
xor    %edx,%edx
div    %rcx
mov    %rdx,%rax
jmp    14ba <calculate_ngtdm+0x39a>
nopl   0x0(%rax)
mov    %r15d,%eax
xor    %edx,%edx
div    %r10d
mov    %edx,%eax
movslq %ebp,%r8
mov    %rax,%rcx
or     %r8,%rcx
shr    $0x20,%rcx
je     14d0 <calculate_ngtdm+0x3b0>
xor    %edx,%edx
div    %r8
jmp    14d4 <calculate_ngtdm+0x3b4>
xor    %edx,%edx
div    %ebp
mov    0x28(%rsp),%rcx
mov    %rax,(%rcx,%r12,8)
mov    0x70(%rsp),%rcx
movslq (%rcx,%r12,4),%rcx
cmp    %rcx,%rax
ja     1450 <calculate_ngtdm+0x330>
movslq -0x4(%r9,%r12,4),%rcx
cmp    %rcx,%rax
jae    1474 <calculate_ngtdm+0x354>
mov    %rcx,%rdx
sub    %rax,%rdx
jmp    1464 <calculate_ngtdm+0x344>
nopl   0x0(%rax,%rax,1)
mov    %r15,%rax
mov    0x50(%rsp),%rcx
or     %rcx,%rax
shr    $0x20,%rax
je     1530 <calculate_ngtdm+0x410>
mov    %r15,%rax
xor    %edx,%edx
div    %rcx
jmp    1537 <calculate_ngtdm+0x417>
nopl   0x0(%rax,%rax,1)
mov    %r15d,%eax
xor    %edx,%edx
div    %ecx
mov    %rax,(%r14)
cmp    0x38(%rsp),%rax
mov    0x18(%rsp),%r8
ja     164e <calculate_ngtdm+0x52e>
mov    0x20(%rsp),%rax
cmpb   $0x0,(%rax,%r15,1)
je     1422 <calculate_ngtdm+0x302>
cmpl   $0x0,0xb0(%rsp)
je     13c0 <calculate_ngtdm+0x2a0>
xorpd  %xmm4,%xmm4
mov    0x30(%rsp),%rax
xorpd  %xmm5,%xmm5
xor    %edx,%edx
jmp    159f <calculate_ngtdm+0x47f>
nopl   0x0(%rax)
mov    0xb8(%rsp),%esi
mov    %r12,%rdi
mov    %r10,%rbx
add    $0x1,%rdx
add    %rdi,%rax
cmp    0x60(%rsp),%rdx
je     13b0 <calculate_ngtdm+0x290>
test   %esi,%esi
je     158d <calculate_ngtdm+0x46d>
mov    %rbx,%r10
mov    %rdi,%r12
xor    %ebp,%ebp
mov    %r15,%rcx
xchg   %ax,%ax
movslq (%rax,%rbp,4),%rdi
mov    (%r14,%rbp,8),%rbx
add    %rdi,%rbx
movslq 0x0(%r13,%rbp,4),%rsi
cmp    %rsi,%rbx
jb     1580 <calculate_ngtdm+0x460>
movslq (%r9,%rbp,4),%rsi
cmp    %rsi,%rbx
ja     1580 <calculate_ngtdm+0x460>
movslq (%r8,%rbp,4),%rsi
imul   %rdi,%rsi
add    %rsi,%rcx
add    $0x1,%rbp
cmp    %rbp,%r11
jne    15b0 <calculate_ngtdm+0x490>
cmp    %r15,%rcx
mov    %r12,%rdi
jne    15f3 <calculate_ngtdm+0x4d3>
mov    0xb8(%rsp),%esi
jmp    158a <calculate_ngtdm+0x46a>
mov    0x20(%rsp),%rsi
cmpb   $0x0,(%rsi,%rcx,1)
mov    %r10,%rbx
je     1616 <calculate_ngtdm+0x4f6>
addsd  %xmm3,%xmm5
mov    0x8(%rsp),%rsi
xorps  %xmm6,%xmm6
cvtsi2sdl (%rsi,%rcx,4),%xmm6
addsd  %xmm6,%xmm4
mov    0xb8(%rsp),%esi
jmp    158d <calculate_ngtdm+0x46d>
mov    0x8(%rsp),%rax
mov    (%rax,%r15,4),%eax
xorps  %xmm6,%xmm6
cvtsi2sd %eax,%xmm6
divsd  %xmm5,%xmm4
subsd  %xmm4,%xmm6
test   %eax,%eax
jg     13d5 <calculate_ngtdm+0x2b5>
mov    %r14,%rdi
call   1030 <free@plt>
xor    %eax,%eax
jmp    165b <calculate_ngtdm+0x53b>
mov    %r14,%rdi
call   1030 <free@plt>
mov    $0x1,%eax
add    $0x78,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define i32 @calculate_ngtdm(i32* %0, i8* %1, i32* %2, i32* %3, i32* %4, i32* %5, i32 %6, i32 %7, double* %8, i32 %9) {
%11 = alloca i32
%12 = alloca i32*
%13 = alloca i8*
%14 = alloca i32*
%15 = alloca i32*
%16 = alloca i32*
%17 = alloca i32*
%18 = alloca i32
%19 = alloca i32
%20 = alloca double*
%21 = alloca i32
%22 = alloca i64
%23 = alloca i64
%24 = alloca i64
%25 = alloca i64*
%26 = alloca i64
%27 = alloca i64
%28 = alloca i32
%29 = alloca double
%30 = alloca double
%31 = alloca double
%32 = alloca i64
%33 = alloca i64
store i32* %0, i32** %12
store i8* %1, i8** %13
store i32* %2, i32** %14
store i32* %3, i32** %15
store i32* %4, i32** %16
store i32* %5, i32** %17
store i32 %6, i32* %18
store i32 %7, i32* %19
store double* %8, double** %20
store i32 %9, i32* %21
%34 = load i32, i32* %19
%35 = sext i32 %34 to i64
%36 = mul i64 8, %35
%37 = trunc i64 %36 to i32
%38 = call i64 @malloc(i32 %37)
%39 = inttoptr i64 %38 to i64*
store i64* %39, i64** %25
%40 = load i32, i32* %21
%41 = mul nsw i32 %40, 3
%42 = sext i32 %41 to i64
store i64 %42, i64* %33
store i32 0, i32* %28
br label %43
43: ; preds = %57, %10
%44 = load i32, i32* %28
%45 = load i32, i32* %21
%46 = icmp slt i32 %44, %45
br i1 %46, label %47, label %60
47: ; preds = %43
%48 = load i32, i32* %28
%49 = add nsw i32 %48, 1
%50 = sitofp i32 %49 to double
%51 = load double*, double** %20
%52 = load i32, i32* %28
%53 = mul nsw i32 %52, 3
%54 = add nsw i32 %53, 2
%55 = sext i32 %54 to i64
%56 = getelementptr inbounds double, double* %51, i64 %55
store double %50, double* %56
br label %57
57: ; preds = %47
%58 = load i32, i32* %28
%59 = add nsw i32 %58, 1
store i32 %59, i32* %28
br label %43
60: ; preds = %43
%61 = load i32*, i32** %14
%62 = getelementptr inbounds i32, i32* %61, i64 0
%63 = load i32, i32* %62
%64 = sext i32 %63 to i64
store i64 %64, i64* %22
%65 = load i32*, i32** %15
%66 = getelementptr inbounds i32, i32* %65, i64 0
%67 = load i32, i32* %66
%68 = load i32*, i32** %16
%69 = getelementptr inbounds i32, i32* %68, i64 0
%70 = load i32, i32* %69
%71 = mul nsw i32 %67, %70
%72 = sext i32 %71 to i64
store i64 %72, i64* %23
store i64 1, i64* %27
br label %73
73: ; preds = %98, %60
%74 = load i64, i64* %27
%75 = load i32, i32* %19
%76 = sext i32 %75 to i64
%77 = icmp ult i64 %74, %76
br i1 %77, label %78, label %101
78: ; preds = %73
%79 = load i32*, i32** %15
%80 = load i64, i64* %27
%81 = getelementptr inbounds i32, i32* %79, i64 %80
%82 = load i32, i32* %81
%83 = load i32*, i32** %16
%84 = load i64, i64* %27
%85 = getelementptr inbounds i32, i32* %83, i64 %84
%86 = load i32, i32* %85
%87 = mul nsw i32 %82, %86
%88 = sext i32 %87 to i64
%89 = load i64, i64* %23
%90 = add i64 %89, %88
store i64 %90, i64* %23
%91 = load i32*, i32** %14
%92 = load i64, i64* %27
%93 = getelementptr inbounds i32, i32* %91, i64 %92
%94 = load i32, i32* %93
%95 = sext i32 %94 to i64
%96 = load i64, i64* %22
%97 = mul i64 %96, %95
store i64 %97, i64* %22
br label %98
98: ; preds = %78
%99 = load i64, i64* %27
%100 = add i64 %99, 1
store i64 %100, i64* %27
br label %73
101: ; preds = %73
br label %102
102: ; preds = %408, %101
%103 = load i64, i64* %23
%104 = load i64, i64* %22
%105 = icmp ult i64 %103, %104
br i1 %105, label %106, label %411
106: ; preds = %102
%107 = load i32, i32* %19
%108 = sub nsw i32 %107, 1
%109 = sext i32 %108 to i64
store i64 %109, i64* %27
br label %110
110: ; preds = %217, %106
%111 = load i64, i64* %27
%112 = icmp ugt i64 %111, 0
br i1 %112, label %113, label %220
113: ; preds = %110
%114 = load i64, i64* %23
%115 = load i32*, i32** %16
%116 = load i64, i64* %27
%117 = sub i64 %116, 1
%118 = getelementptr inbounds i32, i32* %115, i64 %117
%119 = load i32, i32* %118
%120 = sext i32 %119 to i64
%121 = urem i64 %114, %120
%122 = load i32*, i32** %16
%123 = load i64, i64* %27
%124 = getelementptr inbounds i32, i32* %122, i64 %123
%125 = load i32, i32* %124
%126 = sext i32 %125 to i64
%127 = udiv i64 %121, %126
%128 = load i64*, i64** %25
%129 = load i64, i64* %27
%130 = getelementptr inbounds i64, i64* %128, i64 %129
store i64 %127, i64* %130
%131 = load i64*, i64** %25
%132 = load i64, i64* %27
%133 = getelementptr inbounds i64, i64* %131, i64 %132
%134 = load i64, i64* %133
%135 = load i32*, i32** %15
%136 = load i32, i32* %19
%137 = sext i32 %136 to i64
%138 = load i64, i64* %27
%139 = add i64 %137, %138
%140 = getelementptr inbounds i32, i32* %135, i64 %139
%141 = load i32, i32* %140
%142 = sext i32 %141 to i64
%143 = icmp ugt i64 %134, %142
br i1 %143, label %144, label %177
144: ; preds = %113
%145 = load i32*, i32** %14
%146 = load i64, i64* %27
%147 = getelementptr inbounds i32, i32* %145, i64 %146
%148 = load i32, i32* %147
%149 = sext i32 %148 to i64
%150 = load i64*, i64** %25
%151 = load i64, i64* %27
%152 = getelementptr inbounds i64, i64* %150, i64 %151
%153 = load i64, i64* %152
%154 = sub i64 %149, %153
%155 = load i32*, i32** %15
%156 = load i64, i64* %27
%157 = getelementptr inbounds i32, i32* %155, i64 %156
%158 = load i32, i32* %157
%159 = sext i32 %158 to i64
%160 = add i64 %154, %159
%161 = load i32*, i32** %16
%162 = load i64, i64* %27
%163 = getelementptr inbounds i32, i32* %161, i64 %162
%164 = load i32, i32* %163
%165 = sext i32 %164 to i64
%166 = mul i64 %160, %165
%167 = load i64, i64* %23
%168 = add i64 %167, %166
store i64 %168, i64* %23
%169 = load i32*, i32** %15
%170 = load i64, i64* %27
%171 = getelementptr inbounds i32, i32* %169, i64 %170
%172 = load i32, i32* %171
%173 = sext i32 %172 to i64
%174 = load i64*, i64** %25
%175 = load i64, i64* %27
%176 = getelementptr inbounds i64, i64* %174, i64 %175
store i64 %173, i64* %176
br label %216
177: ; preds = %113
%178 = load i64*, i64** %25
%179 = load i64, i64* %27
%180 = getelementptr inbounds i64, i64* %178, i64 %179
%181 = load i64, i64* %180
%182 = load i32*, i32** %15
%183 = load i64, i64* %27
%184 = getelementptr inbounds i32, i32* %182, i64 %183
%185 = load i32, i32* %184
%186 = sext i32 %185 to i64
%187 = icmp ult i64 %181, %186
br i1 %187, label %188, label %215
188: ; preds = %177
%189 = load i32*, i32** %15
%190 = load i64, i64* %27
%191 = getelementptr inbounds i32, i32* %189, i64 %190
%192 = load i32, i32* %191
%193 = sext i32 %192 to i64
%194 = load i64*, i64** %25
%195 = load i64, i64* %27
%196 = getelementptr inbounds i64, i64* %194, i64 %195
%197 = load i64, i64* %196
%198 = sub i64 %193, %197
%199 = load i32*, i32** %16
%200 = load i64, i64* %27
%201 = getelementptr inbounds i32, i32* %199, i64 %200
%202 = load i32, i32* %201
%203 = sext i32 %202 to i64
%204 = mul i64 %198, %203
%205 = load i64, i64* %23
%206 = add i64 %205, %204
store i64 %206, i64* %23
%207 = load i32*, i32** %15
%208 = load i64, i64* %27
%209 = getelementptr inbounds i32, i32* %207, i64 %208
%210 = load i32, i32* %209
%211 = sext i32 %210 to i64
%212 = load i64*, i64** %25
%213 = load i64, i64* %27
%214 = getelementptr inbounds i64, i64* %212, i64 %213
store i64 %211, i64* %214
br label %215
215: ; preds = %188, %177
br label %216
216: ; preds = %215, %144
br label %217
217: ; preds = %216
%218 = load i64, i64* %27
%219 = add i64 %218, -1
store i64 %219, i64* %27
br label %110
220: ; preds = %110
%221 = load i64, i64* %23
%222 = load i32*, i32** %16
%223 = getelementptr inbounds i32, i32* %222, i64 0
%224 = load i32, i32* %223
%225 = sext i32 %224 to i64
%226 = udiv i64 %221, %225
%227 = load i64*, i64** %25
%228 = getelementptr inbounds i64, i64* %227, i64 0
store i64 %226, i64* %228
%229 = load i64*, i64** %25
%230 = getelementptr inbounds i64, i64* %229, i64 0
%231 = load i64, i64* %230
%232 = load i32*, i32** %15
%233 = load i32, i32* %19
%234 = sext i32 %233 to i64
%235 = getelementptr inbounds i32, i32* %232, i64 %234
%236 = load i32, i32* %235
%237 = sext i32 %236 to i64
%238 = icmp ugt i64 %231, %237
br i1 %238, label %239, label %240
239: ; preds = %220
br label %411
240: ; preds = %220
%241 = load i8*, i8** %13
%242 = load i64, i64* %23
%243 = getelementptr inbounds i8, i8* %241, i64 %242
%244 = load i8, i8* %243
%245 = icmp ne i8 %244, 0
br i1 %245, label %246, label %407
246: ; preds = %240
store double 0.000000e+00, double* %29
store double 0.000000e+00, double* %30
store i64 0, i64* %26
br label %247
247: ; preds = %351, %246
%248 = load i64, i64* %26
%249 = load i32, i32* %18
%250 = sext i32 %249 to i64
%251 = icmp ult i64 %248, %250
br i1 %251, label %252, label %354
252: ; preds = %247
%253 = load i64, i64* %23
store i64 %253, i64* %24
store i64 0, i64* %27
br label %254
254: ; preds = %326, %252
%255 = load i64, i64* %27
%256 = load i32, i32* %19
%257 = sext i32 %256 to i64
%258 = icmp ult i64 %255, %257
br i1 %258, label %259, label %329
259: ; preds = %254
%260 = load i64*, i64** %25
%261 = load i64, i64* %27
%262 = getelementptr inbounds i64, i64* %260, i64 %261
%263 = load i64, i64* %262
%264 = load i32*, i32** %17
%265 = load i64, i64* %26
%266 = load i32, i32* %19
%267 = sext i32 %266 to i64
%268 = mul i64 %265, %267
%269 = load i64, i64* %27
%270 = add i64 %268, %269
%271 = getelementptr inbounds i32, i32* %264, i64 %270
%272 = load i32, i32* %271
%273 = sext i32 %272 to i64
%274 = add i64 %263, %273
%275 = load i32*, i32** %15
%276 = load i64, i64* %27
%277 = getelementptr inbounds i32, i32* %275, i64 %276
%278 = load i32, i32* %277
%279 = sext i32 %278 to i64
%280 = icmp ult i64 %274, %279
br i1 %280, label %306, label %281
281: ; preds = %259
%282 = load i64*, i64** %25
%283 = load i64, i64* %27
%284 = getelementptr inbounds i64, i64* %282, i64 %283
%285 = load i64, i64* %284
%286 = load i32*, i32** %17
%287 = load i64, i64* %26
%288 = load i32, i32* %19
%289 = sext i32 %288 to i64
%290 = mul i64 %287, %289
%291 = load i64, i64* %27
%292 = add i64 %290, %291
%293 = getelementptr inbounds i32, i32* %286, i64 %292
%294 = load i32, i32* %293
%295 = sext i32 %294 to i64
%296 = add i64 %285, %295
%297 = load i32*, i32** %15
%298 = load i32, i32* %19
%299 = sext i32 %298 to i64
%300 = load i64, i64* %27
%301 = add i64 %299, %300
%302 = getelementptr inbounds i32, i32* %297, i64 %301
%303 = load i32, i32* %302
%304 = sext i32 %303 to i64
%305 = icmp ugt i64 %296, %304
br i1 %305, label %306, label %308
306: ; preds = %281, %259
%307 = load i64, i64* %23
store i64 %307, i64* %24
br label %329
308: ; preds = %281
%309 = load i32*, i32** %17
%310 = load i64, i64* %26
%311 = load i32, i32* %19
%312 = sext i32 %311 to i64
%313 = mul i64 %310, %312
%314 = load i64, i64* %27
%315 = add i64 %313, %314
%316 = getelementptr inbounds i32, i32* %309, i64 %315
%317 = load i32, i32* %316
%318 = load i32*, i32** %16
%319 = load i64, i64* %27
%320 = getelementptr inbounds i32, i32* %318, i64 %319
%321 = load i32, i32* %320
%322 = mul nsw i32 %317, %321
%323 = sext i32 %322 to i64
%324 = load i64, i64* %24
%325 = add i64 %324, %323
store i64 %325, i64* %24
br label %326
326: ; preds = %308
%327 = load i64, i64* %27
%328 = add i64 %327, 1
store i64 %328, i64* %27
br label %254
329: ; preds = %306, %254
%330 = load i64, i64* %24
%331 = load i64, i64* %23
%332 = icmp ne i64 %330, %331
br i1 %332, label %333, label %350
333: ; preds = %329
%334 = load i8*, i8** %13
%335 = load i64, i64* %24
%336 = getelementptr inbounds i8, i8* %334, i64 %335
%337 = load i8, i8* %336
%338 = sext i8 %337 to i32
%339 = icmp ne i32 %338, 0
br i1 %339, label %340, label %350
340: ; preds = %333
%341 = load double, double* %29
%342 = fadd double %341, 1.000000e+00
store double %342, double* %29
%343 = load i32*, i32** %12
%344 = load i64, i64* %24
%345 = getelementptr inbounds i32, i32* %343, i64 %344
%346 = load i32, i32* %345
%347 = sitofp i32 %346 to double
%348 = load double, double* %30
%349 = fadd double %348, %347
store double %349, double* %30
br label %350
350: ; preds = %340, %333, %329
br label %351
351: ; preds = %350
%352 = load i64, i64* %26
%353 = add i64 %352, 1
store i64 %353, i64* %26
br label %247
354: ; preds = %247
%355 = load double, double* %29
%356 = fcmp oeq double %355, 0.000000e+00
br i1 %356, label %357, label %358
357: ; preds = %354
store double 0.000000e+00, double* %31
br label %368
358: ; preds = %354
%359 = load i32*, i32** %12
%360 = load i64, i64* %23
%361 = getelementptr inbounds i32, i32* %359, i64 %360
%362 = load i32, i32* %361
%363 = sitofp i32 %362 to double
%364 = load double, double* %30
%365 = load double, double* %29
%366 = fdiv double %364, %365
%367 = fsub double %363, %366
store double %367, double* %31
br label %368
368: ; preds = %358, %357
%369 = load double, double* %31
%370 = fcmp olt double %369, 0.000000e+00
br i1 %370, label %371, label %374
371: ; preds = %368
%372 = load double, double* %31
%373 = fmul double %372, -1.000000e+00
store double %373, double* %31
br label %374
374: ; preds = %371, %368
%375 = load i32*, i32** %12
%376 = load i64, i64* %23
%377 = getelementptr inbounds i32, i32* %375, i64 %376
%378 = load i32, i32* %377
%379 = sub nsw i32 %378, 1
%380 = mul nsw i32 %379, 3
%381 = sext i32 %380 to i64
store i64 %381, i64* %32
%382 = load i32*, i32** %12
%383 = load i64, i64* %23
%384 = getelementptr inbounds i32, i32* %382, i64 %383
%385 = load i32, i32* %384
%386 = icmp sle i32 %385, 0
br i1 %386, label %391, label %387
387: ; preds = %374
%388 = load i64, i64* %32
%389 = load i64, i64* %33
%390 = icmp uge i64 %388, %389
br i1 %390, label %391, label %394
391: ; preds = %387, %374
%392 = load i64*, i64** %25
%393 = call i32 @free(i64* %392)
store i32 0, i32* %11
br label %414
394: ; preds = %387
%395 = load double*, double** %20
%396 = load i64, i64* %32
%397 = getelementptr inbounds double, double* %395, i64 %396
%398 = load double, double* %397
%399 = fadd double %398, 1.000000e+00
store double %399, double* %397
%400 = load double, double* %31
%401 = load double*, double** %20
%402 = load i64, i64* %32
%403 = add i64 %402, 1
%404 = getelementptr inbounds double, double* %401, i64 %403
%405 = load double, double* %404
%406 = fadd double %405, %400
store double %406, double* %404
br label %407
407: ; preds = %394, %240
br label %408
408: ; preds = %407
%409 = load i64, i64* %23
%410 = add i64 %409, 1
store i64 %410, i64* %23
br label %102
411: ; preds = %239, %102
%412 = load i64*, i64** %25
%413 = call i32 @free(i64* %412)
store i32 1, i32* %11
br label %414
414: ; preds = %411, %391
%415 = load i32, i32* %11
ret i32 %415
}

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