void interp3_table0_complex_per(
    const double *r_ck, const double *i_ck, const int K1, const int K2,
    const int K3, const double *r_h1, const double *i_h1, const double *r_h2,
    const double *i_h2, const double *r_h3, const double *i_h3, const int J1,
    const int J2, const int J3, const int L1, const int L2, const int L3,
    const double *p_tm, const int M, double *r_fm, double *i_fm) {
  int mm;
  {
    const int ncenter1 = floor(J1 * L1 / 2.);
    r_h1 += ncenter1;
    i_h1 += ncenter1;
  }
  {
    const int ncenter2 = floor(J2 * L2 / 2.);
    r_h2 += ncenter2;
    i_h2 += ncenter2;
  }
  {
    const int ncenter3 = floor(J3 * L3 / 2.);
    r_h3 += ncenter3;
    i_h3 += ncenter3;
  }
  for (mm = 0; mm < M; mm++) {
    int jj1, jj2, jj3;
    const double t3 = p_tm[2 * M];
    const double t2 = p_tm[M];
    const double t1 = *p_tm++;
    double sum3r = 0;
    double sum3i = 0;
    const int koff1 = 1 + floor(t1 - J1 / 2.);
    const int koff2 = 1 + floor(t2 - J2 / 2.);
    int k3 = 1 + floor(t3 - J3 / 2.);
    for (jj3 = 0; jj3 < J3; jj3++, k3++) {
      const double p3 = (t3 - k3) * L3;
      const int n3 = floor(p3 + 0.5);
      const double coef3r = r_h3[n3];
      const double coef3i = i_h3[n3];
      const int k3mod = ((k3) - (K3)*floor((k3) / (double)(K3)));
      double sum2r = 0;
      double sum2i = 0;
      int k2 = koff2;
      for (jj2 = 0; jj2 < J2; jj2++, k2++) {
        const double p2 = (t2 - k2) * L2;
        const int n2 = floor(p2 + 0.5);
        const double coef2r = r_h2[n2];
        const double coef2i = i_h2[n2];
        const int k2mod = ((k2) - (K2)*floor((k2) / (double)(K2)));
        const int k23mod = (k3mod * K2 + k2mod) * K1;
        register double sum1r = 0;
        register double sum1i = 0;
        int k1 = koff1;
        for (jj1 = 0; jj1 < J1; jj1++, k1++) {
          const double p1 = (t1 - k1) * L1;
          const int n1 = floor(p1 + 0.5);
          register const double coef1r = r_h1[n1];
          register const double coef1i = i_h1[n1];
          const int k1mod = ((k1) - (K1)*floor((k1) / (double)(K1)));
          const int kk = k23mod + k1mod;
          sum1r += coef1r * r_ck[kk] - coef1i * i_ck[kk];
          sum1i += coef1r * i_ck[kk] + coef1i * r_ck[kk];
        }
        sum2r += coef2r * sum1r - coef2i * sum1i;
        sum2i += coef2r * sum1i + coef2i * sum1r;
      }
      sum3r += coef3r * sum2r - coef3i * sum2i;
      sum3i += coef3r * sum2i + coef3i * sum2r;
    }
    *r_fm++ = sum3r;
    *i_fm++ = sum3i;
  }
}


