int tpd9(int inverse, const int i[], const double p[], int Nhat,
         const double rawcrd[], double *discrd) {
  if (i[3 + inverse] != 60 || 2 < Nhat) {
    return 1;
  }
  double r, s;
  double u = rawcrd[0];
  double v = rawcrd[1];
  if (i[5]) {
    r = p[0] + p[1] * u + p[2] * v;
    v = p[3] + p[4] * u + p[5] * v;
    u = r;
    p += 6;
  }
  if (inverse)
    p += i[3];
  *discrd = p[0] +
            u * (p[1] +
                 u * (p[4] +
                      u * (p[7] +
                           u * (p[12] +
                                u * (p[17] +
                                     u * (p[24] +
                                          u * (p[31] +
                                               u * (p[40] + u * (p[49])))))))));
  if (Nhat == 1)
    return 0;
  *discrd +=
      v * (p[2] +
           v * (p[6] +
                v * (p[10] +
                     v * (p[16] +
                          v * (p[22] +
                               v * (p[30] +
                                    v * (p[38] +
                                         v * (p[48] + v * (p[58]))))))))) +
      u * (p[5] + v * (p[9] + v * (p[15] + v * (p[21] + v * (p[29] + v * (p[37] + v * (p[47] + v * (p[57]))))))) + u * (p[8] + v * (p[14] + v * (p[20] + v * (p[28] + v * (p[36] + v * (p[46] + v * (p[56])))))) + u * (p[13] + v * (p[19] + v * (p[27] + v * (p[35] + v * (p[45] + v * (p[55]))))) + u * (p[18] + v * (p[26] + v * (p[34] + v * (p[44] + v * (p[54])))) + u * (p[25] + v * (p[33] + v * (p[43] + v * (p[53]))) + u * (p[32] + v * (p[42] + v * (p[52])) + u * (p[41] + v * (p[51]) + u * (p[50])))))))) *
          v;
  if (i[6]) {
    s = u * u + v * v;
    r = sqrt(s);
    *discrd +=
        r * (p[3] + s * (p[11] + s * (p[23] + s * (p[39] + s * (p[59])))));
  }
  return 0;
}


