void all_mu_comp(double *eta, int p, int order, int m, int nknots,
                 double *knots, double *knotsI, double *all_beta,
                 double *all_mu) {
  int i, j, k, supi;
  double *B, *BI, *Ball;
  int *offsets, *offsetsI, *derivs;
  B = (double *)calloc(order * p, sizeof(double));
  BI = (double *)calloc(order * p, sizeof(double));
  Ball = (double *)calloc(m * p, sizeof(double));
  derivs = (int *)calloc(p, sizeof(int));
  for (i = 0; i < p; i++)
    derivs[i] = 0;
  offsets = (int *)calloc(p, sizeof(int));
  offsetsI = (int *)calloc(p, sizeof(int));
  spline_basis(knotsI, (nknots - 1), order, eta, p, derivs, p, offsetsI, BI);
  spline_basis(knots, nknots, order, eta, p, derivs, p, offsets, B);
  for (i = 0; i < order; i++) {
    Ball[offsetsI[0] + i] = BI[i];
  }
  for (j = 1; j < p; j++) {
    if (offsets[j] != 0) {
      for (i = 0; i < order; i++) {
        Ball[m * j + offsets[j] + i - 1] = B[j * order + i];
      }
    } else {
      for (i = 1; i < order; i++) {
        Ball[m * j + offsets[j] + i - 1] = B[j * order + i];
        ;
      }
    }
  }
  for (j = 0; j < p; j++) {
    for (k = 0; k <= j; k++) {
      for (supi = k * m; supi < (k + 1) * m; supi++)
        all_mu[j * (j + 1) / 2 + k] +=
            Ball[supi] * all_beta[j * (j + 1) * m / 2 + supi];
    }
  }
  free(B);
  free(BI);
  free(Ball);
  free(offsets);
  free(offsetsI);
  free(derivs);
}


