void f3_coef(double *f, int l, int m, int n, double a, double b, double c) {
  double Ca[32], Cb[32], Cc[32];
  int iA, iB, iC;
  int I;
  Ca[l] = 1.0;
  for (iA = l - 1; iA >= 0; iA--)
    Ca[iA] = a * (iA + 1) / (l - iA) * Ca[iA + 1];
  Cb[m] = 1.0;
  for (iB = m - 1; iB >= 0; iB--)
    Cb[iB] = b * (iB + 1) / (m - iB) * Cb[iB + 1];
  Cc[n] = 1.0;
  for (iC = n - 1; iC >= 0; iC--)
    Cc[iC] = c * (iC + 1) / (n - iC) * Cc[iC + 1];
  for (I = 0; I <= (l + m + n); I++)
    f[I] = 0.0;
  for (iA = 0; iA < l + 1; iA++)
    for (iB = 0; iB < m + 1; iB++)
      for (iC = 0; iC < n + 1; iC++) {
        I = iA + iB + iC;
        f[I] += Ca[iA] * Cb[iB] * Cc[iC];
      }
  return;
}


x86 assembly (clang, -O0):
<f3_coef>:
push   %rbp
mov    %rsp,%rbp
sub    $0x2c0,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
movsd  %xmm0,-0x20(%rbp)
movsd  %xmm1,-0x28(%rbp)
movsd  %xmm2,-0x30(%rbp)
movslq -0xc(%rbp),%rax
movsd  0xecd(%rip),%xmm0
movsd  %xmm0,-0x130(%rbp,%rax,8)
ff ff
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x334(%rbp)
cmpl   $0x0,-0x334(%rbp)
jl     11b4 <f3_coef+0xb4>
movsd  -0x20(%rbp),%xmm0
mov    -0x334(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm1
mulsd  %xmm1,%xmm0
mov    -0xc(%rbp),%eax
sub    -0x334(%rbp),%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
mov    -0x334(%rbp),%eax
add    $0x1,%eax
cltq
mulsd  -0x130(%rbp,%rax,8),%xmm0
ff ff
movslq -0x334(%rbp),%rax
movsd  %xmm0,-0x130(%rbp,%rax,8)
ff ff
mov    -0x334(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x334(%rbp)
jmp    1148 <f3_coef+0x48>
movslq -0x10(%rbp),%rax
movsd  0xe40(%rip),%xmm0
movsd  %xmm0,-0x230(%rbp,%rax,8)
ff ff
mov    -0x10(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x338(%rbp)
cmpl   $0x0,-0x338(%rbp)
jl     1241 <f3_coef+0x141>
movsd  -0x28(%rbp),%xmm0
mov    -0x338(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm1
mulsd  %xmm1,%xmm0
mov    -0x10(%rbp),%eax
sub    -0x338(%rbp),%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
mov    -0x338(%rbp),%eax
add    $0x1,%eax
cltq
mulsd  -0x230(%rbp,%rax,8),%xmm0
ff ff
movslq -0x338(%rbp),%rax
movsd  %xmm0,-0x230(%rbp,%rax,8)
ff ff
mov    -0x338(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x338(%rbp)
jmp    11d5 <f3_coef+0xd5>
movslq -0x14(%rbp),%rax
movsd  0xdb3(%rip),%xmm0
movsd  %xmm0,-0x330(%rbp,%rax,8)
ff ff
mov    -0x14(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x33c(%rbp)
cmpl   $0x0,-0x33c(%rbp)
jl     12ce <f3_coef+0x1ce>
movsd  -0x30(%rbp),%xmm0
mov    -0x33c(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm1
mulsd  %xmm1,%xmm0
mov    -0x14(%rbp),%eax
sub    -0x33c(%rbp),%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
mov    -0x33c(%rbp),%eax
add    $0x1,%eax
cltq
mulsd  -0x330(%rbp,%rax,8),%xmm0
ff ff
movslq -0x33c(%rbp),%rax
movsd  %xmm0,-0x330(%rbp,%rax,8)
ff ff
mov    -0x33c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x33c(%rbp)
jmp    1262 <f3_coef+0x162>
movl   $0x0,-0x340(%rbp)
mov    -0x340(%rbp),%eax
mov    -0xc(%rbp),%ecx
add    -0x10(%rbp),%ecx
add    -0x14(%rbp),%ecx
cmp    %ecx,%eax
jg     1316 <f3_coef+0x216>
mov    -0x8(%rbp),%rax
movslq -0x340(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x340(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x340(%rbp)
jmp    12d8 <f3_coef+0x1d8>
movl   $0x0,-0x334(%rbp)
mov    -0x334(%rbp),%eax
mov    -0xc(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    141b <f3_coef+0x31b>
movl   $0x0,-0x338(%rbp)
mov    -0x338(%rbp),%eax
mov    -0x10(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    1402 <f3_coef+0x302>
movl   $0x0,-0x33c(%rbp)
mov    -0x33c(%rbp),%eax
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    13e9 <f3_coef+0x2e9>
mov    -0x334(%rbp),%eax
add    -0x338(%rbp),%eax
add    -0x33c(%rbp),%eax
mov    %eax,-0x340(%rbp)
movslq -0x334(%rbp),%rax
movsd  -0x130(%rbp,%rax,8),%xmm0
ff ff
movslq -0x338(%rbp),%rax
mulsd  -0x230(%rbp,%rax,8),%xmm0
ff ff
movslq -0x33c(%rbp),%rax
movsd  -0x330(%rbp,%rax,8),%xmm2
ff ff
mov    -0x8(%rbp),%rax
movslq -0x340(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x33c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x33c(%rbp)
jmp    135c <f3_coef+0x25c>
jmp    13ee <f3_coef+0x2ee>
mov    -0x338(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x338(%rbp)
jmp    133e <f3_coef+0x23e>
jmp    1407 <f3_coef+0x307>
mov    -0x334(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x334(%rbp)
jmp    1320 <f3_coef+0x220>
add    $0x2c0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<f3_coef>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x338,%rsp
mov    %ecx,%ebx
mov    %edx,%r14d
mov    %esi,%ebp
mov    %rdi,0x10(%rsp)
movslq %esi,%rax
movabs $0x3ff0000000000000,%rsi
mov    %rsi,0x230(%rsp,%rax,8)
test   %eax,%eax
mov    %rbx,0x8(%rsp)
jle    124a <f3_coef+0x13a>
mov    %ebp,%edi
cmp    $0x1,%ebp
jne    115e <f3_coef+0x4e>
mov    %rdi,%rax
jmp    1207 <f3_coef+0xf7>
lea    -0x1(%rdi),%rcx
mov    %rbp,%r8
add    $0xffffffff,%ebp
mov    $0x8,%edx
mov    %rcx,%rax
mul    %rdx
seto   %dl
cmp    %ecx,%ebp
jb     18c4 <f3_coef+0x7b4>
shr    $0x20,%rcx
jne    18c4 <f3_coef+0x7b4>
lea    (%rsp,%rbp,8),%rbx
add    $0x230,%rbx
mov    %rbx,%rcx
sub    %rax,%rcx
cmp    %rbx,%rcx
ja     11fc <f3_coef+0xec>
test   %dl,%dl
jne    11fc <f3_coef+0xec>
lea    (%rsp,%rdi,8),%rbx
add    $0x230,%rbx
mov    %rbx,%rcx
sub    %rax,%rcx
cmp    %rbx,%rcx
ja     11fc <f3_coef+0xec>
test   %dl,%dl
jne    11fc <f3_coef+0xec>
lea    0x1(%rbp),%rax
sub    %rdi,%rax
lea    (%rsp,%rax,8),%rax
add    $0x230,%rax
lea    (%rsp,%rdi,8),%rcx
add    $0x238,%rcx
cmp    %rcx,%rax
jae    1683 <f3_coef+0x573>
lea    (%rsp,%rbp,8),%rax
add    $0x238,%rax
lea    0x238(%rsp),%rcx
cmp    %rax,%rcx
jae    1683 <f3_coef+0x573>
mov    %rdi,%rax
mov    %r8,%rbp
mov    0x8(%rsp),%rbx
mov    %ebp,%ecx
sub    %eax,%ecx
add    $0x1,%ecx
xchg   %ax,%ax
xorps  %xmm3,%xmm3
cvtsi2sd %eax,%xmm3
mulsd  %xmm0,%xmm3
xorps  %xmm4,%xmm4
cvtsi2sd %ecx,%xmm4
divsd  %xmm4,%xmm3
mulsd  0x230(%rsp,%rax,8),%xmm3
lea    -0x1(%rax),%rdx
mov    %edx,%edi
movsd  %xmm3,0x230(%rsp,%rdi,8)
add    $0x1,%ecx
cmp    $0x1,%rax
mov    %rdx,%rax
ja     1210 <f3_coef+0x100>
movslq %r14d,%rax
mov    %rsi,0x130(%rsp,%rax,8)
test   %eax,%eax
jle    135a <f3_coef+0x24a>
mov    %r14d,%edi
cmp    $0x1,%r14d
jne    126e <f3_coef+0x15e>
mov    %rdi,%rax
jmp    1318 <f3_coef+0x208>
mov    %rbp,%r8
lea    -0x1(%rdi),%rcx
lea    -0x1(%r14),%ebp
mov    $0x8,%edx
mov    %rcx,%rax
mul    %rdx
seto   %dl
cmp    %ecx,%ebp
jb     18cf <f3_coef+0x7bf>
shr    $0x20,%rcx
jne    18cf <f3_coef+0x7bf>
lea    (%rsp,%rbp,8),%rbx
add    $0x130,%rbx
mov    %rbx,%rcx
sub    %rax,%rcx
cmp    %rbx,%rcx
ja     130d <f3_coef+0x1fd>
test   %dl,%dl
jne    130d <f3_coef+0x1fd>
lea    (%rsp,%rdi,8),%rbx
add    $0x130,%rbx
mov    %rbx,%rcx
sub    %rax,%rcx
cmp    %rbx,%rcx
ja     130d <f3_coef+0x1fd>
test   %dl,%dl
jne    130d <f3_coef+0x1fd>
lea    0x1(%rbp),%rax
sub    %rdi,%rax
lea    (%rsp,%rax,8),%rax
add    $0x130,%rax
lea    (%rsp,%rdi,8),%rcx
add    $0x138,%rcx
cmp    %rcx,%rax
jae    174c <f3_coef+0x63c>
lea    (%rsp,%rbp,8),%rax
add    $0x138,%rax
lea    0x138(%rsp),%rcx
cmp    %rax,%rcx
jae    174c <f3_coef+0x63c>
mov    %rdi,%rax
mov    %r8,%rbp
mov    0x8(%rsp),%rbx
mov    %r14d,%ecx
sub    %eax,%ecx
add    $0x1,%ecx
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mulsd  %xmm1,%xmm0
xorps  %xmm3,%xmm3
cvtsi2sd %ecx,%xmm3
divsd  %xmm3,%xmm0
mulsd  0x130(%rsp,%rax,8),%xmm0
lea    -0x1(%rax),%rdx
mov    %edx,%edi
movsd  %xmm0,0x130(%rsp,%rdi,8)
add    $0x1,%ecx
cmp    $0x1,%rax
mov    %rdx,%rax
ja     1320 <f3_coef+0x210>
movslq %ebx,%rax
mov    %rsi,0x30(%rsp,%rax,8)
test   %eax,%eax
jle    1454 <f3_coef+0x344>
mov    %ebx,%esi
cmp    $0x1,%ebx
jne    1379 <f3_coef+0x269>
mov    %rsi,%rax
jmp    140a <f3_coef+0x2fa>
lea    -0x1(%rsi),%rcx
lea    -0x1(%rbx),%edi
mov    $0x8,%edx
mov    %rcx,%rax
mul    %rdx
seto   %dl
cmp    %ecx,%edi
jb     1371 <f3_coef+0x261>
shr    $0x20,%rcx
jne    1371 <f3_coef+0x261>
mov    %rbp,%rbx
lea    (%rsp,%rdi,8),%rbp
add    $0x30,%rbp
mov    %rbp,%rcx
sub    %rax,%rcx
cmp    %rbp,%rcx
ja     13ff <f3_coef+0x2ef>
test   %dl,%dl
jne    13ff <f3_coef+0x2ef>
lea    (%rsp,%rsi,8),%rbp
add    $0x30,%rbp
mov    %rbp,%rcx
sub    %rax,%rcx
cmp    %rbp,%rcx
ja     13ff <f3_coef+0x2ef>
test   %dl,%dl
jne    13ff <f3_coef+0x2ef>
lea    0x1(%rdi),%rax
sub    %rsi,%rax
lea    (%rsp,%rax,8),%rax
add    $0x30,%rax
lea    (%rsp,%rsi,8),%rcx
add    $0x38,%rcx
cmp    %rcx,%rax
jae    180b <f3_coef+0x6fb>
lea    (%rsp,%rdi,8),%rax
add    $0x38,%rax
lea    0x38(%rsp),%rcx
cmp    %rax,%rcx
jae    180b <f3_coef+0x6fb>
mov    %rsi,%rax
mov    %rbx,%rbp
mov    0x8(%rsp),%rbx
mov    %ebx,%ecx
sub    %eax,%ecx
add    $0x1,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mulsd  %xmm2,%xmm0
xorps  %xmm1,%xmm1
cvtsi2sd %ecx,%xmm1
divsd  %xmm1,%xmm0
mulsd  0x30(%rsp,%rax,8),%xmm0
lea    -0x1(%rax),%rdx
mov    %edx,%esi
movsd  %xmm0,0x30(%rsp,%rsi,8)
add    $0x1,%ecx
cmp    $0x1,%rax
mov    %rdx,%rax
ja     1420 <f3_coef+0x310>
lea    (%r14,%rbp,1),%eax
add    %ebx,%eax
js     1476 <f3_coef+0x366>
lea    (%rbx,%r14,1),%eax
lea    (%rax,%rbp,1),%edx
add    $0x1,%edx
shl    $0x3,%rdx
mov    0x10(%rsp),%rdi
xor    %esi,%esi
call   1030 <memset@plt>
test   %ebp,%ebp
js     1671 <f3_coef+0x561>
lea    0x1(%rbx),%eax
lea    0x1(%r14),%ecx
add    $0x1,%ebp
mov    %eax,%r12d
and    $0xfffffffc,%r12d
lea    -0x4(%r12),%r9
mov    %r9,%r10
shr    $0x2,%r10
add    $0x1,%r10
mov    %r10,%rdx
and    $0xfffffffffffffffe,%rdx
mov    %rdx,0x28(%rsp)
mov    0x10(%rsp),%r13
lea    0x30(%r13),%r11
xor    %esi,%esi
mov    0x8(%rsp),%rbx
mov    %r14,0x20(%rsp)
mov    %rbp,0x18(%rsp)
jmp    14ef <f3_coef+0x3df>
nopw   0x0(%rax,%rax,1)
add    $0x1,%rsi
add    $0x8,%r11
add    $0x8,%r13
mov    0x18(%rsp),%rbp
cmp    %rbp,%rsi
mov    0x20(%rsp),%r14
je     1671 <f3_coef+0x561>
test   %r14d,%r14d
js     14d0 <f3_coef+0x3c0>
mov    %r13,%rdi
mov    %r11,%rbp
xor    %r14d,%r14d
jmp    1511 <f3_coef+0x401>
nop
add    $0x1,%r14
add    $0x8,%rbp
add    $0x8,%rdi
cmp    %rcx,%r14
je     14d0 <f3_coef+0x3c0>
test   %ebx,%ebx
js     1500 <f3_coef+0x3f0>
movsd  0x230(%rsp,%rsi,8),%xmm0
mulsd  0x130(%rsp,%r14,8),%xmm0
cmp    $0x4,%eax
jae    1540 <f3_coef+0x430>
xor    %edx,%edx
jmp    1640 <f3_coef+0x530>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
test   %r9,%r9
je     1662 <f3_coef+0x552>
mov    0x28(%rsp),%r8
xor    %r15d,%r15d
nopl   0x0(%rax)
movapd 0x30(%rsp,%r15,8),%xmm2
mulpd  %xmm1,%xmm2
movupd -0x30(%rbp,%r15,8),%xmm3
addpd  %xmm2,%xmm3
movupd -0x20(%rbp,%r15,8),%xmm2
movupd -0x10(%rbp,%r15,8),%xmm4
movupd 0x0(%rbp,%r15,8),%xmm5
movapd 0x40(%rsp,%r15,8),%xmm6
mulpd  %xmm1,%xmm6
addpd  %xmm2,%xmm6
movupd %xmm3,-0x30(%rbp,%r15,8)
movupd %xmm6,-0x20(%rbp,%r15,8)
movapd 0x50(%rsp,%r15,8),%xmm2
mulpd  %xmm1,%xmm2
addpd  %xmm4,%xmm2
movapd 0x60(%rsp,%r15,8),%xmm3
mulpd  %xmm1,%xmm3
addpd  %xmm5,%xmm3
movupd %xmm2,-0x10(%rbp,%r15,8)
movupd %xmm3,0x0(%rbp,%r15,8)
add    $0x8,%r15
add    $0xfffffffffffffffe,%r8
jne    1560 <f3_coef+0x450>
test   $0x1,%r10b
je     1625 <f3_coef+0x515>
lea    (%r14,%rsi,1),%rdx
add    %r15,%rdx
mov    0x10(%rsp),%rbx
movupd (%rbx,%rdx,8),%xmm2
movupd 0x10(%rbx,%rdx,8),%xmm3
movapd 0x30(%rsp,%r15,8),%xmm4
mulpd  %xmm1,%xmm4
mulpd  0x40(%rsp,%r15,8),%xmm1
addpd  %xmm2,%xmm4
addpd  %xmm3,%xmm1
movupd %xmm4,(%rbx,%rdx,8)
movupd %xmm1,0x10(%rbx,%rdx,8)
mov    0x8(%rsp),%rbx
mov    %r12,%rdx
cmp    %rax,%r12
je     1500 <f3_coef+0x3f0>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  0x30(%rsp,%rdx,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rdi,%rdx,8),%xmm1
movsd  %xmm1,(%rdi,%rdx,8)
add    $0x1,%rdx
cmp    %rdx,%rax
jne    1640 <f3_coef+0x530>
jmp    1500 <f3_coef+0x3f0>
xor    %r15d,%r15d
test   $0x1,%r10b
jne    15e4 <f3_coef+0x4d4>
jmp    1625 <f3_coef+0x515>
add    $0x338,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    %edi,%ebp
and    $0xfffffffe,%ebp
mov    %edi,%eax
and    $0x1,%eax
movd   %r8d,%xmm3
pshufd $0x50,%xmm3,%xmm10
movdqa 0x960(%rip),%xmm4
paddd  %xmm10,%xmm4
movapd %xmm0,%xmm5
unpcklpd %xmm0,%xmm5
movdqa 0x95a(%rip),%xmm8
movdqa 0x961(%rip),%xmm9
mov    %rdi,%rdx
mov    %rbp,%rcx
cs nopw 0x0(%rax,%rax,1)
nop
movq   %rdx,%xmm6
pshufd $0x44,%xmm6,%xmm6
paddq  %xmm8,%xmm6
pshufd $0xfa,%xmm6,%xmm7
movq   %xmm6,%rbx
movd   %ebx,%xmm6
movss  %xmm6,%xmm7
cvtdq2pd %xmm4,%xmm6
mulpd  %xmm5,%xmm6
movdqa %xmm10,%xmm3
psubd  %xmm7,%xmm3
cvtdq2pd %xmm3,%xmm3
divpd  %xmm3,%xmm6
movupd 0x228(%rsp,%rdx,8),%xmm3
mov    %ebx,%ebx
shufpd $0x1,%xmm6,%xmm6
mulpd  %xmm3,%xmm6
movupd %xmm6,0x228(%rsp,%rbx,8)
paddd  %xmm9,%xmm4
add    $0xfffffffffffffffe,%rdx
add    $0xfffffffffffffffe,%rcx
jne    16d0 <f3_coef+0x5c0>
cmp    %rdi,%rbp
mov    %r8,%rbp
mov    0x8(%rsp),%rbx
jne    1207 <f3_coef+0xf7>
jmp    124a <f3_coef+0x13a>
mov    %edi,%ebp
and    $0xfffffffe,%ebp
mov    %edi,%eax
and    $0x1,%eax
movd   %r14d,%xmm0
pshufd $0x50,%xmm0,%xmm6
movdqa 0x898(%rip),%xmm3
paddd  %xmm6,%xmm3
movapd %xmm1,%xmm4
unpcklpd %xmm1,%xmm4
movdqa 0x893(%rip),%xmm8
movdqa 0x89a(%rip),%xmm9
mov    %rdi,%rdx
mov    %rbp,%rcx
nopl   0x0(%rax)
movq   %rdx,%xmm7
pshufd $0x44,%xmm7,%xmm7
paddq  %xmm8,%xmm7
pshufd $0xfa,%xmm7,%xmm5
movq   %xmm7,%rbx
movd   %ebx,%xmm7
movss  %xmm7,%xmm5
cvtdq2pd %xmm3,%xmm7
mulpd  %xmm4,%xmm7
movdqa %xmm6,%xmm0
psubd  %xmm5,%xmm0
cvtdq2pd %xmm0,%xmm0
divpd  %xmm0,%xmm7
movupd 0x128(%rsp,%rdx,8),%xmm0
mov    %ebx,%ebx
shufpd $0x1,%xmm7,%xmm7
mulpd  %xmm0,%xmm7
movupd %xmm7,0x128(%rsp,%rbx,8)
paddd  %xmm9,%xmm3
add    $0xfffffffffffffffe,%rdx
add    $0xfffffffffffffffe,%rcx
jne    1790 <f3_coef+0x680>
cmp    %rdi,%rbp
mov    %r8,%rbp
mov    0x8(%rsp),%rbx
jne    1318 <f3_coef+0x208>
jmp    135a <f3_coef+0x24a>
mov    %esi,%ecx
and    $0xfffffffe,%ecx
mov    %esi,%eax
and    $0x1,%eax
movd   0x8(%rsp),%xmm0
pshufd $0x50,%xmm0,%xmm4
movdqa 0x7d8(%rip),%xmm1
paddd  %xmm4,%xmm1
movapd %xmm2,%xmm3
unpcklpd %xmm2,%xmm3
movdqa 0x7d3(%rip),%xmm8
movdqa 0x7db(%rip),%xmm5
mov    %rsi,%rdx
mov    %rcx,%rdi
nopl   0x0(%rax,%rax,1)
movq   %rdx,%xmm6
pshufd $0x44,%xmm6,%xmm6
paddq  %xmm8,%xmm6
pshufd $0xfa,%xmm6,%xmm7
movq   %xmm6,%rbp
movd   %ebp,%xmm6
movss  %xmm6,%xmm7
cvtdq2pd %xmm1,%xmm6
mulpd  %xmm3,%xmm6
movdqa %xmm4,%xmm0
psubd  %xmm7,%xmm0
cvtdq2pd %xmm0,%xmm0
divpd  %xmm0,%xmm6
movupd 0x28(%rsp,%rdx,8),%xmm0
mov    %ebp,%ebp
shufpd $0x1,%xmm6,%xmm6
mulpd  %xmm0,%xmm6
movupd %xmm6,0x28(%rsp,%rbp,8)
paddd  %xmm5,%xmm1
add    $0xfffffffffffffffe,%rdx
add    $0xfffffffffffffffe,%rdi
jne    1850 <f3_coef+0x740>
cmp    %rsi,%rcx
mov    %rbx,%rbp
mov    0x8(%rsp),%rbx
jne    140a <f3_coef+0x2fa>
jmp    1454 <f3_coef+0x344>
mov    %rdi,%rax
mov    %r8,%rbp
jmp    1207 <f3_coef+0xf7>
mov    %rdi,%rax
mov    %r8,%rbp
jmp    1318 <f3_coef+0x208>

LLVM IR (-O0):
define void @f3_coef(double* %0, i32 %1, i32 %2, i32 %3, double %4, double %5, double %6) {
%8 = alloca double*
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca double
%13 = alloca double
%14 = alloca double
%15 = alloca [32 x double]
%16 = alloca [32 x double]
%17 = alloca [32 x double]
%18 = alloca i32
%19 = alloca i32
%20 = alloca i32
%21 = alloca i32
store double* %0, double** %8
store i32 %1, i32* %9
store i32 %2, i32* %10
store i32 %3, i32* %11
store double %4, double* %12
store double %5, double* %13
store double %6, double* %14
%22 = load i32, i32* %9
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds [32 x double], [32 x double]* %15, i64 0, i64 %23
store double 1.000000e+00, double* %24
%25 = load i32, i32* %9
%26 = sub nsw i32 %25, 1
store i32 %26, i32* %18
br label %27
27: ; preds = %50, %7
%28 = load i32, i32* %18
%29 = icmp sge i32 %28, 0
br i1 %29, label %30, label %53
30: ; preds = %27
%31 = load double, double* %12
%32 = load i32, i32* %18
%33 = add nsw i32 %32, 1
%34 = sitofp i32 %33 to double
%35 = fmul double %31, %34
%36 = load i32, i32* %9
%37 = load i32, i32* %18
%38 = sub nsw i32 %36, %37
%39 = sitofp i32 %38 to double
%40 = fdiv double %35, %39
%41 = load i32, i32* %18
%42 = add nsw i32 %41, 1
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds [32 x double], [32 x double]* %15, i64 0, i64 %43
%45 = load double, double* %44
%46 = fmul double %40, %45
%47 = load i32, i32* %18
%48 = sext i32 %47 to i64
%49 = getelementptr inbounds [32 x double], [32 x double]* %15, i64 0, i64 %48
store double %46, double* %49
br label %50
50: ; preds = %30
%51 = load i32, i32* %18
%52 = add nsw i32 %51, -1
store i32 %52, i32* %18
br label %27
53: ; preds = %27
%54 = load i32, i32* %10
%55 = sext i32 %54 to i64
%56 = getelementptr inbounds [32 x double], [32 x double]* %16, i64 0, i64 %55
store double 1.000000e+00, double* %56
%57 = load i32, i32* %10
%58 = sub nsw i32 %57, 1
store i32 %58, i32* %19
br label %59
59: ; preds = %82, %53
%60 = load i32, i32* %19
%61 = icmp sge i32 %60, 0
br i1 %61, label %62, label %85
62: ; preds = %59
%63 = load double, double* %13
%64 = load i32, i32* %19
%65 = add nsw i32 %64, 1
%66 = sitofp i32 %65 to double
%67 = fmul double %63, %66
%68 = load i32, i32* %10
%69 = load i32, i32* %19
%70 = sub nsw i32 %68, %69
%71 = sitofp i32 %70 to double
%72 = fdiv double %67, %71
%73 = load i32, i32* %19
%74 = add nsw i32 %73, 1
%75 = sext i32 %74 to i64
%76 = getelementptr inbounds [32 x double], [32 x double]* %16, i64 0, i64 %75
%77 = load double, double* %76
%78 = fmul double %72, %77
%79 = load i32, i32* %19
%80 = sext i32 %79 to i64
%81 = getelementptr inbounds [32 x double], [32 x double]* %16, i64 0, i64 %80
store double %78, double* %81
br label %82
82: ; preds = %62
%83 = load i32, i32* %19
%84 = add nsw i32 %83, -1
store i32 %84, i32* %19
br label %59
85: ; preds = %59
%86 = load i32, i32* %11
%87 = sext i32 %86 to i64
%88 = getelementptr inbounds [32 x double], [32 x double]* %17, i64 0, i64 %87
store double 1.000000e+00, double* %88
%89 = load i32, i32* %11
%90 = sub nsw i32 %89, 1
store i32 %90, i32* %20
br label %91
91: ; preds = %114, %85
%92 = load i32, i32* %20
%93 = icmp sge i32 %92, 0
br i1 %93, label %94, label %117
94: ; preds = %91
%95 = load double, double* %14
%96 = load i32, i32* %20
%97 = add nsw i32 %96, 1
%98 = sitofp i32 %97 to double
%99 = fmul double %95, %98
%100 = load i32, i32* %11
%101 = load i32, i32* %20
%102 = sub nsw i32 %100, %101
%103 = sitofp i32 %102 to double
%104 = fdiv double %99, %103
%105 = load i32, i32* %20
%106 = add nsw i32 %105, 1
%107 = sext i32 %106 to i64
%108 = getelementptr inbounds [32 x double], [32 x double]* %17, i64 0, i64 %107
%109 = load double, double* %108
%110 = fmul double %104, %109
%111 = load i32, i32* %20
%112 = sext i32 %111 to i64
%113 = getelementptr inbounds [32 x double], [32 x double]* %17, i64 0, i64 %112
store double %110, double* %113
br label %114
114: ; preds = %94
%115 = load i32, i32* %20
%116 = add nsw i32 %115, -1
store i32 %116, i32* %20
br label %91
117: ; preds = %91
store i32 0, i32* %21
br label %118
118: ; preds = %131, %117
%119 = load i32, i32* %21
%120 = load i32, i32* %9
%121 = load i32, i32* %10
%122 = add nsw i32 %120, %121
%123 = load i32, i32* %11
%124 = add nsw i32 %122, %123
%125 = icmp sle i32 %119, %124
br i1 %125, label %126, label %134
126: ; preds = %118
%127 = load double*, double** %8
%128 = load i32, i32* %21
%129 = sext i32 %128 to i64
%130 = getelementptr inbounds double, double* %127, i64 %129
store double 0.000000e+00, double* %130
br label %131
131: ; preds = %126
%132 = load i32, i32* %21
%133 = add nsw i32 %132, 1
store i32 %133, i32* %21
br label %118
134: ; preds = %118
store i32 0, i32* %18
br label %135
135: ; preds = %185, %134
%136 = load i32, i32* %18
%137 = load i32, i32* %9
%138 = add nsw i32 %137, 1
%139 = icmp slt i32 %136, %138
br i1 %139, label %140, label %188
140: ; preds = %135
store i32 0, i32* %19
br label %141
141: ; preds = %181, %140
%142 = load i32, i32* %19
%143 = load i32, i32* %10
%144 = add nsw i32 %143, 1
%145 = icmp slt i32 %142, %144
br i1 %145, label %146, label %184
146: ; preds = %141
store i32 0, i32* %20
br label %147
147: ; preds = %177, %146
%148 = load i32, i32* %20
%149 = load i32, i32* %11
%150 = add nsw i32 %149, 1
%151 = icmp slt i32 %148, %150
br i1 %151, label %152, label %180
152: ; preds = %147
%153 = load i32, i32* %18
%154 = load i32, i32* %19
%155 = add nsw i32 %153, %154
%156 = load i32, i32* %20
%157 = add nsw i32 %155, %156
store i32 %157, i32* %21
%158 = load i32, i32* %18
%159 = sext i32 %158 to i64
%160 = getelementptr inbounds [32 x double], [32 x double]* %15, i64 0, i64 %159
%161 = load double, double* %160
%162 = load i32, i32* %19
%163 = sext i32 %162 to i64
%164 = getelementptr inbounds [32 x double], [32 x double]* %16, i64 0, i64 %163
%165 = load double, double* %164
%166 = fmul double %161, %165
%167 = load i32, i32* %20
%168 = sext i32 %167 to i64
%169 = getelementptr inbounds [32 x double], [32 x double]* %17, i64 0, i64 %168
%170 = load double, double* %169
%171 = load double*, double** %8
%172 = load i32, i32* %21
%173 = sext i32 %172 to i64
%174 = getelementptr inbounds double, double* %171, i64 %173
%175 = load double, double* %174
%176 = call double @llvm.fmuladd.f64(double %166, double %170, double %175)
store double %176, double* %174
br label %177
177: ; preds = %152
%178 = load i32, i32* %20
%179 = add nsw i32 %178, 1
store i32 %179, i32* %20
br label %147
180: ; preds = %147
br label %181
181: ; preds = %180
%182 = load i32, i32* %19
%183 = add nsw i32 %182, 1
store i32 %183, i32* %19
br label %141
184: ; preds = %141
br label %185
185: ; preds = %184
%186 = load i32, i32* %18
%187 = add nsw i32 %186, 1
store i32 %187, i32* %18
br label %135
188: ; preds = %135
ret void
}

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