#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 signbit(double);
double sqrt(double);
void matrixToHessenbergForm_Householder_double(double *A) {
  double vk[2] = {
      0.0,
      0.0,
  };
  double vk_norm = 0.0;
  double dot_array[3] = {0.0, 0.0, 0.0};
  vk[0] = A[3];
  vk[1] = A[6];
  vk_norm = sqrt(vk[0] * vk[0] + vk[1] * vk[1]);
  vk[0] += signbit(vk[0]) * vk_norm;
  vk_norm = sqrt(vk[0] * vk[0] + vk[1] * vk[1]);
  vk[0] /= vk_norm;
  vk[1] /= vk_norm;
  dot_array[0] = vk[0] * A[3] + vk[1] * A[6];
  dot_array[1] = vk[0] * A[4] + vk[1] * A[7];
  dot_array[2] = vk[0] * A[5] + vk[1] * A[8];
  A[3] -= 2 * vk[0] * dot_array[0];
  A[4] -= 2 * vk[0] * dot_array[1];
  A[5] -= 2 * vk[0] * dot_array[2];
  A[6] -= 2 * vk[1] * dot_array[0];
  A[7] -= 2 * vk[1] * dot_array[1];
  A[8] -= 2 * vk[1] * dot_array[2];
  dot_array[0] = vk[0] * A[1] + vk[1] * A[2];
  dot_array[1] = vk[0] * A[4] + vk[1] * A[5];
  dot_array[2] = vk[0] * A[7] + vk[1] * A[8];
  A[1] -= 2 * vk[0] * dot_array[0];
  A[2] -= 2 * vk[1] * dot_array[0];
  A[4] -= 2 * vk[0] * dot_array[1];
  A[5] -= 2 * vk[1] * dot_array[1];
  A[7] -= 2 * vk[0] * dot_array[2];
  A[8] -= 2 * vk[1] * dot_array[2];
  A[2] = -A[2];
  A[5] = -A[5];
  A[7] = -A[7];
}


