#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
double INFINITY;
double pow(double, int);
double sqrt(double);
void interceptCyl(double s[], double pt0[][3], double norm[][3],
                  double outliner[], double outlinez[], int ix, int jx) {
  int i, j;
  double A, A0, B, B0, C, C0, delr, delz, temp, s1, s2, stemp, res = 1e-6;
  for (j = jx; j--;) {
    stemp = INFINITY;
    C0 = pow(pt0[j][0], 2) + pow(pt0[j][1], 2);
    B0 = 2 * (pt0[j][0] * norm[j][0] + pt0[j][1] * norm[j][1]);
    A0 = pow(norm[j][0], 2) + pow(norm[j][1], 2);
    for (i = ix; i--;) {
      A = 0;
      B = 0;
      C = 0;
      delr = outliner[i + 1] - outliner[i];
      delz = outlinez[i + 1] - outlinez[i];
      if (delz) {
        temp = outliner[i] + (delr / delz) * (pt0[j][2] - outlinez[i]);
        A = A0 - pow((delr / delz) * norm[j][2], 2);
        B = B0 - 2 * temp * (delr / delz) * norm[j][2];
        C = C0 - pow(temp, 2);
      } else if ((norm[j][2] != 0) & (delr != 0)) {
        s1 = (outlinez[i] - pt0[j][2]) / norm[j][2];
        temp = (sqrt(pow(pt0[j][0] + norm[j][0] * s1, 2) +
                     pow(pt0[j][1] + norm[j][1] * s1, 2)) -
                outliner[i]) /
               delr;
        if ((temp > 0) & (temp <= 1) & (s1 < stemp) & (s1 > res)) {
          stemp = s1;
        }
      }
      if (A) {
        temp = B * B - 4 * A * C;
        if (temp >= 0) {
          temp = sqrt(temp);
          s1 = -.5 * (temp + B) / A;
          s2 = .5 * (temp - B) / A;
          temp = (s1 * norm[j][2] + (pt0[j][2] - outlinez[i])) / delz;
          if ((temp > 0) & (temp <= 1) & (s1 < stemp) & (s1 > res)) {
            stemp = s1;
          }
          temp = (s2 * norm[j][2] + (pt0[j][2] - outlinez[i])) / delz;
          if ((temp > 0) & (temp <= 1) & (s2 < stemp) & (s2 > res)) {
            stemp = s2;
          }
        }
      } else if ((B != 0) & (C != 0)) {
        s1 = -B / C;
        temp = (s1 * norm[j][2] + (pt0[j][2] - outlinez[i])) / delz;
        if ((temp > 0) & (temp <= 1) & (s1 < stemp) & (s1 > res)) {
          stemp = s1;
        }
      }
    }
    s[j] = stemp;
  }
}


