#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);
double sqrt(double);
void Householder(int mode, double *v, int p, int l, int m, double *u_p,
                 double *C[], int nC) {
  double vmax, vmax_inv, vtmp, s, b, tj;
  double *cj;
  int i, j;
  if (0 > p || p >= l || l > m)
    return;
  vmax = fabs(v[p]);
  if (mode == 1) {
    for (i = l; i < m; i++)
      vmax = ((vmax) < (fabs(v[i])) ? (fabs(v[i])) : (vmax));
    if (vmax == 0)
      return;
    vmax_inv = 1 / vmax;
    vtmp = v[p] * vmax_inv;
    s = vtmp * vtmp;
    for (i = l; i < m; i++) {
      vtmp = v[i] * vmax_inv;
      s += vtmp * vtmp;
    }
    s = vmax * sqrt(s);
    if (v[p] > 0)
      s = -s;
    *u_p = v[p] - s;
    v[p] = s;
  } else {
    if (vmax == 0)
      return;
  }
  if (nC <= 0 || C == NULL)
    return;
  b = v[p] * (*u_p);
  if (b >= 0)
    return;
  for (j = 0; j < nC; j++) {
    cj = C[j];
    tj = cj[p] * (*u_p);
    for (i = l; i < m; i++)
      tj += cj[i] * v[i];
    if (tj != 0) {
      tj /= b;
      cj[p] += tj * (*u_p);
      for (i = l; i < m; i++)
        cj[i] += tj * v[i];
    }
  }
}


