#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 int32_t;

/* Variables and functions */
float aaf(int *, size_t);
int related(int *gt_types, double *asum, int32_t *N, int32_t *ibs0,
            int32_t *ibs2, int32_t n_samples) {
  int idx, uidx, n_used = 0;
  int32_t j, k;
  float pi = aaf(gt_types, n_samples);
  float numer, val;
  float gtj, gtk;
  float denom = 2.0 * pi * (1.0 - pi);
  if (denom == 0) {
    return 0;
  }
  for (j = 0; j < n_samples; j++) {
    if (gt_types[j] == 3) {
      continue;
    }
    gtj = gt_types[j];
    n_used++;
    for (k = j; k < n_samples; k++) {
      if (gt_types[k] == 3) {
        continue;
      }
      uidx = j + k * n_samples;
      idx = j * n_samples + k;
      gtk = gt_types[k];
      if (j != k) {
        numer = (gtj - 2.0 * pi) * (gtk - 2.0 * pi);
        ibs0[idx] += (gtj != 1 && gtk != 1 && gtj != gtk);
      } else {
        numer = (gtj * gtj) - (1.0 + 2.0 * pi) * gtj + 2.0 * pi * pi;
        asum[idx] += 1;
      }
      val = numer / denom;
      if (gtj == gtk && gtj != 0 && val > 2.5) {
        ibs2[uidx] += 1;
      } else if (val > 2.5) {
        ibs2[idx] += (gtj == gtk && gtk != 1);
      }
      asum[idx] += val;
      N[idx] += 1;
    }
  }
  return n_used;
}