x86 assembly (clang, -O0):
<tpd9>:
push   %rbp
mov    %rsp,%rbp
sub    $0x90,%rsp
mov    %edi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x8(%rbp),%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
cmpl   $0x3c,(%rax,%rcx,4)
jne    1156 <tpd9+0x46>
mov    $0x2,%eax
cmp    -0x1c(%rbp),%eax
jge    1162 <tpd9+0x52>
movl   $0x1,-0x4(%rbp)
jmp    1925 <tpd9+0x815>
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x28(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x10(%rbp),%rax
cmpl   $0x0,0x14(%rax)
je     1214 <tpd9+0x104>
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm1
movsd  -0x48(%rbp),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x18(%rbp),%rax
movsd  0x10(%rax),%xmm0
movsd  -0x50(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x18(%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  0x20(%rax),%xmm1
movsd  -0x48(%rbp),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x18(%rbp),%rax
movsd  0x28(%rax),%xmm0
movsd  -0x50(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x18(%rbp),%rax
add    $0x30,%rax
mov    %rax,-0x18(%rbp)
cmpl   $0x0,-0x8(%rbp)
je     1237 <tpd9+0x127>
mov    -0x10(%rbp),%rax
mov    0xc(%rax),%ecx
mov    -0x18(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x18(%rbp)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm3
movsd  -0x48(%rbp),%xmm2
mov    -0x18(%rbp),%rax
movsd  0x20(%rax),%xmm5
movsd  -0x48(%rbp),%xmm4
mov    -0x18(%rbp),%rax
movsd  0x38(%rax),%xmm7
movsd  -0x48(%rbp),%xmm6
mov    -0x18(%rbp),%rax
movsd  0x60(%rax),%xmm9
movsd  -0x48(%rbp),%xmm8
mov    -0x18(%rbp),%rax
movsd  0x88(%rax),%xmm11
movsd  -0x48(%rbp),%xmm10
mov    -0x18(%rbp),%rax
movsd  0xc0(%rax),%xmm13
movsd  -0x48(%rbp),%xmm12
mov    -0x18(%rbp),%rax
movsd  0xf8(%rax),%xmm15
movsd  -0x48(%rbp),%xmm14
mov    -0x18(%rbp),%rax
movsd  0x140(%rax),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x48(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movsd  0x188(%rax),%xmm1
mulsd  %xmm1,%xmm0
movsd  -0x68(%rbp),%xmm1
addsd  %xmm1,%xmm0
movsd  -0x60(%rbp),%xmm1
mulsd  %xmm0,%xmm14
movsd  -0x58(%rbp),%xmm0
addsd  %xmm15,%xmm14
mulsd  %xmm14,%xmm12
addsd  %xmm13,%xmm12
mulsd  %xmm12,%xmm10
addsd  %xmm11,%xmm10
mulsd  %xmm10,%xmm8
addsd  %xmm9,%xmm8
mulsd  %xmm8,%xmm6
addsd  %xmm7,%xmm6
mulsd  %xmm6,%xmm4
addsd  %xmm5,%xmm4
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
cmpl   $0x1,-0x1c(%rbp)
jne    1361 <tpd9+0x251>
movl   $0x0,-0x4(%rbp)
jmp    1925 <tpd9+0x815>
movsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x10(%rax),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x50(%rbp),%xmm2
mov    -0x18(%rbp),%rax
movsd  0x30(%rax),%xmm3
movsd  -0x50(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movsd  0x50(%rax),%xmm5
movsd  -0x50(%rbp),%xmm4
mov    -0x18(%rbp),%rax
movsd  0x80(%rax),%xmm7
movsd  -0x50(%rbp),%xmm6
mov    -0x18(%rbp),%rax
movsd  0xb0(%rax),%xmm9
movsd  -0x50(%rbp),%xmm8
mov    -0x18(%rbp),%rax
movsd  0xf0(%rax),%xmm11
movsd  -0x50(%rbp),%xmm10
mov    -0x18(%rbp),%rax
movsd  0x130(%rax),%xmm13
movsd  -0x50(%rbp),%xmm12
mov    -0x18(%rbp),%rax
movsd  0x180(%rax),%xmm15
movsd  -0x50(%rbp),%xmm14
mov    -0x18(%rbp),%rax
movsd  0x1d0(%rax),%xmm0
mulsd  %xmm0,%xmm14
movsd  -0x88(%rbp),%xmm0
ff
addsd  %xmm15,%xmm14
mulsd  %xmm14,%xmm12
addsd  %xmm13,%xmm12
mulsd  %xmm12,%xmm10
addsd  %xmm11,%xmm10
mulsd  %xmm10,%xmm8
addsd  %xmm9,%xmm8
mulsd  %xmm8,%xmm6
addsd  %xmm7,%xmm6
mulsd  %xmm6,%xmm4
addsd  %xmm5,%xmm4
mulsd  %xmm4,%xmm1
addsd  %xmm3,%xmm1
mulsd  %xmm1,%xmm2
addsd  %xmm0,%xmm2
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x28(%rax),%xmm0
movsd  -0x50(%rbp),%xmm4
mov    -0x18(%rbp),%rax
movsd  0x48(%rax),%xmm3
movsd  -0x50(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movsd  0x78(%rax),%xmm6
movsd  -0x50(%rbp),%xmm5
mov    -0x18(%rbp),%rax
movsd  0xa8(%rax),%xmm8
movsd  -0x50(%rbp),%xmm7
mov    -0x18(%rbp),%rax
movsd  0xe8(%rax),%xmm10
movsd  -0x50(%rbp),%xmm9
mov    -0x18(%rbp),%rax
movsd  0x128(%rax),%xmm12
movsd  -0x50(%rbp),%xmm11
mov    -0x18(%rbp),%rax
movsd  0x178(%rax),%xmm14
movsd  -0x50(%rbp),%xmm13
mov    -0x18(%rbp),%rax
movsd  0x1c8(%rax),%xmm15
mulsd  %xmm15,%xmm13
addsd  %xmm14,%xmm13
mulsd  %xmm13,%xmm11
addsd  %xmm12,%xmm11
mulsd  %xmm11,%xmm9
addsd  %xmm10,%xmm9
mulsd  %xmm9,%xmm7
addsd  %xmm8,%xmm7
mulsd  %xmm7,%xmm5
addsd  %xmm6,%xmm5
mulsd  %xmm5,%xmm1
addsd  %xmm3,%xmm1
mulsd  %xmm1,%xmm4
addsd  %xmm0,%xmm4
movsd  -0x48(%rbp),%xmm3
mov    -0x18(%rbp),%rax
movsd  0x40(%rax),%xmm0
movsd  -0x50(%rbp),%xmm6
mov    -0x18(%rbp),%rax
movsd  0x70(%rax),%xmm5
movsd  -0x50(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movsd  0xa0(%rax),%xmm8
movsd  -0x50(%rbp),%xmm7
mov    -0x18(%rbp),%rax
movsd  0xe0(%rax),%xmm10
movsd  -0x50(%rbp),%xmm9
mov    -0x18(%rbp),%rax
movsd  0x120(%rax),%xmm12
movsd  -0x50(%rbp),%xmm11
mov    -0x18(%rbp),%rax
movsd  0x170(%rax),%xmm14
movsd  -0x50(%rbp),%xmm13
mov    -0x18(%rbp),%rax
movsd  0x1c0(%rax),%xmm15
mulsd  %xmm15,%xmm13
addsd  %xmm14,%xmm13
mulsd  %xmm13,%xmm11
addsd  %xmm12,%xmm11
mulsd  %xmm11,%xmm9
addsd  %xmm10,%xmm9
mulsd  %xmm9,%xmm7
addsd  %xmm8,%xmm7
mulsd  %xmm7,%xmm1
addsd  %xmm5,%xmm1
mulsd  %xmm1,%xmm6
addsd  %xmm0,%xmm6
movsd  -0x48(%rbp),%xmm5
mov    -0x18(%rbp),%rax
movsd  0x68(%rax),%xmm0
movsd  -0x50(%rbp),%xmm8
mov    -0x18(%rbp),%rax
movsd  0x98(%rax),%xmm7
movsd  -0x50(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movsd  0xd8(%rax),%xmm10
movsd  -0x50(%rbp),%xmm9
mov    -0x18(%rbp),%rax
movsd  0x118(%rax),%xmm12
movsd  -0x50(%rbp),%xmm11
mov    -0x18(%rbp),%rax
movsd  0x168(%rax),%xmm14
movsd  -0x50(%rbp),%xmm13
mov    -0x18(%rbp),%rax
movsd  0x1b8(%rax),%xmm15
mulsd  %xmm15,%xmm13
addsd  %xmm14,%xmm13
mulsd  %xmm13,%xmm11
addsd  %xmm12,%xmm11
mulsd  %xmm11,%xmm9
addsd  %xmm10,%xmm9
mulsd  %xmm9,%xmm1
addsd  %xmm7,%xmm1
mulsd  %xmm1,%xmm8
addsd  %xmm0,%xmm8
movsd  -0x48(%rbp),%xmm7
mov    -0x18(%rbp),%rax
movsd  0x90(%rax),%xmm0
movsd  -0x50(%rbp),%xmm10
mov    -0x18(%rbp),%rax
movsd  0xd0(%rax),%xmm9
movsd  -0x50(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movsd  0x110(%rax),%xmm12
movsd  -0x50(%rbp),%xmm11
mov    -0x18(%rbp),%rax
movsd  0x160(%rax),%xmm14
movsd  -0x50(%rbp),%xmm13
mov    -0x18(%rbp),%rax
movsd  0x1b0(%rax),%xmm15
mulsd  %xmm15,%xmm13
addsd  %xmm14,%xmm13
mulsd  %xmm13,%xmm11
addsd  %xmm12,%xmm11
mulsd  %xmm11,%xmm1
addsd  %xmm9,%xmm1
mulsd  %xmm1,%xmm10
addsd  %xmm0,%xmm10
movsd  -0x48(%rbp),%xmm9
mov    -0x18(%rbp),%rax
movsd  0xc8(%rax),%xmm0
movsd  -0x50(%rbp),%xmm12
mov    -0x18(%rbp),%rax
movsd  0x108(%rax),%xmm11
movsd  -0x50(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movsd  0x158(%rax),%xmm14
movsd  -0x50(%rbp),%xmm13
mov    -0x18(%rbp),%rax
movsd  0x1a8(%rax),%xmm15
mulsd  %xmm15,%xmm13
addsd  %xmm14,%xmm13
mulsd  %xmm13,%xmm1
addsd  %xmm11,%xmm1
mulsd  %xmm1,%xmm12
addsd  %xmm0,%xmm12
movsd  -0x48(%rbp),%xmm11
mov    -0x18(%rbp),%rax
movsd  0x100(%rax),%xmm0
movsd  -0x50(%rbp),%xmm14
mov    -0x18(%rbp),%rax
movsd  0x150(%rax),%xmm13
movsd  -0x50(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movsd  0x1a0(%rax),%xmm15
mulsd  %xmm15,%xmm1
addsd  %xmm13,%xmm1
mulsd  %xmm1,%xmm14
addsd  %xmm0,%xmm14
movsd  -0x48(%rbp),%xmm13
mov    -0x18(%rbp),%rax
movsd  0x148(%rax),%xmm15
movsd  -0x50(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movsd  0x198(%rax),%xmm1
mulsd  %xmm1,%xmm0
movsd  -0x80(%rbp),%xmm1
addsd  %xmm15,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x48(%rbp),%xmm15
mov    -0x18(%rbp),%rax
movsd  0x190(%rax),%xmm0
mulsd  %xmm0,%xmm15
movsd  -0x78(%rbp),%xmm0
addsd  %xmm0,%xmm15
movsd  -0x70(%rbp),%xmm0
mulsd  %xmm15,%xmm13
addsd  %xmm14,%xmm13
mulsd  %xmm13,%xmm11
addsd  %xmm12,%xmm11
mulsd  %xmm11,%xmm9
addsd  %xmm10,%xmm9
mulsd  %xmm9,%xmm7
addsd  %xmm8,%xmm7
mulsd  %xmm7,%xmm5
addsd  %xmm6,%xmm5
mulsd  %xmm5,%xmm3
addsd  %xmm4,%xmm3
mulsd  %xmm3,%xmm1
mulsd  -0x50(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
addsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
mov    -0x10(%rbp),%rax
cmpl   $0x0,0x18(%rax)
je     191e <tpd9+0x80e>
movsd  -0x48(%rbp),%xmm0
movsd  -0x48(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  -0x50(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movsd  0x18(%rax),%xmm1
movsd  -0x40(%rbp),%xmm2
mov    -0x18(%rbp),%rax
movsd  0x58(%rax),%xmm4
movsd  -0x40(%rbp),%xmm3
mov    -0x18(%rbp),%rax
movsd  0xb8(%rax),%xmm6
movsd  -0x40(%rbp),%xmm5
mov    -0x18(%rbp),%rax
movsd  0x138(%rax),%xmm8
movsd  -0x40(%rbp),%xmm7
mov    -0x18(%rbp),%rax
movsd  0x1d8(%rax),%xmm9
mulsd  %xmm9,%xmm7
addsd  %xmm8,%xmm7
mulsd  %xmm7,%xmm5
addsd  %xmm6,%xmm5
mulsd  %xmm5,%xmm3
addsd  %xmm4,%xmm3
mulsd  %xmm3,%xmm2
addsd  %xmm1,%xmm2
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x90,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<tpd9>:
push   %rbp
push   %r14
push   %rbx
sub    $0x10,%rsp
mov    $0x1,%ebp
cmp    $0x2,%ecx
jg     14fb <tpd9+0x3eb>
movslq %edi,%rax
cmpl   $0x3c,0xc(%rsi,%rax,4)
jne    14fb <tpd9+0x3eb>
mov    %r9,%r14
mov    %rdx,%rbx
movupd (%r8),%xmm11
cmpl   $0x0,0x14(%rsi)
je     1187 <tpd9+0x77>
movsd  (%rbx),%xmm0
movsd  0x8(%rbx),%xmm1
movhpd 0x20(%rbx),%xmm1
movsd  0x10(%rbx),%xmm2
movapd %xmm11,%xmm3
unpcklpd %xmm11,%xmm3
movhpd 0x18(%rbx),%xmm0
mulpd  %xmm1,%xmm3
addpd  %xmm3,%xmm0
movhpd 0x28(%rbx),%xmm2
unpckhpd %xmm11,%xmm11
mulpd  %xmm2,%xmm11
addpd  %xmm0,%xmm11
add    $0x30,%rbx
test   %edi,%edi
je     1193 <tpd9+0x83>
movslq 0xc(%rsi),%rax
lea    (%rbx,%rax,8),%rbx
movsd  0x188(%rbx),%xmm0
mulsd  %xmm11,%xmm0
addsd  0x140(%rbx),%xmm0
mulsd  %xmm11,%xmm0
addsd  0xf8(%rbx),%xmm0
mulsd  %xmm11,%xmm0
addsd  0xc0(%rbx),%xmm0
mulsd  %xmm11,%xmm0
addsd  0x88(%rbx),%xmm0
mulsd  %xmm11,%xmm0
addsd  0x60(%rbx),%xmm0
mulsd  %xmm11,%xmm0
addsd  0x38(%rbx),%xmm0
mulsd  %xmm11,%xmm0
addsd  0x20(%rbx),%xmm0
mulsd  %xmm11,%xmm0
addsd  0x8(%rbx),%xmm0
mulsd  %xmm11,%xmm0
addsd  (%rbx),%xmm0
movsd  %xmm0,(%r14)
xor    %ebp,%ebp
cmp    $0x1,%ecx
je     14fb <tpd9+0x3eb>
movapd %xmm11,%xmm1
unpckhpd %xmm11,%xmm1
movsd  0x1d0(%rbx),%xmm9
mulsd  %xmm1,%xmm9
addsd  0x180(%rbx),%xmm9
mulsd  %xmm1,%xmm9
addsd  0x130(%rbx),%xmm9
mulsd  %xmm1,%xmm9
addsd  0xf0(%rbx),%xmm9
mulsd  %xmm1,%xmm9
addsd  0xb0(%rbx),%xmm9
mulsd  %xmm1,%xmm9
addsd  0x80(%rbx),%xmm9
mulsd  %xmm1,%xmm9
addsd  0x50(%rbx),%xmm9
mulsd  %xmm1,%xmm9
addsd  0x30(%rbx),%xmm9
mulsd  %xmm1,%xmm9
addsd  0x10(%rbx),%xmm9
movsd  0x1c8(%rbx),%xmm8
mulsd  %xmm1,%xmm8
addsd  0x178(%rbx),%xmm8
mulsd  %xmm1,%xmm8
addsd  0x128(%rbx),%xmm8
mulsd  %xmm1,%xmm8
addsd  0xe8(%rbx),%xmm8
mulsd  %xmm1,%xmm8
addsd  0xa8(%rbx),%xmm8
mulsd  %xmm1,%xmm8
addsd  0x78(%rbx),%xmm8
mulsd  %xmm1,%xmm8
addsd  0x48(%rbx),%xmm8
mulsd  %xmm1,%xmm8
addsd  0x28(%rbx),%xmm8
movsd  0x1c0(%rbx),%xmm4
mulsd  %xmm1,%xmm4
addsd  0x170(%rbx),%xmm4
mulsd  %xmm1,%xmm4
addsd  0x120(%rbx),%xmm4
mulsd  %xmm1,%xmm4
addsd  0xe0(%rbx),%xmm4
mulsd  %xmm1,%xmm4
addsd  0xa0(%rbx),%xmm4
mulsd  %xmm1,%xmm4
addsd  0x70(%rbx),%xmm4
mulsd  %xmm1,%xmm4
addsd  0x40(%rbx),%xmm4
movsd  0x1b8(%rbx),%xmm5
mulsd  %xmm1,%xmm5
addsd  0x168(%rbx),%xmm5
mulsd  %xmm1,%xmm5
addsd  0x118(%rbx),%xmm5
mulsd  %xmm1,%xmm5
addsd  0xd8(%rbx),%xmm5
mulsd  %xmm1,%xmm5
addsd  0x98(%rbx),%xmm5
mulsd  %xmm1,%xmm5
addsd  0x68(%rbx),%xmm5
movsd  0x1b0(%rbx),%xmm6
mulsd  %xmm1,%xmm6
addsd  0x160(%rbx),%xmm6
mulsd  %xmm1,%xmm6
addsd  0x110(%rbx),%xmm6
mulsd  %xmm1,%xmm6
addsd  0xd0(%rbx),%xmm6
mulsd  %xmm1,%xmm6
addsd  0x90(%rbx),%xmm6
movsd  0x1a8(%rbx),%xmm7
mulsd  %xmm1,%xmm7
addsd  0x158(%rbx),%xmm7
mulsd  %xmm1,%xmm7
addsd  0x108(%rbx),%xmm7
mulsd  %xmm1,%xmm7
addsd  0xc8(%rbx),%xmm7
movsd  0x1a0(%rbx),%xmm2
mulsd  %xmm1,%xmm2
addsd  0x150(%rbx),%xmm2
mulsd  %xmm1,%xmm2
addsd  0x100(%rbx),%xmm2
movsd  0x198(%rbx),%xmm10
mulsd  %xmm1,%xmm10
addsd  0x148(%rbx),%xmm10
movsd  0x190(%rbx),%xmm3
mulsd  %xmm11,%xmm3
addsd  %xmm10,%xmm3
mulsd  %xmm11,%xmm3
addsd  %xmm2,%xmm3
mulsd  %xmm11,%xmm3
addsd  %xmm7,%xmm3
mulsd  %xmm11,%xmm3
addsd  %xmm6,%xmm3
mulsd  %xmm11,%xmm3
addsd  %xmm5,%xmm3
mulsd  %xmm11,%xmm3
addsd  %xmm4,%xmm3
mulsd  %xmm11,%xmm3
addsd  %xmm8,%xmm3
mulsd  %xmm11,%xmm3
mulsd  %xmm1,%xmm3
mulsd  %xmm1,%xmm9
addsd  %xmm3,%xmm9
addsd  %xmm9,%xmm0
movsd  %xmm0,(%r14)
cmpl   $0x0,0x18(%rsi)
je     14fb <tpd9+0x3eb>
mulsd  %xmm1,%xmm1
mulsd  %xmm11,%xmm11
addsd  %xmm1,%xmm11
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm11
jb     14a1 <tpd9+0x391>
xorps  %xmm0,%xmm0
sqrtsd %xmm11,%xmm0
jmp    14b7 <tpd9+0x3a7>
movapd %xmm11,%xmm0
movapd %xmm11,(%rsp)
call   1030 <sqrt@plt>
movapd (%rsp),%xmm11
movsd  0x1d8(%rbx),%xmm1
mulsd  %xmm11,%xmm1
addsd  0x138(%rbx),%xmm1
mulsd  %xmm11,%xmm1
addsd  0xb8(%rbx),%xmm1
mulsd  %xmm11,%xmm1
addsd  0x58(%rbx),%xmm1
mulsd  %xmm11,%xmm1
addsd  0x18(%rbx),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%r14),%xmm1
movsd  %xmm1,(%r14)
mov    %ebp,%eax
add    $0x10,%rsp
pop    %rbx
pop    %r14
pop    %rbp
ret

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

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