#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
double fabs(double);
int mp_covar(int n, double *r, int ldr, int *ipvt, double tol, double *wa) {
  int i, ii, j, jj, k, l;
  int kk, kj, ji, j0, k0, jj0;
  int sing;
  double one = 1.0, temp, tolr, zero = 0.0;
  tolr = tol * fabs(r[0]);
  l = -1;
  for (k = 0; k < n; k++) {
    kk = k * ldr + k;
    if (fabs(r[kk]) <= tolr)
      break;
    r[kk] = one / r[kk];
    for (j = 0; j < k; j++) {
      kj = k * ldr + j;
      temp = r[kk] * r[kj];
      r[kj] = zero;
      k0 = k * ldr;
      j0 = j * ldr;
      for (i = 0; i <= j; i++) {
        r[k0 + i] += (-temp * r[j0 + i]);
      }
    }
    l = k;
  }
  if (l >= 0) {
    for (k = 0; k <= l; k++) {
      k0 = k * ldr;
      for (j = 0; j < k; j++) {
        temp = r[k * ldr + j];
        j0 = j * ldr;
        for (i = 0; i <= j; i++) {
          r[j0 + i] += temp * r[k0 + i];
        }
      }
      temp = r[k0 + k];
      for (i = 0; i <= k; i++) {
        r[k0 + i] *= temp;
      }
    }
  }
  for (j = 0; j < n; j++) {
    jj = ipvt[j];
    sing = (j > l);
    j0 = j * ldr;
    jj0 = jj * ldr;
    for (i = 0; i <= j; i++) {
      ji = j0 + i;
      if (sing)
        r[ji] = zero;
      ii = ipvt[i];
      if (ii > jj)
        r[jj0 + ii] = r[ji];
      if (ii < jj)
        r[ii * ldr + jj] = r[ji];
    }
    wa[jj] = r[j0 + j];
  }
  for (j = 0; j < n; j++) {
    j0 = j * ldr;
    for (i = 0; i < j; i++) {
      r[j0 + i] = r[i * ldr + j];
    }
    r[j0 + j] = wa[j];
  }
  return 0;
}


