#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 */
scalar_t__ fabs(double);
double sqrt(double);
void HouseholderRealSym(double a[], int n, double d[], double e[]) {
  int m, k, j, i;
  double scale, hh, h, g, f;
  for (i = n - 1; i >= 1; i--) {
    m = i - 1;
    h = scale = 0;
    if (m > 0) {
      for (k = 0; k <= m; k++)
        scale += fabs(a[i * n + k]);
      if (scale == 0)
        e[i] = a[i * n + m];
      else {
        for (k = 0; k <= m; k++) {
          a[i * n + k] /= scale;
          h += a[i * n + k] * a[i * n + k];
        }
        f = a[i * n + m];
        g = (f >= 0 ? -sqrt(h) : sqrt(h));
        e[i] = scale * g;
        h -= f * g;
        a[i * n + m] = f - g;
        f = 0;
        for (j = 0; j <= m; j++) {
          a[j * n + i] = a[i * n + j] / h;
          g = 0;
          for (k = 0; k <= j; k++)
            g += a[j * n + k] * a[i * n + k];
          for (k = j + 1; k <= m; k++)
            g += a[k * n + j] * a[i * n + k];
          e[j] = g / h;
          f += e[j] * a[i * n + j];
        }
        hh = f / (h * 2);
        for (j = 0; j <= m; j++) {
          f = a[i * n + j];
          e[j] = g = e[j] - hh * f;
          for (k = 0; k <= j; k++)
            a[j * n + k] -= (f * e[k] + g * a[i * n + k]);
        }
      }
    } else
      e[i] = a[i * n + m];
    d[i] = h;
  }
  d[0] = e[0] = 0;
  for (i = 0; i < n; i++) {
    m = i - 1;
    if (d[i]) {
      for (j = 0; j <= m; j++) {
        g = 0;
        for (k = 0; k <= m; k++)
          g += a[i * n + k] * a[k * n + j];
        for (k = 0; k <= m; k++)
          a[k * n + j] -= g * a[k * n + i];
      }
    }
    d[i] = a[i * n + i];
    a[i * n + i] = 1;
    for (j = 0; j <= m; j++)
      a[j * n + i] = a[i * n + j] = 0;
  }
}


