void apply_householder(float **mat, float *rhs, unsigned int rows,
                       unsigned int cols, float *house, unsigned int row_pos,
                       unsigned int *p) {
  int i, j, k, n;
  float sum;
  float **hhmat;
  float **mat_cpy;
  float *rhs_cpy;
  n = rows - row_pos;
  hhmat = malloc(sizeof(float *) * n);
  for (i = 0; i < n; i++)
    hhmat[i] = malloc(sizeof(float) * n);
  mat_cpy = malloc(sizeof(float *) * rows);
  for (i = 0; i < rows; i++)
    mat_cpy[i] = malloc(sizeof(float) * cols);
  rhs_cpy = malloc(sizeof(float) * rows);
  for (i = 0; i < rows; i++)
    for (j = 0; j < cols; j++)
      mat_cpy[i][j] = mat[i][j];
  for (i = 0; i < rows; i++)
    rhs_cpy[i] = rhs[i];
  for (j = 0; j < n; j++)
    for (i = 0; i < n; i++)
      if (i != j)
        hhmat[i][j] = -2.0 * house[j] * house[i];
      else
        hhmat[i][j] = 1.0 - 2.0 * house[j] * house[i];
  for (k = 0; k < cols; k++)
    for (j = 0; j < n; j++) {
      sum = 0.0;
      for (i = 0; i < n; i++)
        sum += hhmat[j][i] * mat_cpy[i + row_pos][p[k]];
      mat[j + row_pos][p[k]] = sum;
    }
  for (j = 0; j < n; j++) {
    sum = 0.0;
    for (i = 0; i < n; i++)
      sum += hhmat[i][j] * rhs_cpy[i + row_pos];
    rhs[j + row_pos] = sum;
  }
  for (i = 0; i < (rows - row_pos); i++)
    free(hhmat[i]);
  for (i = 0; i < rows; i++)
    free(mat_cpy[i]);
  free(hhmat);
  free(mat_cpy);
  free(rhs_cpy);
}


