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

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
double func(double, double, double, double);
float log10(double);
float pow(int, float);
int pade(int nk, double *k, double *Pk, int npts, double *k_out,
         double *Pk_out) {
  int i;
  double x0, y0, x1, y1, x2, y2;
  int extra_pts = 51;
  x0 = k[nk - 100];
  y0 = Pk[nk - 100];
  x1 = k[nk - 50];
  y1 = Pk[nk - 50];
  x2 = k[nk - 1];
  y2 = Pk[nk - 1];
  double a0 = 0;
  double b1 = 0;
  double b2 = 0;
  b2 = (y1 - y0) + (y1 - y2) * (y0 * x0 - y1 * x1) / (y1 * x1 - y2 * x2);
  b2 /= ((y0 * x0 * x0 - y1 * x1 * x1) + (y2 * x2 * x2 - y1 * x1 * x1) *
                                             (y0 * x0 - y1 * x1) /
                                             (y1 * x1 - y2 * x2));
  b1 = (y2 - y1) + b2 * (y2 * x2 * x2 - y1 * x1 * x1);
  b1 /= (y1 * x1 - y2 * x2);
  a0 = y0 * (1. + b1 * x0 + b2 * x0 * x0);
  float kmin = -7;
  float kmax = log10(k[0]);
  int nk_out = 0;
  float binwidth = (kmax - kmin) / (float)(extra_pts - 1);
  float n = 0.963;
  float A = Pk[0] / pow(k[0], n);
  for (i = 0; i < extra_pts; i++) {
    double new_k = pow(10., kmin + binwidth * i);
    if (i != (extra_pts - 1)) {
      k_out[nk_out] = new_k;
      Pk_out[nk_out] = A * pow(new_k, n);
      nk_out++;
    }
  }
  for (i = 0; i < nk; i++) {
    k_out[nk_out] = k[i];
    Pk_out[nk_out] = Pk[i];
    nk_out++;
  }
  kmin = log10(k[nk - 1]);
  kmax = 4.;
  binwidth = (kmax - kmin) / (float)(extra_pts - 1);
  for (i = 0; i < 51; i++) {
    double new_k = pow(10., kmin + binwidth * i);
    double pade = func(new_k, a0, b1, b2);
    if (i != 0) {
      k_out[nk_out] = new_k;
      Pk_out[nk_out] = pade;
      nk_out++;
    }
  }
  return (0);
}


