void wregress_ch(data, weight, a, ndata, ivar) float **data, *a, *weight;
int ndata, ivar;
{
  int p = ivar + 1, j, i, k;
  float **d, *b, sum, *diag;
  d = matrix(1, p, 1, p);
  b = vector(1, p);
  diag = vector(1, p);
  sum = 0;
  for (k = 1; k <= ndata; k++)
    sum += weight[k] * weight[k];
  d[1][1] = sum;
  for (j = 1; j <= ivar; j++) {
    sum = 0;
    for (k = 1; k <= ndata; k++)
      sum += weight[k] * data[k][j];
    d[1][j + 1] = sum;
  }
  for (i = 1; i <= ivar; i++)
    for (j = i; j <= ivar; j++) {
      sum = 0;
      for (k = 1; k <= ndata; k++)
        sum += data[k][j] * data[k][i] * weight[k] * weight[k];
      d[i + 1][j + 1] = sum;
    }
  sum = 0;
  for (k = 1; k <= ndata; k++)
    sum += weight[k] * data[k][p];
  b[1] = sum;
  for (i = 1; i <= ivar; i++) {
    sum = 0;
    for (k = 1; k <= ndata; k++)
      sum += data[k][p] * data[k][i];
    b[i + 1] = sum;
  }
  choldc_e(d, p, diag);
  cholsl(d, p, diag, b, a);
  free_matrix(d, 1, p, 1, p);
  free_vector(b, 1, p);
  free_vector(diag, 1, p);
}