x86 assembly (clang, -O0):
<apply_householder>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9d,-0x24(%rbp)
mov    -0x14(%rbp),%eax
sub    -0x24(%rbp),%eax
mov    %eax,-0x34(%rbp)
movslq -0x34(%rbp),%rdi
shl    $0x3,%rdi
call   1040 <malloc@plt>
mov    %rax,-0x40(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    1199 <apply_householder+0x79>
movslq -0x34(%rbp),%rdi
shl    $0x2,%rdi
call   1040 <malloc@plt>
mov    %rax,%rdx
mov    -0x40(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1163 <apply_householder+0x43>
mov    -0x14(%rbp),%eax
mov    %eax,%edi
shl    $0x3,%rdi
call   1040 <malloc@plt>
mov    %rax,-0x48(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x14(%rbp),%eax
jae    11e9 <apply_householder+0xc9>
mov    -0x18(%rbp),%eax
mov    %eax,%edi
shl    $0x2,%rdi
call   1040 <malloc@plt>
mov    %rax,%rdx
mov    -0x48(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    11b2 <apply_householder+0x92>
mov    -0x14(%rbp),%eax
mov    %eax,%edi
shl    $0x2,%rdi
call   1040 <malloc@plt>
mov    %rax,-0x50(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x14(%rbp),%eax
jae    126c <apply_householder+0x14c>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x18(%rbp),%eax
jae    1259 <apply_householder+0x139>
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x48(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1215 <apply_householder+0xf5>
jmp    125e <apply_householder+0x13e>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1202 <apply_householder+0xe2>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x14(%rbp),%eax
jae    12a7 <apply_householder+0x187>
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x50(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1273 <apply_householder+0x153>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    13bd <apply_householder+0x29d>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    13aa <apply_householder+0x28a>
mov    -0x28(%rbp),%eax
cmp    -0x2c(%rbp),%eax
je     1329 <apply_householder+0x209>
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xd1e(%rip),%xmm0
mulsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x40(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    1397 <apply_householder+0x277>
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xcc6(%rip),%xmm0
mulsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm1,%xmm0
movsd  0xc86(%rip),%xmm1
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x40(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    139c <apply_householder+0x27c>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    12c1 <apply_householder+0x1a1>
jmp    13af <apply_householder+0x28f>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    12ae <apply_householder+0x18e>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x18(%rbp),%eax
jae    1499 <apply_householder+0x379>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    1486 <apply_householder+0x366>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x38(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    1453 <apply_householder+0x333>
mov    -0x40(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x48(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    -0x24(%rbp),%ecx
mov    %ecx,%ecx
mov    (%rax,%rcx,8),%rax
mov    0x10(%rbp),%rcx
movslq -0x30(%rbp),%rdx
mov    (%rcx,%rdx,4),%ecx
movss  (%rax,%rcx,4),%xmm2
movss  -0x38(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x38(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    13f2 <apply_householder+0x2d2>
movss  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    -0x24(%rbp),%ecx
mov    %ecx,%ecx
mov    (%rax,%rcx,8),%rax
mov    0x10(%rbp),%rcx
movslq -0x30(%rbp),%rdx
mov    (%rcx,%rdx,4),%ecx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    13d7 <apply_householder+0x2b7>
jmp    148b <apply_householder+0x36b>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    13c4 <apply_householder+0x2a4>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    1531 <apply_householder+0x411>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x38(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    150d <apply_householder+0x3ed>
mov    -0x40(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x50(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    -0x24(%rbp),%ecx
mov    %ecx,%ecx
movss  (%rax,%rcx,4),%xmm2
movss  -0x38(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x38(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    14bb <apply_householder+0x39b>
movss  -0x38(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    -0x24(%rbp),%ecx
mov    %ecx,%ecx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    14a0 <apply_householder+0x380>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    -0x24(%rbp),%ecx
cmp    %ecx,%eax
jae    1568 <apply_householder+0x448>
mov    -0x40(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rdi
call   1030 <free@plt>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1538 <apply_householder+0x418>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x14(%rbp),%eax
jae    159a <apply_householder+0x47a>
mov    -0x48(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,8),%rdi
call   1030 <free@plt>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    156f <apply_householder+0x44f>
mov    -0x40(%rbp),%rdi
call   1030 <free@plt>
mov    -0x48(%rbp),%rdi
call   1030 <free@plt>
mov    -0x50(%rbp),%rdi
call   1030 <free@plt>
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<apply_householder>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
mov    %r8,%r12
mov    %ecx,0x14(%rsp)
mov    %rsi,0x28(%rsp)
mov    %rdi,0x38(%rsp)
mov    0x90(%rsp),%rax
mov    %rax,0x50(%rsp)
mov    %edx,0x10(%rsp)
mov    %r9,0x18(%rsp)
sub    %r9d,%edx
mov    %edx,0xc(%rsp)
movslq %edx,%r15
lea    0x0(,%r15,8),%rdi
call   1050 <malloc@plt>
mov    %rax,%r14
test   %r15d,%r15d
jle    11a5 <apply_householder+0x75>
shl    $0x2,%r15
mov    0xc(%rsp),%ebx
xor    %ebp,%ebp
xchg   %ax,%ax
mov    %r15,%rdi
call   1050 <malloc@plt>
mov    %rax,(%r14,%rbp,8)
add    $0x1,%rbp
cmp    %rbp,%rbx
jne    1190 <apply_householder+0x60>
mov    0x10(%rsp),%ebx
mov    %ebx,%ebp
lea    0x0(,%rbp,8),%rdi
call   1050 <malloc@plt>
mov    %rax,%r15
mov    0x14(%rsp),%eax
test   %ebx,%ebx
mov    %rax,0x20(%rsp)
mov    %rbp,0x30(%rsp)
je     141c <apply_householder+0x2ec>
lea    0x0(,%rax,4),%r13
xor    %ebx,%ebx
nopl   0x0(%rax,%rax,1)
mov    %r13,%rdi
call   1050 <malloc@plt>
mov    %rax,(%r15,%rbx,8)
add    $0x1,%rbx
cmp    %rbx,%rbp
jne    11e0 <apply_householder+0xb0>
lea    0x0(,%rbp,4),%rbx
mov    %rbx,%rdi
call   1050 <malloc@plt>
mov    %rax,%r13
cmpl   $0x0,0x10(%rsp)
je     142c <apply_householder+0x2fc>
mov    %rbx,0x40(%rsp)
mov    0x20(%rsp),%rdi
mov    %edi,%r11d
and    $0xfffffff8,%r11d
lea    -0x8(%r11),%r10
mov    %r10,0x48(%rsp)
shr    $0x3,%r10
add    $0x1,%r10
mov    %r10d,%r9d
and    $0x3,%r9d
and    $0xfffffffffffffffc,%r10
mov    %edi,%r8d
and    $0x3,%r8d
mov    %r9,%rsi
shl    $0x5,%rsi
xor    %ecx,%ecx
jmp    126f <apply_householder+0x13f>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%rcx
cmp    0x30(%rsp),%rcx
je     143e <apply_householder+0x30e>
cmpl   $0x0,0x14(%rsp)
je     1260 <apply_householder+0x130>
mov    0x38(%rsp),%rax
mov    (%rax,%rcx,8),%rdx
mov    (%r15,%rcx,8),%rax
cmpl   $0x8,0x14(%rsp)
jae    12a0 <apply_householder+0x170>
xor    %ebx,%ebx
jmp    1390 <apply_householder+0x260>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %rdi,%rbp
lea    (%rdx,%rdi,4),%rdi
cmp    %rdi,%rax
jae    12bf <apply_householder+0x18f>
lea    (%rax,%rbp,4),%rdi
cmp    %rdi,%rdx
jae    12bf <apply_householder+0x18f>
xor    %ebx,%ebx
mov    %rbp,%rdi
jmp    1390 <apply_householder+0x260>
cmpq   $0x18,0x48(%rsp)
jae    12cb <apply_householder+0x19b>
xor    %ebx,%ebx
jmp    132c <apply_householder+0x1fc>
mov    %r10,%rbp
xor    %ebx,%ebx
movups (%rdx,%rbx,4),%xmm0
movups 0x10(%rdx,%rbx,4),%xmm1
movups %xmm0,(%rax,%rbx,4)
movups %xmm1,0x10(%rax,%rbx,4)
movups 0x20(%rdx,%rbx,4),%xmm0
movups 0x30(%rdx,%rbx,4),%xmm1
movups %xmm0,0x20(%rax,%rbx,4)
movups %xmm1,0x30(%rax,%rbx,4)
movups 0x40(%rdx,%rbx,4),%xmm0
movups 0x50(%rdx,%rbx,4),%xmm1
movups %xmm0,0x40(%rax,%rbx,4)
movups %xmm1,0x50(%rax,%rbx,4)
movupd 0x60(%rdx,%rbx,4),%xmm0
movupd 0x70(%rdx,%rbx,4),%xmm1
movupd %xmm0,0x60(%rax,%rbx,4)
movupd %xmm1,0x70(%rax,%rbx,4)
add    $0x20,%rbx
add    $0xfffffffffffffffc,%rbp
jne    12d0 <apply_householder+0x1a0>
test   %r9,%r9
je     1370 <apply_householder+0x240>
lea    (%rax,%rbx,4),%rbp
add    $0x10,%rbp
lea    (%rdx,%rbx,4),%rdi
add    $0x10,%rdi
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movupd -0x10(%rdi,%rbx,1),%xmm0
movupd (%rdi,%rbx,1),%xmm1
movupd %xmm0,-0x10(%rbp,%rbx,1)
movupd %xmm1,0x0(%rbp,%rbx,1)
add    $0x20,%rbx
cmp    %rbx,%rsi
jne    1350 <apply_householder+0x220>
mov    %r11,%rbx
mov    0x20(%rsp),%rdi
cmp    %rdi,%r11
je     1260 <apply_householder+0x130>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %rbx,%rbp
not    %rbp
add    %rdi,%rbp
test   %r8,%r8
je     13c4 <apply_householder+0x294>
mov    %r8,%rdi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movss  (%rdx,%rbx,4),%xmm0
movss  %xmm0,(%rax,%rbx,4)
add    $0x1,%rbx
add    $0xffffffffffffffff,%rdi
jne    13b0 <apply_householder+0x280>
cmp    $0x3,%rbp
mov    0x20(%rsp),%rdi
jb     1260 <apply_householder+0x130>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movss  (%rdx,%rbx,4),%xmm0
movss  %xmm0,(%rax,%rbx,4)
movss  0x4(%rdx,%rbx,4),%xmm0
movss  %xmm0,0x4(%rax,%rbx,4)
movss  0x8(%rdx,%rbx,4),%xmm0
movss  %xmm0,0x8(%rax,%rbx,4)
movss  0xc(%rdx,%rbx,4),%xmm0
movss  %xmm0,0xc(%rax,%rbx,4)
add    $0x4,%rbx
cmp    %rbx,%rdi
jne    13e0 <apply_householder+0x2b0>
jmp    1260 <apply_householder+0x130>
lea    0x0(,%rbp,4),%rdi
call   1050 <malloc@plt>
mov    %rax,%r13
mov    0x18(%rsp),%r11
mov    0xc(%rsp),%ebp
test   %ebp,%ebp
jg     1464 <apply_householder+0x334>
jmp    14dd <apply_householder+0x3ad>
cmpl   $0x0,0x10(%rsp)
mov    0xc(%rsp),%ebp
je     145b <apply_householder+0x32b>
mov    %r13,%rdi
mov    0x28(%rsp),%rsi
mov    0x40(%rsp),%rdx
call   1040 <memcpy@plt>
mov    0x18(%rsp),%r11
test   %ebp,%ebp
jle    14dd <apply_householder+0x3ad>
mov    %ebp,%eax
xor    %ecx,%ecx
movsd  0xb90(%rip),%xmm0
movsd  0xb90(%rip),%xmm1
jmp    1489 <apply_householder+0x359>
nopw   0x0(%rax,%rax,1)
add    $0x1,%rcx
cmp    %rax,%rcx
je     14dd <apply_householder+0x3ad>
xor    %edx,%edx
jmp    14b1 <apply_householder+0x381>
nopl   (%rax)
mulsd  %xmm3,%xmm2
addsd  %xmm1,%xmm2
mov    %rcx,%rsi
cvtsd2ss %xmm2,%xmm2
mov    (%r14,%rsi,8),%rsi
movss  %xmm2,(%rsi,%rcx,4)
add    $0x1,%rdx
cmp    %rdx,%rax
je     1480 <apply_householder+0x350>
movss  (%r12,%rcx,4),%xmm2
xorps  %xmm3,%xmm3
cvtss2sd %xmm2,%xmm3
movaps %xmm3,%xmm2
mulsd  %xmm0,%xmm2
cmp    %rdx,%rcx
je     1490 <apply_householder+0x360>
movss  (%r12,%rdx,4),%xmm3
cvtss2sd %xmm3,%xmm3
mulsd  %xmm3,%xmm2
mov    %rdx,%rsi
jmp    149b <apply_householder+0x36b>
mov    %r11d,%r10d
cmpl   $0x0,0x14(%rsp)
mov    0x30(%rsp),%r12
je     15c8 <apply_householder+0x498>
mov    %ebp,%r9d
mov    %ebp,%edx
and    $0xfffffffe,%edx
lea    0x1(%r11),%esi
xor    %r8d,%r8d
jmp    151f <apply_householder+0x3ef>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
add    $0x1,%r8
cmp    0x20(%rsp),%r8
je     15c8 <apply_householder+0x498>
test   %ebp,%ebp
jle    1510 <apply_householder+0x3e0>
mov    0x50(%rsp),%rax
mov    (%rax,%r8,4),%edi
xor    %ebx,%ebx
jmp    154f <apply_householder+0x41f>
lea    (%rbx,%r11,1),%eax
mov    0x38(%rsp),%rcx
mov    (%rcx,%rax,8),%rax
movss  %xmm0,(%rax,%rdi,4)
add    $0x1,%rbx
cmp    %r9,%rbx
mov    0xc(%rsp),%ebp
je     1510 <apply_householder+0x3e0>
mov    (%r14,%rbx,8),%rcx
cmp    $0x1,%ebp
jne    1560 <apply_householder+0x430>
xorpd  %xmm0,%xmm0
xor    %ebp,%ebp
jmp    15a8 <apply_householder+0x478>
xorpd  %xmm1,%xmm1
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
lea    (%r10,%rbp,1),%eax
mov    (%r15,%rax,8),%rax
movss  (%rcx,%rbp,4),%xmm2
movss  0x4(%rcx,%rbp,4),%xmm0
mulss  (%rax,%rdi,4),%xmm2
lea    (%rsi,%rbp,1),%eax
mov    (%r15,%rax,8),%rax
mulss  (%rax,%rdi,4),%xmm0
addss  %xmm1,%xmm2
addss  %xmm2,%xmm0
add    $0x2,%rbp
movaps %xmm0,%xmm1
cmp    %rbp,%rdx
jne    1570 <apply_householder+0x440>
test   $0x1,%r9b
je     1530 <apply_householder+0x400>
movss  (%rcx,%rbp,4),%xmm1
add    %r11d,%ebp
mov    (%r15,%rbp,8),%rax
mulss  (%rax,%rdi,4),%xmm1
addss  %xmm1,%xmm0
jmp    1530 <apply_householder+0x400>
test   %ebp,%ebp
mov    0x28(%rsp),%r8
jle    167b <apply_householder+0x54b>
mov    %ebp,%ecx
mov    %ebp,%edx
and    $0xfffffffe,%edx
lea    0x1(%r11),%esi
xor    %edi,%edi
jmp    1607 <apply_householder+0x4d7>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
lea    (%rdi,%r11,1),%eax
movss  %xmm0,(%r8,%rax,4)
add    $0x1,%rdi
cmp    %rcx,%rdi
mov    0xc(%rsp),%ebp
je     167b <apply_householder+0x54b>
xorpd  %xmm0,%xmm0
xor    %ebx,%ebx
cmp    $0x1,%ebp
je     165a <apply_householder+0x52a>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    (%r14,%rbx,8),%rax
mov    0x8(%r14,%rbx,8),%rbp
movss  (%rax,%rdi,4),%xmm1
lea    (%r10,%rbx,1),%eax
mulss  0x0(%r13,%rax,4),%xmm1
addss  %xmm0,%xmm1
movss  0x0(%rbp,%rdi,4),%xmm0
lea    (%rsi,%rbx,1),%eax
mulss  0x0(%r13,%rax,4),%xmm0
addss  %xmm1,%xmm0
add    $0x2,%rbx
cmp    %rbx,%rdx
jne    1620 <apply_householder+0x4f0>
test   $0x1,%cl
je     15f0 <apply_householder+0x4c0>
mov    (%r14,%rbx,8),%rax
movss  (%rax,%rdi,4),%xmm1
add    %r11d,%ebx
mulss  0x0(%r13,%rbx,4),%xmm1
addss  %xmm1,%xmm0
jmp    15f0 <apply_householder+0x4c0>
test   %ebp,%ebp
je     16a2 <apply_householder+0x572>
mov    %ebp,%ebx
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    (%r14,%rbp,8),%rdi
call   1030 <free@plt>
add    $0x1,%rbp
cmp    %rbp,%rbx
jne    1690 <apply_householder+0x560>
cmpl   $0x0,0x10(%rsp)
je     16c2 <apply_householder+0x592>
xor    %ebx,%ebx
nopl   0x0(%rax,%rax,1)
mov    (%r15,%rbx,8),%rdi
call   1030 <free@plt>
add    $0x1,%rbx
cmp    %rbx,%r12
jne    16b0 <apply_householder+0x580>
mov    %r14,%rdi
call   1030 <free@plt>
mov    %r15,%rdi
call   1030 <free@plt>
mov    %r13,%rdi
add    $0x58,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

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

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