x86 assembly (clang, -O0):
<pade>:
push   %rbp
mov    %rsp,%rbp
sub    $0xd0,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movl   $0x33,-0x6c(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x64,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x64,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x32,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x32,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x68(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x78(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x80(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x58(%rbp),%xmm0
subsd  -0x48(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm1
subsd  -0x68(%rbp),%xmm1
movsd  -0x48(%rbp),%xmm2
movsd  -0x40(%rbp),%xmm4
movsd  -0x58(%rbp),%xmm3
mulsd  -0x50(%rbp),%xmm3
movq   %xmm3,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm3
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm1
movsd  -0x58(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm4
movsd  -0x68(%rbp),%xmm3
mulsd  -0x60(%rbp),%xmm3
movq   %xmm3,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm3
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x48(%rbp),%xmm1
mulsd  -0x40(%rbp),%xmm1
movsd  -0x40(%rbp),%xmm2
movsd  -0x58(%rbp),%xmm0
mulsd  -0x50(%rbp),%xmm0
mulsd  -0x50(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x68(%rbp),%xmm0
mulsd  -0x60(%rbp),%xmm0
movsd  -0x60(%rbp),%xmm3
movsd  -0x58(%rbp),%xmm2
mulsd  -0x50(%rbp),%xmm2
mulsd  -0x50(%rbp),%xmm2
movq   %xmm2,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
mulsd  %xmm3,%xmm0
addsd  %xmm2,%xmm0
movsd  -0x48(%rbp),%xmm2
movsd  -0x40(%rbp),%xmm4
movsd  -0x58(%rbp),%xmm3
mulsd  -0x50(%rbp),%xmm3
movq   %xmm3,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm3
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm0
movsd  -0x58(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm4
movsd  -0x68(%rbp),%xmm3
mulsd  -0x60(%rbp),%xmm3
movq   %xmm3,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm3
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
divsd  %xmm2,%xmm0
addsd  %xmm0,%xmm1
movsd  -0x88(%rbp),%xmm0
ff
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x68(%rbp),%xmm1
subsd  -0x58(%rbp),%xmm1
movsd  -0x88(%rbp),%xmm0
ff
movsd  -0x68(%rbp),%xmm2
mulsd  -0x60(%rbp),%xmm2
movsd  -0x60(%rbp),%xmm4
movsd  -0x58(%rbp),%xmm3
mulsd  -0x50(%rbp),%xmm3
mulsd  -0x50(%rbp),%xmm3
movq   %xmm3,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm3
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x58(%rbp),%xmm1
movsd  -0x50(%rbp),%xmm2
movsd  -0x68(%rbp),%xmm0
mulsd  -0x60(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x80(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  -0x80(%rbp),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm1,%xmm2
movsd  0xbcd(%rip),%xmm1
addsd  %xmm1,%xmm2
movsd  -0x88(%rbp),%xmm1
ff
mulsd  -0x40(%rbp),%xmm1
movsd  -0x40(%rbp),%xmm3
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
movss  0xb92(%rip),%xmm0
movss  %xmm0,-0x8c(%rbp)
ff
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
call   1050 <log10@plt>
movss  %xmm0,-0x90(%rbp)
ff
movl   $0x0,-0x94(%rbp)
movss  -0x90(%rbp),%xmm0
ff
subss  -0x8c(%rbp),%xmm0
ff
mov    -0x6c(%rbp),%eax
sub    $0x1,%eax
cvtsi2ss %eax,%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x98(%rbp)
ff
movss  0xb39(%rip),%xmm0
movss  %xmm0,-0x9c(%rbp)
ff
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0xc0(%rbp)
ff
mov    -0x10(%rbp),%rax
cvttsd2si (%rax),%edi
movss  -0x9c(%rbp),%xmm0
ff
call   1040 <pow@plt>
movaps %xmm0,%xmm1
movsd  -0xc0(%rbp),%xmm0
ff
cvtss2sd %xmm1,%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0xa0(%rbp)
ff
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x6c(%rbp),%eax
jge    15ec <pade+0x4bc>
movss  -0x8c(%rbp),%xmm1
ff
movss  -0x98(%rbp),%xmm0
ff
cvtsi2ssl -0x34(%rbp),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    $0xa,%edi
call   1040 <pow@plt>
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
mov    -0x34(%rbp),%eax
mov    -0x6c(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
je     15d9 <pade+0x4a9>
movsd  -0xa8(%rbp),%xmm0
ff
mov    -0x28(%rbp),%rax
movslq -0x94(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movss  -0xa0(%rbp),%xmm0
ff
movss  %xmm0,-0xc4(%rbp)
ff
cvttsd2si -0xa8(%rbp),%edi
ff
movss  -0x9c(%rbp),%xmm0
ff
call   1040 <pow@plt>
movaps %xmm0,%xmm1
movss  -0xc4(%rbp),%xmm0
ff
mulss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
mov    -0x30(%rbp),%rax
movslq -0x94(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x94(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x94(%rbp)
jmp    15de <pade+0x4ae>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    151a <pade+0x3ea>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1656 <pade+0x526>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x94(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
movslq -0x94(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x94(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x94(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    15f3 <pade+0x4c3>
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1050 <log10@plt>
movss  %xmm0,-0x8c(%rbp)
ff
movss  0x98b(%rip),%xmm0
movss  %xmm0,-0x90(%rbp)
ff
movss  -0x90(%rbp),%xmm0
ff
subss  -0x8c(%rbp),%xmm0
ff
mov    -0x6c(%rbp),%eax
sub    $0x1,%eax
cvtsi2ss %eax,%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x98(%rbp)
ff
movl   $0x0,-0x34(%rbp)
cmpl   $0x33,-0x34(%rbp)
jge    1772 <pade+0x642>
movss  -0x8c(%rbp),%xmm1
ff
movss  -0x98(%rbp),%xmm0
ff
cvtsi2ssl -0x34(%rbp),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    $0xa,%edi
call   1040 <pow@plt>
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0xb0(%rbp),%xmm0
ff
movsd  -0x78(%rbp),%xmm1
movsd  -0x80(%rbp),%xmm2
movsd  -0x88(%rbp),%xmm3
ff
call   1030 <func@plt>
movsd  %xmm0,-0xb8(%rbp)
ff
cmpl   $0x0,-0x34(%rbp)
je     175f <pade+0x62f>
movsd  -0xb0(%rbp),%xmm0
ff
mov    -0x28(%rbp),%rax
movslq -0x94(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0xb8(%rbp),%xmm0
ff
mov    -0x30(%rbp),%rax
movslq -0x94(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x94(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x94(%rbp)
jmp    1764 <pade+0x634>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    16b2 <pade+0x582>
xor    %eax,%eax
add    $0xd0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<pade>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x68,%rsp
mov    %r9,%r12
mov    %r8,%r15
mov    %rdx,%rbp
mov    %rsi,%rbx
mov    %edi,0x24(%rsp)
movslq %edi,%rax
movsd  -0x320(%rdx,%rax,8),%xmm3
ff ff
movsd  -0x190(%rdx,%rax,8),%xmm4
ff ff
movsd  -0x320(%rsi,%rax,8),%xmm10
fc ff ff
movsd  -0x190(%rsi,%rax,8),%xmm9
fe ff ff
lea    -0x1(%rax),%rcx
mov    %rcx,0x10(%rsp)
movsd  -0x8(%rsi,%rax,8),%xmm8
movsd  -0x8(%rdx,%rax,8),%xmm0
movapd %xmm4,%xmm2
subsd  %xmm0,%xmm2
movapd 0xe63(%rip),%xmm5
xorpd  %xmm4,%xmm5
movapd %xmm8,%xmm1
mulsd  %xmm0,%xmm1
movapd %xmm4,%xmm6
subsd  %xmm4,%xmm0
subsd  %xmm3,%xmm4
mulsd  %xmm9,%xmm5
movsd  %xmm3,0x30(%rsp)
mulsd  %xmm10,%xmm3
movapd %xmm3,%xmm7
addsd  %xmm5,%xmm7
mulsd  %xmm7,%xmm2
mulsd  %xmm9,%xmm6
subsd  %xmm1,%xmm6
divsd  %xmm6,%xmm2
addsd  %xmm4,%xmm2
mulsd  %xmm9,%xmm5
mulsd  %xmm10,%xmm3
addsd  %xmm5,%xmm3
mulsd  %xmm8,%xmm1
addsd  %xmm5,%xmm1
mulsd  %xmm1,%xmm7
divsd  %xmm6,%xmm7
addsd  %xmm3,%xmm7
divsd  %xmm7,%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
divsd  %xmm6,%xmm1
movsd  %xmm1,0x40(%rsp)
mulsd  %xmm10,%xmm1
addsd  0xde7(%rip),%xmm1
movapd %xmm10,%xmm0
movapd %xmm2,0x50(%rsp)
mulsd  %xmm2,%xmm0
mulsd  %xmm10,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x48(%rsp)
movsd  (%rsi),%xmm0
call   1050 <log10@plt>
addss  0xdc0(%rip),%xmm0
divss  0xdbc(%rip),%xmm0
movss  %xmm0,0x28(%rsp)
mov    %rbp,0x38(%rsp)
movsd  0x0(%rbp),%xmm0
movsd  %xmm0,0x18(%rsp)
mov    %rbx,0x8(%rsp)
cvttsd2si (%rbx),%edi
movss  0xd99(%rip),%xmm0
call   1040 <pow@plt>
cvtss2sd %xmm0,%xmm0
movsd  0x18(%rsp),%xmm1
divsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
movss  %xmm0,0x18(%rsp)
xorps  %xmm0,%xmm0
mulss  0x28(%rsp),%xmm0
addss  0xd6c(%rip),%xmm0
mov    $0xa,%edi
call   1040 <pow@plt>
mov    $0xc,%ebp
xor    %r13d,%r13d
mov    %r15,%rbx
mov    %r12,%r14
xorps  %xmm1,%xmm1
cvtss2sd %xmm0,%xmm1
movsd  %xmm1,(%r15,%r13,8)
cvttss2si %xmm0,%edi
movss  0xd37(%rip),%xmm0
call   1040 <pow@plt>
mulss  0x18(%rsp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,(%r12,%r13,8)
add    $0x1,%r13
xorps  %xmm0,%xmm0
cvtsi2ss %r13d,%xmm0
mulss  0x28(%rsp),%xmm0
addss  0xd0c(%rip),%xmm0
mov    $0xa,%edi
call   1040 <pow@plt>
add    $0x2,%rbp
add    $0x8,%r14
add    $0x8,%rbx
cmp    $0x32,%r13d
jne    12d0 <pade+0x1a0>
movsd  0x30(%rsp),%xmm2
mulsd  0x48(%rsp),%xmm2
mov    0x24(%rsp),%ecx
test   %ecx,%ecx
movsd  %xmm2,0x30(%rsp)
jle    136c <pade+0x23c>
mov    %ecx,%eax
cmp    $0x4,%ecx
jae    137e <pade+0x24e>
xor    %r10d,%r10d
mov    0x8(%rsp),%rbx
mov    0x10(%rsp),%r8
mov    0x38(%rsp),%rcx
jmp    1503 <pade+0x3d3>
mov    %r13d,%ebp
mov    0x8(%rsp),%rbx
mov    0x10(%rsp),%r8
jmp    15ac <pade+0x47c>
lea    (%rbx,%rax,8),%rcx
lea    (%r14,%rax,8),%rdi
mov    0x8(%rsp),%r10
lea    (%r10,%rax,8),%rsi
mov    0x38(%rsp),%r8
lea    (%r8,%rax,8),%rdx
cmp    %rdi,%rbx
setb   0x28(%rsp)
cmp    %rcx,%r14
setb   %r9b
cmp    %rsi,%rbx
setb   %r11b
cmp    %r10,%rcx
seta   0x18(%rsp)
cmp    %rdx,%rbx
setb   %bl
cmp    %r8,%rcx
seta   0x24(%rsp)
cmp    %rsi,%r14
setb   %sil
cmp    %r10,%rdi
mov    %r8,%rcx
seta   %r8b
cmp    %rdx,%r14
setb   %r14b
cmp    %rcx,%rdi
seta   %dl
xor    %r10d,%r10d
test   %r9b,0x28(%rsp)
jne    14b4 <pade+0x384>
and    0x18(%rsp),%r11b
jne    14b4 <pade+0x384>
and    0x24(%rsp),%bl
jne    14b4 <pade+0x384>
and    %r8b,%sil
jne    14b4 <pade+0x384>
and    %dl,%r14b
jne    14b4 <pade+0x384>
mov    %eax,%r10d
and    $0xfffffffc,%r10d
lea    -0x4(%r10),%rdx
mov    %rdx,%rsi
shr    $0x2,%rsi
add    $0x1,%rsi
test   %rdx,%rdx
mov    0x8(%rsp),%rbx
mov    0x10(%rsp),%r8
je     14c0 <pade+0x390>
mov    %rsi,%rdi
and    $0xfffffffffffffffe,%rdi
xor    %edx,%edx
nopl   0x0(%rax)
movups (%rbx,%rdx,8),%xmm0
movups 0x10(%rbx,%rdx,8),%xmm1
movups %xmm0,-0x30(%r15,%rbp,4)
movups %xmm1,-0x20(%r15,%rbp,4)
movups (%rcx,%rdx,8),%xmm0
movups 0x10(%rcx,%rdx,8),%xmm1
movups %xmm0,-0x30(%r12,%rbp,4)
movups %xmm1,-0x20(%r12,%rbp,4)
movups 0x20(%rbx,%rdx,8),%xmm0
movups 0x30(%rbx,%rdx,8),%xmm1
movups %xmm0,-0x10(%r15,%rbp,4)
movups %xmm1,(%r15,%rbp,4)
movups 0x20(%rcx,%rdx,8),%xmm0
movups 0x30(%rcx,%rdx,8),%xmm1
movups %xmm0,-0x10(%r12,%rbp,4)
movups %xmm1,(%r12,%rbp,4)
add    $0x8,%rdx
add    $0x10,%rbp
add    $0xfffffffffffffffe,%rdi
jne    1450 <pade+0x320>
jmp    14c2 <pade+0x392>
mov    0x8(%rsp),%rbx
mov    0x10(%rsp),%r8
jmp    1503 <pade+0x3d3>
xor    %edx,%edx
lea    (%r10,%r13,1),%rbp
test   $0x1,%sil
je     14f7 <pade+0x3c7>
add    %rdx,%r13
movups (%rbx,%rdx,8),%xmm0
movups 0x10(%rbx,%rdx,8),%xmm1
movups %xmm0,(%r15,%r13,8)
movups %xmm1,0x10(%r15,%r13,8)
movups (%rcx,%rdx,8),%xmm0
movups 0x10(%rcx,%rdx,8),%xmm1
movups %xmm0,(%r12,%r13,8)
movups %xmm1,0x10(%r12,%r13,8)
mov    %rbp,%r13
cmp    %rax,%r10
je     15ac <pade+0x47c>
mov    %r10,%rdx
not    %rdx
test   $0x1,%al
jne    1517 <pade+0x3e7>
add    %rax,%rdx
jne    153f <pade+0x40f>
jmp    15ac <pade+0x47c>
movsd  (%rbx,%r10,8),%xmm0
movsd  %xmm0,(%r15,%r13,8)
movsd  (%rcx,%r10,8),%xmm0
movsd  %xmm0,(%r12,%r13,8)
add    $0x1,%r13
or     $0x1,%r10
mov    %r13,%rbp
add    %rax,%rdx
je     15ac <pade+0x47c>
lea    (%r15,%r13,8),%rdx
add    $0x8,%rdx
lea    (%r12,%r13,8),%rsi
add    $0x8,%rsi
sub    %r10,%rax
lea    (%rcx,%r10,8),%rdi
add    $0x8,%rdi
lea    (%rbx,%r10,8),%rbp
add    $0x8,%rbp
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  -0x8(%rbp,%rcx,8),%xmm0
movsd  %xmm0,-0x8(%rdx,%rcx,8)
movsd  -0x8(%rdi,%rcx,8),%xmm0
movsd  %xmm0,-0x8(%rsi,%rcx,8)
movsd  0x0(%rbp,%rcx,8),%xmm0
movsd  %xmm0,(%rdx,%rcx,8)
movsd  (%rdi,%rcx,8),%xmm0
movsd  %xmm0,(%rsi,%rcx,8)
add    $0x2,%rcx
cmp    %rcx,%rax
jne    1570 <pade+0x440>
add    %rcx,%r13
mov    %r13,%rbp
movsd  (%rbx,%r8,8),%xmm0
call   1050 <log10@plt>
movaps %xmm0,%xmm1
movss  0xa66(%rip),%xmm2
subss  %xmm0,%xmm2
divss  0xa4e(%rip),%xmm2
xorps  %xmm0,%xmm0
movss  %xmm2,0x8(%rsp)
mulss  %xmm2,%xmm0
movss  %xmm1,0x18(%rsp)
addss  %xmm1,%xmm0
mov    $0xa,%edi
call   1040 <pow@plt>
cvtss2sd %xmm0,%xmm0
movsd  0x30(%rsp),%xmm1
movsd  0x40(%rsp),%xmm2
movaps 0x50(%rsp),%xmm3
call   1030 <func@plt>
movslq %ebp,%rax
lea    (%r12,%rax,8),%r14
lea    (%r15,%rax,8),%r15
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
lea    0x1(%rbx),%rbp
xorps  %xmm0,%xmm0
cvtsi2ss %ebp,%xmm0
mulss  0x8(%rsp),%xmm0
addss  0x18(%rsp),%xmm0
mov    $0xa,%edi
call   1040 <pow@plt>
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,0x28(%rsp)
movsd  0x30(%rsp),%xmm1
movsd  0x40(%rsp),%xmm2
movaps 0x50(%rsp),%xmm3
call   1030 <func@plt>
movsd  0x28(%rsp),%xmm1
movsd  %xmm1,(%r15,%rbx,8)
movsd  %xmm0,(%r14,%rbx,8)
mov    %rbp,%rbx
cmp    $0x32,%ebp
jne    1620 <pade+0x4f0>
xor    %eax,%eax
add    $0x68,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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