void polydplm(__float128 pl1[], __float128 pl2[], int l1, int l2, int m,
              __float128 xi) {
  int k;
  __float128 log2xi = logq(2 * xi);
  if (m == 0) {
    for (k = 1; k <= l1; k++)
      pl1[k - 1] = expq((gammaq(1 + k + l1)) - (gammaq(1 + k)) -
                        (gammaq(1 + l1 - k)) - (gammaq(k)) - (k - 1) * log2xi) /
                   2;
    for (k = 1; k <= l2; k++)
      pl2[k - 1] = expq((gammaq(1 + k + l2)) - (gammaq(1 + k)) -
                        (gammaq(1 + l2 - k)) - (gammaq(k)) - (k - 1) * log2xi) /
                   2;
  } else {
    for (k = 0; k <= l1 - m + 1; k++)
      pl1[k] = 0;
    for (k = 0; k <= l2 - m + 1; k++)
      pl2[k] = 0;
    pl1[l1 + 1 - m] +=
        (l1 - m + 1) * expq((gammaq(2 * l1 + 3)) - (gammaq(l1 + 2)) -
                            (gammaq(l1 - m + 2)) - (l1 + 1 - m) * log2xi);
    for (k = 0; k <= l1 - m; k++) {
      __float128 common =
          expq((gammaq(1 + k + l1 + m)) - (gammaq(1 + k + m)) -
               (gammaq(1 + k)) - (gammaq(1 + l1 - k - m)) - k * log2xi);
      pl1[k] += common *
                (((m) * (m)) + m * (k - l1 - 1) - 2.0 * k * l1 - 2 * k) /
                (m - l1 + k - 1);
      pl1[k + 1] -= common * (l1 + 1) / xi;
    }
    pl2[l2 + 1 - m] +=
        (l2 - m + 1) * expq((gammaq(2 * l2 + 3)) - (gammaq(l2 + 2)) -
                            (gammaq(l2 - m + 2)) - (l2 + 1 - m) * log2xi);
    for (k = 0; k <= l2 - m; k++) {
      __float128 common =
          expq((gammaq(1 + k + l2 + m)) - (gammaq(1 + k + m)) -
               (gammaq(1 + k)) - (gammaq(1 + l2 - k - m)) - k * log2xi);
      pl2[k] += common *
                (((m) * (m)) + m * (k - l2 - 1) - 2.0 * k * l2 - 2 * k) /
                (m - l2 + k - 1);
      pl2[k + 1] -= common * (l2 + 1.) / xi;
    }
  }
}