x86 assembly (clang, -O0):
<matrixToHessenbergForm_Householder_double>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
lea    -0x20(%rbp),%rdi
xor    %esi,%esi
mov    $0x10,%edx
call   1030 <memset@plt>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
lea    -0x40(%rbp),%rdi
xor    %esi,%esi
mov    $0x18,%edx
call   1030 <memset@plt>
mov    -0x8(%rbp),%rax
movsd  0x18(%rax),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x8(%rbp),%rax
movsd  0x30(%rax),%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  -0x20(%rbp),%xmm0
movsd  -0x20(%rbp),%xmm2
movsd  -0x18(%rbp),%xmm1
mulsd  -0x18(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1050 <sqrt@plt>
movsd  %xmm0,-0x28(%rbp)
movsd  -0x20(%rbp),%xmm0
call   1040 <signbit@plt>
movsd  -0x28(%rbp),%xmm2
movsd  -0x20(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x20(%rbp),%xmm0
movsd  -0x20(%rbp),%xmm2
movsd  -0x18(%rbp),%xmm1
mulsd  -0x18(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1050 <sqrt@plt>
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm1
movsd  -0x20(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x28(%rbp),%xmm1
movsd  -0x18(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movsd  0x18(%rax),%xmm2
movsd  -0x18(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mulsd  0x30(%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movsd  0x20(%rax),%xmm2
movsd  -0x18(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mulsd  0x38(%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movsd  0x28(%rax),%xmm2
movsd  -0x18(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mulsd  0x40(%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  0xd6a(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm2
mov    -0x8(%rbp),%rax
movsd  0x18(%rax),%xmm1
movq   %xmm0,%rcx
movabs $0x8000000000000000,%rdx
xor    %rdx,%rcx
movq   %rcx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x18(%rax)
movsd  0xd2b(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm2
mov    -0x8(%rbp),%rax
movsd  0x20(%rax),%xmm1
movq   %xmm0,%rcx
movabs $0x8000000000000000,%rdx
xor    %rdx,%rcx
movq   %rcx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x20(%rax)
movsd  0xcec(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
mov    -0x8(%rbp),%rax
movsd  0x28(%rax),%xmm1
movq   %xmm0,%rcx
movabs $0x8000000000000000,%rdx
xor    %rdx,%rcx
movq   %rcx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x28(%rax)
movsd  0xcad(%rip),%xmm0
mulsd  -0x18(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm2
mov    -0x8(%rbp),%rax
movsd  0x30(%rax),%xmm1
movq   %xmm0,%rcx
movabs $0x8000000000000000,%rdx
xor    %rdx,%rcx
movq   %rcx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x30(%rax)
movsd  0xc6e(%rip),%xmm0
mulsd  -0x18(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm2
mov    -0x8(%rbp),%rax
movsd  0x38(%rax),%xmm1
movq   %xmm0,%rcx
movabs $0x8000000000000000,%rdx
xor    %rdx,%rcx
movq   %rcx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x38(%rax)
movsd  0xc2f(%rip),%xmm0
mulsd  -0x18(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
mov    -0x8(%rbp),%rax
movsd  0x40(%rax),%xmm1
movq   %xmm0,%rcx
movabs $0x8000000000000000,%rdx
xor    %rdx,%rcx
movq   %rcx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x40(%rax)
movsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movsd  0x8(%rax),%xmm2
movsd  -0x18(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mulsd  0x10(%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movsd  0x20(%rax),%xmm2
movsd  -0x18(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mulsd  0x28(%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movsd  0x38(%rax),%xmm2
movsd  -0x18(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mulsd  0x40(%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  0xb75(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm2
mov    -0x8(%rbp),%rax
movsd  0x8(%rax),%xmm1
movq   %xmm0,%rcx
movabs $0x8000000000000000,%rdx
xor    %rdx,%rcx
movq   %rcx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x8(%rax)
movsd  0xb36(%rip),%xmm0
mulsd  -0x18(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm2
mov    -0x8(%rbp),%rax
movsd  0x10(%rax),%xmm1
movq   %xmm0,%rcx
movabs $0x8000000000000000,%rdx
xor    %rdx,%rcx
movq   %rcx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x10(%rax)
movsd  0xaf7(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm2
mov    -0x8(%rbp),%rax
movsd  0x20(%rax),%xmm1
movq   %xmm0,%rcx
movabs $0x8000000000000000,%rdx
xor    %rdx,%rcx
movq   %rcx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x20(%rax)
movsd  0xab8(%rip),%xmm0
mulsd  -0x18(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm2
mov    -0x8(%rbp),%rax
movsd  0x28(%rax),%xmm1
movq   %xmm0,%rcx
movabs $0x8000000000000000,%rdx
xor    %rdx,%rcx
movq   %rcx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x28(%rax)
movsd  0xa79(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
mov    -0x8(%rbp),%rax
movsd  0x38(%rax),%xmm1
movq   %xmm0,%rcx
movabs $0x8000000000000000,%rdx
xor    %rdx,%rcx
movq   %rcx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x38(%rax)
movsd  0xa3a(%rip),%xmm0
mulsd  -0x18(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
mov    -0x8(%rbp),%rax
movsd  0x40(%rax),%xmm1
movq   %xmm0,%rcx
movabs $0x8000000000000000,%rdx
xor    %rdx,%rcx
movq   %rcx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x40(%rax)
mov    -0x8(%rbp),%rax
movsd  0x10(%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x8(%rbp),%rax
movsd  %xmm0,0x10(%rax)
mov    -0x8(%rbp),%rax
movsd  0x28(%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x8(%rbp),%rax
movsd  %xmm0,0x28(%rax)
mov    -0x8(%rbp),%rax
movsd  0x38(%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x8(%rbp),%rax
movsd  %xmm0,0x38(%rax)
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<matrixToHessenbergForm_Householder_double>:
push   %rbx
sub    $0x40,%rsp
mov    %rdi,%rbx
movsd  0x18(%rdi),%xmm2
movsd  0x30(%rdi),%xmm0
movapd %xmm0,0x30(%rsp)
movapd %xmm0,%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm2,0x10(%rsp)
movapd %xmm2,%xmm0
mulsd  %xmm2,%xmm0
movsd  %xmm1,0x28(%rsp)
addsd  %xmm1,%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     1168 <matrixToHessenbergForm_Householder_double+0x48>
sqrtsd %xmm0,%xmm0
jmp    116d <matrixToHessenbergForm_Householder_double+0x4d>
call   1040 <sqrt@plt>
movsd  %xmm0,0x8(%rsp)
movsd  0x10(%rsp),%xmm0
call   1030 <signbit@plt>
movapd %xmm0,%xmm1
mulsd  0x8(%rsp),%xmm1
addsd  0x10(%rsp),%xmm1
movapd %xmm1,%xmm0
mulsd  %xmm1,%xmm0
addsd  0x28(%rsp),%xmm0
ucomisd 0xe64(%rip),%xmm0
jb     11ac <matrixToHessenbergForm_Householder_double+0x8c>
sqrtsd %xmm0,%xmm0
jmp    11bd <matrixToHessenbergForm_Householder_double+0x9d>
movapd %xmm1,0x10(%rsp)
call   1040 <sqrt@plt>
movapd 0x10(%rsp),%xmm1
divsd  %xmm0,%xmm1
movapd 0x30(%rsp),%xmm3
divsd  %xmm0,%xmm3
movsd  0x30(%rbx),%xmm9
movsd  0x20(%rbx),%xmm4
movsd  0x38(%rbx),%xmm8
movapd %xmm3,%xmm0
mulsd  %xmm8,%xmm0
movapd %xmm1,%xmm6
mulsd  %xmm4,%xmm6
addsd  %xmm0,%xmm6
movsd  0x28(%rbx),%xmm7
movsd  0x40(%rbx),%xmm11
movapd %xmm3,%xmm0
mulsd  %xmm11,%xmm0
movapd %xmm1,%xmm2
mulsd  %xmm7,%xmm2
addsd  %xmm0,%xmm2
movsd  0xde6(%rip),%xmm13
movapd %xmm1,%xmm5
mulsd  %xmm13,%xmm5
movapd %xmm5,%xmm0
mulsd  %xmm6,%xmm0
addsd  %xmm4,%xmm0
movapd %xmm5,%xmm12
mulsd  %xmm2,%xmm12
addsd  %xmm7,%xmm12
mulsd  %xmm3,%xmm13
movapd %xmm3,%xmm4
unpcklpd %xmm3,%xmm4
movapd %xmm9,%xmm7
unpcklpd %xmm12,%xmm7
mulpd  %xmm4,%xmm7
movapd %xmm1,%xmm4
unpcklpd %xmm1,%xmm4
movsd  0x18(%rbx),%xmm10
unpcklpd %xmm0,%xmm10
mulpd  %xmm10,%xmm4
addpd  %xmm7,%xmm4
movapd %xmm13,%xmm0
mulsd  %xmm4,%xmm0
addsd  %xmm9,%xmm0
movsd  %xmm0,0x30(%rbx)
mulsd  %xmm13,%xmm6
addsd  %xmm8,%xmm6
mulsd  %xmm13,%xmm2
addsd  %xmm11,%xmm2
movsd  0x8(%rbx),%xmm9
movsd  0x10(%rbx),%xmm8
movapd %xmm3,%xmm0
mulsd  %xmm8,%xmm0
movapd %xmm1,%xmm7
mulsd  %xmm9,%xmm7
addsd  %xmm0,%xmm7
mulsd  %xmm2,%xmm3
mulsd  %xmm6,%xmm1
addsd  %xmm3,%xmm1
movapd %xmm5,%xmm0
mulsd  %xmm7,%xmm0
addsd  %xmm9,%xmm0
movsd  %xmm0,0x8(%rbx)
mulsd  %xmm13,%xmm7
addsd  %xmm8,%xmm7
movapd %xmm5,%xmm0
unpcklpd %xmm5,%xmm0
mulpd  %xmm4,%xmm0
addpd  %xmm10,%xmm0
movupd %xmm0,0x18(%rbx)
unpckhpd %xmm4,%xmm4
mulsd  %xmm13,%xmm4
addsd  %xmm12,%xmm4
mulsd  %xmm1,%xmm5
addsd  %xmm6,%xmm5
mulsd  %xmm13,%xmm1
addsd  %xmm2,%xmm1
movsd  %xmm1,0x40(%rbx)
movapd 0xce8(%rip),%xmm0
xorpd  %xmm0,%xmm7
movlpd %xmm7,0x10(%rbx)
xorpd  %xmm0,%xmm4
movlpd %xmm4,0x28(%rbx)
xorpd  %xmm0,%xmm5
movlpd %xmm5,0x38(%rbx)
add    $0x40,%rsp
pop    %rbx
ret

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

LLVM IR (-O2):
define void @matrixToHessenbergForm_Householder_double(double* nocapture %0) {
%2 = getelementptr inbounds double, double* %0, i64 3
%3 = load double, double* %2
%4 = getelementptr inbounds double, double* %0, i64 6
%5 = load double, double* %4
%6 = fmul double %5, %5
%7 = tail call double @llvm.fmuladd.f64(double %3, double %3, double %6)
%8 = tail call double @sqrt(double %7)
%9 = tail call double @signbit(double %3)
%10 = tail call double @llvm.fmuladd.f64(double %9, double %8, double %3)
%11 = tail call double @llvm.fmuladd.f64(double %10, double %10, double %6)
%12 = tail call double @sqrt(double %11)
%13 = fdiv double %10, %12
%14 = fdiv double %5, %12
%15 = load double, double* %2
%16 = load double, double* %4
%17 = getelementptr inbounds double, double* %0, i64 4
%18 = load double, double* %17
%19 = getelementptr inbounds double, double* %0, i64 7
%20 = load double, double* %19
%21 = fmul double %14, %20
%22 = tail call double @llvm.fmuladd.f64(double %13, double %18, double %21)
%23 = getelementptr inbounds double, double* %0, i64 5
%24 = load double, double* %23
%25 = getelementptr inbounds double, double* %0, i64 8
%26 = load double, double* %25
%27 = fmul double %14, %26
%28 = tail call double @llvm.fmuladd.f64(double %13, double %24, double %27)
%29 = fmul double %13, -2.000000e+00
%30 = tail call double @llvm.fmuladd.f64(double %29, double %22, double %18)
%31 = tail call double @llvm.fmuladd.f64(double %29, double %28, double %24)
%32 = fmul double %14, -2.000000e+00
%33 = insertelement <2 x double> poison, double %14, i64 0
%34 = shufflevector <2 x double> %33, <2 x double> poison, <2 x i32> zeroinitializer
%35 = insertelement <2 x double> poison, double %16, i64 0
%36 = insertelement <2 x double> %35, double %31, i64 1
%37 = fmul <2 x double> %34, %36
%38 = insertelement <2 x double> poison, double %13, i64 0
%39 = shufflevector <2 x double> %38, <2 x double> poison, <2 x i32> zeroinitializer
%40 = insertelement <2 x double> poison, double %15, i64 0
%41 = insertelement <2 x double> %40, double %30, i64 1
%42 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %39, <2 x double> %41, <2 x double> %37)
%43 = extractelement <2 x double> %42, i64 0
%44 = tail call double @llvm.fmuladd.f64(double %32, double %43, double %16)
store double %44, double* %4
%45 = tail call double @llvm.fmuladd.f64(double %32, double %22, double %20)
%46 = tail call double @llvm.fmuladd.f64(double %32, double %28, double %26)
%47 = getelementptr inbounds double, double* %0, i64 1
%48 = load double, double* %47
%49 = getelementptr inbounds double, double* %0, i64 2
%50 = load double, double* %49
%51 = fmul double %14, %50
%52 = tail call double @llvm.fmuladd.f64(double %13, double %48, double %51)
%53 = fmul double %14, %46
%54 = tail call double @llvm.fmuladd.f64(double %13, double %45, double %53)
%55 = tail call double @llvm.fmuladd.f64(double %29, double %52, double %48)
store double %55, double* %47
%56 = tail call double @llvm.fmuladd.f64(double %32, double %52, double %50)
%57 = insertelement <2 x double> poison, double %29, i64 0
%58 = shufflevector <2 x double> %57, <2 x double> poison, <2 x i32> zeroinitializer
%59 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %58, <2 x double> %42, <2 x double> %41)
%60 = bitcast double* %2 to <2 x double>*
store <2 x double> %59, <2 x double>* %60
%61 = extractelement <2 x double> %42, i64 1
%62 = tail call double @llvm.fmuladd.f64(double %32, double %61, double %31)
%63 = tail call double @llvm.fmuladd.f64(double %29, double %54, double %45)
%64 = tail call double @llvm.fmuladd.f64(double %32, double %54, double %46)
store double %64, double* %25
%65 = fneg double %56
store double %65, double* %49
%66 = fneg double %62
store double %66, double* %23
%67 = fneg double %63
store double %67, double* %19
ret void
}