void interp2_table0_complex_per_adj(double *r_ck, 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,
                                    const double *r_fm, const 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;
  }
  (void)memset((void *)r_ck, 0, K1 * K2 * sizeof(*r_ck));
  (void)memset((void *)i_ck, 0, K1 * K2 * sizeof(*i_ck));
  for (mm = 0; mm < M; mm++) {
    int jj1, jj2;
    const double t2 = p_tm[M];
    const double t1 = *p_tm++;
    const double fmr = *r_fm++;
    const double fmi = *i_fm++;
    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;
      const double v2r = coef2r * fmr + coef2i * fmi;
      const double v2i = coef2r * fmi - coef2i * fmr;
      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;
        r_ck[kk] += coef1r * v2r + coef1i * v2i;
        i_ck[kk] += coef1r * v2i - coef1i * v2r;
      }
    }
  }
}


x86 assembly (clang, -O0):
<interp2_table0_complex_per_adj>:
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  0xe89(%rip),%xmm1
divsd  %xmm1,%xmm0
call   1040 <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  0xe3c(%rip),%xmm1
divsd  %xmm1,%xmm0
call   1040 <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)
mov    -0x8(%rbp),%rdi
mov    -0x14(%rbp),%eax
imul   -0x18(%rbp),%eax
movslq %eax,%rdx
shl    $0x3,%rdx
xor    %esi,%esi
call   1030 <memset@plt>
mov    -0x10(%rbp),%rdi
mov    -0x14(%rbp),%eax
imul   -0x18(%rbp),%eax
movslq %eax,%rdx
shl    $0x3,%rdx
xor    %esi,%esi
call   1030 <memset@plt>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    0x48(%rbp),%eax
jge    163f <interp2_table0_complex_per_adj+0x51f>
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)
mov    0x50(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,0x50(%rbp)
movsd  (%rax),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    0x58(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,0x58(%rbp)
movsd  (%rax),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x50(%rbp),%xmm0
cvtsi2sdl 0x20(%rbp),%xmm1
movsd  0xd51(%rip),%xmm2
divsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
call   1040 <floor@plt>
movaps %xmm0,%xmm1
movsd  0xd41(%rip),%xmm0
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x64(%rbp)
movsd  -0x48(%rbp),%xmm0
cvtsi2sdl 0x28(%rbp),%xmm1
movsd  0xd1c(%rip),%xmm2
divsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
call   1040 <floor@plt>
movaps %xmm0,%xmm1
movsd  0xd0c(%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    162c <interp2_table0_complex_per_adj+0x50c>
movsd  -0x48(%rbp),%xmm0
cvtsi2sdl -0x68(%rbp),%xmm1
subsd  %xmm1,%xmm0
cvtsi2sdl 0x38(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  0xcd2(%rip),%xmm0
addsd  -0x70(%rbp),%xmm0
call   1040 <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   1040 <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)
movsd  -0x80(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm2
movsd  -0x88(%rbp),%xmm1
mulsd  -0x60(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x80(%rbp),%xmm0
movsd  -0x60(%rbp),%xmm2
movsd  -0x88(%rbp),%xmm1
ff
mulsd  -0x58(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%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    1610 <interp2_table0_complex_per_adj+0x4f0>
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  0xb76(%rip),%xmm0
addsd  -0xb0(%rbp),%xmm0
ff
call   1040 <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   1040 <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
movsd  -0x98(%rbp),%xmm2
ff
movsd  -0xc8(%rbp),%xmm1
ff
mulsd  -0xa0(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movslq -0xd0(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0xc0(%rbp),%xmm0
ff
movsd  -0xa0(%rbp),%xmm2
ff
movsd  -0xc8(%rbp),%xmm1
ff
mulsd  -0x98(%rbp),%xmm1
ff
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0xd0(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
mov    -0xa4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xa4(%rbp)
jmp    1464 <interp2_table0_complex_per_adj+0x344>
jmp    1615 <interp2_table0_complex_per_adj+0x4f5>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x68(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x68(%rbp)
jmp    130e <interp2_table0_complex_per_adj+0x1ee>
jmp    1631 <interp2_table0_complex_per_adj+0x511>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1237 <interp2_table0_complex_per_adj+0x117>
add    $0xf0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<interp2_table0_complex_per_adj>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x118,%rsp
mov    %r9,0x8(%rsp)
mov    %r8,0x10(%rsp)
mov    %ecx,%ebx
mov    %edx,%r15d
mov    %rsi,%r14
mov    %rdi,%rbp
mov    0x178(%rsp),%r12d
mov    0x170(%rsp),%eax
imul   0x160(%rsp),%eax
cvtsi2sd %eax,%xmm0
mulsd  0xe97(%rip),%xmm0
call   1040 <floor@plt>
movsd  %xmm0,0x28(%rsp)
mov    %r12d,%eax
imul   0x168(%rsp),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mulsd  0xe72(%rip),%xmm0
call   1040 <floor@plt>
movsd  %xmm0,0x20(%rsp)
mov    %ebx,%r13d
mov    %r15d,0x1c(%rsp)
imul   %r15d,%ebx
movslq %ebx,%rbx
shl    $0x3,%rbx
xor    %r12d,%r12d
mov    %rbp,0xa0(%rsp)
mov    %rbp,%rdi
xor    %esi,%esi
mov    %rbx,%rdx
call   1030 <memset@plt>
mov    %r14,0xa8(%rsp)
mov    %r14,%rdi
xor    %esi,%esi
mov    %rbx,%rdx
call   1030 <memset@plt>
movsd  0xe1f(%rip),%xmm0
cmpl   $0x0,0x188(%rsp)
jle    1606 <interp2_table0_complex_per_adj+0x4e6>
mov    0x198(%rsp),%rax
mov    %rax,0x48(%rsp)
mov    0x190(%rsp),%rax
mov    %rax,0x40(%rsp)
mov    0x180(%rsp),%rdx
cvttsd2si 0x28(%rsp),%eax
cltq
mov    0x10(%rsp),%rcx
lea    (%rcx,%rax,8),%rcx
mov    %rcx,0xb8(%rsp)
cvttsd2si 0x20(%rsp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
mov    0x150(%rsp),%rsi
add    %rcx,%rsi
mov    %rsi,0x70(%rsp)
add    0x158(%rsp),%rcx
mov    %rcx,0x78(%rsp)
mov    0x8(%rsp),%rcx
lea    (%rcx,%rax,8),%rax
mov    %rax,0x98(%rsp)
cvtsi2sdl 0x160(%rsp),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x68(%rsp)
xorps  %xmm1,%xmm1
cvtsi2sdl 0x168(%rsp),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x60(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sdl 0x178(%rsp),%xmm0
movsd  %xmm0,0x90(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sd %r13d,%xmm0
movapd 0xd5d(%rip),%xmm1
movapd %xmm0,0xe0(%rsp)
xorpd  %xmm1,%xmm0
movapd %xmm0,0xd0(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sdl 0x170(%rsp),%xmm0
movsd  %xmm0,0x58(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sdl 0x1c(%rsp),%xmm0
movapd %xmm0,0xf0(%rsp)
xorpd  %xmm0,%xmm1
movapd %xmm1,0x100(%rsp)
mov    0x98(%rsp),%rbx
mov    0xa8(%rsp),%rbp
mov    0xa0(%rsp),%r15
jmp    134c <interp2_table0_complex_per_adj+0x22c>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    0x50(%rsp),%rdx
add    $0x8,%rdx
addq   $0x8,0x40(%rsp)
addq   $0x8,0x48(%rsp)
mov    0x30(%rsp),%r12d
add    $0x1,%r12d
cmp    0x188(%rsp),%r12d
je     1606 <interp2_table0_complex_per_adj+0x4e6>
mov    %r12d,0x30(%rsp)
mov    %rdx,0x50(%rsp)
movsd  (%rdx),%xmm0
movsd  %xmm0,0x20(%rsp)
subsd  0x68(%rsp),%xmm0
call   1040 <floor@plt>
addsd  0xc95(%rip),%xmm0
cmpl   $0x0,0x168(%rsp)
jle    1320 <interp2_table0_complex_per_adj+0x200>
mov    0x50(%rsp),%rax
mov    0x188(%rsp),%ecx
movsd  (%rax,%rcx,8),%xmm1
mov    0x40(%rsp),%rax
movsd  (%rax),%xmm2
movsd  %xmm2,0x88(%rsp)
mov    0x48(%rsp),%rax
movsd  (%rax),%xmm2
movsd  %xmm2,0x80(%rsp)
cvttsd2si %xmm0,%eax
mov    %eax,0x34(%rsp)
movsd  %xmm1,0xb0(%rsp)
movapd %xmm1,%xmm0
subsd  0x60(%rsp),%xmm0
call   1040 <floor@plt>
movsd  0xc26(%rip),%xmm1
addsd  0xc26(%rip),%xmm0
cvttsd2si %xmm0,%eax
xor    %ecx,%ecx
jmp    1413 <interp2_table0_complex_per_adj+0x2f3>
nopw   0x0(%rax,%rax,1)
mov    0x38(%rsp),%ecx
add    $0x1,%ecx
mov    0x3c(%rsp),%eax
add    $0x1,%eax
cmp    0x168(%rsp),%ecx
movsd  0xbf3(%rip),%xmm1
je     1320 <interp2_table0_complex_per_adj+0x200>
mov    %ecx,0x38(%rsp)
mov    %eax,0x3c(%rsp)
xorps  %xmm2,%xmm2
cvtsi2sd %eax,%xmm2
movsd  %xmm2,0x8(%rsp)
movsd  0xb0(%rsp),%xmm0
subsd  %xmm2,%xmm0
mulsd  0x90(%rsp),%xmm0
addsd  %xmm1,%xmm0
call   1040 <floor@plt>
movsd  %xmm0,0x10(%rsp)
movsd  0x8(%rsp),%xmm0
divsd  0xe0(%rsp),%xmm0
call   1040 <floor@plt>
cmpl   $0x0,0x160(%rsp)
movsd  0x58(%rsp),%xmm5
movsd  0x20(%rsp),%xmm6
jle    13f0 <interp2_table0_complex_per_adj+0x2d0>
mulsd  0xd0(%rsp),%xmm0
addsd  0x8(%rsp),%xmm0
cvttsd2si %xmm0,%r13d
imul   0x1c(%rsp),%r13d
cvttsd2si 0x10(%rsp),%eax
cltq
mov    0x70(%rsp),%rcx
movsd  (%rcx,%rax,8),%xmm7
mov    0x78(%rsp),%rcx
movsd  (%rcx,%rax,8),%xmm0
movsd  0x80(%rsp),%xmm8
movapd %xmm8,%xmm1
mulsd  %xmm0,%xmm1
movapd %xmm7,%xmm3
movsd  0x88(%rsp),%xmm2
mulsd  %xmm2,%xmm3
addsd  %xmm1,%xmm3
movsd  %xmm3,0xc0(%rsp)
mulsd  %xmm2,%xmm0
mulsd  %xmm8,%xmm7
subsd  %xmm0,%xmm7
movsd  %xmm7,0xc8(%rsp)
mov    0x160(%rsp),%r14d
mov    0x34(%rsp),%r12d
cs nopw 0x0(%rax,%rax,1)
nop
xorps  %xmm1,%xmm1
cvtsi2sd %r12d,%xmm1
movsd  %xmm1,0x8(%rsp)
movapd %xmm6,%xmm0
subsd  %xmm1,%xmm0
mulsd  %xmm5,%xmm0
addsd  0xace(%rip),%xmm0
call   1040 <floor@plt>
cvttsd2si %xmm0,%eax
cltq
mov    0xb8(%rsp),%rcx
movsd  (%rcx,%rax,8),%xmm0
movsd  %xmm0,0x10(%rsp)
movsd  (%rbx,%rax,8),%xmm0
movsd  %xmm0,0x28(%rsp)
movsd  0x8(%rsp),%xmm0
divsd  0xf0(%rsp),%xmm0
call   1040 <floor@plt>
movsd  0x20(%rsp),%xmm6
movsd  0x58(%rsp),%xmm5
movapd 0x100(%rsp),%xmm9
mulsd  %xmm9,%xmm0
addsd  0x8(%rsp),%xmm0
cvttsd2si %xmm0,%eax
add    %r13d,%eax
movsd  0xc8(%rsp),%xmm2
movapd %xmm2,%xmm0
movsd  0x28(%rsp),%xmm7
mulsd  %xmm7,%xmm0
movsd  0x10(%rsp),%xmm4
movapd %xmm4,%xmm1
movsd  0xc0(%rsp),%xmm3
mulsd  %xmm3,%xmm1
addsd  %xmm0,%xmm1
cltq
addsd  (%r15,%rax,8),%xmm1
movsd  %xmm1,(%r15,%rax,8)
mulsd  %xmm3,%xmm7
movapd %xmm4,%xmm0
mulsd  %xmm2,%xmm0
subsd  %xmm7,%xmm0
addsd  0x0(%rbp,%rax,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rax,8)
add    $0x1,%r12d
add    $0xffffffff,%r14d
jne    1510 <interp2_table0_complex_per_adj+0x3f0>
jmp    13f0 <interp2_table0_complex_per_adj+0x2d0>
add    $0x118,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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