x86 assembly (clang, -O0):
<wregress_ch>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%esi
mov    -0x24(%rbp),%ecx
mov    $0x1,%edx
mov    %edx,%edi
mov    $0x0,%al
call   1050 <matrix@plt>
cltq
mov    %rax,-0x38(%rbp)
mov    -0x24(%rbp),%esi
mov    $0x1,%edi
mov    $0x0,%al
call   1060 <vector@plt>
cltq
mov    %rax,-0x40(%rbp)
mov    -0x24(%rbp),%esi
mov    $0x1,%edi
mov    $0x0,%al
call   1060 <vector@plt>
cltq
mov    %rax,-0x50(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x44(%rbp)
movl   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jg     121d <wregress_ch+0xbd>
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x44(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    11d7 <wregress_ch+0x77>
movss  -0x44(%rbp),%xmm0
mov    -0x38(%rbp),%rax
mov    0x8(%rax),%rax
movss  %xmm0,0x4(%rax)
movl   $0x1,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x20(%rbp),%eax
jg     12c8 <wregress_ch+0x168>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x44(%rbp)
movl   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jg     129f <wregress_ch+0x13f>
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x44(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1251 <wregress_ch+0xf1>
movss  -0x44(%rbp),%xmm0
mov    -0x38(%rbp),%rax
mov    0x8(%rax),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1236 <wregress_ch+0xd6>
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jg     13b1 <wregress_ch+0x251>
mov    -0x2c(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x20(%rbp),%eax
jg     139e <wregress_ch+0x23e>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x44(%rbp)
movl   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jg     136c <wregress_ch+0x20c>
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x44(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    12fc <wregress_ch+0x19c>
movss  -0x44(%rbp),%xmm0
mov    -0x38(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x28(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    12e1 <wregress_ch+0x181>
jmp    13a3 <wregress_ch+0x243>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    12cf <wregress_ch+0x16f>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x44(%rbp)
movl   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jg     140e <wregress_ch+0x2ae>
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x44(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    13c0 <wregress_ch+0x260>
movss  -0x44(%rbp),%xmm0
mov    -0x40(%rbp),%rax
movss  %xmm0,0x4(%rax)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jg     14b9 <wregress_ch+0x359>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x44(%rbp)
movl   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jg     1494 <wregress_ch+0x334>
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x44(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    143e <wregress_ch+0x2de>
movss  -0x44(%rbp),%xmm0
mov    -0x40(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1423 <wregress_ch+0x2c3>
mov    -0x38(%rbp),%rdi
mov    -0x24(%rbp),%esi
mov    -0x50(%rbp),%rdx
mov    $0x0,%al
call   1070 <choldc_e@plt>
mov    -0x38(%rbp),%rdi
mov    -0x24(%rbp),%esi
mov    -0x50(%rbp),%rdx
mov    -0x40(%rbp),%rcx
mov    -0x18(%rbp),%r8
mov    $0x0,%al
call   1080 <cholsl@plt>
mov    -0x38(%rbp),%rdi
mov    -0x24(%rbp),%edx
mov    -0x24(%rbp),%r8d
mov    $0x1,%ecx
mov    %ecx,%esi
mov    $0x0,%al
call   1030 <free_matrix@plt>
mov    -0x40(%rbp),%rdi
mov    -0x24(%rbp),%edx
mov    $0x1,%esi
mov    $0x0,%al
call   1040 <free_vector@plt>
mov    -0x50(%rbp),%rdi
mov    -0x24(%rbp),%edx
mov    $0x1,%esi
mov    $0x0,%al
call   1040 <free_vector@plt>
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<wregress_ch>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %r8d,%r12d
mov    %ecx,%r15d
mov    %rdx,0x30(%rsp)
mov    %rsi,%r14
mov    %rdi,%rbx
lea    0x1(%r12),%ebp
mov    $0x1,%edi
mov    %ebp,%esi
mov    $0x1,%edx
mov    %ebp,%ecx
xor    %eax,%eax
call   1050 <matrix@plt>
cltq
mov    %rax,0x10(%rsp)
mov    $0x1,%edi
mov    %ebp,%esi
xor    %eax,%eax
call   1060 <vector@plt>
mov    %eax,0x18(%rsp)
mov    $0x1,%edi
mov    %ebp,0xc(%rsp)
mov    %ebp,%esi
xor    %eax,%eax
call   1060 <vector@plt>
mov    %eax,0x1c(%rsp)
mov    %r15d,%eax
test   %r15d,%r15d
mov    %rax,0x28(%rsp)
jle    11f6 <wregress_ch+0x96>
lea    -0x1(%rax),%rdx
mov    %eax,%ecx
and    $0x3,%ecx
cmp    $0x3,%rdx
jae    11fe <wregress_ch+0x9e>
xorps  %xmm0,%xmm0
mov    $0x1,%edx
test   %rcx,%rcx
jne    125e <wregress_ch+0xfe>
jmp    1286 <wregress_ch+0x126>
xorps  %xmm0,%xmm0
jmp    1286 <wregress_ch+0x126>
mov    %eax,%esi
and    $0xfffffffc,%esi
xorps  %xmm0,%xmm0
xor    %edx,%edx
nopl   0x0(%rax,%rax,1)
movss  0x4(%r14,%rdx,4),%xmm1
movss  0x8(%r14,%rdx,4),%xmm2
mulss  %xmm1,%xmm1
addss  %xmm0,%xmm1
mulss  %xmm2,%xmm2
addss  %xmm1,%xmm2
movss  0xc(%r14,%rdx,4),%xmm1
mulss  %xmm1,%xmm1
addss  %xmm2,%xmm1
movss  0x10(%r14,%rdx,4),%xmm0
mulss  %xmm0,%xmm0
addss  %xmm1,%xmm0
add    $0x4,%rdx
cmp    %rdx,%rsi
jne    1210 <wregress_ch+0xb0>
add    $0x1,%rdx
test   %rcx,%rcx
je     1286 <wregress_ch+0x126>
lea    (%r14,%rdx,4),%rdx
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movss  (%rdx,%rsi,4),%xmm1
mulss  %xmm1,%xmm1
addss  %xmm1,%xmm0
add    $0x1,%rsi
cmp    %rsi,%rcx
jne    1270 <wregress_ch+0x110>
mov    0x10(%rsp),%rax
mov    0x8(%rax),%r9
movss  %xmm0,0x4(%r9)
mov    0xc(%rsp),%r10d
mov    %r12,0x20(%rsp)
test   %r12d,%r12d
jle    14a3 <wregress_ch+0x343>
lea    0x1(%r15),%r13d
lea    -0x1(%r13),%r8
add    $0xfffffffffffffffe,%r13
mov    %r8d,%edi
and    $0x3,%edi
mov    %r8,%r11
and    $0xfffffffffffffffc,%r11
mov    $0x1,%ecx
jmp    12e4 <wregress_ch+0x184>
nopl   0x0(%rax,%rax,1)
movss  %xmm0,0x4(%r9,%rcx,4)
add    $0x1,%rcx
cmp    %r10,%rcx
je     13b2 <wregress_ch+0x252>
xorps  %xmm0,%xmm0
test   %r15d,%r15d
jle    12d0 <wregress_ch+0x170>
cmp    $0x3,%r13
jae    1300 <wregress_ch+0x1a0>
mov    $0x1,%esi
test   %rdi,%rdi
jne    137b <wregress_ch+0x21b>
jmp    12d0 <wregress_ch+0x170>
xchg   %ax,%ax
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movss  0x4(%r14,%rsi,4),%xmm1
movss  0x8(%r14,%rsi,4),%xmm2
mov    0x8(%rbx,%rsi,8),%rbp
mov    0x10(%rbx,%rsi,8),%rdx
mulss  0x0(%rbp,%rcx,4),%xmm1
addss  %xmm0,%xmm1
mulss  (%rdx,%rcx,4),%xmm2
addss  %xmm1,%xmm2
movss  0xc(%r14,%rsi,4),%xmm1
mov    0x18(%rbx,%rsi,8),%rdx
mulss  (%rdx,%rcx,4),%xmm1
movss  0x10(%r14,%rsi,4),%xmm0
mov    0x20(%rbx,%rsi,8),%rdx
mulss  (%rdx,%rcx,4),%xmm0
addss  %xmm2,%xmm1
addss  %xmm1,%xmm0
add    $0x4,%rsi
cmp    %rsi,%r11
jne    1310 <wregress_ch+0x1b0>
add    $0x1,%rsi
test   %rdi,%rdi
je     12d0 <wregress_ch+0x170>
lea    (%rbx,%rsi,8),%rbp
lea    (%r14,%rsi,4),%rsi
xor    %r12d,%r12d
cs nopw 0x0(%rax,%rax,1)
movss  (%rsi,%r12,4),%xmm1
mov    0x0(%rbp,%r12,8),%rdx
mulss  (%rdx,%rcx,4),%xmm1
addss  %xmm1,%xmm0
add    $0x1,%r12
cmp    %r12,%rdi
jne    1390 <wregress_ch+0x230>
jmp    12d0 <wregress_ch+0x170>
cmpl   $0x0,0x20(%rsp)
jle    14a3 <wregress_ch+0x343>
mov    %r8,%rdx
and    $0xfffffffffffffffe,%rdx
neg    %rdx
mov    $0x1,%ebp
jmp    13dc <wregress_ch+0x27c>
xchg   %ax,%ax
mov    %r9,%rbp
cmp    %r10,%r9
je     14a3 <wregress_ch+0x343>
lea    0x1(%rbp),%r9
mov    0x10(%rsp),%rax
mov    0x8(%rax,%rbp,8),%r11
mov    %rbp,%rcx
jmp    1400 <wregress_ch+0x2a0>
nop
movss  %xmm0,0x4(%r11,%rcx,4)
add    $0x1,%rcx
cmp    %r10,%rcx
je     13d0 <wregress_ch+0x270>
xorps  %xmm0,%xmm0
test   %r15d,%r15d
jle    13f0 <wregress_ch+0x290>
mov    $0x1,%esi
test   %r13,%r13
je     1474 <wregress_ch+0x314>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    (%rbx,%rsi,8),%rax
mov    0x8(%rbx,%rsi,8),%rdi
movss  (%rax,%rcx,4),%xmm1
mulss  (%rax,%rbp,4),%xmm1
movss  (%r14,%rsi,4),%xmm2
movss  0x4(%r14,%rsi,4),%xmm3
mulss  %xmm2,%xmm1
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
movss  (%rdi,%rcx,4),%xmm0
mulss  (%rdi,%rbp,4),%xmm0
mulss  %xmm3,%xmm0
mulss  %xmm3,%xmm0
addss  %xmm1,%xmm0
lea    (%rdx,%rsi,1),%rax
add    $0x2,%rax
add    $0x2,%rsi
cmp    $0x1,%rax
jne    1420 <wregress_ch+0x2c0>
test   $0x1,%r8b
je     13f0 <wregress_ch+0x290>
mov    (%rbx,%rsi,8),%rax
movss  (%rax,%rcx,4),%xmm1
mulss  (%rax,%rbp,4),%xmm1
movss  (%r14,%rsi,4),%xmm2
mulss  %xmm2,%xmm1
mulss  %xmm2,%xmm1
addss  %xmm1,%xmm0
jmp    13f0 <wregress_ch+0x290>
movslq 0x18(%rsp),%r13
movslq 0xc(%rsp),%rsi
test   %r15d,%r15d
jle    14de <wregress_ch+0x37e>
mov    0x28(%rsp),%rcx
lea    -0x1(%rcx),%rdx
mov    %ecx,%r8d
and    $0x3,%r8d
cmp    $0x3,%rdx
jae    14e6 <wregress_ch+0x386>
xorps  %xmm0,%xmm0
mov    $0x1,%edx
test   %r8,%r8
jne    1556 <wregress_ch+0x3f6>
jmp    157b <wregress_ch+0x41b>
xorps  %xmm0,%xmm0
jmp    157b <wregress_ch+0x41b>
and    $0xfffffffc,%ecx
xorps  %xmm0,%xmm0
xor    %edx,%edx
xchg   %ax,%ax
movss  0x4(%r14,%rdx,4),%xmm1
movss  0x8(%r14,%rdx,4),%xmm2
mov    0x8(%rbx,%rdx,8),%rax
mov    0x10(%rbx,%rdx,8),%rdi
mulss  (%rax,%rsi,4),%xmm1
addss  %xmm0,%xmm1
mulss  (%rdi,%rsi,4),%xmm2
addss  %xmm1,%xmm2
movss  0xc(%r14,%rdx,4),%xmm1
mov    0x18(%rbx,%rdx,8),%rax
mulss  (%rax,%rsi,4),%xmm1
movss  0x10(%r14,%rdx,4),%xmm0
mov    0x20(%rbx,%rdx,8),%rax
mulss  (%rax,%rsi,4),%xmm0
addss  %xmm2,%xmm1
addss  %xmm1,%xmm0
add    $0x4,%rdx
cmp    %rdx,%rcx
jne    14f0 <wregress_ch+0x390>
add    $0x1,%rdx
test   %r8,%r8
je     157b <wregress_ch+0x41b>
lea    (%rbx,%rdx,8),%rdi
lea    (%r14,%rdx,4),%rdx
xor    %ebp,%ebp
movss  (%rdx,%rbp,4),%xmm1
mov    (%rdi,%rbp,8),%rax
mulss  (%rax,%rsi,4),%xmm1
addss  %xmm1,%xmm0
add    $0x1,%rbp
cmp    %rbp,%r8
jne    1560 <wregress_ch+0x400>
movslq 0x1c(%rsp),%r14
movss  %xmm0,0x4(%r13)
cmpl   $0x0,0x20(%rsp)
jle    1681 <wregress_ch+0x521>
lea    0x1(%r15),%r8d
lea    -0x1(%r8),%r9
add    $0xfffffffffffffffe,%r8
mov    %r9d,%edi
and    $0x3,%edi
and    $0xfffffffffffffffc,%r9
mov    $0x1,%ecx
jmp    15c4 <wregress_ch+0x464>
xchg   %ax,%ax
movss  %xmm0,0x4(%r13,%rcx,4)
add    $0x1,%rcx
cmp    %r10,%rcx
je     1681 <wregress_ch+0x521>
xorps  %xmm0,%xmm0
test   %r15d,%r15d
jle    15b0 <wregress_ch+0x450>
cmp    $0x3,%r8
jae    15e0 <wregress_ch+0x480>
mov    $0x1,%ebp
test   %rdi,%rdi
jne    1652 <wregress_ch+0x4f2>
jmp    15b0 <wregress_ch+0x450>
xchg   %ax,%ax
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    0x8(%rbx,%rbp,8),%rax
mov    0x10(%rbx,%rbp,8),%rdx
movss  (%rax,%rsi,4),%xmm1
mulss  (%rax,%rcx,4),%xmm1
addss  %xmm0,%xmm1
movss  (%rdx,%rsi,4),%xmm2
mulss  (%rdx,%rcx,4),%xmm2
addss  %xmm1,%xmm2
mov    0x18(%rbx,%rbp,8),%rax
movss  (%rax,%rsi,4),%xmm1
mulss  (%rax,%rcx,4),%xmm1
mov    0x20(%rbx,%rbp,8),%rax
movss  (%rax,%rsi,4),%xmm0
mulss  (%rax,%rcx,4),%xmm0
addss  %xmm2,%xmm1
addss  %xmm1,%xmm0
add    $0x4,%rbp
cmp    %rbp,%r9
jne    15f0 <wregress_ch+0x490>
add    $0x1,%rbp
test   %rdi,%rdi
je     15b0 <wregress_ch+0x450>
lea    (%rbx,%rbp,8),%rbp
xor    %eax,%eax
nopl   0x0(%rax,%rax,1)
mov    0x0(%rbp,%rax,8),%rdx
movss  (%rdx,%rsi,4),%xmm1
mulss  (%rdx,%rcx,4),%xmm1
addss  %xmm1,%xmm0
add    $0x1,%rax
cmp    %rax,%rdi
jne    1660 <wregress_ch+0x500>
jmp    15b0 <wregress_ch+0x450>
mov    0x10(%rsp),%rbp
mov    %rbp,%rdi
mov    0xc(%rsp),%ebx
mov    %ebx,%esi
mov    %r14,%rdx
xor    %eax,%eax
call   1070 <choldc_e@plt>
mov    %rbp,%rdi
mov    %ebx,%esi
mov    %r14,%rdx
mov    %r13,%rcx
mov    0x30(%rsp),%r8
xor    %eax,%eax
call   1080 <cholsl@plt>
mov    %rbp,%rdi
mov    $0x1,%esi
mov    %ebx,%edx
mov    $0x1,%ecx
mov    %ebx,%r8d
xor    %eax,%eax
call   1030 <free_matrix@plt>
mov    %r13,%rdi
mov    $0x1,%esi
mov    %ebx,%edx
xor    %eax,%eax
call   1040 <free_vector@plt>
mov    %r14,%rdi
mov    $0x1,%esi
mov    %ebx,%edx
xor    %eax,%eax
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1040 <free_vector@plt>

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

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