void interp2_table0_complex_per(const double *r_ck, const double *i_ck,
                                const int K1, const int K2, const double *r_h1,
                                const double *i_h1, const double *r_h2,
                                const double *i_h2, const int J1, const int J2,
                                const int L1, const int L2, 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;
  }
  for (mm = 0; mm < M; mm++) {
    int jj1, jj2;
    const double t2 = p_tm[M];
    const double t1 = *p_tm++;
    double sum2r = 0;
    double sum2i = 0;
    const int koff1 = 1 + floor(t1 - J1 / 2.);
    int k2 = 1 + floor(t2 - J2 / 2.);
    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 k12mod = 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 = k12mod + 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;
    }
    *r_fm++ = sum2r;
    *i_fm++ = sum2i;
  }
}


x86 assembly (clang, -O0):
<interp2_table0_complex_per>:
push   %rbp
mov    %rsp,%rbp
sub    $0xf0,%rsp
mov    0x58(%rbp),%rax
mov    0x50(%rbp),%rax
mov    0x48(%rbp),%eax
mov    0x40(%rbp),%rax
mov    0x38(%rbp),%eax
mov    0x30(%rbp),%eax
mov    0x28(%rbp),%eax
mov    0x20(%rbp),%eax
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    %r8,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
mov    0x20(%rbp),%eax
imul   0x30(%rbp),%eax
cvtsi2sd %eax,%xmm0
movsd  0xe99(%rip),%xmm1
divsd  %xmm1,%xmm0
call   1030 <floor@plt>
cvttsd2si %xmm0,%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%ecx
mov    -0x20(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x20(%rbp)
mov    -0x30(%rbp),%ecx
mov    -0x28(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x28(%rbp)
mov    0x28(%rbp),%eax
imul   0x38(%rbp),%eax
cvtsi2sd %eax,%xmm0
movsd  0xe4c(%rip),%xmm1
divsd  %xmm1,%xmm0
call   1030 <floor@plt>
cvttsd2si %xmm0,%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%ecx
mov    0x10(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,0x10(%rbp)
mov    -0x34(%rbp),%ecx
mov    0x18(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,0x18(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    0x48(%rbp),%eax
jge    163f <interp2_table0_complex_per+0x52f>
mov    0x40(%rbp),%rax
movslq 0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,0x40(%rbp)
movsd  (%rax),%xmm0
movsd  %xmm0,-0x50(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x50(%rbp),%xmm0
cvtsi2sdl 0x20(%rbp),%xmm1
movsd  0xdb3(%rip),%xmm2
divsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
call   1030 <floor@plt>
movaps %xmm0,%xmm1
movsd  0xda3(%rip),%xmm0
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x64(%rbp)
movsd  -0x48(%rbp),%xmm0
cvtsi2sdl 0x28(%rbp),%xmm1
movsd  0xd7e(%rip),%xmm2
divsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
call   1030 <floor@plt>
movaps %xmm0,%xmm1
movsd  0xd6e(%rip),%xmm0
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x68(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    0x28(%rbp),%eax
jge    1601 <interp2_table0_complex_per+0x4f1>
movsd  -0x48(%rbp),%xmm0
cvtsi2sdl -0x68(%rbp),%xmm1
subsd  %xmm1,%xmm0
cvtsi2sdl 0x38(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  0xd34(%rip),%xmm0
addsd  -0x70(%rbp),%xmm0
call   1030 <floor@plt>
cvttsd2si %xmm0,%eax
mov    %eax,-0x74(%rbp)
mov    0x10(%rbp),%rax
movslq -0x74(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    0x18(%rbp),%rax
movslq -0x74(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
cvtsi2sdl -0x68(%rbp),%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
cvtsi2sdl -0x18(%rbp),%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
cvtsi2sdl -0x68(%rbp),%xmm0
cvtsi2sdl -0x18(%rbp),%xmm1
divsd  %xmm1,%xmm0
call   1030 <floor@plt>
movsd  -0xe0(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0xd8(%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,-0x8c(%rbp)
mov    -0x8c(%rbp),%eax
imul   -0x14(%rbp),%eax
mov    %eax,-0x90(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
mov    -0x64(%rbp),%eax
mov    %eax,-0xa4(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    0x20(%rbp),%eax
jge    1575 <interp2_table0_complex_per+0x465>
movsd  -0x50(%rbp),%xmm0
cvtsi2sdl -0xa4(%rbp),%xmm1
ff
subsd  %xmm1,%xmm0
cvtsi2sdl 0x30(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  0xc27(%rip),%xmm0
addsd  -0xb0(%rbp),%xmm0
ff
call   1030 <floor@plt>
cvttsd2si %xmm0,%eax
mov    %eax,-0xb4(%rbp)
mov    -0x20(%rbp),%rax
movslq -0xb4(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xc0(%rbp)
ff
mov    -0x28(%rbp),%rax
movslq -0xb4(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
cvtsi2sdl -0xa4(%rbp),%xmm0
ff
movsd  %xmm0,-0xf0(%rbp)
ff
cvtsi2sdl -0x14(%rbp),%xmm0
movsd  %xmm0,-0xe8(%rbp)
ff
cvtsi2sdl -0xa4(%rbp),%xmm0
ff
cvtsi2sdl -0x14(%rbp),%xmm1
divsd  %xmm1,%xmm0
call   1030 <floor@plt>
movsd  -0xf0(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0xe8(%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,-0xcc(%rbp)
mov    -0x90(%rbp),%eax
add    -0xcc(%rbp),%eax
mov    %eax,-0xd0(%rbp)
movsd  -0xc0(%rbp),%xmm0
ff
mov    -0x8(%rbp),%rax
movslq -0xd0(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0xc8(%rbp),%xmm1
ff
mov    -0x10(%rbp),%rax
movslq -0xd0(%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  -0x98(%rbp),%xmm0
ff
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0xc0(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movslq -0xd0(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0xc8(%rbp),%xmm1
ff
mov    -0x8(%rbp),%rax
movslq -0xd0(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
addsd  -0xa0(%rbp),%xmm0
ff
movsd  %xmm0,-0xa0(%rbp)
ff
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
mov    -0xa4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xa4(%rbp)
jmp    13b3 <interp2_table0_complex_per+0x2a3>
movsd  -0x80(%rbp),%xmm0
movsd  -0x98(%rbp),%xmm2
ff
movsd  -0x88(%rbp),%xmm1
ff
mulsd  -0xa0(%rbp),%xmm1
ff
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
addsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x80(%rbp),%xmm0
movsd  -0xa0(%rbp),%xmm2
ff
movsd  -0x88(%rbp),%xmm1
ff
mulsd  -0x98(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
addsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x68(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x68(%rbp)
jmp    12ac <interp2_table0_complex_per+0x19c>
movsd  -0x58(%rbp),%xmm0
mov    0x50(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,0x50(%rbp)
movsd  %xmm0,(%rax)
movsd  -0x60(%rbp),%xmm0
mov    0x58(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,0x58(%rbp)
movsd  %xmm0,(%rax)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11f5 <interp2_table0_complex_per+0xe5>
add    $0xf0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<interp2_table0_complex_per>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x128,%rsp
mov    %r9,%r15
mov    %r8,(%rsp)
mov    %ecx,%ebx
mov    %edx,0x8(%rsp)
mov    %rsi,%rbp
mov    %rdi,%r13
mov    0x188(%rsp),%r14d
mov    0x180(%rsp),%r12d
mov    %r12d,%eax
imul   0x170(%rsp),%eax
cvtsi2sd %eax,%xmm0
mulsd  0xea5(%rip),%xmm0
call   1030 <floor@plt>
movsd  %xmm0,0x20(%rsp)
mov    %r14d,%eax
imul   0x178(%rsp),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mulsd  0xe80(%rip),%xmm0
call   1030 <floor@plt>
movsd  0xe73(%rip),%xmm1
cmpl   $0x0,0x198(%rsp)
jle    15b0 <interp2_table0_complex_per+0x4a0>
mov    %r15,%rdi
mov    %rbp,%rdx
mov    0x1a8(%rsp),%r15
mov    %r14d,%esi
mov    0x1a0(%rsp),%r14
mov    0x190(%rsp),%rbp
cvttsd2si 0x20(%rsp),%eax
cltq
mov    (%rsp),%rcx
lea    (%rcx,%rax,8),%rcx
mov    %rcx,0x90(%rsp)
cvttsd2si %xmm0,%ecx
lea    (%rdi,%rax,8),%rax
mov    %rax,0x88(%rsp)
movslq %ecx,%rax
shl    $0x3,%rax
mov    0x160(%rsp),%rcx
add    %rax,%rcx
mov    %rcx,0x78(%rsp)
add    0x168(%rsp),%rax
mov    %rax,0x80(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sdl 0x170(%rsp),%xmm0
cvtsi2sdl 0x178(%rsp),%xmm2
mulsd  %xmm1,%xmm0
movsd  %xmm0,0x50(%rsp)
mulsd  %xmm1,%xmm2
movsd  %xmm2,0x48(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sd %esi,%xmm0
movsd  %xmm0,0x70(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sd %ebx,%xmm0
movapd 0xdc0(%rip),%xmm1
movapd %xmm0,0xc0(%rsp)
xorpd  %xmm1,%xmm0
movapd %xmm0,0xb0(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sd %r12d,%xmm0
movsd  %xmm0,0x38(%rsp)
mov    %rbp,%rbx
xorps  %xmm0,%xmm0
cvtsi2sdl 0x8(%rsp),%xmm0
movapd %xmm0,0x100(%rsp)
xorpd  %xmm0,%xmm1
movapd %xmm1,0x110(%rsp)
xor    %ebp,%ebp
mov    %rdx,0x98(%rsp)
mov    %r13,0xa8(%rsp)
mov    0x88(%rsp),%r12
jmp    12fd <interp2_table0_complex_per+0x1ed>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
xorpd  %xmm0,%xmm0
mov    0x58(%rsp),%rbx
add    $0x8,%rbx
mov    0x60(%rsp),%r14
movhpd %xmm0,(%r14)
add    $0x8,%r14
mov    0x68(%rsp),%r15
movlpd %xmm0,(%r15)
add    $0x8,%r15
mov    0xc(%rsp),%ebp
add    $0x1,%ebp
cmp    0x198(%rsp),%ebp
je     15b0 <interp2_table0_complex_per+0x4a0>
movsd  (%rbx),%xmm0
movsd  %xmm0,0x40(%rsp)
subsd  0x50(%rsp),%xmm0
call   1030 <floor@plt>
movsd  0xcee(%rip),%xmm1
addsd  %xmm1,%xmm0
cmpl   $0x0,0x178(%rsp)
mov    %r15,0x68(%rsp)
mov    %r14,0x60(%rsp)
mov    %rbx,0x58(%rsp)
mov    %ebp,0xc(%rsp)
jle    12c0 <interp2_table0_complex_per+0x1b0>
mov    0x198(%rsp),%eax
movsd  (%rbx,%rax,8),%xmm2
cvttsd2si %xmm0,%eax
mov    %eax,0x10(%rsp)
movsd  %xmm2,0xa0(%rsp)
movapd %xmm2,%xmm0
subsd  0x48(%rsp),%xmm0
call   1030 <floor@plt>
addsd  0xc99(%rip),%xmm0
cvttsd2si %xmm0,%eax
xorpd  %xmm0,%xmm0
xor    %ecx,%ecx
movsd  0xc7f(%rip),%xmm1
jmp    13ff <interp2_table0_complex_per+0x2ef>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movapd 0xd0(%rsp),%xmm2
movapd %xmm2,%xmm0
movapd 0xc6b(%rip),%xmm1
xorpd  %xmm1,%xmm0
unpcklpd %xmm3,%xmm2
movapd %xmm3,%xmm1
shufpd $0x1,%xmm0,%xmm1
mulpd  %xmm2,%xmm1
movapd 0xe0(%rsp),%xmm2
unpcklpd %xmm2,%xmm2
mulpd  %xmm3,%xmm2
addpd  %xmm1,%xmm2
movapd 0xf0(%rsp),%xmm0
addpd  %xmm2,%xmm0
mov    0x14(%rsp),%ecx
add    $0x1,%ecx
mov    0x18(%rsp),%eax
add    $0x1,%eax
cmp    0x178(%rsp),%ecx
movsd  0xc07(%rip),%xmm1
je     12c4 <interp2_table0_complex_per+0x1b4>
mov    %ecx,0x14(%rsp)
movapd %xmm0,0xf0(%rsp)
mov    %eax,0x18(%rsp)
xorps  %xmm2,%xmm2
cvtsi2sd %eax,%xmm2
movsd  %xmm2,(%rsp)
movsd  0xa0(%rsp),%xmm0
subsd  %xmm2,%xmm0
mulsd  0x70(%rsp),%xmm0
addsd  %xmm1,%xmm0
call   1030 <floor@plt>
movsd  %xmm0,0x20(%rsp)
movsd  (%rsp),%xmm0
divsd  0xc0(%rsp),%xmm0
call   1030 <floor@plt>
cvttsd2si 0x20(%rsp),%eax
cltq
mov    0x78(%rsp),%rcx
movsd  (%rcx,%rax,8),%xmm1
movaps %xmm1,0xe0(%rsp)
mov    0x80(%rsp),%rcx
movsd  (%rcx,%rax,8),%xmm1
movaps %xmm1,0xd0(%rsp)
xorpd  %xmm3,%xmm3
cmpl   $0x0,0x170(%rsp)
mov    0x98(%rsp),%r15
mov    0x90(%rsp),%r14
movsd  0x38(%rsp),%xmm1
movsd  0x40(%rsp),%xmm2
jle    1390 <interp2_table0_complex_per+0x280>
mulsd  0xb0(%rsp),%xmm0
addsd  (%rsp),%xmm0
cvttsd2si %xmm0,%edi
imul   0x8(%rsp),%edi
mov    %edi,0x1c(%rsp)
xorpd  %xmm3,%xmm3
mov    0x170(%rsp),%r13d
mov    0x10(%rsp),%ebx
nopl   0x0(%rax)
movapd %xmm3,0x20(%rsp)
xorps  %xmm3,%xmm3
cvtsi2sd %ebx,%xmm3
movsd  %xmm3,(%rsp)
movapd %xmm2,%xmm0
subsd  %xmm3,%xmm0
mulsd  %xmm1,%xmm0
addsd  0xafa(%rip),%xmm0
call   1030 <floor@plt>
cvttsd2si %xmm0,%eax
movslq %eax,%rbp
movsd  (%rsp),%xmm0
divsd  0x100(%rsp),%xmm0
call   1030 <floor@plt>
movapd 0x110(%rsp),%xmm4
mulsd  %xmm4,%xmm0
addsd  (%rsp),%xmm0
cvttsd2si %xmm0,%eax
add    0x1c(%rsp),%eax
cltq
movsd  (%r12,%rbp,8),%xmm0
movapd %xmm0,%xmm1
xorpd  0xabd(%rip),%xmm1
movsd  (%r15,%rax,8),%xmm2
unpcklpd %xmm2,%xmm0
mov    0xa8(%rsp),%rdi
movsd  (%rdi,%rax,8),%xmm3
unpcklpd %xmm3,%xmm2
unpcklpd %xmm1,%xmm3
movsd  0x38(%rsp),%xmm1
mulpd  %xmm0,%xmm3
movsd  (%r14,%rbp,8),%xmm0
unpcklpd %xmm0,%xmm0
mulpd  %xmm0,%xmm2
addpd  %xmm3,%xmm2
movapd 0x20(%rsp),%xmm3
addpd  %xmm2,%xmm3
movsd  0x40(%rsp),%xmm2
add    $0x1,%ebx
add    $0xffffffff,%r13d
jne    14e0 <interp2_table0_complex_per+0x3d0>
jmp    1390 <interp2_table0_complex_per+0x280>
add    $0x128,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @interp2_table0_complex_per(double* nocapture %0, double* nocapture %1, i32 %2, i32 %3, double* nocapture %4, double* nocapture %5, double* nocapture %6, double* nocapture %7, i32 %8, i32 %9, i32 %10, i32 %11, double* nocapture %12, i32 %13, double* nocapture writeonly %14, double* nocapture writeonly %15) {
%17 = mul nsw i32 %10, %8
%18 = sitofp i32 %17 to double
%19 = fmul double %18, 5.000000e-01
%20 = tail call double @llvm.floor.f64(double %19)
%21 = fptosi double %20 to i32
%22 = sext i32 %21 to i64
%23 = getelementptr inbounds double, double* %4, i64 %22
%24 = getelementptr inbounds double, double* %5, i64 %22
%25 = mul nsw i32 %11, %9
%26 = sitofp i32 %25 to double
%27 = fmul double %26, 5.000000e-01
%28 = tail call double @llvm.floor.f64(double %27)
%29 = fptosi double %28 to i32
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds double, double* %6, i64 %30
%32 = getelementptr inbounds double, double* %7, i64 %30
%33 = icmp sgt i32 %13, 0
br i1 %33, label %34, label %150
34: ; preds = %16
%35 = zext i32 %13 to i64
%36 = sitofp i32 %8 to double
%37 = fmul double %36, 5.000000e-01
%38 = sitofp i32 %9 to double
%39 = fmul double %38, 5.000000e-01
%40 = icmp sgt i32 %9, 0
%41 = sitofp i32 %11 to double
%42 = sitofp i32 %3 to double
%43 = fneg double %42
%44 = icmp sgt i32 %8, 0
%45 = sitofp i32 %10 to double
%46 = sitofp i32 %2 to double
%47 = fneg double %46
br label %48
48: ; preds = %34, %142
%49 = phi double* [ %12, %34 ], [ %55, %142 ]
%50 = phi double* [ %14, %34 ], [ %144, %142 ]
%51 = phi double* [ %15, %34 ], [ %146, %142 ]
%52 = phi i32 [ 0, %34 ], [ %148, %142 ]
%53 = getelementptr inbounds double, double* %49, i64 %35
%54 = load double, double* %53
%55 = getelementptr inbounds double, double* %49, i64 1
%56 = load double, double* %49
%57 = fsub double %56, %37
%58 = tail call double @llvm.floor.f64(double %57)
%59 = fadd double %58, 1.000000e+00
%60 = fptosi double %59 to i32
br i1 %40, label %61, label %142
61: ; preds = %48
%62 = fsub double %54, %39
%63 = tail call double @llvm.floor.f64(double %62)
%64 = fadd double %63, 1.000000e+00
%65 = fptosi double %64 to i32
br label %66
66: ; preds = %61, %126
%67 = phi i32 [ %139, %126 ], [ 0, %61 ]
%68 = phi i32 [ %140, %126 ], [ %65, %61 ]
%69 = phi <2 x double> [ %138, %126 ], [ zeroinitializer, %61 ]
%70 = sitofp i32 %68 to double
%71 = fsub double %54, %70
%72 = fmul double %71, %41
%73 = fadd double %72, 5.000000e-01
%74 = tail call double @llvm.floor.f64(double %73)
%75 = fptosi double %74 to i32
%76 = sext i32 %75 to i64
%77 = getelementptr inbounds double, double* %31, i64 %76
%78 = load double, double* %77
%79 = getelementptr inbounds double, double* %32, i64 %76
%80 = load double, double* %79
%81 = fdiv double %70, %42
%82 = tail call double @llvm.floor.f64(double %81)
%83 = tail call double @llvm.fmuladd.f64(double %43, double %82, double %70)
%84 = fptosi double %83 to i32
%85 = mul nsw i32 %84, %2
br i1 %44, label %86, label %126
86: ; preds = %66, %86
%87 = phi i32 [ %124, %86 ], [ %60, %66 ]
%88 = phi i32 [ %123, %86 ], [ 0, %66 ]
%89 = phi <2 x double> [ %122, %86 ], [ zeroinitializer, %66 ]
%90 = sitofp i32 %87 to double
%91 = fsub double %56, %90
%92 = fmul double %91, %45
%93 = fadd double %92, 5.000000e-01
%94 = tail call double @llvm.floor.f64(double %93)
%95 = fptosi double %94 to i32
%96 = sext i32 %95 to i64
%97 = getelementptr inbounds double, double* %23, i64 %96
%98 = load double, double* %97
%99 = getelementptr inbounds double, double* %24, i64 %96
%100 = load double, double* %99
%101 = fdiv double %90, %46
%102 = tail call double @llvm.floor.f64(double %101)
%103 = tail call double @llvm.fmuladd.f64(double %47, double %102, double %90)
%104 = fptosi double %103 to i32
%105 = add nsw i32 %85, %104
%106 = sext i32 %105 to i64
%107 = getelementptr inbounds double, double* %0, i64 %106
%108 = load double, double* %107
%109 = getelementptr inbounds double, double* %1, i64 %106
%110 = load double, double* %109
%111 = fneg double %100
%112 = insertelement <2 x double> poison, double %100, i64 0
%113 = insertelement <2 x double> %112, double %110, i64 1
%114 = insertelement <2 x double> poison, double %108, i64 0
%115 = insertelement <2 x double> %114, double %111, i64 1
%116 = fmul <2 x double> %113, %115
%117 = insertelement <2 x double> poison, double %98, i64 0
%118 = shufflevector <2 x double> %117, <2 x double> poison, <2 x i32> zeroinitializer
%119 = insertelement <2 x double> poison, double %110, i64 0
%120 = insertelement <2 x double> %119, double %108, i64 1
%121 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %118, <2 x double> %120, <2 x double> %116)
%122 = fadd <2 x double> %89, %121
%123 = add nuw nsw i32 %88, 1
%124 = add nsw i32 %87, 1
%125 = icmp eq i32 %123, %8
br i1 %125, label %126, label %86
126: ; preds = %86, %66
%127 = phi <2 x double> [ zeroinitializer, %66 ], [ %122, %86 ]
%128 = fneg double %80
%129 = insertelement <2 x double> poison, double %80, i64 0
%130 = shufflevector <2 x double> %129, <2 x double> %127, <2 x i32> <i32 0, i32 2>
%131 = extractelement <2 x double> %127, i64 1
%132 = insertelement <2 x double> poison, double %131, i64 0
%133 = insertelement <2 x double> %132, double %128, i64 1
%134 = fmul <2 x double> %130, %133
%135 = insertelement <2 x double> poison, double %78, i64 0
%136 = shufflevector <2 x double> %135, <2 x double> poison, <2 x i32> zeroinitializer
%137 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %136, <2 x double> %127, <2 x double> %134)
%138 = fadd <2 x double> %69, %137
%139 = add nuw nsw i32 %67, 1
%140 = add nsw i32 %68, 1
%141 = icmp eq i32 %139, %9
br i1 %141, label %142, label %66
142: ; preds = %126, %48
%143 = phi <2 x double> [ zeroinitializer, %48 ], [ %138, %126 ]
%144 = getelementptr inbounds double, double* %50, i64 1
%145 = extractelement <2 x double> %143, i64 1
store double %145, double* %50
%146 = getelementptr inbounds double, double* %51, i64 1
%147 = extractelement <2 x double> %143, i64 0
store double %147, double* %51
%148 = add nuw nsw i32 %52, 1
%149 = icmp eq i32 %148, %13
br i1 %149, label %150, label %48
150: ; preds = %142, %16
ret void
}