x86 assembly (clang, -O0):
<mp_covar>:
push   %rbp
mov    %rsp,%rbp
push   %rax
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
movsd  %xmm0,-0x28(%rbp)
mov    %r8,-0x30(%rbp)
movsd  0xeec(%rip),%xmm0
movsd  %xmm0,-0x70(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x28(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm1
movaps 0xeb8(%rip),%xmm2
pand   %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movl   $0xffffffff,-0x48(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    12c4 <mp_covar+0x1c4>
mov    -0x44(%rbp),%eax
imul   -0x14(%rbp),%eax
add    -0x44(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movaps 0xe70(%rip),%xmm0
pand   %xmm0,%xmm1
movsd  -0x80(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jb     11a8 <mp_covar+0xa8>
jmp    12c4 <mp_covar+0x1c4>
movsd  -0x70(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x4c(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x44(%rbp),%eax
jge    12b0 <mp_covar+0x1b0>
mov    -0x44(%rbp),%eax
imul   -0x14(%rbp),%eax
add    -0x3c(%rbp),%eax
mov    %eax,-0x50(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x50(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x88(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movslq -0x50(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x44(%rbp),%eax
imul   -0x14(%rbp),%eax
mov    %eax,-0x5c(%rbp)
mov    -0x3c(%rbp),%eax
imul   -0x14(%rbp),%eax
mov    %eax,-0x58(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jg     129d <mp_covar+0x19d>
movsd  -0x78(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x10(%rbp),%rax
mov    -0x58(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x10(%rbp),%rax
mov    -0x5c(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1236 <mp_covar+0x136>
jmp    12a2 <mp_covar+0x1a2>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    11ce <mp_covar+0xce>
mov    -0x44(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1163 <mp_covar+0x63>
cmpl   $0x0,-0x48(%rbp)
jl     13f9 <mp_covar+0x2f9>
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x48(%rbp),%eax
jg     13f4 <mp_covar+0x2f4>
mov    -0x44(%rbp),%eax
imul   -0x14(%rbp),%eax
mov    %eax,-0x5c(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x44(%rbp),%eax
jge    138d <mp_covar+0x28d>
mov    -0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x3c(%rbp),%eax
imul   -0x14(%rbp),%eax
mov    %eax,-0x58(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jg     137a <mp_covar+0x27a>
movsd  -0x78(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x5c(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x10(%rbp),%rax
mov    -0x58(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    132a <mp_covar+0x22a>
jmp    137f <mp_covar+0x27f>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    12f2 <mp_covar+0x1f2>
mov    -0x10(%rbp),%rax
mov    -0x5c(%rbp),%ecx
add    -0x44(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x44(%rbp),%eax
jg     13e1 <mp_covar+0x2e1>
movsd  -0x78(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x5c(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    13ab <mp_covar+0x2ab>
jmp    13e6 <mp_covar+0x2e6>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    12d5 <mp_covar+0x1d5>
jmp    13f9 <mp_covar+0x2f9>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1522 <mp_covar+0x422>
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x40(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x48(%rbp),%eax
setg   %al
and    $0x1,%al
movzbl %al,%eax
mov    %eax,-0x64(%rbp)
mov    -0x3c(%rbp),%eax
imul   -0x14(%rbp),%eax
mov    %eax,-0x58(%rbp)
mov    -0x40(%rbp),%eax
imul   -0x14(%rbp),%eax
mov    %eax,-0x60(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jg     14f5 <mp_covar+0x3f5>
mov    -0x58(%rbp),%eax
add    -0x34(%rbp),%eax
mov    %eax,-0x54(%rbp)
cmpl   $0x0,-0x64(%rbp)
je     147a <mp_covar+0x37a>
movsd  -0x88(%rbp),%xmm0
ff
mov    -0x10(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x40(%rbp),%eax
jle    14b3 <mp_covar+0x3b3>
mov    -0x10(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x60(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
cmp    -0x40(%rbp),%eax
jge    14e2 <mp_covar+0x3e2>
mov    -0x10(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x40(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    14e7 <mp_covar+0x3e7>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1446 <mp_covar+0x346>
mov    -0x10(%rbp),%rax
mov    -0x58(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1400 <mp_covar+0x300>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    15b5 <mp_covar+0x4b5>
mov    -0x3c(%rbp),%eax
imul   -0x14(%rbp),%eax
mov    %eax,-0x58(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jge    1588 <mp_covar+0x488>
mov    -0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x58(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1546 <mp_covar+0x446>
mov    -0x30(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x58(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1529 <mp_covar+0x429>
xor    %eax,%eax
add    $0x8,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<mp_covar>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x10,%rsp
mov    %r8,-0x10(%rsp)
test   %edi,%edi
jle    1a64 <mp_covar+0x964>
movsd  (%rsi),%xmm10
andpd  0xed7(%rip),%xmm10
mulsd  %xmm0,%xmm10
movslq %edx,%r10
mov    %rdi,-0x38(%rsp)
mov    %edi,%eax
mov    %rax,-0x30(%rsp)
movsd  (%rsi),%xmm5
movapd 0xeb7(%rip),%xmm0
andpd  %xmm5,%xmm0
mov    $0xffffffff,%eax
ucomisd %xmm0,%xmm10
mov    %rsi,-0x68(%rsp)
mov    %r10,-0x70(%rsp)
mov    %edx,-0x7c(%rsp)
mov    %rcx,-0x60(%rsp)
jae    14cc <mp_covar+0x3cc>
mov    -0x38(%rsp),%rax
add    $0xffffffff,%eax
mov    %rax,-0x50(%rsp)
lea    0x8(%rsi),%r8
lea    0x0(,%r10,8),%r12
lea    0x8(,%r10,8),%rax
mov    %rax,-0x18(%rsp)
lea    0x30(%rsi),%rbx
xor    %edi,%edi
movsd  0xe7b(%rip),%xmm9
movapd 0xe52(%rip),%xmm8
movapd 0xe5a(%rip),%xmm3
mov    %r8,%rbp
mov    %rbx,-0x28(%rsp)
mov    %rsi,%r11
xor    %eax,%eax
mov    %r8,-0x20(%rsp)
nopl   0x0(%rax,%rax,1)
mov    %rax,%r9
movapd %xmm9,%xmm4
divsd  %xmm5,%xmm4
mov    %r11,-0x78(%rsp)
movsd  %xmm4,(%r11)
test   %rax,%rax
je     1470 <mp_covar+0x370>
lea    (%rsi,%rdi,8),%rax
mov    %rax,-0x40(%rsp)
lea    0x1(%rdi),%rax
mov    %rax,-0x48(%rsp)
mov    $0x1,%edx
mov    %rsi,%r13
mov    -0x28(%rsp),%r11
xor    %r14d,%r14d
mov    %rdi,-0x58(%rsp)
cs nopw 0x0(%rax,%rax,1)
xorpd  %xmm3,%xmm4
lea    (%r14,%rdi,1),%rax
mulsd  (%rsi,%rax,8),%xmm4
movq   $0x0,(%rsi,%rax,8)
mov    %r14,%rcx
imul   %r10,%rcx
cmp    $0x4,%rdx
jae    1250 <mp_covar+0x150>
xor    %r15d,%r15d
jmp    13b0 <mp_covar+0x2b0>
nopw   0x0(%rax,%rax,1)
mov    -0x18(%rsp),%rax
imul   %r14,%rax
add    %r8,%rax
cmp    -0x40(%rsp),%rax
jbe    1290 <mp_covar+0x190>
lea    (%rsi,%rcx,8),%rax
mov    -0x48(%rsp),%rdi
add    %r14,%rdi
lea    (%rsi,%rdi,8),%rdi
cmp    %rdi,%rax
jae    1290 <mp_covar+0x190>
xor    %r15d,%r15d
mov    -0x58(%rsp),%rdi
jmp    13b0 <mp_covar+0x2b0>
cs nopw 0x0(%rax,%rax,1)
nop
mov    %rdx,%r15
movabs $0x7ffffffffffffffc,%rax
ff ff 7f
and    %rax,%r15
lea    -0x4(%r15),%rdi
mov    %rdi,%rax
shr    $0x2,%rax
add    $0x1,%rax
movapd %xmm4,%xmm5
unpcklpd %xmm4,%xmm5
test   %rdi,%rdi
je     144d <mp_covar+0x34d>
mov    %rdx,%r8
and    $0xfffffffffffffffc,%r8
add    $0xfffffffffffffffc,%r8
shr    $0x2,%r8
add    $0x1,%r8
and    $0xfffffffffffffffe,%r8
xor    %r10d,%r10d
mov    -0x58(%rsp),%rdi
nop
movupd -0x30(%r11,%r10,8),%xmm6
movupd -0x20(%r11,%r10,8),%xmm7
mulpd  %xmm5,%xmm6
movupd -0x30(%rbx,%r10,8),%xmm2
addpd  %xmm6,%xmm2
movupd -0x20(%rbx,%r10,8),%xmm6
movupd -0x10(%rbx,%r10,8),%xmm0
movupd (%rbx,%r10,8),%xmm1
mulpd  %xmm5,%xmm7
addpd  %xmm6,%xmm7
movupd %xmm2,-0x30(%rbx,%r10,8)
movupd %xmm7,-0x20(%rbx,%r10,8)
movupd -0x10(%r11,%r10,8),%xmm2
movupd (%r11,%r10,8),%xmm6
mulpd  %xmm5,%xmm2
addpd  %xmm0,%xmm2
mulpd  %xmm5,%xmm6
addpd  %xmm1,%xmm6
movupd %xmm2,-0x10(%rbx,%r10,8)
movupd %xmm6,(%rbx,%r10,8)
add    $0x8,%r10
add    $0xfffffffffffffffe,%r8
jne    12e0 <mp_covar+0x1e0>
test   $0x1,%al
je     13a0 <mp_covar+0x2a0>
lea    (%r10,%rcx,1),%rax
mov    -0x68(%rsp),%rsi
movupd (%rsi,%rax,8),%xmm0
movupd 0x10(%rsi,%rax,8),%xmm1
add    %rdi,%r10
movupd (%rsi,%r10,8),%xmm2
movupd 0x10(%rsi,%r10,8),%xmm6
mulpd  %xmm5,%xmm0
addpd  %xmm2,%xmm0
mulpd  %xmm5,%xmm1
addpd  %xmm6,%xmm1
movupd %xmm0,(%rsi,%r10,8)
movupd %xmm1,0x10(%rsi,%r10,8)
cmp    %r15,%rdx
mov    -0x70(%rsp),%r10
mov    -0x20(%rsp),%r8
je     1427 <mp_covar+0x327>
nop
mov    %r15,%rax
test   $0x1,%dl
je     13de <mp_covar+0x2de>
add    %r15,%rcx
lea    (%r15,%rdi,1),%rax
mov    -0x68(%rsp),%rsi
movsd  (%rsi,%rcx,8),%xmm0
mulsd  %xmm4,%xmm0
addsd  (%rsi,%rax,8),%xmm0
movsd  %xmm0,(%rsi,%rax,8)
mov    %r15,%rax
or     $0x1,%rax
cmp    %r15,%r14
je     1427 <mp_covar+0x327>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  0x0(%r13,%rax,8),%xmm0
mulsd  %xmm4,%xmm0
addsd  -0x8(%rbp,%rax,8),%xmm0
movsd  %xmm0,-0x8(%rbp,%rax,8)
movsd  0x8(%r13,%rax,8),%xmm0
mulsd  %xmm4,%xmm0
addsd  0x0(%rbp,%rax,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rax,8)
add    $0x2,%rax
cmp    %rax,%rdx
jne    13f0 <mp_covar+0x2f0>
add    $0x1,%r14
cmp    %r9,%r14
mov    -0x68(%rsp),%rsi
je     1470 <mp_covar+0x370>
add    $0x1,%rdx
mov    -0x78(%rsp),%rax
movsd  (%rax),%xmm4
add    %r12,%r11
add    %r12,%r13
jmp    1220 <mp_covar+0x120>
xor    %r10d,%r10d
mov    -0x58(%rsp),%rdi
test   $0x1,%al
jne    135f <mp_covar+0x25f>
jmp    13a0 <mp_covar+0x2a0>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
lea    0x1(%r9),%rax
cmp    -0x30(%rsp),%rax
je     14be <mp_covar+0x3be>
mov    %rax,%rdi
imul   %r10,%rdi
mov    %eax,%ecx
mov    -0x7c(%rsp),%edx
imul   %edx,%ecx
movslq %ecx,%rcx
add    %rax,%rcx
lea    (%rsi,%rcx,8),%r11
movsd  (%rsi,%rcx,8),%xmm5
movapd %xmm5,%xmm4
andpd  %xmm8,%xmm4
add    %r12,%rbx
add    %r12,%rbp
ucomisd %xmm4,%xmm10
jb     11d0 <mp_covar+0xd0>
mov    %r9d,%eax
mov    -0x60(%rsp),%rcx
jmp    14cc <mp_covar+0x3cc>
mov    -0x60(%rsp),%rcx
mov    -0x7c(%rsp),%edx
mov    -0x50(%rsp),%rax
mov    %rax,-0x50(%rsp)
test   %eax,%eax
js     18ee <mp_covar+0x7ee>
mov    -0x50(%rsp),%rax
add    $0x1,%eax
mov    %rax,(%rsp)
lea    0x8(%rsi),%rax
mov    %rax,-0x48(%rsp)
lea    0x0(,%r10,8),%rax
mov    %rax,-0x40(%rsp)
lea    0x8(,%r10,8),%rax
mov    %rax,-0x28(%rsp)
lea    0x30(%rsi),%rdi
mov    $0x1,%r9d
xor    %r15d,%r15d
mov    %rsi,%r11
mov    %rdi,-0x8(%rsp)
jmp    153d <mp_covar+0x43d>
nop
add    $0x1,%r15
add    $0x1,%r9
mov    -0x40(%rsp),%rax
add    %rax,%rdi
add    %rax,%r11
cmp    (%rsp),%r15
je     18ee <mp_covar+0x7ee>
mov    %r9,0x8(%rsp)
mov    %r15,%rax
imul   %r10,%rax
mov    %rax,-0x78(%rsp)
test   %r15,%r15
je     17c0 <mp_covar+0x6c0>
mov    -0x78(%rsp),%rax
lea    (%rsi,%rax,8),%rcx
mov    %rcx,-0x20(%rsp)
add    $0x1,%rax
mov    %rax,-0x18(%rsp)
mov    $0x1,%r14d
mov    $0xfffffffffffffffd,%r8
mov    -0x48(%rsp),%r13
mov    -0x8(%rsp),%rdx
xor    %r12d,%r12d
mov    %r15,-0x58(%rsp)
jmp    15b6 <mp_covar+0x4b6>
nop
add    $0x1,%r14
mov    -0x40(%rsp),%rax
add    %rax,%rdx
add    $0x1,%r8
add    %rax,%r13
mov    -0x58(%rsp),%r15
cmp    %r15,%r12
mov    -0x68(%rsp),%rsi
je     17c0 <mp_covar+0x6c0>
mov    %r12,%r9
lea    -0x3(%r12),%r15
add    $0x1,%r12
mov    -0x78(%rsp),%rax
lea    (%r9,%rax,1),%rbp
mov    %r9,%rax
imul   %r10,%rax
mov    %r9,%rbx
movsd  (%rsi,%rbp,8),%xmm0
sub    $0x3,%rbx
jb     1613 <mp_covar+0x513>
mov    -0x18(%rsp),%rcx
lea    (%rcx,%r9,1),%rbp
lea    (%rsi,%rbp,8),%rbp
lea    (%rsi,%rax,8),%rcx
cmp    %rbp,%rcx
jae    16a0 <mp_covar+0x5a0>
mov    -0x28(%rsp),%rcx
imul   %r9,%rcx
add    -0x48(%rsp),%rcx
cmp    -0x20(%rsp),%rcx
jbe    16a0 <mp_covar+0x5a0>
xor    %r15d,%r15d
mov    %r15,%rcx
test   $0x1,%r14b
je     1649 <mp_covar+0x549>
mov    -0x78(%rsp),%rcx
add    %r15,%rcx
add    %r15,%rax
mov    -0x68(%rsp),%rsi
movsd  (%rsi,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rsi,%rax,8),%xmm1
movsd  %xmm1,(%rsi,%rax,8)
mov    %r15,%rcx
or     $0x1,%rcx
cmp    %r15,%r9
je     1590 <mp_covar+0x490>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  (%r11,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  -0x8(%r13,%rcx,8),%xmm1
movsd  %xmm1,-0x8(%r13,%rcx,8)
movsd  0x8(%r11,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  0x0(%r13,%rcx,8),%xmm1
movsd  %xmm1,0x0(%r13,%rcx,8)
add    $0x2,%rcx
cmp    %rcx,%r14
jne    1660 <mp_covar+0x560>
jmp    1590 <mp_covar+0x490>
nop
shr    $0x2,%r15
add    $0x1,%r15
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
mov    %r15,%rcx
cmp    $0x4,%rbx
jae    16c1 <mp_covar+0x5c1>
xor    %r10d,%r10d
jmp    175b <mp_covar+0x65b>
mov    %r8,%rbx
shr    $0x2,%rbx
add    $0x1,%rbx
and    $0xfffffffffffffffe,%rbx
xor    %r10d,%r10d
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movupd -0x30(%rdi,%r10,8),%xmm2
movupd -0x20(%rdi,%r10,8),%xmm3
mulpd  %xmm1,%xmm2
movupd -0x30(%rdx,%r10,8),%xmm4
addpd  %xmm2,%xmm4
movupd -0x20(%rdx,%r10,8),%xmm2
movupd -0x10(%rdx,%r10,8),%xmm5
movupd (%rdx,%r10,8),%xmm6
mulpd  %xmm1,%xmm3
addpd  %xmm2,%xmm3
movupd %xmm4,-0x30(%rdx,%r10,8)
movupd %xmm3,-0x20(%rdx,%r10,8)
movupd -0x10(%rdi,%r10,8),%xmm2
movupd (%rdi,%r10,8),%xmm3
mulpd  %xmm1,%xmm2
addpd  %xmm5,%xmm2
mulpd  %xmm1,%xmm3
addpd  %xmm6,%xmm3
movupd %xmm2,-0x10(%rdx,%r10,8)
movupd %xmm3,(%rdx,%r10,8)
add    $0x8,%r10
add    $0xfffffffffffffffe,%rbx
jne    16e0 <mp_covar+0x5e0>
mov    %r12,%r15
and    $0xfffffffffffffffc,%r15
test   $0x1,%cl
je     17ac <mp_covar+0x6ac>
mov    -0x78(%rsp),%rcx
add    %r10,%rcx
mov    -0x68(%rsp),%rsi
movupd (%rsi,%rcx,8),%xmm2
movupd 0x10(%rsi,%rcx,8),%xmm3
add    %rax,%r10
movupd (%rsi,%r10,8),%xmm4
movupd 0x10(%rsi,%r10,8),%xmm5
mulpd  %xmm1,%xmm2
addpd  %xmm4,%xmm2
mulpd  %xmm1,%xmm3
addpd  %xmm5,%xmm3
movupd %xmm2,(%rsi,%r10,8)
movupd %xmm3,0x10(%rsi,%r10,8)
cmp    %r15,%r12
mov    -0x70(%rsp),%r10
je     1590 <mp_covar+0x490>
jmp    1616 <mp_covar+0x516>
nop
mov    -0x78(%rsp),%rax
add    %r15,%rax
movsd  (%rsi,%rax,8),%xmm0
mov    0x8(%rsp),%r9
cmp    $0x4,%r9
jae    17f0 <mp_covar+0x6f0>
xor    %eax,%eax
mov    -0x60(%rsp),%rcx
mov    -0x7c(%rsp),%edx
jmp    18d0 <mp_covar+0x7d0>
nopl   0x0(%rax,%rax,1)
mov    %r9,%rax
movabs $0x7ffffffffffffffc,%rcx
ff ff 7f
and    %rcx,%rax
lea    -0x4(%rax),%rdx
mov    %rdx,%r8
shr    $0x2,%r8
add    $0x1,%r8
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
test   %rdx,%rdx
je     1890 <mp_covar+0x790>
mov    %r9,%rbp
and    $0xfffffffffffffffc,%rbp
add    $0xfffffffffffffffc,%rbp
shr    $0x2,%rbp
add    $0x1,%rbp
and    $0xfffffffffffffffe,%rbp
xor    %ebx,%ebx
mov    -0x60(%rsp),%rcx
mov    -0x7c(%rsp),%edx
xchg   %ax,%ax
movupd -0x30(%rdi,%rbx,8),%xmm2
movupd -0x20(%rdi,%rbx,8),%xmm3
movupd -0x10(%rdi,%rbx,8),%xmm4
movupd (%rdi,%rbx,8),%xmm5
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,-0x30(%rdi,%rbx,8)
movupd %xmm3,-0x20(%rdi,%rbx,8)
mulpd  %xmm1,%xmm4
mulpd  %xmm1,%xmm5
movupd %xmm4,-0x10(%rdi,%rbx,8)
movupd %xmm5,(%rdi,%rbx,8)
add    $0x8,%rbx
add    $0xfffffffffffffffe,%rbp
jne    1840 <mp_covar+0x740>
test   $0x1,%r8b
jne    18a1 <mp_covar+0x7a1>
jmp    18e4 <mp_covar+0x7e4>
xor    %ebx,%ebx
mov    -0x60(%rsp),%rcx
mov    -0x7c(%rsp),%edx
test   $0x1,%r8b
je     18e4 <mp_covar+0x7e4>
add    -0x78(%rsp),%rbx
movupd (%rsi,%rbx,8),%xmm2
movupd 0x10(%rsi,%rbx,8),%xmm3
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm3
movupd %xmm2,(%rsi,%rbx,8)
movupd %xmm3,0x10(%rsi,%rbx,8)
jmp    18e4 <mp_covar+0x7e4>
cs nopw 0x0(%rax,%rax,1)
movsd  (%r11,%rax,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%r11,%rax,8)
add    $0x1,%rax
cmp    %rax,%r9
jne    18d0 <mp_covar+0x7d0>
jmp    1520 <mp_covar+0x420>
cmpl   $0x0,-0x38(%rsp)
jle    1a64 <mp_covar+0x964>
movslq -0x50(%rsp),%r11
mov    -0x70(%rsp),%rax
lea    0x0(,%rax,8),%r8
mov    $0x1,%r14d
xor    %r15d,%r15d
mov    %rsi,%r12
jmp    1945 <mp_covar+0x845>
nopl   0x0(%rax)
add    %r15,%r9
movsd  (%rsi,%r9,8),%xmm0
mov    -0x10(%rsp),%rax
movsd  %xmm0,(%rax,%rbp,8)
add    $0x1,%r15
add    $0x1,%r14
add    %r8,%r12
cmp    -0x30(%rsp),%r15
je     19a7 <mp_covar+0x8a7>
movslq (%rcx,%r15,4),%rbp
mov    %r15,%r9
mov    -0x70(%rsp),%r10
imul   %r10,%r9
imul   %ebp,%r10d
xor    %ebx,%ebx
jmp    1969 <mp_covar+0x869>
nopl   (%rax)
add    $0x1,%rbx
cmp    %rbx,%r14
je     1920 <mp_covar+0x820>
cmp    %r11,%r15
jle    1976 <mp_covar+0x876>
movq   $0x0,(%r12,%rbx,8)
mov    (%rcx,%rbx,4),%eax
cmp    %ebp,%eax
jle    198f <mp_covar+0x88f>
movsd  (%r12,%rbx,8),%xmm0
lea    (%rax,%r10,1),%edi
movslq %edi,%rdi
movsd  %xmm0,(%rsi,%rdi,8)
cmp    %ebp,%eax
jge    1960 <mp_covar+0x860>
movsd  (%r12,%rbx,8),%xmm0
imul   %edx,%eax
add    %ebp,%eax
cltq
movsd  %xmm0,(%rsi,%rax,8)
jmp    1960 <mp_covar+0x860>
cmpl   $0x0,-0x38(%rsp)
jle    1a64 <mp_covar+0x964>
mov    -0x70(%rsp),%rcx
shl    $0x4,%rcx
lea    0x8(%rsi),%rdx
xor    %ebx,%ebx
mov    %rsi,%r9
jmp    19f5 <mp_covar+0x8f5>
cs nopw 0x0(%rax,%rax,1)
mov    -0x10(%rsp),%rax
movsd  (%rax,%rbx,8),%xmm0
add    %rbx,%r10
movsd  %xmm0,(%rsi,%r10,8)
add    $0x1,%rbx
add    $0x8,%r9
add    %r8,%rdx
cmp    -0x30(%rsp),%rbx
je     1a64 <mp_covar+0x964>
mov    %rbx,%r10
imul   -0x70(%rsp),%r10
test   %rbx,%rbx
je     19d0 <mp_covar+0x8d0>
cmp    $0x1,%rbx
jne    1a10 <mp_covar+0x910>
xor    %eax,%eax
jmp    1a41 <mp_covar+0x941>
nopl   (%rax)
mov    %rbx,%rbp
and    $0xfffffffffffffffe,%rbp
mov    %r9,%rdi
xor    %eax,%eax
nopl   0x0(%rax)
movsd  (%rdi),%xmm0
movsd  %xmm0,-0x8(%rdx,%rax,8)
movsd  (%rdi,%r8,1),%xmm0
movsd  %xmm0,(%rdx,%rax,8)
add    $0x2,%rax
add    %rcx,%rdi
cmp    %rax,%rbp
jne    1a20 <mp_covar+0x920>
test   $0x1,%bl
je     19d0 <mp_covar+0x8d0>
mov    %rax,%rdi
imul   -0x70(%rsp),%rdi
add    %rbx,%rdi
movsd  (%rsi,%rdi,8),%xmm0
add    %r10,%rax
movsd  %xmm0,(%rsi,%rax,8)
jmp    19d0 <mp_covar+0x8d0>
xor    %eax,%eax
add    $0x10,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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