x86 assembly (clang, -O0):
<all_mu_comp>:
push   %rbp
mov    %rsp,%rbp
sub    $0x90,%rsp
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8d,-0x18(%rbp)
mov    %r9,-0x20(%rbp)
mov    -0x10(%rbp),%eax
imul   -0xc(%rbp),%eax
movslq %eax,%rdi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %rax,-0x38(%rbp)
mov    -0x10(%rbp),%eax
imul   -0xc(%rbp),%eax
movslq %eax,%rdi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %rax,-0x40(%rbp)
mov    -0x14(%rbp),%eax
imul   -0xc(%rbp),%eax
movslq %eax,%rdi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %rax,-0x48(%rbp)
movslq -0xc(%rbp),%rdi
mov    $0x4,%esi
call   1050 <calloc@plt>
mov    %rax,-0x60(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11e6 <all_mu_comp+0xb6>
mov    -0x60(%rbp),%rax
movslq -0x24(%rbp),%rcx
movl   $0x0,(%rax,%rcx,4)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11bd <all_mu_comp+0x8d>
movslq -0xc(%rbp),%rdi
mov    $0x4,%esi
call   1050 <calloc@plt>
mov    %rax,-0x50(%rbp)
movslq -0xc(%rbp),%rdi
mov    $0x4,%esi
call   1050 <calloc@plt>
mov    %rax,-0x58(%rbp)
mov    0x10(%rbp),%rdi
mov    -0x18(%rbp),%esi
sub    $0x1,%esi
mov    -0x10(%rbp),%edx
mov    -0x8(%rbp),%rcx
mov    -0xc(%rbp),%r8d
mov    -0x60(%rbp),%r9
mov    -0xc(%rbp),%r11d
mov    -0x58(%rbp),%r10
mov    -0x40(%rbp),%rax
mov    %r11d,(%rsp)
mov    %r10,0x8(%rsp)
mov    %rax,0x10(%rsp)
mov    $0x0,%al
call   1040 <spline_basis@plt>
mov    -0x20(%rbp),%rdi
mov    -0x18(%rbp),%esi
mov    -0x10(%rbp),%edx
mov    -0x8(%rbp),%rcx
mov    -0xc(%rbp),%r8d
mov    -0x60(%rbp),%r9
mov    -0xc(%rbp),%r11d
mov    -0x50(%rbp),%r10
mov    -0x38(%rbp),%rax
mov    %r11d,(%rsp)
mov    %r10,0x8(%rsp)
mov    %rax,0x10(%rsp)
mov    $0x0,%al
call   1040 <spline_basis@plt>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    12be <all_mu_comp+0x18e>
mov    -0x40(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x48(%rbp),%rax
mov    -0x58(%rbp),%rcx
mov    (%rcx),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1282 <all_mu_comp+0x152>
movl   $0x1,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    13b6 <all_mu_comp+0x286>
mov    -0x50(%rbp),%rax
movslq -0x28(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
je     1343 <all_mu_comp+0x213>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    133e <all_mu_comp+0x20e>
mov    -0x38(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x10(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x48(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   -0x28(%rbp),%ecx
mov    -0x50(%rbp),%rdx
movslq -0x28(%rbp),%rsi
add    (%rdx,%rsi,4),%ecx
add    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    12ea <all_mu_comp+0x1ba>
jmp    13a3 <all_mu_comp+0x273>
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    139e <all_mu_comp+0x26e>
mov    -0x38(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x10(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x48(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   -0x28(%rbp),%ecx
mov    -0x50(%rbp),%rdx
movslq -0x28(%rbp),%rsi
add    (%rdx,%rsi,4),%ecx
add    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    134a <all_mu_comp+0x21a>
jmp    13a3 <all_mu_comp+0x273>
jmp    13a8 <all_mu_comp+0x278>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    12c5 <all_mu_comp+0x195>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    14a7 <all_mu_comp+0x377>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x28(%rbp),%eax
jg     1494 <all_mu_comp+0x364>
mov    -0x2c(%rbp),%eax
imul   -0x14(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
imul   -0x14(%rbp),%ecx
cmp    %ecx,%eax
jge    1481 <all_mu_comp+0x351>
mov    -0x48(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
mov    %rax,-0x70(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
imul   -0x14(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x70(%rbp),%rax
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    0x20(%rbp),%rax
mov    %rax,-0x68(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x68(%rbp),%rax
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    13e6 <all_mu_comp+0x2b6>
jmp    1486 <all_mu_comp+0x356>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    13d0 <all_mu_comp+0x2a0>
jmp    1499 <all_mu_comp+0x369>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    13bd <all_mu_comp+0x28d>
mov    -0x38(%rbp),%rdi
call   1030 <free@plt>
mov    -0x40(%rbp),%rdi
call   1030 <free@plt>
mov    -0x48(%rbp),%rdi
call   1030 <free@plt>
mov    -0x50(%rbp),%rdi
call   1030 <free@plt>
mov    -0x58(%rbp),%rdi
call   1030 <free@plt>
mov    -0x60(%rbp),%rdi
call   1030 <free@plt>
add    $0x90,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<all_mu_comp>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xd8,%rsp
mov    %r9,0x38(%rsp)
mov    %r8,0x18(%rsp)
mov    %ecx,%r12d
mov    %edx,(%rsp)
mov    %esi,%ebx
mov    %rdi,0x8(%rsp)
mov    0x120(%rsp),%r15
mov    0x118(%rsp),%rax
mov    %rax,0x70(%rsp)
mov    0x110(%rsp),%rax
mov    %rax,0x50(%rsp)
mov    %edx,%eax
imul   %esi,%eax
movslq %eax,%rbp
mov    $0x8,%esi
mov    %rbp,%rdi
call   1050 <calloc@plt>
mov    %rax,0x10(%rsp)
mov    $0x8,%esi
mov    %rbp,%rdi
call   1050 <calloc@plt>
mov    %rax,%rbp
mov    %r12d,%eax
imul   %ebx,%eax
movslq %eax,%rdi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %rax,%r13
mov    %ebx,0x4(%rsp)
movslq %ebx,%r14
mov    $0x4,%esi
mov    %r14,%rdi
call   1050 <calloc@plt>
mov    %rax,0x40(%rsp)
mov    $0x4,%esi
mov    %r14,%rdi
call   1050 <calloc@plt>
mov    %rax,0x28(%rsp)
mov    $0x4,%esi
mov    %r14,%rdi
call   1050 <calloc@plt>
mov    %rax,%r10
mov    %rax,0x58(%rsp)
mov    0x18(%rsp),%rax
lea    -0x1(%rax),%esi
sub    $0x8,%rsp
mov    0x58(%rsp),%rdi
mov    0x8(%rsp),%ebx
mov    %ebx,%edx
mov    0x10(%rsp),%rcx
mov    %r14d,%r8d
mov    0x48(%rsp),%r9
xor    %eax,%eax
push   %rbp
push   %r10
push   %r14
call   1040 <spline_basis@plt>
add    $0x18,%rsp
mov    0x40(%rsp),%rdi
mov    0x20(%rsp),%rsi
mov    %ebx,%edx
mov    0x10(%rsp),%rcx
mov    %r14d,%r8d
mov    0x48(%rsp),%r9
xor    %eax,%eax
push   0x18(%rsp)
push   0x38(%rsp)
push   %r14
call   1040 <spline_basis@plt>
add    $0x20,%rsp
mov    %ebx,%r14d
test   %ebx,%ebx
jle    1296 <all_mu_comp+0x156>
mov    0x58(%rsp),%rax
movslq (%rax),%rax
lea    0x0(,%rax,8),%rdi
add    %r13,%rdi
lea    0x0(,%r14,8),%rdx
mov    %rbp,%rsi
call   1060 <memcpy@plt>
mov    %r12d,%eax
mov    %rax,0x18(%rsp)
cmpl   $0x2,0x4(%rsp)
mov    0x28(%rsp),%r9
jl     1752 <all_mu_comp+0x612>
movslq (%rsp),%rcx
mov    0x4(%rsp),%eax
mov    %rax,0x50(%rsp)
lea    -0x2(%r14),%rax
mov    %rax,0xa8(%rsp)
lea    -0x1(%r14),%rax
lea    -0x1(%r12),%edx
mov    %r14d,%esi
and    $0xfffffffc,%esi
mov    %rsi,0xc0(%rsp)
add    $0xfffffffffffffffc,%rsi
mov    %rsi,0xb8(%rsp)
mov    %rsi,%rdi
shr    $0x2,%rdi
add    $0x1,%rdi
lea    -0x5(%r14),%rsi
mov    %rsi,0x98(%rsp)
mov    %rsi,%rbx
shr    $0x2,%rbx
add    $0x1,%rbx
mov    %rdi,%rsi
mov    %rdi,0xb0(%rsp)
and    $0xfffffffffffffffe,%rdi
mov    %rdi,0xa0(%rsp)
mov    %r14,%rsi
neg    %rsi
mov    %rsi,0x60(%rsp)
mov    %rax,0x48(%rsp)
and    $0xfffffffffffffffc,%rax
mov    %rax,0x88(%rsp)
add    $0x1,%rax
mov    %rax,0x80(%rsp)
mov    %rbx,%rax
mov    %rbx,0x90(%rsp)
and    $0xfffffffffffffffe,%rbx
mov    %rbx,0x78(%rsp)
lea    0x3(%r12),%eax
mov    %rax,0x8(%rsp)
lea    0x0(,%rcx,8),%rax
mov    %rax,0xd0(%rsp)
mov    0x10(%rsp),%rax
lea    0x30(%rax,%rcx,8),%r11
mov    %rcx,0x68(%rsp)
lea    (%rax,%rcx,8),%rsi
add    $0x8,%rsi
mov    $0x1,%r8d
xor    %ebx,%ebx
mov    0x18(%rsp),%r10
mov    %rdx,0xc8(%rsp)
jmp    13ed <all_mu_comp+0x2ad>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r8
mov    0x38(%rsp),%rbx
add    $0x1,%ebx
mov    0x8(%rsp),%rax
add    %r12d,%eax
mov    %rax,0x8(%rsp)
mov    0xd0(%rsp),%rax
add    %rax,%r11
add    %rax,%rsi
mov    0x18(%rsp),%rax
add    %rax,%rdx
add    %rax,%r10
cmp    0x50(%rsp),%r8
je     1752 <all_mu_comp+0x612>
mov    %rbx,0x38(%rsp)
imul   %r12d,%ebx
mov    (%r9,%r8,4),%edi
test   %edi,%edi
je     14c0 <all_mu_comp+0x380>
cmpl   $0x0,(%rsp)
jle    13b0 <all_mu_comp+0x270>
mov    %r8,%rax
imul   0x68(%rsp),%rax
mov    %rax,0x30(%rsp)
mov    %r8d,%eax
imul   %r12d,%eax
add    %edi,%eax
add    $0xffffffff,%eax
mov    %rax,0x20(%rsp)
cmpl   $0x4,(%rsp)
jae    1534 <all_mu_comp+0x3f4>
xor    %eax,%eax
mov    %rax,%rbx
test   $0x1,%r14b
je     1466 <all_mu_comp+0x326>
mov    0x30(%rsp),%rbx
add    %rax,%rbx
mov    0x10(%rsp),%rcx
movsd  (%rcx,%rbx,8),%xmm0
mov    0x20(%rsp),%rcx
add    %eax,%ecx
movslq %ecx,%rcx
movsd  %xmm0,0x0(%r13,%rcx,8)
mov    %rax,%rbx
or     $0x1,%rbx
not    %rax
cmp    0x60(%rsp),%rax
je     13b0 <all_mu_comp+0x270>
lea    (%rdi,%rdx,1),%eax
lea    (%r10,%rdi,1),%ecx
nopl   0x0(%rax,%rax,1)
movsd  -0x8(%rsi,%rbx,8),%xmm0
lea    (%rax,%rbx,1),%edi
movslq %edi,%rdi
movsd  %xmm0,0x0(%r13,%rdi,8)
movsd  (%rsi,%rbx,8),%xmm0
lea    (%rcx,%rbx,1),%edi
movslq %edi,%rdi
movsd  %xmm0,0x0(%r13,%rdi,8)
add    $0x2,%rbx
cmp    %rbx,%r14
jne    1480 <all_mu_comp+0x340>
jmp    13b0 <all_mu_comp+0x270>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
cmpl   $0x2,(%rsp)
jl     13b0 <all_mu_comp+0x270>
mov    %r8,%rax
imul   0x68(%rsp),%rax
mov    %rax,0x30(%rsp)
mov    %r8d,%edi
imul   %r12d,%edi
add    $0xffffffff,%edi
mov    $0x1,%eax
cmpq   $0x4,0x48(%rsp)
jb     16e6 <all_mu_comp+0x5a6>
mov    %edi,0x20(%rsp)
add    %r12d,%ebx
mov    0xa8(%rsp),%rdi
lea    (%rbx,%rdi,1),%ecx
cmp    %ebx,%ecx
jl     1626 <all_mu_comp+0x4e6>
mov    %rdi,%rcx
shr    $0x20,%rcx
mov    0x20(%rsp),%edi
jne    16e6 <all_mu_comp+0x5a6>
cmpq   $0x4,0x98(%rsp)
jae    162f <all_mu_comp+0x4ef>
xor    %ebx,%ebx
jmp    1686 <all_mu_comp+0x546>
add    0xc8(%rsp),%ebx
add    %edi,%ebx
mov    0x48(%rsp),%rcx
lea    (%rbx,%rcx,1),%eax
cmp    %ebx,%eax
jl     1431 <all_mu_comp+0x2f1>
mov    %rcx,%rax
shr    $0x20,%rax
mov    $0x0,%eax
jne    1433 <all_mu_comp+0x2f3>
cmpq   $0x0,0xb8(%rsp)
je     15cd <all_mu_comp+0x48d>
mov    0x8(%rsp),%rax
lea    (%rdi,%rax,1),%r9d
mov    0xa0(%rsp),%rax
xor    %ebx,%ebx
nopl   (%rax)
movups -0x30(%r11,%rbx,8),%xmm0
movups -0x20(%r11,%rbx,8),%xmm1
lea    -0x4(%r9),%ecx
movslq %ecx,%rcx
movups %xmm0,0x0(%r13,%rcx,8)
movups %xmm1,0x10(%r13,%rcx,8)
movupd -0x10(%r11,%rbx,8),%xmm0
movupd (%r11,%rbx,8),%xmm1
movslq %r9d,%r9
movupd %xmm0,0x0(%r13,%r9,8)
movupd %xmm1,0x10(%r13,%r9,8)
add    $0x8,%rbx
add    $0x8,%r9d
add    $0xfffffffffffffffe,%rax
jne    1580 <all_mu_comp+0x440>
jmp    15cf <all_mu_comp+0x48f>
xor    %ebx,%ebx
testb  $0x1,0xb0(%rsp)
01
mov    0x28(%rsp),%r9
je     160d <all_mu_comp+0x4cd>
mov    %rbx,%rax
add    0x30(%rsp),%rax
mov    0x10(%rsp),%rcx
movupd (%rcx,%rax,8),%xmm0
movupd 0x10(%rcx,%rax,8),%xmm1
mov    0x20(%rsp),%rax
add    %ebx,%eax
cltq
movupd %xmm0,0x0(%r13,%rax,8)
movupd %xmm1,0x10(%r13,%rax,8)
mov    0xc0(%rsp),%rcx
mov    %rcx,%rax
cmp    %r14,%rcx
je     13b0 <all_mu_comp+0x270>
jmp    1433 <all_mu_comp+0x2f3>
mov    0x20(%rsp),%edi
jmp    16e6 <all_mu_comp+0x5a6>
mov    0x78(%rsp),%rax
mov    %r10d,%ecx
xor    %ebx,%ebx
nopl   0x0(%rax)
movups (%rsi,%rbx,8),%xmm0
movups 0x10(%rsi,%rbx,8),%xmm1
movslq %ecx,%rcx
movups %xmm0,0x0(%r13,%rcx,8)
movups %xmm1,0x10(%r13,%rcx,8)
movupd 0x20(%rsi,%rbx,8),%xmm0
movupd 0x30(%rsi,%rbx,8),%xmm1
lea    0x4(%rcx),%r9d
movslq %r9d,%rdi
movupd %xmm0,0x0(%r13,%rdi,8)
movupd %xmm1,0x10(%r13,%rdi,8)
add    $0x8,%rbx
add    $0x8,%ecx
add    $0xfffffffffffffffe,%rax
jne    1640 <all_mu_comp+0x500>
testb  $0x1,0x90(%rsp)
01
mov    0x20(%rsp),%edi
je     16c6 <all_mu_comp+0x586>
mov    %rbx,%rax
or     $0x1,%rax
or     $0x1,%ebx
add    0x30(%rsp),%rax
mov    0x10(%rsp),%rcx
movupd (%rcx,%rax,8),%xmm0
movupd 0x10(%rcx,%rax,8),%xmm1
add    %edi,%ebx
movslq %ebx,%rax
movupd %xmm0,0x0(%r13,%rax,8)
movupd %xmm1,0x10(%r13,%rax,8)
mov    0x80(%rsp),%rax
mov    0x88(%rsp),%rcx
cmp    %rcx,0x48(%rsp)
mov    0x28(%rsp),%r9
je     13b0 <all_mu_comp+0x270>
mov    %rax,%rbx
test   $0x1,%r14b
jne    1711 <all_mu_comp+0x5d1>
mov    0x30(%rsp),%rbx
add    %rax,%rbx
mov    0x10(%rsp),%rcx
movsd  (%rcx,%rbx,8),%xmm0
add    %eax,%edi
movslq %edi,%rcx
movsd  %xmm0,0x0(%r13,%rcx,8)
lea    0x1(%rax),%rbx
not    %rax
cmp    0x60(%rsp),%rax
je     13b0 <all_mu_comp+0x270>
nop
movsd  -0x8(%rsi,%rbx,8),%xmm0
lea    (%rdx,%rbx,1),%eax
cltq
movsd  %xmm0,0x0(%r13,%rax,8)
movsd  (%rsi,%rbx,8),%xmm0
lea    (%r10,%rbx,1),%eax
cltq
movsd  %xmm0,0x0(%r13,%rax,8)
add    $0x2,%rbx
cmp    %rbx,%r14
jne    1720 <all_mu_comp+0x5e0>
jmp    13b0 <all_mu_comp+0x270>
cmpl   $0x0,0x4(%rsp)
mov    0x58(%rsp),%rbx
jle    187a <all_mu_comp+0x73a>
mov    0x70(%rsp),%rax
add    $0x8,%rax
mov    %rax,0x38(%rsp)
xor    %r11d,%r11d
mov    $0x1,%r14d
jmp    1794 <all_mu_comp+0x654>
nopl   0x0(%rax,%rax,1)
add    $0x1,%r14
mov    0x8(%rsp),%r11
cmp    0x4(%rsp),%r11d
je     187a <all_mu_comp+0x73a>
lea    0x1(%r11),%eax
mov    %rax,0x8(%rsp)
imul   %eax,%r11d
mov    %r11d,%eax
imul   %r12d,%eax
mov    %eax,%ecx
shr    $0x1f,%ecx
add    %eax,%ecx
sar    %ecx
shr    %r11d
movslq %ecx,%r10
mov    0x38(%rsp),%rax
lea    (%rax,%r10,8),%rsi
mov    0x18(%rsp),%rax
xor    %r9d,%r9d
xor    %ecx,%ecx
jmp    17dd <all_mu_comp+0x69d>
nopl   0x0(%rax)
add    %r12d,%ecx
add    0x18(%rsp),%rax
cmp    %r14,%r9
je     1780 <all_mu_comp+0x640>
mov    %r9,%rdi
mov    %edi,%edx
imul   %r12d,%edx
add    $0x1,%r9
mov    %r9d,%r8d
imul   %r12d,%r8d
cmp    %r8d,%edx
jge    17d0 <all_mu_comp+0x690>
add    %r11,%rdi
movslq %ecx,%rdx
movsd  (%r15,%rdi,8),%xmm0
test   $0x1,%r12b
je     1836 <all_mu_comp+0x6f6>
movsd  0x0(%r13,%rdx,8),%xmm1
mov    %rbp,%r8
lea    (%rdx,%r10,1),%rbp
mov    0x70(%rsp),%rbx
mulsd  (%rbx,%rbp,8),%xmm1
mov    0x58(%rsp),%rbx
mov    %r8,%rbp
addsd  %xmm1,%xmm0
movsd  %xmm0,(%r15,%rdi,8)
add    $0x1,%rdx
cmp    $0x1,%r12d
je     17d0 <all_mu_comp+0x690>
nopl   0x0(%rax)
movsd  0x0(%r13,%rdx,8),%xmm1
mulsd  -0x8(%rsi,%rdx,8),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%r15,%rdi,8)
movsd  0x8(%r13,%rdx,8),%xmm0
mulsd  (%rsi,%rdx,8),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%r15,%rdi,8)
add    $0x2,%rdx
cmp    %edx,%eax
jne    1840 <all_mu_comp+0x700>
jmp    17d0 <all_mu_comp+0x690>
mov    0x10(%rsp),%rdi
call   1030 <free@plt>
mov    %rbp,%rdi
call   1030 <free@plt>
mov    %r13,%rdi
call   1030 <free@plt>
mov    0x28(%rsp),%rdi
call   1030 <free@plt>
mov    %rbx,%rdi
call   1030 <free@plt>
mov    0x40(%rsp),%rdi
add    $0xd8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

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

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