x86 assembly (clang, -O0):
<interceptCyl>:
push   %rbp
mov    %rsp,%rbp
sub    $0xd0,%rsp
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9d,-0x2c(%rbp)
movsd  0xeb2(%rip),%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
mov    0x10(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    %eax,%ecx
add    $0xffffffff,%ecx
mov    %ecx,-0x34(%rbp)
cmp    $0x0,%eax
je     19da <interceptCyl+0x8ba>
mov    0x2e71(%rip),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movsd  (%rax),%xmm0
mov    $0x2,%edi
call   1030 <pow@plt>
movsd  %xmm0,-0xb0(%rbp)
ff
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movsd  0x8(%rax),%xmm0
mov    $0x2,%edi
call   1030 <pow@plt>
movaps %xmm0,%xmm1
movsd  -0xb0(%rbp),%xmm0
ff
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movsd  (%rax),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movsd  (%rax),%xmm2
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movsd  0x8(%rax),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
mulsd  0x8(%rax),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  0xdd0(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movsd  (%rax),%xmm0
mov    $0x2,%edi
call   1030 <pow@plt>
movsd  %xmm0,-0xa8(%rbp)
ff
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movsd  0x8(%rax),%xmm0
mov    $0x2,%edi
call   1030 <pow@plt>
movaps %xmm0,%xmm1
movsd  -0xa8(%rbp),%xmm0
ff
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    %eax,%ecx
add    $0xffffffff,%ecx
mov    %ecx,-0x30(%rbp)
cmp    $0x0,%eax
je     19c0 <interceptCyl+0x8a0>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x30(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x28(%rbp),%rax
mov    -0x30(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x78(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    132f <interceptCyl+0x20f>
jp     132f <interceptCyl+0x20f>
jmp    144a <interceptCyl+0x32a>
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x70(%rbp),%xmm0
divsd  -0x78(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movsd  0x10(%rax),%xmm2
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0xc0(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
divsd  -0x78(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
mulsd  0x10(%rax),%xmm0
mov    $0x2,%edi
call   1030 <pow@plt>
movaps %xmm0,%xmm1
movsd  -0xc0(%rbp),%xmm0
ff
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x58(%rbp),%xmm1
movsd  0xc3e(%rip),%xmm0
mulsd  -0x80(%rbp),%xmm0
movsd  -0x70(%rbp),%xmm2
divsd  -0x78(%rbp),%xmm2
mulsd  %xmm2,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movsd  0x10(%rax),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0x80(%rbp),%xmm0
mov    $0x2,%edi
call   1030 <pow@plt>
movaps %xmm0,%xmm1
movsd  -0xb8(%rbp),%xmm0
ff
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
jmp    1612 <interceptCyl+0x4f2>
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movsd  0x10(%rax),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
setne  %al
setp   %cl
or     %cl,%al
and    $0x1,%al
movzbl %al,%eax
movsd  -0x70(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
setne  %cl
setp   %dl
or     %dl,%cl
and    $0x1,%cl
movzbl %cl,%ecx
and    %ecx,%eax
cmp    $0x0,%eax
je     160d <interceptCyl+0x4ed>
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
subsd  0x10(%rax),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
divsd  0x10(%rax),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movsd  (%rax),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movsd  (%rax),%xmm0
movsd  -0x88(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    $0x2,%edi
call   1030 <pow@plt>
movsd  %xmm0,-0xc8(%rbp)
ff
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movsd  0x8(%rax),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movsd  0x8(%rax),%xmm0
movsd  -0x88(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    $0x2,%edi
call   1030 <pow@plt>
movaps %xmm0,%xmm1
movsd  -0xc8(%rbp),%xmm0
ff
addsd  %xmm1,%xmm0
call   1040 <sqrt@plt>
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
divsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x80(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
movsd  0xa6b(%rip),%xmm0
ucomisd -0x80(%rbp),%xmm0
setae  %cl
and    $0x1,%cl
movzbl %cl,%ecx
and    %ecx,%eax
movsd  -0x88(%rbp),%xmm1
ff
movsd  -0x98(%rbp),%xmm0
ff
ucomisd %xmm1,%xmm0
seta   %cl
and    $0x1,%cl
movzbl %cl,%ecx
and    %ecx,%eax
movsd  -0x88(%rbp),%xmm0
ff
ucomisd -0xa0(%rbp),%xmm0
ff
seta   %cl
and    $0x1,%cl
movzbl %cl,%ecx
and    %ecx,%eax
cmp    $0x0,%eax
je     1608 <interceptCyl+0x4e8>
movsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
jmp    160d <interceptCyl+0x4ed>
jmp    1612 <interceptCyl+0x4f2>
movsd  -0x40(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    162f <interceptCyl+0x50f>
jp     162f <interceptCyl+0x50f>
jmp    187c <interceptCyl+0x75c>
movsd  -0x50(%rbp),%xmm0
movsd  -0x50(%rbp),%xmm2
movsd  0x9d7(%rip),%xmm1
mulsd  -0x40(%rbp),%xmm1
mulsd  -0x60(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x80(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     1877 <interceptCyl+0x757>
movsd  -0x80(%rbp),%xmm0
call   1040 <sqrt@plt>
movsd  %xmm0,-0x80(%rbp)
movsd  -0x80(%rbp),%xmm1
addsd  -0x50(%rbp),%xmm1
movsd  0x986(%rip),%xmm0
mulsd  %xmm1,%xmm0
divsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x80(%rbp),%xmm1
subsd  -0x50(%rbp),%xmm1
movsd  0x95b(%rip),%xmm0
mulsd  %xmm1,%xmm0
divsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x88(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movsd  0x10(%rax),%xmm2
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movsd  0x10(%rax),%xmm1
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
divsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x80(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
movsd  0x8cf(%rip),%xmm0
ucomisd -0x80(%rbp),%xmm0
setae  %cl
and    $0x1,%cl
movzbl %cl,%ecx
and    %ecx,%eax
movsd  -0x88(%rbp),%xmm1
ff
movsd  -0x98(%rbp),%xmm0
ff
ucomisd %xmm1,%xmm0
seta   %cl
and    $0x1,%cl
movzbl %cl,%ecx
and    %ecx,%eax
movsd  -0x88(%rbp),%xmm0
ff
ucomisd -0xa0(%rbp),%xmm0
ff
seta   %cl
and    $0x1,%cl
movzbl %cl,%ecx
and    %ecx,%eax
cmp    $0x0,%eax
je     17a4 <interceptCyl+0x684>
movsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x90(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movsd  0x10(%rax),%xmm2
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movsd  0x10(%rax),%xmm1
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
divsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x80(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
movsd  0x801(%rip),%xmm0
ucomisd -0x80(%rbp),%xmm0
setae  %cl
and    $0x1,%cl
movzbl %cl,%ecx
and    %ecx,%eax
movsd  -0x90(%rbp),%xmm1
ff
movsd  -0x98(%rbp),%xmm0
ff
ucomisd %xmm1,%xmm0
seta   %cl
and    $0x1,%cl
movzbl %cl,%ecx
and    %ecx,%eax
movsd  -0x90(%rbp),%xmm0
ff
ucomisd -0xa0(%rbp),%xmm0
ff
seta   %cl
and    $0x1,%cl
movzbl %cl,%ecx
and    %ecx,%eax
cmp    $0x0,%eax
je     1872 <interceptCyl+0x752>
movsd  -0x90(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
jmp    1877 <interceptCyl+0x757>
jmp    19bb <interceptCyl+0x89b>
movsd  -0x50(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
setne  %al
setp   %cl
or     %cl,%al
and    $0x1,%al
movzbl %al,%eax
movsd  -0x60(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
setne  %cl
setp   %dl
or     %dl,%cl
and    $0x1,%cl
movzbl %cl,%ecx
and    %ecx,%eax
cmp    $0x0,%eax
je     19b6 <interceptCyl+0x896>
movsd  -0x50(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
divsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x88(%rbp),%xmm0
ff
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movsd  0x10(%rax),%xmm2
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movsd  0x10(%rax),%xmm1
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
divsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x80(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
movsd  0x6c2(%rip),%xmm0
ucomisd -0x80(%rbp),%xmm0
setae  %cl
and    $0x1,%cl
movzbl %cl,%ecx
and    %ecx,%eax
movsd  -0x88(%rbp),%xmm1
ff
movsd  -0x98(%rbp),%xmm0
ff
ucomisd %xmm1,%xmm0
seta   %cl
and    $0x1,%cl
movzbl %cl,%ecx
and    %ecx,%eax
movsd  -0x88(%rbp),%xmm0
ff
ucomisd -0xa0(%rbp),%xmm0
ff
seta   %cl
and    $0x1,%cl
movzbl %cl,%ecx
and    %ecx,%eax
cmp    $0x0,%eax
je     19b1 <interceptCyl+0x891>
movsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
jmp    19b6 <interceptCyl+0x896>
jmp    19bb <interceptCyl+0x89b>
jmp    129e <interceptCyl+0x17e>
movsd  -0x98(%rbp),%xmm0
ff
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    115c <interceptCyl+0x3c>
add    $0xd0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<interceptCyl>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xb8,%rsp
mov    %r9d,0xc(%rsp)
mov    %rdi,0x88(%rsp)
mov    0xf0(%rsp),%eax
test   %eax,%eax
je     1733 <interceptCyl+0x613>
mov    %rdx,%r12
mov    %rsi,%r13
movslq 0xc(%rsp),%rdx
movslq %eax,%rbx
lea    (%rcx,%rdx,8),%r15
add    $0xfffffffffffffff8,%r15
lea    (%r8,%rdx,8),%r14
add    $0xfffffffffffffff8,%r14
neg    %edx
mov    %rdx,0x70(%rsp)
mov    %r12,0x80(%rsp)
mov    %rsi,0x78(%rsp)
jmp    11c0 <interceptCyl+0xa0>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    0x88(%rsp),%rax
mov    0x90(%rsp),%rbx
movapd 0x20(%rsp),%xmm0
movsd  %xmm0,(%rax,%rbx,8)
test   %ebx,%ebx
mov    0x80(%rsp),%r12
mov    0x78(%rsp),%r13
je     1733 <interceptCyl+0x613>
mov    0x2e21(%rip),%rax
movsd  (%rax),%xmm0
movaps %xmm0,0x20(%rsp)
lea    -0x8(,%rbx,8),%rax
ff
add    $0xffffffffffffffff,%rbx
lea    (%rax,%rax,2),%rbp
movsd  0x0(%r13,%rbp,1),%xmm0
mov    $0x2,%edi
call   1030 <pow@plt>
movsd  %xmm0,0x50(%rsp)
movsd  0x8(%r13,%rbp,1),%xmm0
mov    $0x2,%edi
call   1030 <pow@plt>
movsd  %xmm0,0x30(%rsp)
movsd  (%r12,%rbp,1),%xmm0
movsd  0x8(%r13,%rbp,1),%xmm1
mulsd  0x8(%r12,%rbp,1),%xmm1
movsd  %xmm1,0x10(%rsp)
movsd  0x0(%r13,%rbp,1),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x48(%rsp)
mov    $0x2,%edi
call   1030 <pow@plt>
movsd  %xmm0,0x40(%rsp)
movsd  0x8(%r12,%rbp,1),%xmm0
mov    $0x2,%edi
call   1030 <pow@plt>
cmpl   $0x0,0xc(%rsp)
mov    %rbx,0x90(%rsp)
je     1190 <interceptCyl+0x70>
mov    %r13,%rax
add    %rbp,%rax
mov    %rax,0xb0(%rsp)
lea    (%r12,%rbp,1),%rax
mov    %rax,0xa8(%rsp)
lea    0x8(%r13,%rbp,1),%rax
mov    %rax,0xa0(%rsp)
movsd  0x50(%rsp),%xmm1
addsd  0x30(%rsp),%xmm1
movsd  %xmm1,0x50(%rsp)
lea    0x8(%r12,%rbp,1),%rax
mov    %rax,0x98(%rsp)
movsd  0x48(%rsp),%xmm5
addsd  0x10(%rsp),%xmm5
addsd  %xmm5,%xmm5
movsd  0x40(%rsp),%xmm4
addsd  %xmm0,%xmm4
lea    0x0(,%rbx,8),%rax
lea    (%rax,%rax,2),%rax
add    %rax,%r12
add    $0x10,%r12
add    %rax,%r13
add    $0x10,%r13
xor    %ebp,%ebp
mov    0x70(%rsp),%rbx
xorpd  %xmm6,%xmm6
movsd  %xmm5,0x48(%rsp)
movsd  %xmm4,0x40(%rsp)
jmp    1344 <interceptCyl+0x224>
movapd %xmm2,%xmm3
xorpd  %xmm4,%xmm4
cmpltsd %xmm0,%xmm4
andpd  %xmm4,%xmm3
andnpd %xmm2,%xmm4
orpd   %xmm3,%xmm4
movsd  0xce7(%rip),%xmm0
cmpltsd %xmm1,%xmm0
andpd  %xmm0,%xmm4
andnpd %xmm2,%xmm0
orpd   %xmm4,%xmm0
movapd %xmm0,0x20(%rsp)
add    $0xffffffffffffffff,%rbp
cmp    %ebp,%ebx
je     1190 <interceptCyl+0x70>
movsd  (%r15,%rbp,8),%xmm0
movsd  0x8(%r15,%rbp,8),%xmm3
subsd  %xmm0,%xmm3
movsd  (%r14,%rbp,8),%xmm2
movsd  0x8(%r14,%rbp,8),%xmm1
subsd  %xmm2,%xmm1
ucomisd %xmm6,%xmm1
jne    1372 <interceptCyl+0x252>
jnp    1460 <interceptCyl+0x340>
movapd %xmm1,0x60(%rsp)
divsd  %xmm1,%xmm3
movsd  %xmm3,0x10(%rsp)
movsd  0x0(%r13),%xmm1
subsd  %xmm2,%xmm1
mulsd  %xmm3,%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0x58(%rsp)
movsd  (%r12),%xmm0
mulsd  %xmm3,%xmm0
mov    $0x2,%edi
call   1030 <pow@plt>
movsd  0x40(%rsp),%xmm1
subsd  %xmm0,%xmm1
movapd %xmm1,0x30(%rsp)
movsd  0x58(%rsp),%xmm0
movapd %xmm0,%xmm1
mulsd  0xc40(%rip),%xmm1
mulsd  0x10(%rsp),%xmm1
mulsd  (%r12),%xmm1
addsd  0x48(%rsp),%xmm1
movapd %xmm1,0x10(%rsp)
mov    $0x2,%edi
call   1030 <pow@plt>
movapd 0x30(%rsp),%xmm4
xorpd  %xmm6,%xmm6
movsd  0x50(%rsp),%xmm1
subsd  %xmm0,%xmm1
ucomisd %xmm6,%xmm4
jne    1412 <interceptCyl+0x2f2>
jnp    1510 <interceptCyl+0x3f0>
movapd %xmm4,%xmm2
mulsd  0xbfa(%rip),%xmm2
mulsd  %xmm1,%xmm2
movapd 0x10(%rsp),%xmm1
movapd %xmm1,%xmm0
mulsd  %xmm1,%xmm0
addsd  %xmm2,%xmm0
ucomisd %xmm6,%xmm0
jb     1338 <interceptCyl+0x218>
jb     15c3 <interceptCyl+0x4a3>
movapd %xmm1,%xmm2
sqrtsd %xmm0,%xmm0
jmp    15d8 <interceptCyl+0x4b8>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
ucomisd %xmm6,%xmm3
jne    146c <interceptCyl+0x34c>
jnp    1338 <interceptCyl+0x218>
movsd  (%r12),%xmm0
ucomisd %xmm6,%xmm0
jne    147e <interceptCyl+0x35e>
jnp    1338 <interceptCyl+0x218>
movsd  %xmm3,0x10(%rsp)
subsd  0x0(%r13),%xmm2
divsd  %xmm0,%xmm2
movapd %xmm2,0x30(%rsp)
mov    0xa8(%rsp),%rax
movsd  (%rax),%xmm0
mulsd  %xmm2,%xmm0
mov    0xb0(%rsp),%rax
addsd  (%rax),%xmm0
mov    $0x2,%edi
call   1030 <pow@plt>
movsd  %xmm0,0x58(%rsp)
mov    0x98(%rsp),%rax
movsd  (%rax),%xmm0
mulsd  0x30(%rsp),%xmm0
mov    0xa0(%rsp),%rax
addsd  (%rax),%xmm0
mov    $0x2,%edi
call   1030 <pow@plt>
xorpd  %xmm6,%xmm6
addsd  0x58(%rsp),%xmm0
ucomisd %xmm6,%xmm0
jb     16b1 <interceptCyl+0x591>
sqrtsd %xmm0,%xmm0
jmp    16ba <interceptCyl+0x59a>
cs nopw 0x0(%rax,%rax,1)
nop
movapd 0x10(%rsp),%xmm4
ucomisd %xmm6,%xmm4
jne    1522 <interceptCyl+0x402>
jnp    1338 <interceptCyl+0x218>
ucomisd %xmm6,%xmm1
jne    152e <interceptCyl+0x40e>
jnp    1338 <interceptCyl+0x218>
xorpd  0xaea(%rip),%xmm4
divsd  %xmm1,%xmm4
movsd  0x0(%r13),%xmm0
subsd  (%r14,%rbp,8),%xmm0
movsd  (%r12),%xmm1
mulsd  %xmm4,%xmm1
addsd  %xmm0,%xmm1
divsd  0x60(%rsp),%xmm1
xorpd  %xmm0,%xmm0
cmpltsd %xmm1,%xmm0
cmplesd 0xa94(%rip),%xmm1
movapd %xmm1,%xmm2
andpd  %xmm4,%xmm2
movapd 0x20(%rsp),%xmm3
andnpd %xmm3,%xmm1
orpd   %xmm2,%xmm1
andpd  %xmm0,%xmm1
andnpd %xmm3,%xmm0
orpd   %xmm1,%xmm0
movsd  0xa72(%rip),%xmm1
cmpltsd %xmm4,%xmm1
cmpltsd %xmm3,%xmm4
andpd  %xmm4,%xmm0
andnpd %xmm3,%xmm4
orpd   %xmm0,%xmm4
andpd  %xmm1,%xmm4
andnpd %xmm3,%xmm1
orpd   %xmm4,%xmm1
movapd %xmm1,0x20(%rsp)
jmp    1338 <interceptCyl+0x218>
call   1040 <sqrt@plt>
movapd 0x30(%rsp),%xmm4
movapd 0x10(%rsp),%xmm2
xorpd  %xmm6,%xmm6
movapd 0x60(%rsp),%xmm3
movapd %xmm0,%xmm1
subsd  %xmm2,%xmm1
addsd  %xmm2,%xmm0
unpcklpd %xmm0,%xmm1
mulpd  0xa3a(%rip),%xmm1
unpcklpd %xmm4,%xmm4
movsd  0x0(%r13),%xmm2
subsd  (%r14,%rbp,8),%xmm2
divpd  %xmm4,%xmm1
movsd  (%r12),%xmm0
unpcklpd %xmm0,%xmm0
unpcklpd %xmm2,%xmm2
mulpd  %xmm1,%xmm0
addpd  %xmm2,%xmm0
unpcklpd %xmm3,%xmm3
divpd  %xmm3,%xmm0
movapd %xmm0,%xmm2
cmplepd 0xa0b(%rip),%xmm2
movmskpd %xmm2,%eax
movapd %xmm1,%xmm3
unpckhpd %xmm1,%xmm3
movapd %xmm3,%xmm2
test   $0x2,%al
jne    164f <interceptCyl+0x52f>
movapd 0x20(%rsp),%xmm2
movapd %xmm0,%xmm4
unpckhpd %xmm0,%xmm4
xorpd  %xmm5,%xmm5
cmpltsd %xmm4,%xmm5
andpd  %xmm5,%xmm2
movapd 0x20(%rsp),%xmm4
andnpd %xmm4,%xmm5
orpd   %xmm2,%xmm5
movsd  0x98e(%rip),%xmm2
cmpltsd %xmm3,%xmm2
cmpltsd %xmm4,%xmm3
andpd  %xmm3,%xmm5
andnpd %xmm4,%xmm3
orpd   %xmm5,%xmm3
andpd  %xmm2,%xmm3
andnpd %xmm4,%xmm2
orpd   %xmm3,%xmm2
test   $0x1,%al
je     1300 <interceptCyl+0x1e0>
movapd %xmm1,%xmm3
minsd  %xmm2,%xmm3
jmp    1304 <interceptCyl+0x1e4>
call   1040 <sqrt@plt>
xorpd  %xmm6,%xmm6
movapd 0x30(%rsp),%xmm7
movsd  0x10(%rsp),%xmm1
subsd  (%r15,%rbp,8),%xmm0
divsd  %xmm1,%xmm0
xorpd  %xmm1,%xmm1
cmpltsd %xmm0,%xmm1
cmplesd 0x91e(%rip),%xmm0
movapd %xmm0,%xmm2
andpd  %xmm7,%xmm2
movapd 0x20(%rsp),%xmm3
andnpd %xmm3,%xmm0
orpd   %xmm2,%xmm0
andpd  %xmm1,%xmm0
andnpd %xmm3,%xmm1
orpd   %xmm0,%xmm1
movsd  0x8fc(%rip),%xmm0
cmpltsd %xmm7,%xmm0
cmpltsd %xmm3,%xmm7
andpd  %xmm7,%xmm1
andnpd %xmm3,%xmm7
orpd   %xmm1,%xmm7
andpd  %xmm0,%xmm7
andnpd %xmm3,%xmm0
orpd   %xmm7,%xmm0
jmp    1332 <interceptCyl+0x212>
add    $0xb8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
@INFINITY = global double 0.000000e+00
define void @interceptCyl(double* nocapture writeonly %0, [3 x double]* nocapture %1, [3 x double]* nocapture %2, double* nocapture %3, double* nocapture %4, i32 %5, i32 %6) {
%8 = icmp eq i32 %6, 0
br i1 %8, label %182, label %9
9: ; preds = %7
%10 = icmp eq i32 %5, 0
%11 = sext i32 %5 to i64
%12 = sext i32 %6 to i64
br label %13
13: ; preds = %9, %177
%14 = phi i64 [ %12, %9 ], [ %15, %177 ]
%15 = add nsw i64 %14, -1
%16 = load double, double* @INFINITY
%17 = getelementptr inbounds [3 x double], [3 x double]* %1, i64 %15, i64 0
%18 = load double, double* %17
%19 = tail call double @pow(double %18, i32 2)
%20 = getelementptr inbounds [3 x double], [3 x double]* %1, i64 %15, i64 1
%21 = load double, double* %20
%22 = tail call double @pow(double %21, i32 2)
%23 = fadd double %19, %22
%24 = load double, double* %17
%25 = getelementptr inbounds [3 x double], [3 x double]* %2, i64 %15, i64 0
%26 = load double, double* %25
%27 = load double, double* %20
%28 = getelementptr inbounds [3 x double], [3 x double]* %2, i64 %15, i64 1
%29 = load double, double* %28
%30 = fmul double %27, %29
%31 = tail call double @llvm.fmuladd.f64(double %24, double %26, double %30)
%32 = fmul double %31, 2.000000e+00
%33 = tail call double @pow(double %26, i32 2)
%34 = load double, double* %28
%35 = tail call double @pow(double %34, i32 2)
%36 = fadd double %33, %35
br i1 %10, label %177, label %37
37: ; preds = %13
%38 = getelementptr inbounds [3 x double], [3 x double]* %2, i64 %15, i64 2
%39 = getelementptr inbounds [3 x double], [3 x double]* %1, i64 %15, i64 2
br label %40
40: ; preds = %37, %173
%41 = phi i64 [ %11, %37 ], [ %43, %173 ]
%42 = phi double [ %16, %37 ], [ %174, %173 ]
%43 = add nsw i64 %41, -1
%44 = getelementptr inbounds double, double* %3, i64 %41
%45 = load double, double* %44
%46 = getelementptr inbounds double, double* %3, i64 %43
%47 = load double, double* %46
%48 = fsub double %45, %47
%49 = getelementptr inbounds double, double* %4, i64 %41
%50 = load double, double* %49
%51 = getelementptr inbounds double, double* %4, i64 %43
%52 = load double, double* %51
%53 = fsub double %50, %52
%54 = fcmp une double %53, 0.000000e+00
br i1 %54, label %85, label %55
55: ; preds = %40
%56 = load double, double* %38
%57 = fcmp une double %56, 0.000000e+00
%58 = fcmp une double %48, 0.000000e+00
%59 = and i1 %58, %57
br i1 %59, label %60, label %173
60: ; preds = %55
%61 = load double, double* %39
%62 = fsub double %52, %61
%63 = fdiv double %62, %56
%64 = load double, double* %17
%65 = load double, double* %25
%66 = tail call double @llvm.fmuladd.f64(double %65, double %63, double %64)
%67 = tail call double @pow(double %66, i32 2)
%68 = load double, double* %20
%69 = load double, double* %28
%70 = tail call double @llvm.fmuladd.f64(double %69, double %63, double %68)
%71 = tail call double @pow(double %70, i32 2)
%72 = fadd double %67, %71
%73 = tail call double @sqrt(double %72)
%74 = load double, double* %46
%75 = fsub double %73, %74
%76 = fdiv double %75, %48
%77 = fcmp ogt double %76, 0.000000e+00
%78 = fcmp ole double %76, 1.000000e+00
%79 = and i1 %77, %78
%80 = fcmp olt double %63, %42
%81 = and i1 %80, %79
%82 = fcmp ogt double %63, 0x3EB0C6F7A0B5ED8D
%83 = and i1 %82, %81
%84 = select i1 %83, double %63, double %42
br label %173
85: ; preds = %40
%86 = fdiv double %48, %53
%87 = load double, double* %39
%88 = fsub double %87, %52
%89 = tail call double @llvm.fmuladd.f64(double %86, double %88, double %47)
%90 = load double, double* %38
%91 = fmul double %86, %90
%92 = tail call double @pow(double %91, i32 2)
%93 = fsub double %36, %92
%94 = load double, double* %38
%95 = fmul double %89, -2.000000e+00
%96 = fmul double %86, %95
%97 = tail call double @llvm.fmuladd.f64(double %96, double %94, double %32)
%98 = tail call double @pow(double %89, i32 2)
%99 = fsub double %23, %98
%100 = fcmp une double %93, 0.000000e+00
br i1 %100, label %101, label %152
101: ; preds = %85
%102 = fmul double %93, -4.000000e+00
%103 = fmul double %102, %99
%104 = tail call double @llvm.fmuladd.f64(double %97, double %97, double %103)
%105 = fcmp ult double %104, 0.000000e+00
br i1 %105, label %173, label %106
106: ; preds = %101
%107 = tail call double @sqrt(double %104)
%108 = insertelement <2 x double> poison, double %107, i64 0
%109 = shufflevector <2 x double> %108, <2 x double> poison, <2 x i32> zeroinitializer
%110 = insertelement <2 x double> poison, double %97, i64 0
%111 = shufflevector <2 x double> %110, <2 x double> poison, <2 x i32> zeroinitializer
%112 = fsub <2 x double> %109, %111
%113 = fadd <2 x double> %109, %111
%114 = shufflevector <2 x double> %112, <2 x double> %113, <2 x i32> <i32 0, i32 3>
%115 = fmul <2 x double> %114, <double 5.000000e-01, double -5.000000e-01>
%116 = insertelement <2 x double> poison, double %93, i64 0
%117 = shufflevector <2 x double> %116, <2 x double> poison, <2 x i32> zeroinitializer
%118 = fdiv <2 x double> %115, %117
%119 = load double, double* %38
%120 = load double, double* %39
%121 = load double, double* %51
%122 = fsub double %120, %121
%123 = insertelement <2 x double> poison, double %119, i64 0
%124 = shufflevector <2 x double> %123, <2 x double> poison, <2 x i32> zeroinitializer
%125 = insertelement <2 x double> poison, double %122, i64 0
%126 = shufflevector <2 x double> %125, <2 x double> poison, <2 x i32> zeroinitializer
%127 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %118, <2 x double> %124, <2 x double> %126)
%128 = insertelement <2 x double> poison, double %53, i64 0
%129 = shufflevector <2 x double> %128, <2 x double> poison, <2 x i32> zeroinitializer
%130 = fdiv <2 x double> %127, %129
%131 = extractelement <2 x double> %130, i64 1
%132 = fcmp ogt double %131, 0.000000e+00
%133 = fcmp ole <2 x double> %130, <double 1.000000e+00, double 1.000000e+00>
%134 = extractelement <2 x i1> %133, i64 1
%135 = and i1 %132, %134
%136 = extractelement <2 x double> %118, i64 1
%137 = fcmp olt double %136, %42
%138 = and i1 %137, %135
%139 = fcmp ogt double %136, 0x3EB0C6F7A0B5ED8D
%140 = and i1 %139, %138
%141 = select i1 %140, double %136, double %42
%142 = extractelement <2 x double> %130, i64 0
%143 = fcmp ogt double %142, 0.000000e+00
%144 = extractelement <2 x i1> %133, i64 0
%145 = and i1 %143, %144
%146 = extractelement <2 x double> %118, i64 0
%147 = fcmp olt double %146, %141
%148 = and i1 %145, %147
%149 = fcmp ogt double %146, 0x3EB0C6F7A0B5ED8D
%150 = and i1 %149, %148
%151 = select i1 %150, double %146, double %141
br label %173
152: ; preds = %85
%153 = fcmp une double %97, 0.000000e+00
%154 = fcmp une double %99, 0.000000e+00
%155 = and i1 %153, %154
br i1 %155, label %156, label %173
156: ; preds = %152
%157 = fneg double %97
%158 = fdiv double %157, %99
%159 = load double, double* %38
%160 = load double, double* %39
%161 = load double, double* %51
%162 = fsub double %160, %161
%163 = tail call double @llvm.fmuladd.f64(double %158, double %159, double %162)
%164 = fdiv double %163, %53
%165 = fcmp ogt double %164, 0.000000e+00
%166 = fcmp ole double %164, 1.000000e+00
%167 = and i1 %165, %166
%168 = fcmp olt double %158, %42
%169 = and i1 %168, %167
%170 = fcmp ogt double %158, 0x3EB0C6F7A0B5ED8D
%171 = and i1 %170, %169
%172 = select i1 %171, double %158, double %42
br label %173
173: ; preds = %55, %60, %106, %156, %152, %101
%174 = phi double [ %42, %101 ], [ %42, %152 ], [ %172, %156 ], [ %151, %106 ], [ %84, %60 ], [ %42, %55 ]
%175 = trunc i64 %43 to i32
%176 = icmp eq i32 %175, 0
br i1 %176, label %177, label %40
177: ; preds = %173, %13
%178 = phi double [ %16, %13 ], [ %174, %173 ]
%179 = getelementptr inbounds double, double* %0, i64 %15
store double %178, double* %179
%180 = trunc i64 %15 to i32
%181 = icmp eq i32 %180, 0
br i1 %181, label %182, label %13
182: ; preds = %177, %7
ret void
}