x86 assembly (clang, -O0):
<HouseholderRealSym>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x30(%rbp)
cmpl   $0x1,-0x30(%rbp)
jl     162a <HouseholderRealSym+0x50a>
mov    -0x30(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x24(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
cmpl   $0x0,-0x24(%rbp)
jle    15e7 <HouseholderRealSym+0x4c7>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x24(%rbp),%eax
jg     11b8 <HouseholderRealSym+0x98>
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1030 <fabs@plt>
cvtsi2sd %rax,%xmm0
addsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1174 <HouseholderRealSym+0x54>
movsd  -0x38(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    11f8 <HouseholderRealSym+0xd8>
jp     11f8 <HouseholderRealSym+0xd8>
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    15e2 <HouseholderRealSym+0x4c2>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x24(%rbp),%eax
jg     127b <HouseholderRealSym+0x15b>
movsd  -0x38(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    11ff <HouseholderRealSym+0xdf>
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x58(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     12d3 <HouseholderRealSym+0x1b3>
movsd  -0x48(%rbp),%xmm0
call   1040 <sqrt@plt>
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x60(%rbp)
jmp    12e2 <HouseholderRealSym+0x1c2>
movsd  -0x48(%rbp),%xmm0
call   1040 <sqrt@plt>
movsd  %xmm0,-0x60(%rbp)
movsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x38(%rbp),%xmm0
mulsd  -0x50(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x58(%rbp),%xmm0
movsd  -0x50(%rbp),%xmm2
movsd  -0x48(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x58(%rbp),%xmm0
subsd  -0x50(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jg     14c4 <HouseholderRealSym+0x3a4>
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  -0x48(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jg     1409 <HouseholderRealSym+0x2e9>
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    13b1 <HouseholderRealSym+0x291>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x24(%rbp),%eax
jg     146a <HouseholderRealSym+0x34a>
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1412 <HouseholderRealSym+0x2f2>
movsd  -0x50(%rbp),%xmm0
divsd  -0x48(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x58(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1365 <HouseholderRealSym+0x245>
movsd  -0x58(%rbp),%xmm0
movsd  0xb37(%rip),%xmm1
mulsd  -0x48(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jg     15dd <HouseholderRealSym+0x4bd>
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x40(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jg     15ca <HouseholderRealSym+0x4aa>
movsd  -0x58(%rbp),%xmm1
mov    -0x20(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x50(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    155c <HouseholderRealSym+0x43c>
jmp    15cf <HouseholderRealSym+0x4af>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    14e6 <HouseholderRealSym+0x3c6>
jmp    15e2 <HouseholderRealSym+0x4c2>
jmp    160a <HouseholderRealSym+0x4ea>
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x48(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x30(%rbp)
jmp    1140 <HouseholderRealSym+0x20>
mov    -0x20(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
mov    -0x18(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1830 <HouseholderRealSym+0x710>
mov    -0x30(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1681 <HouseholderRealSym+0x561>
jp     1681 <HouseholderRealSym+0x561>
jmp    1789 <HouseholderRealSym+0x669>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jg     1784 <HouseholderRealSym+0x664>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x24(%rbp),%eax
jg     16fb <HouseholderRealSym+0x5db>
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    16a3 <HouseholderRealSym+0x583>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x24(%rbp),%eax
jg     1771 <HouseholderRealSym+0x651>
movsd  -0x50(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movq   %xmm0,%rdx
movabs $0x8000000000000000,%rsi
xor    %rsi,%rdx
movq   %rdx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1702 <HouseholderRealSym+0x5e2>
jmp    1776 <HouseholderRealSym+0x656>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1688 <HouseholderRealSym+0x568>
jmp    1789 <HouseholderRealSym+0x669>
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  0x83b(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jg     181d <HouseholderRealSym+0x6fd>
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    17d1 <HouseholderRealSym+0x6b1>
jmp    1822 <HouseholderRealSym+0x702>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1647 <HouseholderRealSym+0x527>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<HouseholderRealSym>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xc8,%rsp
mov    %rcx,%rbx
mov    %rdx,0x30(%rsp)
mov    %rdi,%r9
mov    %esi,%edi
mov    %rsi,0x48(%rsp)
cmp    $0x2,%esi
mov    %r9,0x10(%rsp)
mov    %rdi,0x18(%rsp)
jl     1916 <HouseholderRealSym+0x7f6>
mov    0x48(%rsp),%rax
lea    -0x1(%rax),%ecx
lea    0x8(%r9),%rax
mov    %rax,0x58(%rsp)
lea    0x0(,%rdi,8),%rax
mov    %rax,0x50(%rsp)
lea    0x8(,%rdi,8),%rdx
mov    %rdx,0xc0(%rsp)
lea    -0x1(%rdi),%rdx
imul   %rdi,%rdx
lea    0x1(%rdx),%rsi
mov    %rsi,0xa0(%rsp)
mov    %rcx,%rsi
neg    %rsi
mov    %rsi,0x98(%rsp)
lea    (%r9,%rdx,8),%r14
add    $0x18,%r14
neg    %rax
mov    %rax,0xb0(%rsp)
lea    (%r9,%rdx,8),%r12
add    $0x8,%r12
lea    (%r9,%rdx,8),%r13
lea    0x10(%r9),%rax
mov    %rax,0x90(%rsp)
xor    %esi,%esi
mov    %rdx,0xa8(%rsp)
mov    %rdx,0x60(%rsp)
mov    %rdi,%rdx
mov    %rcx,%rax
mov    %rcx,0xb8(%rsp)
mov    %rcx,0x20(%rsp)
jmp    1261 <HouseholderRealSym+0x141>
cs nopw 0x0(%rax,%rax,1)
nop
mov    0x28(%rsp),%rax
add    %r15,%rax
mov    0x10(%rsp),%r9
movsd  (%r9,%rax,8),%xmm0
mov    0x40(%rsp),%rdx
movsd  %xmm0,(%rbx,%rdx,8)
xorpd  %xmm9,%xmm9
mov    0x18(%rsp),%rdi
mov    0x68(%rsp),%r14
mov    0x30(%rsp),%rax
movsd  %xmm9,(%rax,%rdx,8)
addq   $0xffffffffffffffff,0x20(%rsp)
add    $0x1,%rsi
mov    0xb0(%rsp),%rax
add    %rax,%r14
sub    %rdi,0x60(%rsp)
add    %rax,%r12
add    %rax,%r13
cmp    $0x2,%rcx
jle    1916 <HouseholderRealSym+0x7f6>
lea    -0x1(%rdx),%rcx
mov    %rcx,%rax
imul   %rdi,%rax
mov    %rax,0x28(%rsp)
cmp    $0x2,%rdx
je     18f9 <HouseholderRealSym+0x7d9>
mov    %rcx,0x40(%rsp)
mov    %rdx,0x70(%rsp)
mov    %rsi,0x78(%rsp)
not    %rsi
mov    0xb8(%rsp),%rax
mov    %rsi,0x38(%rsp)
add    %rsi,%rax
cmp    $0x3,%rax
jae    12b0 <HouseholderRealSym+0x190>
xorpd  %xmm4,%xmm4
xor    %ebp,%ebp
jmp    134c <HouseholderRealSym+0x22c>
xchg   %ax,%ax
mov    0x20(%rsp),%r15
and    $0xfffffffffffffffc,%r15
xorps  %xmm4,%xmm4
xor    %ebp,%ebp
xchg   %ax,%ax
movsd  %xmm4,0x8(%rsp)
movsd  -0x18(%r14,%rbp,8),%xmm0
call   1030 <fabs@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %rax,%xmm0
addsd  0x8(%rsp),%xmm0
movsd  %xmm0,0x8(%rsp)
movsd  -0x10(%r14,%rbp,8),%xmm0
call   1030 <fabs@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %rax,%xmm0
addsd  0x8(%rsp),%xmm0
movsd  %xmm0,0x8(%rsp)
movsd  -0x8(%r14,%rbp,8),%xmm0
call   1030 <fabs@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %rax,%xmm0
addsd  0x8(%rsp),%xmm0
movsd  %xmm0,0x8(%rsp)
movsd  (%r14,%rbp,8),%xmm0
call   1030 <fabs@plt>
xorps  %xmm4,%xmm4
cvtsi2sd %rax,%xmm4
addsd  0x8(%rsp),%xmm4
add    $0x4,%rbp
cmp    %rbp,%r15
jne    12c0 <HouseholderRealSym+0x1a0>
mov    %r14,0x68(%rsp)
testb  $0x3,0x20(%rsp)
je     13ad <HouseholderRealSym+0x28d>
mov    0x20(%rsp),%r15
and    $0x3,%r15d
add    0x60(%rsp),%rbp
mov    0x10(%rsp),%rax
lea    (%rax,%rbp,8),%rbp
xor    %r14d,%r14d
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  %xmm4,0x8(%rsp)
movsd  0x0(%rbp,%r14,8),%xmm0
call   1030 <fabs@plt>
movsd  0x8(%rsp),%xmm4
xorps  %xmm0,%xmm0
cvtsi2sd %rax,%xmm0
addsd  %xmm0,%xmm4
add    $0x1,%r14
cmp    %r14,%r15
jne    1380 <HouseholderRealSym+0x260>
mov    0x70(%rsp),%rcx
lea    -0x2(%rcx),%r15
ucomisd 0xc52(%rip),%xmm4
mov    0x78(%rsp),%rsi
jne    13cb <HouseholderRealSym+0x2ab>
jnp    1200 <HouseholderRealSym+0xe0>
mov    0x38(%rsp),%rax
cmp    0x98(%rsp),%rax
mov    0x10(%rsp),%r9
mov    0x18(%rsp),%rdi
jne    13ed <HouseholderRealSym+0x2cd>
xorpd  %xmm9,%xmm9
xor    %eax,%eax
jmp    1443 <HouseholderRealSym+0x323>
mov    0x20(%rsp),%rcx
and    $0xfffffffffffffffe,%rcx
xorpd  %xmm0,%xmm0
xor    %eax,%eax
nopl   0x0(%rax)
movsd  -0x8(%r12,%rax,8),%xmm1
movsd  (%r12,%rax,8),%xmm9
divsd  %xmm4,%xmm1
movsd  %xmm1,-0x8(%r12,%rax,8)
mulsd  %xmm1,%xmm1
addsd  %xmm0,%xmm1
divsd  %xmm4,%xmm9
movsd  %xmm9,(%r12,%rax,8)
mulsd  %xmm9,%xmm9
addsd  %xmm1,%xmm9
add    $0x2,%rax
movapd %xmm9,%xmm0
cmp    %rax,%rcx
jne    1400 <HouseholderRealSym+0x2e0>
mov    %rsi,%r14
imul   %rdi,%r14
testb  $0x1,0x20(%rsp)
je     146f <HouseholderRealSym+0x34f>
add    0x28(%rsp),%rax
movsd  (%r9,%rax,8),%xmm0
divsd  %xmm4,%xmm0
movsd  %xmm0,(%r9,%rax,8)
mulsd  %xmm0,%xmm0
addsd  %xmm0,%xmm9
mov    0xa8(%rsp),%rbp
sub    %r14,%rbp
mov    0x28(%rsp),%rax
lea    (%rax,%r15,1),%rcx
movsd  (%r9,%rcx,8),%xmm1
ucomisd 0xb7e(%rip),%xmm9
jb     149e <HouseholderRealSym+0x37e>
xorps  %xmm0,%xmm0
sqrtsd %xmm9,%xmm0
jmp    14ee <HouseholderRealSym+0x3ce>
movapd %xmm9,%xmm0
movsd  %xmm9,0x38(%rsp)
movsd  %xmm4,0x8(%rsp)
mov    %rcx,0x88(%rsp)
movsd  %xmm1,0x80(%rsp)
call   1040 <sqrt@plt>
movsd  0x80(%rsp),%xmm1
mov    0x88(%rsp),%rcx
movsd  0x8(%rsp),%xmm4
movsd  0x38(%rsp),%xmm9
mov    0x18(%rsp),%rdi
mov    0x10(%rsp),%r9
lea    (%r9,%rbp,8),%rax
mov    %rax,0x8(%rsp)
mov    0xa0(%rsp),%rax
sub    %r14,%rax
mov    %rax,0x38(%rsp)
xorpd  %xmm8,%xmm8
xorpd  %xmm2,%xmm2
cmpnlesd %xmm1,%xmm2
movapd %xmm2,%xmm3
andpd  %xmm0,%xmm3
xorpd  0xadb(%rip),%xmm0
andnpd %xmm0,%xmm2
orpd   %xmm3,%xmm2
mulsd  %xmm2,%xmm4
mov    0x40(%rsp),%rax
movsd  %xmm4,(%rbx,%rax,8)
movapd %xmm1,%xmm0
mulsd  %xmm2,%xmm0
subsd  %xmm0,%xmm9
subsd  %xmm2,%xmm1
movsd  %xmm1,(%r9,%rcx,8)
mov    $0x1,%r8d
mov    0x58(%rsp),%rdx
xor    %ecx,%ecx
mov    0x50(%rsp),%r14
jmp    159a <HouseholderRealSym+0x47a>
cs nopw 0x0(%rax,%rax,1)
divsd  %xmm9,%xmm0
movsd  %xmm0,(%rbx,%rcx,8)
mulsd  (%r9,%r10,8),%xmm0
addsd  %xmm0,%xmm8
add    $0x1,%r8
add    %r14,%rdx
mov    %r11,%rcx
cmp    0x20(%rsp),%r11
je     1680 <HouseholderRealSym+0x560>
lea    0x1(%rcx),%r11
mov    0x28(%rsp),%rax
lea    (%rcx,%rax,1),%r10
movsd  (%r9,%r10,8),%xmm0
divsd  %xmm9,%xmm0
mov    %rcx,%rsi
imul   %rdi,%rsi
mov    0x40(%rsp),%rax
add    %rsi,%rax
movsd  %xmm0,(%r9,%rax,8)
test   %rcx,%rcx
je     1670 <HouseholderRealSym+0x550>
mov    %r8,%rax
and    $0xfffffffffffffffe,%rax
xorpd  %xmm2,%xmm2
xor    %ebp,%ebp
nopl   (%rax)
movsd  -0x8(%rdx,%rbp,8),%xmm3
movsd  (%rdx,%rbp,8),%xmm0
mulsd  -0x8(%r12,%rbp,8),%xmm3
mulsd  (%r12,%rbp,8),%xmm0
addsd  %xmm2,%xmm3
addsd  %xmm3,%xmm0
add    $0x2,%rbp
movapd %xmm0,%xmm2
cmp    %rbp,%rax
jne    15e0 <HouseholderRealSym+0x4c0>
test   $0x1,%r11b
je     162b <HouseholderRealSym+0x50b>
add    %rbp,%rsi
movsd  (%r9,%rsi,8),%xmm2
add    0x28(%rsp),%rbp
mulsd  (%r9,%rbp,8),%xmm2
addsd  %xmm2,%xmm0
cmp    %r15,%rcx
jge    1570 <HouseholderRealSym+0x450>
movslq %r8d,%rsi
mov    %rdi,%rax
imul   %rsi,%rax
add    %rcx,%rax
add    $0xffffffffffffffff,%rsi
lea    (%r9,%rax,8),%rax
nopl   0x0(%rax)
movsd  (%rax),%xmm2
mulsd  0x8(%r13,%rsi,8),%xmm2
addsd  %xmm2,%xmm0
add    $0x1,%rsi
add    %r14,%rax
cmp    %r15,%rsi
jl     1650 <HouseholderRealSym+0x530>
jmp    1570 <HouseholderRealSym+0x450>
xor    %ebp,%ebp
xorpd  %xmm0,%xmm0
test   $0x1,%r11b
jne    1613 <HouseholderRealSym+0x4f3>
jmp    162b <HouseholderRealSym+0x50b>
xchg   %ax,%ax
movapd %xmm9,%xmm0
addsd  %xmm9,%xmm0
xorpd  0x96d(%rip),%xmm8
divsd  %xmm0,%xmm8
mov    $0x1,%edx
mov    0x58(%rsp),%rsi
mov    0x90(%rsp),%r10
xor    %r8d,%r8d
jmp    16c5 <HouseholderRealSym+0x5a5>
nop
add    $0x1,%rdx
add    %r14,%r10
add    %r14,%rsi
cmp    0x20(%rsp),%r8
je     18e0 <HouseholderRealSym+0x7c0>
mov    %r8,%rax
add    $0x1,%r8
mov    %rax,%r11
imul   %rdi,%r11
mov    0x28(%rsp),%rcx
add    %rax,%rcx
movsd  (%r9,%rcx,8),%xmm0
movapd %xmm8,%xmm2
mulsd  %xmm0,%xmm2
addsd  (%rbx,%rax,8),%xmm2
movsd  %xmm2,(%rbx,%rax,8)
cmp    $0x3,%rax
jae    1710 <HouseholderRealSym+0x5f0>
xor    %r15d,%r15d
jmp    1830 <HouseholderRealSym+0x710>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    %r8,%rdi
lea    (%r9,%r11,8),%r14
mov    0xc0(%rsp),%rbp
imul   %rax,%rbp
add    0x58(%rsp),%rbp
mov    0x38(%rsp),%rcx
add    %rax,%rcx
lea    (%r9,%rcx,8),%r15
lea    (%rbx,%rax,8),%rcx
add    $0x8,%rcx
cmp    %rcx,%r14
setb   %r8b
cmp    %rbx,%rbp
seta   %r9b
cmp    %r15,%r14
setb   %cl
cmp    0x8(%rsp),%rbp
seta   %bpl
test   %r9b,%r8b
jne    1817 <HouseholderRealSym+0x6f7>
mov    $0x0,%r15d
and    %bpl,%cl
mov    0x10(%rsp),%r9
mov    0x18(%rsp),%rcx
mov    0x50(%rsp),%r14
mov    %rdi,%r8
mov    %rcx,%rdi
jne    1830 <HouseholderRealSym+0x710>
mov    %rdx,%rcx
and    $0xfffffffffffffffc,%rcx
mov    %r8,%r15
and    $0xfffffffffffffffc,%r15
movapd %xmm2,%xmm3
unpcklpd %xmm2,%xmm3
movapd %xmm0,%xmm4
unpcklpd %xmm0,%xmm4
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
movupd (%rbx,%rbp,8),%xmm5
movupd 0x10(%rbx,%rbp,8),%xmm6
movupd 0x0(%r13,%rbp,8),%xmm7
movupd 0x10(%r13,%rbp,8),%xmm1
mulpd  %xmm3,%xmm7
mulpd  %xmm3,%xmm1
mulpd  %xmm4,%xmm5
addpd  %xmm7,%xmm5
mulpd  %xmm4,%xmm6
addpd  %xmm1,%xmm6
movupd -0x10(%r10,%rbp,8),%xmm1
subpd  %xmm5,%xmm1
movupd (%r10,%rbp,8),%xmm5
subpd  %xmm6,%xmm5
movupd %xmm1,-0x10(%r10,%rbp,8)
movupd %xmm5,(%r10,%rbp,8)
add    $0x4,%rbp
cmp    %rbp,%rcx
jne    17b0 <HouseholderRealSym+0x690>
cmp    %r15,%r8
je     16b0 <HouseholderRealSym+0x590>
jmp    1830 <HouseholderRealSym+0x710>
xor    %r15d,%r15d
mov    0x10(%rsp),%r9
mov    0x18(%rsp),%rcx
mov    0x50(%rsp),%r14
mov    %rdi,%r8
mov    %rcx,%rdi
nop
mov    %r15,%rcx
test   $0x1,%dl
je     1872 <HouseholderRealSym+0x752>
mov    0x28(%rsp),%rcx
add    %r15,%rcx
movsd  (%r9,%rcx,8),%xmm1
mulsd  %xmm2,%xmm1
movsd  (%rbx,%r15,8),%xmm3
mulsd  %xmm0,%xmm3
addsd  %xmm1,%xmm3
add    %r15,%r11
movsd  (%r9,%r11,8),%xmm1
subsd  %xmm3,%xmm1
movsd  %xmm1,(%r9,%r11,8)
mov    %r15,%rcx
or     $0x1,%rcx
cmp    %r15,%rax
je     16b0 <HouseholderRealSym+0x590>
nopl   0x0(%rax,%rax,1)
movsd  0x0(%r13,%rcx,8),%xmm1
mulsd  %xmm2,%xmm1
movsd  (%rbx,%rcx,8),%xmm3
mulsd  %xmm0,%xmm3
addsd  %xmm1,%xmm3
movsd  -0x8(%rsi,%rcx,8),%xmm1
subsd  %xmm3,%xmm1
movsd  (%rsi,%rcx,8),%xmm3
movsd  %xmm1,-0x8(%rsi,%rcx,8)
movsd  0x8(%r13,%rcx,8),%xmm1
mulsd  %xmm2,%xmm1
movsd  0x8(%rbx,%rcx,8),%xmm4
mulsd  %xmm0,%xmm4
addsd  %xmm1,%xmm4
subsd  %xmm4,%xmm3
movsd  %xmm3,(%rsi,%rcx,8)
add    $0x2,%rcx
cmp    %rcx,%rdx
jne    1880 <HouseholderRealSym+0x760>
jmp    16b0 <HouseholderRealSym+0x590>
nopl   (%rax)
mov    0x68(%rsp),%r14
mov    0x78(%rsp),%rsi
mov    0x70(%rsp),%rcx
mov    0x40(%rsp),%rdx
jmp    122c <HouseholderRealSym+0x10c>
mov    0x28(%rsp),%rax
movsd  (%r9,%rax,8),%xmm0
movsd  %xmm0,(%rbx,%rcx,8)
mov    0x30(%rsp),%rax
movq   $0x0,(%rax,%rcx,8)
movq   $0x0,(%rbx)
mov    0x30(%rsp),%rax
movq   $0x0,(%rax)
cmpl   $0x0,0x48(%rsp)
jle    1bbb <HouseholderRealSym+0xa9b>
mov    0x10(%rsp),%r12
lea    0x8(%r12),%rax
mov    0x18(%rsp),%rcx
lea    0x0(,%rcx,8),%r14
mov    %rcx,%rdx
shl    $0x4,%rdx
lea    (%r12,%rcx,8),%r9
xorpd  %xmm0,%xmm0
xor    %r11d,%r11d
movapd 0x69b(%rip),%xmm1
mov    %r9,%r13
xorpd  %xmm2,%xmm2
nopl   0x0(%rax)
mov    %r11,%rdi
and    $0xfffffffffffffffe,%rdi
mov    %r11,%r10
imul   0x18(%rsp),%r10
ucomisd %xmm0,%xmm2
jne    198c <HouseholderRealSym+0x86c>
jnp    1ae0 <HouseholderRealSym+0x9c0>
test   %r11,%r11
je     1ae0 <HouseholderRealSym+0x9c0>
mov    %r9,%rbp
mov    0x10(%rsp),%rsi
xor    %r15d,%r15d
jmp    19c5 <HouseholderRealSym+0x8a5>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%r15
add    $0x8,%rsi
add    $0x8,%rbp
cmp    %r11,%r15
je     1ae0 <HouseholderRealSym+0x9c0>
cmp    $0x1,%r11
jne    19e0 <HouseholderRealSym+0x8c0>
xorpd  %xmm2,%xmm2
xor    %ebx,%ebx
jmp    1a1d <HouseholderRealSym+0x8fd>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
xorpd  %xmm3,%xmm3
mov    %rsi,%rcx
xor    %ebx,%ebx
nopl   0x0(%rax)
movsd  -0x8(%rax,%rbx,8),%xmm4
mulsd  (%rcx),%xmm4
movsd  (%rax,%rbx,8),%xmm2
addsd  %xmm3,%xmm4
mulsd  (%rcx,%r14,1),%xmm2
addsd  %xmm4,%xmm2
add    $0x2,%rbx
add    %rdx,%rcx
movapd %xmm2,%xmm3
cmp    %rbx,%rdi
jne    19f0 <HouseholderRealSym+0x8d0>
test   $0x1,%r11b
je     1a45 <HouseholderRealSym+0x925>
lea    (%rbx,%r10,1),%rcx
mov    0x10(%rsp),%r8
movsd  (%r8,%rcx,8),%xmm3
imul   0x18(%rsp),%rbx
add    %r15,%rbx
mulsd  (%r8,%rbx,8),%xmm3
addsd  %xmm3,%xmm2
xorpd  %xmm1,%xmm2
cmp    $0x1,%r11
jne    1a60 <HouseholderRealSym+0x940>
xor    %ebx,%ebx
jmp    1aa7 <HouseholderRealSym+0x987>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
xor    %ecx,%ecx
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%r12,%rcx,1),%xmm3
mulsd  %xmm2,%xmm3
addsd  (%rsi,%rcx,1),%xmm3
movsd  %xmm3,(%rsi,%rcx,1)
movsd  0x0(%r13,%rcx,1),%xmm3
mulsd  %xmm2,%xmm3
addsd  0x0(%rbp,%rcx,1),%xmm3
movsd  %xmm3,0x0(%rbp,%rcx,1)
add    $0x2,%rbx
add    %rdx,%rcx
cmp    %rbx,%rdi
jne    1a70 <HouseholderRealSym+0x950>
test   $0x1,%r11b
je     19b0 <HouseholderRealSym+0x890>
imul   0x18(%rsp),%rbx
lea    (%rbx,%r11,1),%rcx
mov    0x10(%rsp),%r8
mulsd  (%r8,%rcx,8),%xmm2
add    %r15,%rbx
addsd  (%r8,%rbx,8),%xmm2
movsd  %xmm2,(%r8,%rbx,8)
jmp    19b0 <HouseholderRealSym+0x890>
nopw   0x0(%rax,%rax,1)
mov    %r11d,%ecx
imul   0x48(%rsp),%ecx
movslq %ecx,%rcx
add    %r11,%rcx
mov    0x10(%rsp),%rbp
movsd  0x0(%rbp,%rcx,8),%xmm2
mov    0x30(%rsp),%rsi
movsd  %xmm2,(%rsi,%r11,8)
movabs $0x3ff0000000000000,%rsi
mov    %rsi,0x0(%rbp,%rcx,8)
test   %r11,%r11
je     1b95 <HouseholderRealSym+0xa75>
cmp    $0x1,%r11
jne    1b30 <HouseholderRealSym+0xa10>
xor    %ecx,%ecx
jmp    1b6e <HouseholderRealSym+0xa4e>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
xor    %esi,%esi
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movq   $0x0,-0x8(%rax,%rcx,8)
movq   $0x0,(%r12,%rsi,1)
movq   $0x0,(%rax,%rcx,8)
movq   $0x0,0x0(%r13,%rsi,1)
add    $0x2,%rcx
add    %rdx,%rsi
cmp    %rcx,%rdi
jne    1b40 <HouseholderRealSym+0xa20>
test   $0x1,%r11b
je     1b95 <HouseholderRealSym+0xa75>
add    %rcx,%r10
mov    0x10(%rsp),%rsi
movq   $0x0,(%rsi,%r10,8)
imul   0x18(%rsp),%rcx
add    %r11,%rcx
movq   $0x0,(%rsi,%rcx,8)
add    $0x1,%r11
cmp    0x18(%rsp),%r11
je     1bbb <HouseholderRealSym+0xa9b>
mov    0x30(%rsp),%rcx
movsd  (%rcx,%r11,8),%xmm2
add    %r14,%rax
add    $0x8,%r13
add    $0x8,%r12
jmp    1970 <HouseholderRealSym+0x850>
add    $0xc8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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