x86 assembly (clang, -O0):
<interp3_table0_complex_per>:
push   %rbp
mov    %rsp,%rbp
sub    $0x160,%rsp
mov    0x80(%rbp),%rax
mov    0x78(%rbp),%rax
mov    0x70(%rbp),%eax
mov    0x68(%rbp),%rax
mov    0x60(%rbp),%eax
mov    0x58(%rbp),%eax
mov    0x50(%rbp),%eax
mov    0x48(%rbp),%eax
mov    0x40(%rbp),%eax
mov    0x38(%rbp),%eax
mov    0x30(%rbp),%rax
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8d,-0x1c(%rbp)
mov    %r9,-0x28(%rbp)
mov    0x38(%rbp),%eax
imul   0x50(%rbp),%eax
cvtsi2sd %eax,%xmm0
movsd  0xe84(%rip),%xmm1
divsd  %xmm1,%xmm0
call   1030 <floor@plt>
cvttsd2si %xmm0,%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%ecx
mov    -0x28(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x28(%rbp)
mov    -0x30(%rbp),%ecx
mov    0x10(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,0x10(%rbp)
mov    0x40(%rbp),%eax
imul   0x58(%rbp),%eax
cvtsi2sd %eax,%xmm0
movsd  0xe37(%rip),%xmm1
divsd  %xmm1,%xmm0
call   1030 <floor@plt>
cvttsd2si %xmm0,%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%ecx
mov    0x18(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,0x18(%rbp)
mov    -0x34(%rbp),%ecx
mov    0x20(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,0x20(%rbp)
mov    0x48(%rbp),%eax
imul   0x60(%rbp),%eax
cvtsi2sd %eax,%xmm0
movsd  0xdea(%rip),%xmm1
divsd  %xmm1,%xmm0
call   1030 <floor@plt>
cvttsd2si %xmm0,%eax
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%ecx
mov    0x28(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,0x28(%rbp)
mov    -0x38(%rbp),%ecx
mov    0x30(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,0x30(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    0x70(%rbp),%eax
jge    18cb <interp3_table0_complex_per+0x7bb>
mov    0x68(%rbp),%rax
mov    0x70(%rbp),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    0x68(%rbp),%rax
movslq 0x70(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    0x68(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,0x68(%rbp)
movsd  (%rax),%xmm0
movsd  %xmm0,-0x60(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x68(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x60(%rbp),%xmm0
cvtsi2sdl 0x38(%rbp),%xmm1
movsd  0xd3a(%rip),%xmm2
divsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
call   1030 <floor@plt>
movaps %xmm0,%xmm1
movsd  0xd2a(%rip),%xmm0
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x74(%rbp)
movsd  -0x58(%rbp),%xmm0
cvtsi2sdl 0x40(%rbp),%xmm1
movsd  0xd05(%rip),%xmm2
divsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
call   1030 <floor@plt>
movaps %xmm0,%xmm1
movsd  0xcf5(%rip),%xmm0
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x78(%rbp)
movsd  -0x50(%rbp),%xmm0
cvtsi2sdl 0x48(%rbp),%xmm1
movsd  0xcd0(%rip),%xmm2
divsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
call   1030 <floor@plt>
movaps %xmm0,%xmm1
movsd  0xcc0(%rip),%xmm0
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x7c(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    0x48(%rbp),%eax
jge    1887 <interp3_table0_complex_per+0x777>
movsd  -0x50(%rbp),%xmm0
cvtsi2sdl -0x7c(%rbp),%xmm1
subsd  %xmm1,%xmm0
cvtsi2sdl 0x60(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  0xc83(%rip),%xmm0
addsd  -0x88(%rbp),%xmm0
ff
call   1030 <floor@plt>
cvttsd2si %xmm0,%eax
mov    %eax,-0x8c(%rbp)
mov    0x28(%rbp),%rax
movslq -0x8c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
mov    0x30(%rbp),%rax
movslq -0x8c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
cvtsi2sdl -0x7c(%rbp),%xmm0
movsd  %xmm0,-0x138(%rbp)
ff
cvtsi2sdl -0x1c(%rbp),%xmm0
movsd  %xmm0,-0x130(%rbp)
ff
cvtsi2sdl -0x7c(%rbp),%xmm0
cvtsi2sdl -0x1c(%rbp),%xmm1
divsd  %xmm1,%xmm0
call   1030 <floor@plt>
movsd  -0x138(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0x130(%rbp),%xmm0
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0xa4(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
mov    -0x78(%rbp),%eax
mov    %eax,-0xbc(%rbp)
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    0x40(%rbp),%eax
jge    17f5 <interp3_table0_complex_per+0x6e5>
movsd  -0x58(%rbp),%xmm0
cvtsi2sdl -0xbc(%rbp),%xmm1
ff
subsd  %xmm1,%xmm0
cvtsi2sdl 0x58(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
movsd  0xb77(%rip),%xmm0
addsd  -0xc8(%rbp),%xmm0
ff
call   1030 <floor@plt>
cvttsd2si %xmm0,%eax
mov    %eax,-0xcc(%rbp)
mov    0x18(%rbp),%rax
movslq -0xcc(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
mov    0x20(%rbp),%rax
movslq -0xcc(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
cvtsi2sdl -0xbc(%rbp),%xmm0
ff
movsd  %xmm0,-0x148(%rbp)
ff
cvtsi2sdl -0x18(%rbp),%xmm0
movsd  %xmm0,-0x140(%rbp)
ff
cvtsi2sdl -0xbc(%rbp),%xmm0
ff
cvtsi2sdl -0x18(%rbp),%xmm1
divsd  %xmm1,%xmm0
call   1030 <floor@plt>
movsd  -0x148(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0x140(%rbp),%xmm0
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0xe4(%rbp)
mov    -0xa4(%rbp),%eax
imul   -0x18(%rbp),%eax
add    -0xe4(%rbp),%eax
imul   -0x14(%rbp),%eax
mov    %eax,-0xe8(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xf0(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xf8(%rbp)
ff
mov    -0x74(%rbp),%eax
mov    %eax,-0xfc(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    0x38(%rbp),%eax
jge    1751 <interp3_table0_complex_per+0x641>
movsd  -0x60(%rbp),%xmm0
cvtsi2sdl -0xfc(%rbp),%xmm1
ff
subsd  %xmm1,%xmm0
cvtsi2sdl 0x50(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x108(%rbp)
ff
movsd  0xa4b(%rip),%xmm0
addsd  -0x108(%rbp),%xmm0
ff
call   1030 <floor@plt>
cvttsd2si %xmm0,%eax
mov    %eax,-0x10c(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x10c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x118(%rbp)
ff
mov    0x10(%rbp),%rax
movslq -0x10c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x120(%rbp)
ff
cvtsi2sdl -0xfc(%rbp),%xmm0
ff
movsd  %xmm0,-0x158(%rbp)
ff
cvtsi2sdl -0x14(%rbp),%xmm0
movsd  %xmm0,-0x150(%rbp)
ff
cvtsi2sdl -0xfc(%rbp),%xmm0
ff
cvtsi2sdl -0x14(%rbp),%xmm1
divsd  %xmm1,%xmm0
call   1030 <floor@plt>
movsd  -0x158(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0x150(%rbp),%xmm0
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x124(%rbp)
mov    -0xe8(%rbp),%eax
add    -0x124(%rbp),%eax
mov    %eax,-0x128(%rbp)
movsd  -0x118(%rbp),%xmm0
ff
mov    -0x8(%rbp),%rax
movslq -0x128(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x120(%rbp),%xmm1
ff
mov    -0x10(%rbp),%rax
movslq -0x128(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
addsd  -0xf0(%rbp),%xmm0
ff
movsd  %xmm0,-0xf0(%rbp)
ff
movsd  -0x118(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movslq -0x128(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x120(%rbp),%xmm1
ff
mov    -0x8(%rbp),%rax
movslq -0x128(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
addsd  -0xf8(%rbp),%xmm0
ff
movsd  %xmm0,-0xf8(%rbp)
ff
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0xfc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xfc(%rbp)
jmp    158f <interp3_table0_complex_per+0x47f>
movsd  -0xd8(%rbp),%xmm0
ff
movsd  -0xf0(%rbp),%xmm2
ff
movsd  -0xe0(%rbp),%xmm1
ff
mulsd  -0xf8(%rbp),%xmm1
ff
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
addsd  -0xb0(%rbp),%xmm0
ff
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0xd8(%rbp),%xmm0
ff
movsd  -0xf8(%rbp),%xmm2
ff
movsd  -0xe0(%rbp),%xmm1
ff
mulsd  -0xf0(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
addsd  -0xb8(%rbp),%xmm0
ff
movsd  %xmm0,-0xb8(%rbp)
ff
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
mov    -0xbc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xbc(%rbp)
jmp    1463 <interp3_table0_complex_per+0x353>
movsd  -0x98(%rbp),%xmm0
ff
movsd  -0xb0(%rbp),%xmm2
ff
movsd  -0xa0(%rbp),%xmm1
ff
mulsd  -0xb8(%rbp),%xmm1
ff
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
addsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x98(%rbp),%xmm0
ff
movsd  -0xb8(%rbp),%xmm2
ff
movsd  -0xa0(%rbp),%xmm1
ff
mulsd  -0xb0(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
addsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
mov    -0x7c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x7c(%rbp)
jmp    135a <interp3_table0_complex_per+0x24a>
movsd  -0x68(%rbp),%xmm0
mov    0x78(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,0x78(%rbp)
movsd  %xmm0,(%rax)
movsd  -0x70(%rbp),%xmm0
mov    0x80(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,0x80(%rbp)
movsd  %xmm0,(%rax)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1257 <interp3_table0_complex_per+0x147>
add    $0x160,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<interp3_table0_complex_per>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x1c8,%rsp
mov    %r9,(%rsp)
mov    %r8d,%r14d
mov    %ecx,0x8(%rsp)
mov    %edx,0xc(%rsp)
mov    %rsi,0xf0(%rsp)
mov    %rdi,%r12
mov    0x260(%rsp),%ebx
mov    0x250(%rsp),%ebp
mov    0x248(%rsp),%r13d
mov    0x240(%rsp),%r15d
mov    %r15d,%eax
imul   0x228(%rsp),%eax
cvtsi2sd %eax,%xmm0
mulsd  0xe90(%rip),%xmm0
call   1030 <floor@plt>
movsd  %xmm0,0x40(%rsp)
mov    %r13d,%eax
imul   0x230(%rsp),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mulsd  0xe6b(%rip),%xmm0
call   1030 <floor@plt>
movsd  %xmm0,0x10(%rsp)
imul   0x238(%rsp),%ebp
xorps  %xmm0,%xmm0
cvtsi2sd %ebp,%xmm0
mulsd  0xe49(%rip),%xmm0
call   1030 <floor@plt>
movsd  0xe3c(%rip),%xmm1
test   %ebx,%ebx
jle    1805 <interp3_table0_complex_per+0x6f5>
mov    %r15d,%edx
mov    0x270(%rsp),%r15
mov    %r13d,%esi
mov    0x268(%rsp),%r13
cvttsd2si 0x40(%rsp),%eax
cltq
mov    (%rsp),%rcx
lea    (%rcx,%rax,8),%rcx
mov    %rcx,0xe0(%rsp)
shl    $0x3,%rax
add    0x200(%rsp),%rax
mov    %rax,0xe8(%rsp)
cvttsd2si 0x10(%rsp),%eax
cltq
shl    $0x3,%rax
mov    0x208(%rsp),%rcx
add    %rax,%rcx
mov    %rcx,0xd0(%rsp)
add    0x210(%rsp),%rax
mov    %rax,0xd8(%rsp)
cvttsd2si %xmm0,%eax
cltq
shl    $0x3,%rax
mov    0x218(%rsp),%rcx
add    %rax,%rcx
mov    %rcx,0xb0(%rsp)
add    0x220(%rsp),%rax
mov    %rax,0xb8(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sdl 0x250(%rsp),%xmm0
movsd  %xmm0,0xa8(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sd %r14d,%xmm0
mov    0x258(%rsp),%rbp
mov    0x260(%rsp),%eax
mov    %rax,%rcx
add    %eax,%eax
cvtsi2sd %esi,%xmm2
movsd  %xmm2,0xc8(%rsp)
cltq
mov    %rax,0x70(%rsp)
mov    %ecx,%eax
mov    %rax,0x88(%rsp)
xorps  %xmm2,%xmm2
cvtsi2sdl 0x228(%rsp),%xmm2
mulsd  %xmm1,%xmm2
movsd  %xmm2,0x80(%rsp)
xorps  %xmm2,%xmm2
cvtsi2sdl 0x230(%rsp),%xmm2
mulsd  %xmm1,%xmm2
movsd  %xmm2,0x78(%rsp)
xorps  %xmm2,%xmm2
cvtsi2sdl 0x238(%rsp),%xmm2
mulsd  %xmm1,%xmm2
movsd  %xmm2,0x68(%rsp)
movapd 0xd0a(%rip),%xmm1
movapd %xmm0,0x110(%rsp)
xorpd  %xmm1,%xmm0
movapd %xmm0,0x100(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sdl 0x8(%rsp),%xmm0
movapd %xmm0,0x160(%rsp)
xorpd  %xmm1,%xmm0
movapd %xmm0,0x150(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sd %edx,%xmm0
movsd  %xmm0,0x60(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sdl 0xc(%rsp),%xmm0
movapd %xmm0,0x1a0(%rsp)
xorpd  %xmm0,%xmm1
movapd %xmm1,0x1b0(%rsp)
xor    %ebx,%ebx
mov    %r12,0xf8(%rsp)
mov    0xf0(%rsp),%r14
mov    0xe0(%rsp),%r12
jmp    13d9 <interp3_table0_complex_per+0x2c9>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
xorpd  %xmm0,%xmm0
mov    0x90(%rsp),%rbp
add    $0x8,%rbp
mov    0x98(%rsp),%r13
movhpd %xmm0,0x0(%r13)
add    $0x8,%r13
mov    0xa0(%rsp),%r15
movlpd %xmm0,(%r15)
add    $0x8,%r15
mov    0x1c(%rsp),%ebx
add    $0x1,%ebx
mov    0x260(%rsp),%eax
cmp    %eax,%ebx
je     1805 <interp3_table0_complex_per+0x6f5>
mov    0x88(%rsp),%rax
movsd  0x0(%rbp,%rax,8),%xmm0
movsd  %xmm0,0x58(%rsp)
movsd  0x0(%rbp),%xmm0
movsd  %xmm0,0x10(%rsp)
subsd  0x80(%rsp),%xmm0
call   1030 <floor@plt>
movsd  %xmm0,(%rsp)
movsd  0x58(%rsp),%xmm0
subsd  0x78(%rsp),%xmm0
call   1030 <floor@plt>
cmpl   $0x0,0x238(%rsp)
mov    %r15,0xa0(%rsp)
mov    %r13,0x98(%rsp)
mov    %rbp,0x90(%rsp)
mov    %ebx,0x1c(%rsp)
jle    1390 <interp3_table0_complex_per+0x280>
mov    0x70(%rsp),%rax
movsd  0x0(%rbp,%rax,8),%xmm2
movsd  0xbaf(%rip),%xmm3
movsd  (%rsp),%xmm1
addsd  %xmm3,%xmm1
cvttsd2si %xmm1,%eax
mov    %eax,0x2c(%rsp)
addsd  %xmm3,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,0x20(%rsp)
movsd  %xmm2,0xc0(%rsp)
movapd %xmm2,%xmm0
subsd  0x68(%rsp),%xmm0
call   1030 <floor@plt>
addsd  0xb72(%rip),%xmm0
cvttsd2si %xmm0,%eax
xorpd  %xmm0,%xmm0
xor    %ecx,%ecx
movsd  0xb58(%rip),%xmm1
jmp    151f <interp3_table0_complex_per+0x40f>
nopw   0x0(%rax,%rax,1)
movapd 0x120(%rsp),%xmm2
movapd %xmm2,%xmm0
movapd 0xb4b(%rip),%xmm1
xorpd  %xmm1,%xmm0
unpcklpd %xmm3,%xmm2
movapd %xmm3,%xmm1
shufpd $0x1,%xmm0,%xmm1
mulpd  %xmm2,%xmm1
movapd 0x130(%rsp),%xmm2
unpcklpd %xmm2,%xmm2
mulpd  %xmm3,%xmm2
addpd  %xmm1,%xmm2
movapd 0x140(%rsp),%xmm0
addpd  %xmm2,%xmm0
mov    0x24(%rsp),%ecx
add    $0x1,%ecx
mov    0x28(%rsp),%eax
add    $0x1,%eax
cmp    0x238(%rsp),%ecx
movsd  0xae7(%rip),%xmm1
je     1394 <interp3_table0_complex_per+0x284>
mov    %ecx,0x24(%rsp)
movapd %xmm0,0x140(%rsp)
mov    %eax,0x28(%rsp)
xorps  %xmm2,%xmm2
cvtsi2sd %eax,%xmm2
movsd  0xc0(%rsp),%xmm0
movsd  %xmm2,(%rsp)
subsd  %xmm2,%xmm0
mulsd  0xa8(%rsp),%xmm0
addsd  %xmm1,%xmm0
call   1030 <floor@plt>
cvttsd2si %xmm0,%eax
cltq
mov    0xb0(%rsp),%rcx
movsd  (%rcx,%rax,8),%xmm0
movaps %xmm0,0x130(%rsp)
mov    0xb8(%rsp),%rcx
movsd  (%rcx,%rax,8),%xmm0
movapd %xmm0,0x120(%rsp)
xorpd  %xmm3,%xmm3
cmpl   $0x0,0x230(%rsp)
jle    14b0 <interp3_table0_complex_per+0x3a0>
movsd  (%rsp),%xmm0
divsd  0x110(%rsp),%xmm0
call   1030 <floor@plt>
mulsd  0x100(%rsp),%xmm0
addsd  (%rsp),%xmm0
cvttsd2si %xmm0,%eax
imul   0x8(%rsp),%eax
mov    %eax,0x30(%rsp)
xorpd  %xmm3,%xmm3
xor    %eax,%eax
mov    0x20(%rsp),%ecx
jmp    163f <interp3_table0_complex_per+0x52f>
nopl   0x0(%rax,%rax,1)
movapd 0x170(%rsp),%xmm2
movapd %xmm2,%xmm0
xorpd  %xmm4,%xmm0
unpcklpd %xmm3,%xmm2
movapd %xmm3,%xmm1
shufpd $0x1,%xmm0,%xmm1
mulpd  %xmm2,%xmm1
movapd 0x180(%rsp),%xmm0
unpcklpd %xmm0,%xmm0
mulpd  %xmm3,%xmm0
addpd  %xmm1,%xmm0
movapd 0x190(%rsp),%xmm3
addpd  %xmm0,%xmm3
mov    0x38(%rsp),%eax
add    $0x1,%eax
mov    0x34(%rsp),%ecx
add    $0x1,%ecx
cmp    0x230(%rsp),%eax
je     14b0 <interp3_table0_complex_per+0x3a0>
mov    %eax,0x38(%rsp)
movapd %xmm3,0x190(%rsp)
mov    %ecx,0x34(%rsp)
xorps  %xmm1,%xmm1
cvtsi2sd %ecx,%xmm1
movsd  %xmm1,(%rsp)
movsd  0x58(%rsp),%xmm0
subsd  %xmm1,%xmm0
mulsd  0xc8(%rsp),%xmm0
movsd  0x989(%rip),%xmm1
addsd  %xmm1,%xmm0
call   1030 <floor@plt>
movsd  %xmm0,0x40(%rsp)
movsd  (%rsp),%xmm0
divsd  0x160(%rsp),%xmm0
call   1030 <floor@plt>
cvttsd2si 0x40(%rsp),%eax
cltq
mov    0xd0(%rsp),%rcx
movsd  (%rcx,%rax,8),%xmm1
movaps %xmm1,0x180(%rsp)
mov    0xd8(%rsp),%rcx
movsd  (%rcx,%rax,8),%xmm1
movaps %xmm1,0x170(%rsp)
xorpd  %xmm3,%xmm3
cmpl   $0x0,0x228(%rsp)
mov    0xe8(%rsp),%r13
movsd  0x60(%rsp),%xmm2
movapd 0x923(%rip),%xmm4
movsd  0x10(%rsp),%xmm1
jle    15e0 <interp3_table0_complex_per+0x4d0>
mulsd  0x150(%rsp),%xmm0
addsd  (%rsp),%xmm0
cvttsd2si %xmm0,%edi
add    0x30(%rsp),%edi
imul   0xc(%rsp),%edi
mov    %edi,0x3c(%rsp)
xorpd  %xmm3,%xmm3
mov    0x228(%rsp),%r15d
mov    0x2c(%rsp),%ebx
nopl   0x0(%rax,%rax,1)
movapd %xmm3,0x40(%rsp)
xorps  %xmm3,%xmm3
cvtsi2sd %ebx,%xmm3
movsd  %xmm3,(%rsp)
movapd %xmm1,%xmm0
subsd  %xmm3,%xmm0
mulsd  %xmm2,%xmm0
addsd  0x8aa(%rip),%xmm0
call   1030 <floor@plt>
cvttsd2si %xmm0,%eax
movslq %eax,%rbp
movsd  (%rsp),%xmm0
divsd  0x1a0(%rsp),%xmm0
call   1030 <floor@plt>
movapd 0x1b0(%rsp),%xmm5
movapd 0x88a(%rip),%xmm4
mulsd  %xmm5,%xmm0
addsd  (%rsp),%xmm0
cvttsd2si %xmm0,%eax
add    0x3c(%rsp),%eax
cltq
movsd  0x0(%r13,%rbp,8),%xmm0
movapd %xmm0,%xmm1
xorpd  %xmm4,%xmm1
movsd  (%r14,%rax,8),%xmm2
unpcklpd %xmm2,%xmm0
mov    0xf8(%rsp),%rdi
movsd  (%rdi,%rax,8),%xmm3
unpcklpd %xmm3,%xmm2
unpcklpd %xmm1,%xmm3
movsd  0x10(%rsp),%xmm1
mulpd  %xmm0,%xmm3
movsd  (%r12,%rbp,8),%xmm0
unpcklpd %xmm0,%xmm0
mulpd  %xmm0,%xmm2
addpd  %xmm3,%xmm2
movapd 0x40(%rsp),%xmm3
addpd  %xmm2,%xmm3
movsd  0x60(%rsp),%xmm2
add    $0x1,%ebx
add    $0xffffffff,%r15d
jne    1730 <interp3_table0_complex_per+0x620>
jmp    15e0 <interp3_table0_complex_per+0x4d0>
add    $0x1c8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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