x86 assembly (clang, -O0):
<related>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
mov    %r8,-0x30(%rbp)
mov    %r9,-0x38(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x10(%rbp),%rdi
mov    -0x38(%rbp),%rsi
call   1030 <aaf@plt>
movss  %xmm0,-0x5c(%rbp)
movss  -0x5c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xeab(%rip),%xmm0
mulsd  %xmm1,%xmm0
movss  -0x5c(%rbp),%xmm1
cvtss2sd %xmm1,%xmm2
movsd  0xe8e(%rip),%xmm1
subsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x70(%rbp)
movss  -0x70(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jne    11a6 <related+0x96>
jp     11a6 <related+0x96>
movl   $0x0,-0x4(%rbp)
jmp    14f7 <related+0x3e7>
movq   $0x0,-0x50(%rbp)
mov    -0x50(%rbp),%rax
cmp    -0x38(%rbp),%rax
jae    14f1 <related+0x3e1>
mov    -0x10(%rbp),%rax
mov    -0x50(%rbp),%rcx
cmpl   $0x3,(%rax,%rcx,4)
jne    11d3 <related+0xc3>
jmp    14e0 <related+0x3d0>
mov    -0x10(%rbp),%rax
mov    -0x50(%rbp),%rcx
cvtsi2ssl (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x68(%rbp)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
mov    -0x50(%rbp),%rax
mov    %rax,-0x58(%rbp)
mov    -0x58(%rbp),%rax
cmp    -0x38(%rbp),%rax
jae    14db <related+0x3cb>
mov    -0x10(%rbp),%rax
mov    -0x58(%rbp),%rcx
cmpl   $0x3,(%rax,%rcx,4)
jne    121b <related+0x10b>
jmp    14ca <related+0x3ba>
mov    -0x50(%rbp),%rax
mov    -0x58(%rbp),%rcx
imul   -0x38(%rbp),%rcx
add    %rcx,%rax
mov    %eax,-0x40(%rbp)
mov    -0x50(%rbp),%rax
imul   -0x38(%rbp),%rax
add    -0x58(%rbp),%rax
mov    %eax,-0x3c(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x58(%rbp),%rcx
cvtsi2ssl (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x6c(%rbp)
mov    -0x50(%rbp),%rax
cmp    -0x58(%rbp),%rax
je     131f <related+0x20f>
movss  -0x68(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movss  -0x5c(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
addsd  %xmm1,%xmm1
subsd  %xmm1,%xmm0
movss  -0x6c(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
movss  -0x5c(%rbp),%xmm2
cvtss2sd %xmm2,%xmm2
addsd  %xmm2,%xmm2
subsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x60(%rbp)
movss  -0x68(%rbp),%xmm0
xor    %eax,%eax
movss  0xd6a(%rip),%xmm1
ucomiss %xmm1,%xmm0
mov    %al,-0x71(%rbp)
jne    12c5 <related+0x1b5>
jp     12c5 <related+0x1b5>
jmp    12ff <related+0x1ef>
movss  -0x6c(%rbp),%xmm0
xor    %eax,%eax
movss  0xd44(%rip),%xmm1
ucomiss %xmm1,%xmm0
mov    %al,-0x71(%rbp)
jne    12eb <related+0x1db>
jp     12eb <related+0x1db>
jmp    12ff <related+0x1ef>
movss  -0x68(%rbp),%xmm0
ucomiss -0x6c(%rbp),%xmm0
setne  %al
setp   %cl
or     %cl,%al
mov    %al,-0x71(%rbp)
mov    -0x71(%rbp),%al
and    $0x1,%al
movzbl %al,%eax
movslq %eax,%rdx
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
add    (%rax,%rcx,8),%rdx
mov    %rdx,(%rax,%rcx,8)
jmp    13b7 <related+0x2a7>
movss  -0x68(%rbp),%xmm0
mulss  -0x68(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movss  -0x5c(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
addsd  %xmm1,%xmm1
movsd  0xcbe(%rip),%xmm2
addsd  %xmm2,%xmm1
movss  -0x68(%rbp),%xmm2
cvtss2sd %xmm2,%xmm2
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movss  -0x5c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm2
movsd  0xc89(%rip),%xmm0
mulsd  %xmm2,%xmm0
movss  -0x5c(%rbp),%xmm2
cvtss2sd %xmm2,%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x60(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  0xc53(%rip),%xmm0
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movss  -0x60(%rbp),%xmm0
divss  -0x70(%rbp),%xmm0
movss  %xmm0,-0x64(%rbp)
movss  -0x68(%rbp),%xmm0
ucomiss -0x6c(%rbp),%xmm0
jne    142b <related+0x31b>
jp     142b <related+0x31b>
movss  -0x68(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jne    13f7 <related+0x2e7>
jp     13f7 <related+0x2e7>
jmp    142b <related+0x31b>
movss  -0x64(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xc08(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    142b <related+0x31b>
mov    -0x30(%rbp),%rax
movslq -0x40(%rbp),%rcx
mov    (%rax,%rcx,8),%rdx
add    $0x1,%rdx
mov    %rdx,(%rax,%rcx,8)
jmp    149b <related+0x38b>
movss  -0x64(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xbd4(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1496 <related+0x386>
movss  -0x68(%rbp),%xmm0
xor    %eax,%eax
ucomiss -0x6c(%rbp),%xmm0
mov    %al,-0x72(%rbp)
jne    147b <related+0x36b>
jp     147b <related+0x36b>
movss  -0x6c(%rbp),%xmm0
movss  0xbab(%rip),%xmm1
ucomiss %xmm1,%xmm0
setne  %al
setp   %cl
or     %cl,%al
mov    %al,-0x72(%rbp)
mov    -0x72(%rbp),%al
and    $0x1,%al
movzbl %al,%eax
movslq %eax,%rdx
mov    -0x30(%rbp),%rax
movslq -0x3c(%rbp),%rcx
add    (%rax,%rcx,8),%rdx
mov    %rdx,(%rax,%rcx,8)
jmp    149b <related+0x38b>
movss  -0x64(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mov    (%rax,%rcx,8),%rdx
add    $0x1,%rdx
mov    %rdx,(%rax,%rcx,8)
mov    -0x58(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x58(%rbp)
jmp    11f6 <related+0xe6>
jmp    14e0 <related+0x3d0>
mov    -0x50(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x50(%rbp)
jmp    11ae <related+0x9e>
mov    -0x44(%rbp),%eax
mov    %eax,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<related>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
mov    %r9,%r13
mov    %r8,0x38(%rsp)
mov    %rcx,0x30(%rsp)
mov    %rdx,0x40(%rsp)
mov    %rsi,%r12
mov    %rdi,%rbx
mov    %r9,%rsi
call   1030 <aaf@plt>
cvtss2sd %xmm0,%xmm9
movaps %xmm9,%xmm1
addsd  %xmm9,%xmm1
movsd  0xeab(%rip),%xmm10
movapd %xmm10,%xmm0
subsd  %xmm9,%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm3
xor    %eax,%eax
xorps  %xmm0,%xmm0
ucomiss %xmm0,%xmm3
jne    1177 <related+0x67>
jnp    13b6 <related+0x2a6>
test   %r13,%r13
je     13b6 <related+0x2a6>
movsd  0xe77(%rip),%xmm8
addsd  %xmm1,%xmm8
xorpd  0xe79(%rip),%xmm8
movabs $0x100000000,%rax
lea    0x1(%r13),%rcx
mov    %r13,%rdx
shl    $0x20,%rdx
add    %rax,%rdx
lea    0x4(%rbx),%rdi
xor    %r10d,%r10d
mov    $0x1,%r11d
mulsd  %xmm1,%xmm9
movss  0xe57(%rip),%xmm5
xor    %esi,%esi
xor    %eax,%eax
xor    %r8d,%r8d
mov    %rbx,0x20(%rsp)
mov    %rcx,0x18(%rsp)
mov    %rdx,0x10(%rsp)
jmp    1223 <related+0x113>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    0x20(%rsp),%rbx
mov    0x18(%rsp),%rcx
mov    0x10(%rsp),%rdx
mov    0xc(%rsp),%eax
mov    0x28(%rsp),%r10
add    $0x1,%r8
add    %rcx,%rsi
add    %rdx,%r10
add    $0x4,%rdi
add    $0x1,%r11
cmp    %r13,%r8
je     13b6 <related+0x2a6>
mov    (%rbx,%r8,4),%r14d
cmp    $0x3,%r14d
je     1208 <related+0xf8>
xorps  %xmm6,%xmm6
cvtsi2ss %r14d,%xmm6
add    $0x1,%eax
mov    %eax,0xc(%rsp)
xorps  %xmm7,%xmm7
cvtss2sd %xmm6,%xmm7
movaps %xmm6,%xmm0
mulss  %xmm6,%xmm0
cvtss2sd %xmm0,%xmm0
movapd %xmm8,%xmm4
mulsd  %xmm7,%xmm4
subsd  %xmm1,%xmm7
addsd  %xmm0,%xmm4
addsd  %xmm9,%xmm4
mov    %r10,0x28(%rsp)
xor    %r9d,%r9d
mov    %r14d,%ecx
mov    %r8,%r15
movabs $0x100000000,%rbx
cmp    $0x3,%ecx
je     138f <related+0x27f>
cs nopw 0x0(%rax,%rax,1)
nop
xorps  %xmm0,%xmm0
cvtsi2ss %ecx,%xmm0
test   %r9,%r9
je     12f0 <related+0x1e0>
cmp    $0x1,%r14d
setne  %al
xorps  %xmm2,%xmm2
cvtss2sd %xmm0,%xmm2
subsd  %xmm1,%xmm2
mulsd  %xmm7,%xmm2
cmp    $0x1,%ecx
setne  %bl
and    %al,%bl
ucomiss %xmm0,%xmm6
setp   %al
setne  %dl
or     %al,%dl
and    %bl,%dl
movzbl %dl,%eax
mov    %r10,%rdx
sar    $0x1d,%rdx
mov    0x30(%rsp),%rbx
add    %rax,(%rbx,%rdx,1)
movabs $0x100000000,%rbx
jmp    130c <related+0x1fc>
cs nopw 0x0(%rax,%rax,1)
mov    %r10,%rax
sar    $0x1d,%rax
movsd  (%r12,%rax,1),%xmm2
addsd  %xmm10,%xmm2
movsd  %xmm2,(%r12,%rax,1)
movapd %xmm4,%xmm2
cvtsd2ss %xmm2,%xmm2
divss  %xmm3,%xmm2
test   %r14d,%r14d
je     1340 <related+0x230>
ucomiss %xmm0,%xmm6
jne    1340 <related+0x230>
jp     1340 <related+0x230>
ucomiss %xmm5,%xmm2
jbe    1340 <related+0x230>
mov    %r15,%rbp
imul   %r13,%rbp
add    %r8,%rbp
mov    $0x1,%ecx
jmp    135f <related+0x24f>
cs nopw 0x0(%rax,%rax,1)
ucomiss %xmm5,%xmm2
jbe    136b <related+0x25b>
ucomiss %xmm0,%xmm6
lea    (%rsi,%r9,1),%rbp
setnp  %al
sete   %dl
and    %al,%dl
cmp    $0x1,%ecx
setne  %al
and    %dl,%al
movzbl %al,%ecx
movslq %ebp,%rax
mov    0x38(%rsp),%rdx
add    %rcx,(%rdx,%rax,8)
xorps  %xmm0,%xmm0
cvtss2sd %xmm2,%xmm0
mov    %r10,%rax
sar    $0x1d,%rax
addsd  (%r12,%rax,1),%xmm0
movsd  %xmm0,(%r12,%rax,1)
mov    0x40(%rsp),%rcx
addq   $0x1,(%rcx,%rax,1)
lea    (%r11,%r9,1),%rax
cmp    %r13,%rax
jae    11f0 <related+0xe0>
mov    (%rdi,%r9,4),%ecx
add    $0x1,%r9
add    %rbx,%r10
add    $0x1,%r15
cmp    $0x3,%ecx
jne    1290 <related+0x180>
jmp    138f <related+0x27f>
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define i32 @related(i32* %0, double* nocapture %1, i64* nocapture %2, i64* nocapture %3, i64* nocapture %4, i64 %5) {
%7 = tail call float @aaf(i32* %0, i64 %5)
%8 = fpext float %7 to double
%9 = fmul double %8, 2.000000e+00
%10 = fsub double 1.000000e+00, %8
%11 = fmul double %9, %10
%12 = fptrunc double %11 to float
%13 = fcmp une float %12, 0.000000e+00
%14 = icmp ne i64 %5, 0
%15 = and i1 %13, %14
br i1 %15, label %16, label %108
16: ; preds = %6
%17 = tail call double @llvm.fmuladd.f64(double %8, double 2.000000e+00, double 1.000000e+00)
%18 = fneg double %17
br label %19
19: ; preds = %16, %104
%20 = phi i32 [ 0, %16 ], [ %105, %104 ]
%21 = phi i64 [ 0, %16 ], [ %106, %104 ]
%22 = getelementptr inbounds i32, i32* %0, i64 %21
%23 = load i32, i32* %22
%24 = icmp eq i32 %23, 3
br i1 %24, label %104, label %25
25: ; preds = %19
%26 = sitofp i32 %23 to float
%27 = add nsw i32 %20, 1
%28 = mul i64 %21, %5
%29 = fpext float %26 to double
%30 = tail call double @llvm.fmuladd.f64(double %8, double -2.000000e+00, double %29)
%31 = icmp ne i32 %23, 1
%32 = fmul float %26, %26
%33 = fpext float %32 to double
%34 = tail call double @llvm.fmuladd.f64(double %18, double %29, double %33)
%35 = tail call double @llvm.fmuladd.f64(double %9, double %8, double %34)
%36 = icmp ne i32 %23, 0
br label %37
37: ; preds = %101, %25
%38 = phi i32 [ %23, %25 ], [ %103, %101 ]
%39 = phi i64 [ %21, %25 ], [ %99, %101 ]
%40 = icmp eq i32 %38, 3
br i1 %40, label %98, label %41
41: ; preds = %37
%42 = mul i64 %39, %5
%43 = add i64 %42, %21
%44 = add i64 %39, %28
%45 = sitofp i32 %38 to float
%46 = icmp eq i64 %21, %39
br i1 %46, label %61, label %47
47: ; preds = %41
%48 = fpext float %45 to double
%49 = tail call double @llvm.fmuladd.f64(double %8, double -2.000000e+00, double %48)
%50 = fmul double %30, %49
%51 = icmp ne i32 %38, 1
%52 = select i1 %31, i1 %51, i1 false
%53 = fcmp une float %26, %45
%54 = select i1 %52, i1 %53, i1 false
%55 = zext i1 %54 to i64
%56 = shl i64 %44, 32
%57 = ashr exact i64 %56, 32
%58 = getelementptr inbounds i64, i64* %3, i64 %57
%59 = load i64, i64* %58
%60 = add i64 %59, %55
store i64 %60, i64* %58
br label %67
61: ; preds = %41
%62 = shl i64 %44, 32
%63 = ashr exact i64 %62, 32
%64 = getelementptr inbounds double, double* %1, i64 %63
%65 = load double, double* %64
%66 = fadd double %65, 1.000000e+00
store double %66, double* %64
br label %67
67: ; preds = %61, %47
%68 = phi double [ %50, %47 ], [ %35, %61 ]
%69 = fptrunc double %68 to float
%70 = fdiv float %69, %12
%71 = fcmp oeq float %26, %45
%72 = and i1 %36, %71
%73 = fcmp ogt float %70, 2.500000e+00
%74 = select i1 %72, i1 %73, i1 false
br i1 %74, label %80, label %75
75: ; preds = %67
br i1 %73, label %76, label %88
76: ; preds = %75
%77 = icmp ne i32 %38, 1
%78 = and i1 %77, %71
%79 = zext i1 %78 to i64
br label %80
80: ; preds = %67, %76
%81 = phi i64 [ %44, %76 ], [ %43, %67 ]
%82 = phi i64 [ %79, %76 ], [ 1, %67 ]
%83 = shl i64 %81, 32
%84 = ashr exact i64 %83, 32
%85 = getelementptr inbounds i64, i64* %4, i64 %84
%86 = load i64, i64* %85
%87 = add i64 %86, %82
store i64 %87, i64* %85
br label %88
88: ; preds = %80, %75
%89 = fpext float %70 to double
%90 = shl i64 %44, 32
%91 = ashr exact i64 %90, 32
%92 = getelementptr inbounds double, double* %1, i64 %91
%93 = load double, double* %92
%94 = fadd double %93, %89
store double %94, double* %92
%95 = getelementptr inbounds i64, i64* %2, i64 %91
%96 = load i64, i64* %95
%97 = add i64 %96, 1
store i64 %97, i64* %95
br label %98
98: ; preds = %37, %88
%99 = add nuw i64 %39, 1
%100 = icmp ult i64 %99, %5
br i1 %100, label %101, label %104
101: ; preds = %98
%102 = getelementptr inbounds i32, i32* %0, i64 %99
%103 = load i32, i32* %102
br label %37
104: ; preds = %98, %19
%105 = phi i32 [ %20, %19 ], [ %27, %98 ]
%106 = add nuw i64 %21, 1
%107 = icmp eq i64 %106, %5
br i1 %107, label %108, label %19
108: ; preds = %104, %6
%109 = phi i32 [ 0, %6 ], [ %105, %104 ]
ret i32 %109
}