x86 assembly (clang, -O0):
<polydplm>:
push   %rbp
mov    %rsp,%rbp
sub    $0x250,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8d,-0x1c(%rbp)
movaps %xmm0,-0x30(%rbp)
movaps -0x30(%rbp),%xmm1
movaps %xmm1,%xmm0
call   1ac0 <__addtf3>
mov    $0x1,%al
call   1040 <logq@plt>
mov    %eax,%edi
call   5800 <__floatsitf>
movaps %xmm0,-0x50(%rbp)
cmpl   $0x0,-0x1c(%rbp)
jne    137c <polydplm+0x24c>
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     1275 <polydplm+0x145>
mov    -0x34(%rbp),%edx
mov    -0x14(%rbp),%eax
mov    %eax,%ecx
mov    %edx,%eax
lea    0x1(%rax,%rcx,1),%edi
xor    %eax,%eax
mov    %al,-0x85(%rbp)
call   1030 <gammaq@plt>
mov    %eax,%ecx
mov    -0x85(%rbp),%al
mov    %ecx,-0x90(%rbp)
mov    -0x34(%rbp),%edi
add    $0x1,%edi
call   1030 <gammaq@plt>
mov    -0x90(%rbp),%ecx
mov    %eax,%edx
mov    -0x85(%rbp),%al
sub    %edx,%ecx
mov    %ecx,-0x8c(%rbp)
mov    -0x14(%rbp),%edi
mov    -0x34(%rbp),%ecx
sub    %ecx,%edi
add    $0x1,%edi
call   1030 <gammaq@plt>
mov    -0x8c(%rbp),%ecx
mov    %eax,%edx
mov    -0x85(%rbp),%al
sub    %edx,%ecx
mov    %ecx,-0x84(%rbp)
mov    -0x34(%rbp),%edi
call   1030 <gammaq@plt>
mov    -0x84(%rbp),%edi
sub    %eax,%edi
call   5800 <__floatsitf>
movaps %xmm0,-0x80(%rbp)
mov    -0x34(%rbp),%edi
add    $0xffffffff,%edi
call   5800 <__floatsitf>
movaps -0x50(%rbp),%xmm1
call   39f0 <__multf3>
movaps %xmm0,%xmm1
movaps -0x80(%rbp),%xmm0
call   44c0 <__subtf3>
mov    $0x1,%al
call   1050 <expq@plt>
mov    %eax,%edi
mov    %edi,%eax
shr    $0x1f,%eax
add    %eax,%edi
sar    %edi
call   5800 <__floatsitf>
mov    -0x8(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    $0xffffffff,%ecx
movslq %ecx,%rcx
shl    $0x4,%rcx
mov    %rax,%rdx
add    %rcx,%rdx
movaps %xmm0,(%rax,%rcx,1)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1180 <polydplm+0x50>
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     1377 <polydplm+0x247>
mov    -0x34(%rbp),%edx
mov    -0x18(%rbp),%eax
mov    %eax,%ecx
mov    %edx,%eax
lea    0x1(%rax,%rcx,1),%edi
xor    %eax,%eax
mov    %al,-0xa5(%rbp)
call   1030 <gammaq@plt>
mov    %eax,%ecx
mov    -0xa5(%rbp),%al
mov    %ecx,-0xb0(%rbp)
mov    -0x34(%rbp),%edi
add    $0x1,%edi
call   1030 <gammaq@plt>
mov    -0xb0(%rbp),%ecx
mov    %eax,%edx
mov    -0xa5(%rbp),%al
sub    %edx,%ecx
mov    %ecx,-0xac(%rbp)
mov    -0x18(%rbp),%edi
mov    -0x34(%rbp),%ecx
sub    %ecx,%edi
add    $0x1,%edi
call   1030 <gammaq@plt>
mov    -0xac(%rbp),%ecx
mov    %eax,%edx
mov    -0xa5(%rbp),%al
sub    %edx,%ecx
mov    %ecx,-0xa4(%rbp)
mov    -0x34(%rbp),%edi
call   1030 <gammaq@plt>
mov    -0xa4(%rbp),%edi
sub    %eax,%edi
call   5800 <__floatsitf>
movaps %xmm0,-0xa0(%rbp)
mov    -0x34(%rbp),%edi
add    $0xffffffff,%edi
call   5800 <__floatsitf>
movaps -0x50(%rbp),%xmm1
call   39f0 <__multf3>
movaps %xmm0,%xmm1
movaps -0xa0(%rbp),%xmm0
call   44c0 <__subtf3>
mov    $0x1,%al
call   1050 <expq@plt>
mov    %eax,%edi
mov    %edi,%eax
shr    $0x1f,%eax
add    %eax,%edi
sar    %edi
call   5800 <__floatsitf>
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    $0xffffffff,%ecx
movslq %ecx,%rcx
shl    $0x4,%rcx
mov    %rax,%rdx
add    %rcx,%rdx
movaps %xmm0,(%rax,%rcx,1)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    127c <polydplm+0x14c>
jmp    1aa8 <polydplm+0x978>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jg     13c0 <polydplm+0x290>
xor    %eax,%eax
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
shl    $0x4,%rcx
mov    %rax,%rdx
add    %rcx,%rdx
xorps  %xmm0,%xmm0
movaps %xmm0,(%rax,%rcx,1)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1383 <polydplm+0x253>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x18(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jg     1404 <polydplm+0x2d4>
xor    %eax,%eax
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
shl    $0x4,%rcx
mov    %rax,%rdx
add    %rcx,%rdx
xorps  %xmm0,%xmm0
movaps %xmm0,(%rax,%rcx,1)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    13c7 <polydplm+0x297>
mov    -0x14(%rbp),%ecx
mov    -0x1c(%rbp),%edx
mov    %ecx,%eax
sub    %edx,%eax
add    $0x1,%eax
mov    %eax,-0xc4(%rbp)
mov    %ecx,%eax
lea    0x3(%rax,%rax,1),%edi
xor    %eax,%eax
mov    %al,-0xe5(%rbp)
call   1030 <gammaq@plt>
mov    %eax,%ecx
mov    -0xe5(%rbp),%al
mov    %ecx,-0xec(%rbp)
mov    -0x14(%rbp),%edi
add    $0x2,%edi
call   1030 <gammaq@plt>
mov    -0xec(%rbp),%ecx
mov    %eax,%edx
mov    -0xe5(%rbp),%al
sub    %edx,%ecx
mov    %ecx,-0xe4(%rbp)
mov    -0x14(%rbp),%edi
mov    -0x1c(%rbp),%ecx
sub    %ecx,%edi
add    $0x2,%edi
call   1030 <gammaq@plt>
mov    -0xe4(%rbp),%edi
sub    %eax,%edi
call   5800 <__floatsitf>
movaps %xmm0,-0xe0(%rbp)
mov    -0x14(%rbp),%edi
mov    -0x1c(%rbp),%eax
sub    %eax,%edi
add    $0x1,%edi
call   5800 <__floatsitf>
movaps -0x50(%rbp),%xmm1
call   39f0 <__multf3>
movaps %xmm0,%xmm1
movaps -0xe0(%rbp),%xmm0
call   44c0 <__subtf3>
mov    $0x1,%al
call   1050 <expq@plt>
mov    -0xc4(%rbp),%edi
imul   %eax,%edi
call   5800 <__floatsitf>
movaps %xmm0,%xmm1
mov    -0x8(%rbp),%rax
mov    %rax,-0xb8(%rbp)
mov    -0x14(%rbp),%ecx
mov    -0x1c(%rbp),%edx
sub    %edx,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
shl    $0x4,%rcx
mov    %rcx,-0xc0(%rbp)
mov    %rax,%rdx
add    %rcx,%rdx
movaps (%rax,%rcx,1),%xmm0
call   1ac0 <__addtf3>
mov    -0xc0(%rbp),%rcx
mov    -0xb8(%rbp),%rax
movaps %xmm0,(%rax,%rcx,1)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
cmp    %ecx,%eax
jg     174e <polydplm+0x61e>
mov    -0x34(%rbp),%edx
mov    -0x14(%rbp),%eax
add    %eax,%edx
mov    -0x1c(%rbp),%ecx
mov    %ecx,%eax
mov    %edx,%ecx
lea    0x1(%rax,%rcx,1),%edi
xor    %eax,%eax
mov    %al,-0x175(%rbp)
call   1030 <gammaq@plt>
mov    %eax,%ecx
mov    -0x175(%rbp),%al
mov    %ecx,-0x180(%rbp)
mov    -0x34(%rbp),%esi
mov    -0x1c(%rbp),%ecx
mov    %ecx,%edx
mov    %esi,%ecx
lea    0x1(%rcx,%rdx,1),%edi
call   1030 <gammaq@plt>
mov    -0x180(%rbp),%ecx
mov    %eax,%edx
mov    -0x175(%rbp),%al
sub    %edx,%ecx
mov    %ecx,-0x17c(%rbp)
mov    -0x34(%rbp),%edi
add    $0x1,%edi
call   1030 <gammaq@plt>
mov    -0x17c(%rbp),%ecx
mov    %eax,%edx
mov    -0x175(%rbp),%al
sub    %edx,%ecx
mov    %ecx,-0x174(%rbp)
mov    -0x14(%rbp),%edi
mov    -0x34(%rbp),%ecx
sub    %ecx,%edi
mov    -0x1c(%rbp),%ecx
sub    %ecx,%edi
add    $0x1,%edi
call   1030 <gammaq@plt>
mov    -0x174(%rbp),%edi
sub    %eax,%edi
call   5800 <__floatsitf>
movaps %xmm0,-0x170(%rbp)
mov    -0x34(%rbp),%edi
call   5800 <__floatsitf>
movaps -0x50(%rbp),%xmm1
call   39f0 <__multf3>
movaps %xmm0,%xmm1
movaps -0x170(%rbp),%xmm0
call   44c0 <__subtf3>
mov    $0x1,%al
call   1050 <expq@plt>
mov    %eax,%edi
call   5800 <__floatsitf>
movaps %xmm0,-0x60(%rbp)
movaps -0x60(%rbp),%xmm0
movaps %xmm0,-0x160(%rbp)
mov    -0x1c(%rbp),%esi
mov    %esi,-0x14c(%rbp)
mov    %esi,%edx
imul   %edx,%edx
movslq -0x34(%rbp),%rax
mov    %rax,-0x128(%rbp)
mov    %eax,-0x144(%rbp)
mov    -0x14(%rbp),%ecx
mov    %ecx,-0x148(%rbp)
mov    %ecx,%edi
not    %edi
add    %eax,%edi
imul   %edi,%esi
add    %esi,%edx
cvtsi2sd %edx,%xmm0
cvtsi2sd %eax,%xmm1
addsd  %xmm1,%xmm1
cvtsi2sd %ecx,%xmm2
mulsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
add    %eax,%eax
cvtsi2sd %eax,%xmm1
subsd  %xmm1,%xmm0
call   58a0 <__extenddftf2>
movaps %xmm0,%xmm1
movaps -0x160(%rbp),%xmm0
call   39f0 <__multf3>
mov    -0x14c(%rbp),%edx
mov    -0x148(%rbp),%ecx
mov    -0x144(%rbp),%eax
movaps %xmm0,-0x140(%rbp)
sub    %ecx,%edx
mov    %eax,%ecx
mov    %edx,%eax
lea    -0x1(%rax,%rcx,1),%edi
call   5800 <__floatsitf>
movaps %xmm0,%xmm1
movaps -0x140(%rbp),%xmm0
call   2ea0 <__divtf3>
mov    -0x128(%rbp),%rcx
movaps %xmm0,%xmm1
mov    -0x8(%rbp),%rax
mov    %rax,-0x118(%rbp)
shl    $0x4,%rcx
mov    %rcx,-0x120(%rbp)
movaps (%rax,%rcx,1),%xmm0
call   1ac0 <__addtf3>
mov    -0x120(%rbp),%rcx
mov    -0x118(%rbp),%rax
movaps %xmm0,(%rax,%rcx,1)
movaps -0x60(%rbp),%xmm0
movaps %xmm0,-0x110(%rbp)
mov    -0x14(%rbp),%edi
add    $0x1,%edi
call   5800 <__floatsitf>
movaps %xmm0,%xmm1
movaps -0x110(%rbp),%xmm0
call   39f0 <__multf3>
movaps -0x30(%rbp),%xmm1
call   2ea0 <__divtf3>
movaps %xmm0,%xmm1
mov    -0x8(%rbp),%rcx
mov    %rcx,-0x100(%rbp)
movslq -0x34(%rbp),%rax
shl    $0x4,%rax
mov    %rax,-0xf8(%rbp)
lea    0x10(%rax,%rcx,1),%rdx
movaps 0x10(%rax,%rcx,1),%xmm0
call   44c0 <__subtf3>
mov    -0x100(%rbp),%rcx
mov    -0xf8(%rbp),%rax
movaps %xmm0,0x10(%rax,%rcx,1)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1509 <polydplm+0x3d9>
mov    -0x18(%rbp),%ecx
mov    -0x1c(%rbp),%edx
mov    %ecx,%eax
sub    %edx,%eax
add    $0x1,%eax
mov    %eax,-0x194(%rbp)
mov    %ecx,%eax
lea    0x3(%rax,%rax,1),%edi
xor    %eax,%eax
mov    %al,-0x1b5(%rbp)
call   1030 <gammaq@plt>
mov    %eax,%ecx
mov    -0x1b5(%rbp),%al
mov    %ecx,-0x1bc(%rbp)
mov    -0x18(%rbp),%edi
add    $0x2,%edi
call   1030 <gammaq@plt>
mov    -0x1bc(%rbp),%ecx
mov    %eax,%edx
mov    -0x1b5(%rbp),%al
sub    %edx,%ecx
mov    %ecx,-0x1b4(%rbp)
mov    -0x18(%rbp),%edi
mov    -0x1c(%rbp),%ecx
sub    %ecx,%edi
add    $0x2,%edi
call   1030 <gammaq@plt>
mov    -0x1b4(%rbp),%edi
sub    %eax,%edi
call   5800 <__floatsitf>
movaps %xmm0,-0x1b0(%rbp)
mov    -0x18(%rbp),%edi
mov    -0x1c(%rbp),%eax
sub    %eax,%edi
add    $0x1,%edi
call   5800 <__floatsitf>
movaps -0x50(%rbp),%xmm1
call   39f0 <__multf3>
movaps %xmm0,%xmm1
movaps -0x1b0(%rbp),%xmm0
call   44c0 <__subtf3>
mov    $0x1,%al
call   1050 <expq@plt>
mov    -0x194(%rbp),%edi
imul   %eax,%edi
call   5800 <__floatsitf>
movaps %xmm0,%xmm1
mov    -0x10(%rbp),%rax
mov    %rax,-0x188(%rbp)
mov    -0x18(%rbp),%ecx
mov    -0x1c(%rbp),%edx
sub    %edx,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
shl    $0x4,%rcx
mov    %rcx,-0x190(%rbp)
mov    %rax,%rdx
add    %rcx,%rdx
movaps (%rax,%rcx,1),%xmm0
call   1ac0 <__addtf3>
mov    -0x190(%rbp),%rcx
mov    -0x188(%rbp),%rax
movaps %xmm0,(%rax,%rcx,1)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x18(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
cmp    %ecx,%eax
jg     1aa3 <polydplm+0x973>
mov    -0x34(%rbp),%edx
mov    -0x18(%rbp),%eax
add    %eax,%edx
mov    -0x1c(%rbp),%ecx
mov    %ecx,%eax
mov    %edx,%ecx
lea    0x1(%rax,%rcx,1),%edi
xor    %eax,%eax
mov    %al,-0x245(%rbp)
call   1030 <gammaq@plt>
mov    %eax,%ecx
mov    -0x245(%rbp),%al
mov    %ecx,-0x250(%rbp)
mov    -0x34(%rbp),%esi
mov    -0x1c(%rbp),%ecx
mov    %ecx,%edx
mov    %esi,%ecx
lea    0x1(%rcx,%rdx,1),%edi
call   1030 <gammaq@plt>
mov    -0x250(%rbp),%ecx
mov    %eax,%edx
mov    -0x245(%rbp),%al
sub    %edx,%ecx
mov    %ecx,-0x24c(%rbp)
mov    -0x34(%rbp),%edi
add    $0x1,%edi
call   1030 <gammaq@plt>
mov    -0x24c(%rbp),%ecx
mov    %eax,%edx
mov    -0x245(%rbp),%al
sub    %edx,%ecx
mov    %ecx,-0x244(%rbp)
mov    -0x18(%rbp),%edi
mov    -0x34(%rbp),%ecx
sub    %ecx,%edi
mov    -0x1c(%rbp),%ecx
sub    %ecx,%edi
add    $0x1,%edi
call   1030 <gammaq@plt>
mov    -0x244(%rbp),%edi
sub    %eax,%edi
call   5800 <__floatsitf>
movaps %xmm0,-0x240(%rbp)
mov    -0x34(%rbp),%edi
call   5800 <__floatsitf>
movaps -0x50(%rbp),%xmm1
call   39f0 <__multf3>
movaps %xmm0,%xmm1
movaps -0x240(%rbp),%xmm0
call   44c0 <__subtf3>
mov    $0x1,%al
call   1050 <expq@plt>
mov    %eax,%edi
call   5800 <__floatsitf>
movaps %xmm0,-0x70(%rbp)
movaps -0x70(%rbp),%xmm0
movaps %xmm0,-0x230(%rbp)
mov    -0x1c(%rbp),%esi
mov    %esi,-0x21c(%rbp)
mov    %esi,%edx
imul   %edx,%edx
movslq -0x34(%rbp),%rax
mov    %rax,-0x1f8(%rbp)
mov    %eax,-0x214(%rbp)
mov    -0x18(%rbp),%ecx
mov    %ecx,-0x218(%rbp)
mov    %ecx,%edi
not    %edi
add    %eax,%edi
imul   %edi,%esi
add    %esi,%edx
cvtsi2sd %edx,%xmm0
cvtsi2sd %eax,%xmm1
addsd  %xmm1,%xmm1
cvtsi2sd %ecx,%xmm2
mulsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
add    %eax,%eax
cvtsi2sd %eax,%xmm1
subsd  %xmm1,%xmm0
call   58a0 <__extenddftf2>
movaps %xmm0,%xmm1
movaps -0x230(%rbp),%xmm0
call   39f0 <__multf3>
mov    -0x21c(%rbp),%edx
mov    -0x218(%rbp),%ecx
mov    -0x214(%rbp),%eax
movaps %xmm0,-0x210(%rbp)
sub    %ecx,%edx
mov    %eax,%ecx
mov    %edx,%eax
lea    -0x1(%rax,%rcx,1),%edi
call   5800 <__floatsitf>
movaps %xmm0,%xmm1
movaps -0x210(%rbp),%xmm0
call   2ea0 <__divtf3>
mov    -0x1f8(%rbp),%rcx
movaps %xmm0,%xmm1
mov    -0x10(%rbp),%rax
mov    %rax,-0x1e8(%rbp)
shl    $0x4,%rcx
mov    %rcx,-0x1f0(%rbp)
movaps (%rax,%rcx,1),%xmm0
call   1ac0 <__addtf3>
mov    -0x1f0(%rbp),%rcx
mov    -0x1e8(%rbp),%rax
movaps %xmm0,(%rax,%rcx,1)
movaps -0x70(%rbp),%xmm0
movaps %xmm0,-0x1e0(%rbp)
cvtsi2sdl -0x18(%rbp),%xmm0
movsd  0x45cb(%rip),%xmm1
addsd  %xmm1,%xmm0
call   58a0 <__extenddftf2>
movaps %xmm0,%xmm1
movaps -0x1e0(%rbp),%xmm0
call   39f0 <__multf3>
movaps -0x30(%rbp),%xmm1
call   2ea0 <__divtf3>
movaps %xmm0,%xmm1
mov    -0x10(%rbp),%rcx
mov    %rcx,-0x1d0(%rbp)
movslq -0x34(%rbp),%rax
shl    $0x4,%rax
mov    %rax,-0x1c8(%rbp)
lea    0x10(%rax,%rcx,1),%rdx
movaps 0x10(%rax,%rcx,1),%xmm0
call   44c0 <__subtf3>
mov    -0x1d0(%rbp),%rcx
mov    -0x1c8(%rbp),%rax
movaps %xmm0,0x10(%rax,%rcx,1)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1853 <polydplm+0x723>
jmp    1aa8 <polydplm+0x978>
add    $0x250,%rsp
pop    %rbp
ret
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)

x86 assembly (clang, -O2):
<polydplm>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xd8,%rsp
mov    %r8d,%r12d
mov    %ecx,%ebp
mov    %edx,%r13d
mov    %rsi,%r15
mov    %rdi,%rbx
movapd %xmm0,0x70(%rsp)
movapd %xmm0,%xmm1
call   18e0 <__addtf3>
mov    $0x1,%al
call   1050 <logq@plt>
mov    %eax,%edi
call   5620 <__floatsitf>
test   %r12d,%r12d
mov    %r15,0x60(%rsp)
mov    %rbx,0x20(%rsp)
movapd %xmm0,0x80(%rsp)
mov    %rbp,0x40(%rsp)
je     1718 <polydplm+0x5d8>
mov    %r13d,%r14d
sub    %r12d,%r14d
lea    0x1(%r14),%eax
mov    %eax,0x10(%rsp)
cmp    $0xffffffff,%r14d
jl     11ce <polydplm+0x8e>
xor    %edx,%edx
mov    0x10(%rsp),%eax
test   %eax,%eax
cmovg  %eax,%edx
shl    $0x4,%rdx
add    $0x10,%rdx
mov    %rbx,%rdi
xor    %esi,%esi
call   1040 <memset@plt>
mov    %ebp,%eax
sub    %r12d,%eax
lea    0x1(%rax),%ebp
mov    %rax,0xc8(%rsp)
cmp    $0xffffffff,%eax
jl     11fc <polydplm+0xbc>
xor    %edx,%edx
test   %ebp,%ebp
cmovg  %ebp,%edx
shl    $0x4,%rdx
add    $0x10,%rdx
mov    %r15,%rdi
xor    %esi,%esi
call   1040 <memset@plt>
mov    %ebp,0x6c(%rsp)
lea    0x3(%r13,%r13,1),%edi
xor    %eax,%eax
call   1030 <gammaq@plt>
mov    %eax,%ebx
lea    0x2(%r13),%edi
xor    %eax,%eax
call   1030 <gammaq@plt>
mov    %eax,%ebp
lea    0x2(%r14),%edi
xor    %eax,%eax
call   1030 <gammaq@plt>
add    %ebp,%eax
sub    %eax,%ebx
mov    %ebx,%edi
call   5620 <__floatsitf>
movapd %xmm0,0x90(%rsp)
lea    0x1(%r13),%ebp
mov    %ebp,0xc(%rsp)
sub    %r12d,%ebp
mov    %ebp,%edi
call   5620 <__floatsitf>
movaps 0x80(%rsp),%xmm1
call   3810 <__multf3>
movapd %xmm0,%xmm1
movaps 0x90(%rsp),%xmm0
call   42e0 <__subtf3>
mov    $0x1,%al
call   1060 <expq@plt>
mov    0x10(%rsp),%edi
imul   %eax,%edi
call   5620 <__floatsitf>
movaps %xmm0,%xmm1
mov    %ebp,%r15d
movslq %ebp,%rbp
shl    $0x4,%rbp
mov    0x20(%rsp),%rbx
movapd (%rbx,%rbp,1),%xmm0
call   18e0 <__addtf3>
movapd %xmm0,(%rbx,%rbp,1)
test   %r14d,%r14d
mov    %r12,0x90(%rsp)
mov    %r12,%rbp
js     148c <polydplm+0x34c>
xorps  %xmm0,%xmm0
cvtsi2sd %r13d,%xmm0
movsd  %xmm0,0xb0(%rsp)
mov    0xc(%rsp),%edi
call   5620 <__floatsitf>
movapd %xmm0,0xa0(%rsp)
mov    %r15d,%eax
mov    %rax,0x58(%rsp)
mov    %r13d,%eax
not    %eax
add    %ebp,%eax
mov    %rax,0x50(%rsp)
imul   %ebp,%eax
mov    %eax,0xc(%rsp)
lea    0x1(%rbp),%eax
mov    %rax,0x48(%rsp)
lea    0x1(%rbp,%r13,1),%eax
mov    %rax,0xd0(%rsp)
xor    %r14d,%r14d
xor    %ebx,%ebx
nopw   0x0(%rax,%rax,1)
lea    0x1(%rbx),%r12
mov    %r12,0x28(%rsp)
mov    0x48(%rsp),%rax
lea    (%rax,%rbx,1),%r13d
mov    0xd0(%rsp),%rax
lea    (%rax,%rbx,1),%edi
xor    %eax,%eax
call   1030 <gammaq@plt>
mov    %eax,%ebp
mov    %r13d,%edi
xor    %eax,%eax
call   1030 <gammaq@plt>
mov    %eax,%r13d
mov    %r12d,%edi
xor    %eax,%eax
call   1030 <gammaq@plt>
mov    %eax,%r12d
mov    %r15d,%edi
xor    %eax,%eax
call   1030 <gammaq@plt>
add    %r13d,%r12d
mov    0x20(%rsp),%r13
add    %eax,%r12d
sub    %r12d,%ebp
mov    %ebp,%edi
call   5620 <__floatsitf>
movapd %xmm0,0x10(%rsp)
mov    %ebx,%edi
call   5620 <__floatsitf>
movaps 0x80(%rsp),%xmm1
call   3810 <__multf3>
movapd %xmm0,%xmm1
movaps 0x10(%rsp),%xmm0
call   42e0 <__subtf3>
mov    $0x1,%al
call   1060 <expq@plt>
mov    %eax,%edi
call   5620 <__floatsitf>
movaps %xmm0,0x10(%rsp)
mov    0x50(%rsp),%rax
lea    (%rax,%rbx,1),%ebp
mov    0xc(%rsp),%r12d
xorps  %xmm1,%xmm1
cvtsi2sd %r12d,%xmm1
xorps  %xmm0,%xmm0
cvtsi2sd %ebx,%xmm0
mulsd  0x4c32(%rip),%xmm0
mulsd  0xb0(%rsp),%xmm0
addsd  %xmm1,%xmm0
xorps  %xmm1,%xmm1
cvtsi2sd %r14d,%xmm1
subsd  %xmm1,%xmm0
call   56c0 <__extenddftf2>
movaps 0x10(%rsp),%xmm1
call   3810 <__multf3>
movapd %xmm0,0x30(%rsp)
mov    %ebp,%edi
mov    0x90(%rsp),%rbp
call   5620 <__floatsitf>
movapd %xmm0,%xmm1
movaps 0x30(%rsp),%xmm0
call   2cc0 <__divtf3>
movaps %xmm0,%xmm1
movaps 0x0(%r13,%r14,8),%xmm0
movaps 0x10(%r13,%r14,8),%xmm2
movaps %xmm2,0x30(%rsp)
call   18e0 <__addtf3>
movaps %xmm0,0x0(%r13,%r14,8)
movaps 0xa0(%rsp),%xmm0
movaps 0x10(%rsp),%xmm1
call   3810 <__multf3>
movaps 0x70(%rsp),%xmm1
call   2cc0 <__divtf3>
movaps %xmm0,%xmm1
movapd 0x30(%rsp),%xmm0
call   42e0 <__subtf3>
movapd %xmm0,0x10(%r13,%r14,8)
add    $0x2,%r14
add    %ebp,%r12d
mov    %r12d,0xc(%rsp)
add    $0xffffffff,%r15d
mov    0x28(%rsp),%rax
mov    %rax,%rbx
cmp    %rax,0x58(%rsp)
jne    1310 <polydplm+0x1d0>
mov    %rbp,%r12
mov    0x40(%rsp),%r14
lea    (%r14,%r14,1),%edi
add    $0x3,%edi
xor    %eax,%eax
call   1030 <gammaq@plt>
mov    %eax,%ebx
lea    0x2(%r14),%edi
xor    %eax,%eax
call   1030 <gammaq@plt>
mov    %eax,%ebp
mov    0xc8(%rsp),%r15
lea    0x2(%r15),%edi
xor    %eax,%eax
call   1030 <gammaq@plt>
add    %ebp,%eax
sub    %eax,%ebx
mov    %ebx,%edi
call   5620 <__floatsitf>
movapd %xmm0,0x10(%rsp)
mov    0x6c(%rsp),%ebx
mov    %ebx,%edi
call   5620 <__floatsitf>
movaps 0x80(%rsp),%xmm1
call   3810 <__multf3>
movapd %xmm0,%xmm1
movapd 0x10(%rsp),%xmm0
call   42e0 <__subtf3>
mov    $0x1,%al
call   1060 <expq@plt>
imul   %ebx,%eax
mov    %eax,%edi
call   5620 <__floatsitf>
movapd %xmm0,%xmm1
mov    %ebx,%r13d
movslq %ebx,%rbx
shl    $0x4,%rbx
mov    0x60(%rsp),%rbp
movapd 0x0(%rbp,%rbx,1),%xmm0
call   18e0 <__addtf3>
movapd %xmm0,0x0(%rbp,%rbx,1)
test   %r15d,%r15d
js     18cb <polydplm+0x78b>
xorps  %xmm1,%xmm1
cvtsi2sd %r14d,%xmm1
movsd  0x4abe(%rip),%xmm0
movsd  %xmm1,0x20(%rsp)
addsd  %xmm1,%xmm0
call   56c0 <__extenddftf2>
movapd %xmm0,0xb0(%rsp)
mov    %r13d,%eax
mov    %rax,0xa0(%rsp)
mov    %r14d,%eax
not    %eax
add    %r12d,%eax
mov    %rax,0x58(%rsp)
imul   %r12d,%eax
mov    %eax,0xc(%rsp)
lea    0x1(%r12),%eax
mov    %rax,0x50(%rsp)
lea    (%r12,%r14,1),%eax
add    $0x1,%eax
mov    %rax,0x48(%rsp)
xor    %r14d,%r14d
xor    %ebx,%ebx
nopl   (%rax)
lea    0x1(%rbx),%rbp
mov    %rbp,0x28(%rsp)
mov    0x50(%rsp),%rax
lea    (%rax,%rbx,1),%r15d
mov    0x48(%rsp),%rax
lea    (%rax,%rbx,1),%edi
xor    %eax,%eax
call   1030 <gammaq@plt>
mov    %eax,%r12d
mov    %r15d,%edi
xor    %eax,%eax
call   1030 <gammaq@plt>
mov    %eax,%r15d
mov    %ebp,%edi
xor    %eax,%eax
call   1030 <gammaq@plt>
mov    %eax,%ebp
mov    %r13d,%edi
xor    %eax,%eax
call   1030 <gammaq@plt>
add    %r15d,%ebp
add    %eax,%ebp
sub    %ebp,%r12d
mov    %r12d,%edi
mov    0x60(%rsp),%r12
call   5620 <__floatsitf>
movapd %xmm0,0x10(%rsp)
mov    %ebx,%edi
call   5620 <__floatsitf>
movaps 0x80(%rsp),%xmm1
call   3810 <__multf3>
movapd %xmm0,%xmm1
movaps 0x10(%rsp),%xmm0
call   42e0 <__subtf3>
mov    $0x1,%al
call   1060 <expq@plt>
mov    %eax,%edi
call   5620 <__floatsitf>
movaps %xmm0,0x10(%rsp)
mov    0x58(%rsp),%rax
lea    (%rax,%rbx,1),%r15d
mov    0xc(%rsp),%ebp
xorps  %xmm1,%xmm1
cvtsi2sd %ebp,%xmm1
xorps  %xmm0,%xmm0
cvtsi2sd %ebx,%xmm0
mulsd  0x49a7(%rip),%xmm0
mulsd  0x20(%rsp),%xmm0
addsd  %xmm1,%xmm0
xorps  %xmm1,%xmm1
cvtsi2sd %r14d,%xmm1
subsd  %xmm1,%xmm0
call   56c0 <__extenddftf2>
movaps 0x10(%rsp),%xmm1
call   3810 <__multf3>
movapd %xmm0,0x30(%rsp)
mov    %r15d,%edi
call   5620 <__floatsitf>
movapd %xmm0,%xmm1
movaps 0x30(%rsp),%xmm0
call   2cc0 <__divtf3>
movaps %xmm0,%xmm1
movaps (%r12,%r14,8),%xmm0
movaps 0x10(%r12,%r14,8),%xmm2
movaps %xmm2,0x30(%rsp)
call   18e0 <__addtf3>
movaps %xmm0,(%r12,%r14,8)
movapd 0xb0(%rsp),%xmm0
movaps 0x10(%rsp),%xmm1
call   3810 <__multf3>
movaps 0x70(%rsp),%xmm1
call   2cc0 <__divtf3>
movapd %xmm0,%xmm1
movapd 0x30(%rsp),%xmm0
call   42e0 <__subtf3>
movapd %xmm0,0x10(%r12,%r14,8)
add    $0x2,%r14
add    0x90(%rsp),%ebp
mov    %ebp,0xc(%rsp)
add    $0xffffffff,%r13d
mov    0x28(%rsp),%rax
mov    %rax,%rbx
cmp    %rax,0xa0(%rsp)
jne    15a0 <polydplm+0x460>
jmp    18cb <polydplm+0x78b>
test   %r13d,%r13d
jle    17ed <polydplm+0x6ad>
mov    %r13d,%eax
shl    $0x4,%rax
mov    %rax,0x70(%rsp)
mov    $0x2,%ebp
xor    %r15d,%r15d
mov    %r13d,%r14d
mov    %r13,0x28(%rsp)
nopl   (%rax)
lea    0x0(%r13,%rbp,1),%edi
xor    %eax,%eax
call   1030 <gammaq@plt>
mov    %eax,%ebx
mov    %ebp,%edi
xor    %eax,%eax
call   1030 <gammaq@plt>
mov    %eax,0x10(%rsp)
lea    -0x1(%rbp),%r12d
mov    %r14d,%edi
xor    %eax,%eax
call   1030 <gammaq@plt>
mov    %eax,%r13d
mov    %r12d,%edi
xor    %eax,%eax
call   1030 <gammaq@plt>
add    0x10(%rsp),%r13d
add    %eax,%r13d
sub    %r13d,%ebx
mov    0x28(%rsp),%r13
mov    %ebx,%edi
call   5620 <__floatsitf>
movapd %xmm0,0x10(%rsp)
lea    -0x2(%rbp),%edi
call   5620 <__floatsitf>
movaps 0x80(%rsp),%xmm1
call   3810 <__multf3>
movapd %xmm0,%xmm1
movapd 0x10(%rsp),%xmm0
call   42e0 <__subtf3>
mov    $0x1,%al
call   1060 <expq@plt>
mov    %eax,%edi
shr    $0x1f,%edi
add    %eax,%edi
sar    %edi
call   5620 <__floatsitf>
mov    0x20(%rsp),%rax
movapd %xmm0,(%rax,%r15,1)
add    $0x1,%ebp
add    $0xffffffff,%r14d
add    $0x10,%r15
cmp    %r15,0x70(%rsp)
jne    1740 <polydplm+0x600>
mov    0x40(%rsp),%r12
test   %r12d,%r12d
jle    18cb <polydplm+0x78b>
mov    %r12d,%eax
shl    $0x4,%rax
mov    %rax,0x70(%rsp)
mov    $0x2,%ebx
xor    %r15d,%r15d
mov    %r12d,%r14d
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
lea    (%r12,%rbx,1),%edi
xor    %eax,%eax
call   1030 <gammaq@plt>
mov    %eax,%r13d
mov    %ebx,%edi
xor    %eax,%eax
call   1030 <gammaq@plt>
mov    %eax,0x10(%rsp)
lea    -0x1(%rbx),%r12d
mov    %r14d,%edi
xor    %eax,%eax
call   1030 <gammaq@plt>
mov    %eax,%ebp
mov    %r12d,%edi
mov    0x40(%rsp),%r12
xor    %eax,%eax
call   1030 <gammaq@plt>
add    0x10(%rsp),%ebp
add    %eax,%ebp
sub    %ebp,%r13d
mov    %r13d,%edi
call   5620 <__floatsitf>
movapd %xmm0,0x10(%rsp)
lea    -0x2(%rbx),%edi
call   5620 <__floatsitf>
movaps 0x80(%rsp),%xmm1
call   3810 <__multf3>
movapd %xmm0,%xmm1
movapd 0x10(%rsp),%xmm0
call   42e0 <__subtf3>
mov    $0x1,%al
call   1060 <expq@plt>
mov    %eax,%edi
shr    $0x1f,%edi
add    %eax,%edi
sar    %edi
call   5620 <__floatsitf>
mov    0x60(%rsp),%rax
movapd %xmm0,(%rax,%r15,1)
add    $0x1,%ebx
add    $0xffffffff,%r14d
add    $0x10,%r15
cmp    %r15,0x70(%rsp)
jne    1820 <polydplm+0x6e0>
add    $0xd8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
nopl   (%rax)

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

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