x86 assembly (clang, -O0):
<Householder>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8d,-0x1c(%rbp)
mov    %r9,-0x28(%rbp)
xor    %eax,%eax
cmp    -0x14(%rbp),%eax
jg     115a <Householder+0x4a>
mov    -0x14(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    115a <Householder+0x4a>
mov    -0x18(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jle    115f <Householder+0x4f>
jmp    14ed <Householder+0x3dd>
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xe8d(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
cmpl   $0x1,-0x4(%rbp)
jne    133c <Householder+0x22c>
mov    -0x18(%rbp),%eax
mov    %eax,-0x64(%rbp)
mov    -0x64(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1203 <Householder+0xf3>
movsd  -0x30(%rbp),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xe4f(%rip),%xmm2
pand   %xmm2,%xmm0
ucomisd %xmm1,%xmm0
jbe    11e1 <Householder+0xd1>
mov    -0x10(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xe2d(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
jmp    11eb <Householder+0xdb>
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x64(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
jmp    118c <Householder+0x7c>
movsd  -0x30(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1220 <Householder+0x110>
jp     1220 <Householder+0x110>
jmp    14ed <Householder+0x3dd>
movsd  0xde8(%rip),%xmm0
divsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
mulsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x18(%rbp),%eax
mov    %eax,-0x64(%rbp)
mov    -0x64(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    12ab <Householder+0x19b>
mov    -0x10(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x64(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
jmp    125e <Householder+0x14e>
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x48(%rbp),%xmm0
call   1030 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0x78(%rbp),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    130b <Householder+0x1fb>
movsd  -0x48(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x48(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x48(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    135e <Householder+0x24e>
movsd  -0x30(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1359 <Householder+0x249>
jp     1359 <Householder+0x249>
jmp    14ed <Householder+0x3dd>
jmp    135e <Householder+0x24e>
cmpl   $0x0,0x18(%rbp)
jle    1373 <Householder+0x263>
cmpq   $0x0,0x10(%rbp)
jne    1378 <Householder+0x268>
jmp    14ed <Householder+0x3dd>
mov    -0x10(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
mulsd  (%rax),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     13a9 <Householder+0x299>
jmp    14ed <Householder+0x3dd>
movl   $0x0,-0x68(%rbp)
mov    -0x68(%rbp),%eax
cmp    0x18(%rbp),%eax
jge    14ed <Householder+0x3dd>
mov    0x10(%rbp),%rax
movslq -0x68(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x60(%rbp)
mov    -0x60(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
mulsd  (%rax),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x18(%rbp),%eax
mov    %eax,-0x64(%rbp)
mov    -0x64(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1432 <Householder+0x322>
mov    -0x60(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x58(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x64(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
jmp    13ec <Householder+0x2dc>
movsd  -0x58(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    144f <Householder+0x33f>
jp     144f <Householder+0x33f>
jmp    14da <Householder+0x3ca>
movsd  -0x50(%rbp),%xmm1
movsd  -0x58(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x58(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm2
mov    -0x60(%rbp),%rax
movslq -0x14(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%eax
mov    %eax,-0x64(%rbp)
mov    -0x64(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    14d5 <Householder+0x3c5>
movsd  -0x58(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x60(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x64(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
jmp    148f <Householder+0x37f>
jmp    14da <Householder+0x3ca>
jmp    14df <Householder+0x3cf>
mov    -0x68(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x68(%rbp)
jmp    13b0 <Householder+0x2a0>
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<Householder>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x78,%rsp
test   %edx,%edx
js     16ec <Householder+0x5dc>
mov    %ecx,%ebx
cmp    %ecx,%edx
jge    16ec <Householder+0x5dc>
mov    %r8d,%r15d
cmp    %r8d,%ebx
jg     16ec <Householder+0x5dc>
mov    %rsi,%rbp
mov    %edx,%r12d
movsd  (%rsi,%r12,8),%xmm0
cmp    $0x1,%edi
jne    12c7 <Householder+0x1b7>
movapd 0xea7(%rip),%xmm6
andpd  %xmm0,%xmm6
cmp    %r15d,%ebx
jge    12e0 <Householder+0x1d0>
movslq %ebx,%rsi
movslq %r15d,%rax
mov    %rax,%rdx
sub    %rsi,%rdx
mov    %rsi,%rcx
not    %rcx
add    %rax,%rcx
mov    %rdx,%rdi
and    $0x3,%rdi
je     11ac <Householder+0x9c>
movapd 0xe74(%rip),%xmm1
movapd %xmm6,%xmm2
movsd  0x0(%rbp,%rsi,8),%xmm6
andpd  %xmm1,%xmm6
maxsd  %xmm2,%xmm6
add    $0x1,%rsi
movapd %xmm6,%xmm2
add    $0xffffffffffffffff,%rdi
jne    1190 <Householder+0x80>
movapd %xmm6,%xmm1
cmp    $0x3,%rcx
jb     1205 <Householder+0xf5>
movapd 0xe42(%rip),%xmm2
xchg   %ax,%ax
movsd  0x0(%rbp,%rsi,8),%xmm3
movsd  0x8(%rbp,%rsi,8),%xmm4
movsd  0x10(%rbp,%rsi,8),%xmm5
movsd  0x18(%rbp,%rsi,8),%xmm1
andpd  %xmm2,%xmm3
maxsd  %xmm6,%xmm3
andpd  %xmm2,%xmm4
maxsd  %xmm3,%xmm4
andpd  %xmm2,%xmm5
maxsd  %xmm4,%xmm5
andpd  %xmm2,%xmm1
maxsd  %xmm5,%xmm1
add    $0x4,%rsi
movapd %xmm1,%xmm6
cmp    %rsi,%rax
jne    11c0 <Householder+0xb0>
xorpd  %xmm2,%xmm2
ucomisd %xmm2,%xmm1
jne    1215 <Householder+0x105>
jnp    16ec <Householder+0x5dc>
movsd  0xe03(%rip),%xmm2
divsd  %xmm1,%xmm2
mulsd  %xmm2,%xmm0
mulsd  %xmm0,%xmm0
cmp    %r15d,%ebx
jge    12c1 <Householder+0x1b1>
movslq %ebx,%rsi
and    $0x3,%rdx
je     125c <Householder+0x14c>
nopl   0x0(%rax,%rax,1)
movsd  0x0(%rbp,%rsi,8),%xmm3
mulsd  %xmm2,%xmm3
mulsd  %xmm3,%xmm3
addsd  %xmm3,%xmm0
add    $0x1,%rsi
add    $0xffffffffffffffff,%rdx
jne    1240 <Householder+0x130>
cmp    $0x3,%rcx
jb     12c1 <Householder+0x1b1>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  0x0(%rbp,%rsi,8),%xmm3
mulsd  %xmm2,%xmm3
mulsd  %xmm3,%xmm3
addsd  %xmm0,%xmm3
movsd  0x8(%rbp,%rsi,8),%xmm0
mulsd  %xmm2,%xmm0
mulsd  %xmm0,%xmm0
addsd  %xmm3,%xmm0
movsd  0x10(%rbp,%rsi,8),%xmm3
mulsd  %xmm2,%xmm3
mulsd  %xmm3,%xmm3
addsd  %xmm0,%xmm3
movsd  0x18(%rbp,%rsi,8),%xmm0
mulsd  %xmm2,%xmm0
mulsd  %xmm0,%xmm0
addsd  %xmm3,%xmm0
add    $0x4,%rsi
cmp    %rsi,%rax
jne    1270 <Householder+0x160>
movapd %xmm1,%xmm6
jmp    1304 <Householder+0x1f4>
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    136b <Householder+0x25b>
jnp    16ec <Householder+0x5dc>
jmp    136b <Householder+0x25b>
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    12f0 <Householder+0x1e0>
jnp    16ec <Householder+0x5dc>
movsd  0xd28(%rip),%xmm1
divsd  %xmm6,%xmm1
mulsd  %xmm1,%xmm0
mulsd  %xmm0,%xmm0
xorpd  %xmm2,%xmm2
ucomisd %xmm2,%xmm0
jb     1314 <Householder+0x204>
sqrtsd %xmm0,%xmm0
jmp    132f <Householder+0x21f>
mov    %r9,%r14
movapd %xmm6,0x20(%rsp)
call   1030 <sqrt@plt>
xorpd  %xmm2,%xmm2
movapd 0x20(%rsp),%xmm6
mov    %r14,%r9
mulsd  %xmm0,%xmm6
movsd  0x0(%rbp,%r12,8),%xmm0
movapd 0xcce(%rip),%xmm1
xorpd  %xmm6,%xmm1
cmpltsd %xmm0,%xmm2
andpd  %xmm2,%xmm1
andnpd %xmm6,%xmm2
orpd   %xmm1,%xmm2
subsd  %xmm2,%xmm0
movsd  %xmm0,(%r9)
movlpd %xmm2,0x0(%rbp,%r12,8)
movapd %xmm2,%xmm0
mov    0xb0(%rsp),%r10
test   %r10,%r10
je     16ec <Householder+0x5dc>
mov    0xb8(%rsp),%eax
test   %eax,%eax
jle    16ec <Householder+0x5dc>
movsd  (%r9),%xmm3
mulsd  %xmm3,%xmm0
xorpd  %xmm8,%xmm8
ucomisd %xmm8,%xmm0
jae    16ec <Householder+0x5dc>
movslq %ebx,%r14
mov    %eax,%r8d
movslq %r15d,%rsi
lea    0x0(,%r14,8),%rdi
add    %rbp,%rdi
lea    0x0(,%rsi,8),%rax
add    %rbp,%rax
mov    %rax,0x68(%rsp)
mov    %rsi,%rdx
sub    %r14,%rdx
mov    %r14,%r11
not    %r11
add    %rsi,%r11
lea    -0x4(%rdx),%rcx
mov    %rcx,0x60(%rsp)
shr    $0x2,%rcx
add    $0x1,%rcx
mov    %edx,%eax
and    $0x3,%eax
mov    %rax,0x20(%rsp)
mov    %rdx,0x18(%rsp)
mov    %rdx,%rax
mov    %r8,%rdx
and    $0xfffffffffffffffc,%rax
mov    %rax,0x50(%rsp)
add    %r14,%rax
mov    %rax,0x48(%rsp)
mov    %rcx,%rax
mov    %rcx,0x58(%rsp)
and    $0xfffffffffffffffe,%rcx
mov    %rcx,0x40(%rsp)
mov    %rsi,%rax
neg    %rax
mov    %rax,0x70(%rsp)
lea    0x6(%r14),%rax
mov    %rax,0x38(%rsp)
xor    %r13d,%r13d
mov    %r11,0x10(%rsp)
mov    %r8,0x8(%rsp)
mov    (%r10,%r13,8),%r8
movsd  (%r8,%r12,8),%xmm4
movapd %xmm4,%xmm2
mulsd  %xmm3,%xmm2
cmp    %r15d,%ebx
jge    1540 <Householder+0x430>
mov    0x20(%rsp),%rax
mov    %r14,%rcx
test   %rax,%rax
je     148a <Householder+0x37a>
nopl   0x0(%rax,%rax,1)
movsd  (%r8,%rcx,8),%xmm5
mulsd  0x0(%rbp,%rcx,8),%xmm5
addsd  %xmm5,%xmm2
add    $0x1,%rcx
add    $0xffffffffffffffff,%rax
jne    1470 <Householder+0x360>
cmp    $0x3,%r11
jb     14dc <Householder+0x3cc>
movsd  (%r8,%rcx,8),%xmm5
movsd  0x8(%r8,%rcx,8),%xmm6
mulsd  0x0(%rbp,%rcx,8),%xmm5
mulsd  0x8(%rbp,%rcx,8),%xmm6
addsd  %xmm2,%xmm5
movsd  0x10(%r8,%rcx,8),%xmm7
mulsd  0x10(%rbp,%rcx,8),%xmm7
addsd  %xmm5,%xmm6
movsd  0x18(%r8,%rcx,8),%xmm2
mulsd  0x18(%rbp,%rcx,8),%xmm2
addsd  %xmm6,%xmm7
addsd  %xmm7,%xmm2
add    $0x4,%rcx
cmp    %rcx,%rsi
jne    1490 <Householder+0x380>
ucomisd %xmm8,%xmm2
jne    14e5 <Householder+0x3d5>
jnp    155b <Householder+0x44b>
divsd  %xmm0,%xmm2
mulsd  %xmm2,%xmm3
addsd  %xmm4,%xmm3
movsd  %xmm3,(%r8,%r12,8)
cmp    %r15d,%ebx
jge    155b <Householder+0x44b>
mov    %r14,%rax
cmpq   $0x4,0x18(%rsp)
jb     1669 <Householder+0x559>
lea    (%r8,%r14,8),%rax
cmp    0x68(%rsp),%rax
jae    1526 <Householder+0x416>
lea    (%r8,%rsi,8),%rcx
mov    %r14,%rax
cmp    %rcx,%rdi
jb     1669 <Householder+0x559>
movapd %xmm2,%xmm3
unpcklpd %xmm2,%xmm3
cmpq   $0x4,0x60(%rsp)
jae    1572 <Householder+0x462>
mov    %r10,%rax
xor    %ecx,%ecx
jmp    1607 <Householder+0x4f7>
ucomisd %xmm8,%xmm2
jne    1549 <Householder+0x439>
jnp    155b <Householder+0x44b>
divsd  %xmm0,%xmm2
mulsd  %xmm3,%xmm2
addsd  %xmm4,%xmm2
movsd  %xmm2,(%r8,%r12,8)
add    $0x1,%r13
cmp    %rdx,%r13
je     16ec <Householder+0x5dc>
movsd  (%r9),%xmm3
jmp    1440 <Householder+0x330>
mov    %r10,%rax
mov    0x38(%rsp),%rcx
lea    (%r8,%rcx,8),%r10
mov    0x40(%rsp),%r11
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nop
movupd (%rdi,%rcx,8),%xmm4
movupd 0x10(%rdi,%rcx,8),%xmm5
mulpd  %xmm3,%xmm4
movupd -0x30(%r10,%rcx,8),%xmm6
addpd  %xmm4,%xmm6
movupd -0x20(%r10,%rcx,8),%xmm4
movupd -0x10(%r10,%rcx,8),%xmm7
movupd (%r10,%rcx,8),%xmm1
mulpd  %xmm3,%xmm5
addpd  %xmm4,%xmm5
movupd %xmm6,-0x30(%r10,%rcx,8)
movupd %xmm5,-0x20(%r10,%rcx,8)
movupd 0x20(%rdi,%rcx,8),%xmm4
movupd 0x30(%rdi,%rcx,8),%xmm5
mulpd  %xmm3,%xmm4
addpd  %xmm7,%xmm4
mulpd  %xmm3,%xmm5
addpd  %xmm1,%xmm5
movupd %xmm4,-0x10(%r10,%rcx,8)
movupd %xmm5,(%r10,%rcx,8)
add    $0x8,%rcx
add    $0xfffffffffffffffe,%r11
jne    1590 <Householder+0x480>
testb  $0x1,0x58(%rsp)
mov    %rax,%r10
je     164a <Householder+0x53a>
add    %r14,%rcx
movupd 0x0(%rbp,%rcx,8),%xmm1
movupd 0x10(%rbp,%rcx,8),%xmm4
movupd (%r8,%rcx,8),%xmm5
movupd 0x10(%r8,%rcx,8),%xmm6
mulpd  %xmm3,%xmm1
addpd  %xmm5,%xmm1
mulpd  %xmm3,%xmm4
addpd  %xmm6,%xmm4
movupd %xmm1,(%r8,%rcx,8)
movupd %xmm4,0x10(%r8,%rcx,8)
mov    0x48(%rsp),%rax
mov    0x50(%rsp),%rcx
cmp    %rcx,0x18(%rsp)
mov    0x10(%rsp),%r11
mov    0x8(%rsp),%rdx
je     155b <Householder+0x44b>
mov    %r15d,%edx
sub    %eax,%edx
mov    %rax,%rcx
test   $0x1,%dl
je     1690 <Householder+0x580>
movsd  0x0(%rbp,%rax,8),%xmm1
mulsd  %xmm2,%xmm1
addsd  (%r8,%rax,8),%xmm1
movsd  %xmm1,(%r8,%rax,8)
lea    0x1(%rax),%rcx
not    %rax
cmp    0x70(%rsp),%rax
mov    0x10(%rsp),%r11
mov    0x8(%rsp),%rdx
je     155b <Householder+0x44b>
nopl   0x0(%rax,%rax,1)
movsd  0x0(%rbp,%rcx,8),%xmm1
mulsd  %xmm2,%xmm1
addsd  (%r8,%rcx,8),%xmm1
movsd  %xmm1,(%r8,%rcx,8)
movsd  0x8(%rbp,%rcx,8),%xmm1
mulsd  %xmm2,%xmm1
addsd  0x8(%r8,%rcx,8),%xmm1
movsd  %xmm1,0x8(%r8,%rcx,8)
add    $0x2,%rcx
cmp    %rcx,%rsi
jne    16b0 <Householder+0x5a0>
jmp    155b <Householder+0x44b>
add    $0x78,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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