#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 */
int /*<<< orphan*/ free(double *);
double house(double *, int, double *);
scalar_t__ malloc(int);
int /*<<< orphan*/ mmult(double *, double *, double *, int, int, int);
int /*<<< orphan*/ scale(double *, int, int, double);
void hessenberg(double *A, int N) {
  int k, i, j, t, u;
  double *x, *v, *AT, *w;
  double beta;
  x = (double *)malloc(sizeof(double) * N);
  v = (double *)malloc(sizeof(double) * N);
  AT = (double *)malloc(sizeof(double) * N * N);
  w = (double *)malloc(sizeof(double) * N);
  for (k = 0; k < N - 2; ++k) {
    for (i = k + 1; i < N; ++i) {
      x[i - k - 1] = A[i * N + k];
    }
    beta = house(x, N - k - 1, v);
    for (i = k + 1; i < N; i++) {
      t = i * N;
      u = 0;
      for (j = k; j < N; j++) {
        AT[u + i - k - 1] = A[j + t];
        u += (N - k - 1);
      }
    }
    mmult(AT, v, w, N - k, N - k - 1, 1);
    scale(w, N - k, 1, beta);
    mmult(v, w, AT, N - k - 1, 1, N - k);
    for (i = k + 1; i < N; i++) {
      t = i * N;
      for (j = k; j < N; j++) {
        A[t + j] -= AT[(i - k - 1) * (N - k) + j - k];
      }
    }
    for (i = 0; i < N; i++) {
      t = i * N;
      u = i * (N - k - 1);
      for (j = k + 1; j < N; j++) {
        AT[u + j - k - 1] = A[t + j];
      }
    }
    mmult(AT, v, w, N, N - k - 1, 1);
    scale(w, N, 1, beta);
    mmult(w, v, AT, N, 1, N - k - 1);
    for (i = 0; i < N; i++) {
      t = i * N;
      u = i * (N - k - 1);
      for (j = k + 1; j < N; j++) {
        A[t + j] -= AT[u + j - k - 1];
      }
    }
  }
  free(x);
  free(v);
  free(AT);
  free(w);
}


x86 assembly (clang, -O0):
<hessenberg>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
movslq -0xc(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   1050 <malloc@plt>
mov    %rax,-0x28(%rbp)
movslq -0xc(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   1050 <malloc@plt>
mov    %rax,-0x30(%rbp)
movslq -0xc(%rbp),%rax
shl    $0x3,%rax
movslq -0xc(%rbp),%rcx
imul   %rcx,%rax
mov    %eax,%edi
call   1050 <malloc@plt>
mov    %rax,-0x38(%rbp)
movslq -0xc(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   1050 <malloc@plt>
mov    %rax,-0x40(%rbp)
movl   $0x0,-0x10(%rbp)
mov    -0x10(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    $0x2,%ecx
cmp    %ecx,%eax
jge    155d <hessenberg+0x40d>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1219 <hessenberg+0xc9>
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x10(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    -0x10(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    11d4 <hessenberg+0x84>
mov    -0x28(%rbp),%rdi
mov    -0xc(%rbp),%esi
sub    -0x10(%rbp),%esi
sub    $0x1,%esi
mov    -0x30(%rbp),%rdx
call   1070 <house@plt>
movsd  %xmm0,-0x48(%rbp)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12c6 <hessenberg+0x176>
mov    -0x14(%rbp),%eax
imul   -0xc(%rbp),%eax
mov    %eax,-0x1c(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x10(%rbp),%eax
mov    %eax,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12b3 <hessenberg+0x163>
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x38(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    -0x14(%rbp),%ecx
sub    -0x10(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xc(%rbp),%eax
sub    -0x10(%rbp),%eax
sub    $0x1,%eax
add    -0x20(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1260 <hessenberg+0x110>
jmp    12b8 <hessenberg+0x168>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    123d <hessenberg+0xed>
mov    -0x38(%rbp),%rdi
mov    -0x30(%rbp),%rsi
mov    -0x40(%rbp),%rdx
mov    -0xc(%rbp),%ecx
sub    -0x10(%rbp),%ecx
mov    -0xc(%rbp),%r8d
sub    -0x10(%rbp),%r8d
sub    $0x1,%r8d
mov    $0x1,%r9d
call   1040 <mmult@plt>
mov    -0x40(%rbp),%rdi
mov    -0xc(%rbp),%esi
sub    -0x10(%rbp),%esi
movsd  -0x48(%rbp),%xmm0
mov    $0x1,%edx
call   1060 <scale@plt>
mov    -0x30(%rbp),%rdi
mov    -0x40(%rbp),%rsi
mov    -0x38(%rbp),%rdx
mov    -0xc(%rbp),%ecx
sub    -0x10(%rbp),%ecx
sub    $0x1,%ecx
mov    -0xc(%rbp),%r9d
sub    -0x10(%rbp),%r9d
mov    $0x1,%r8d
call   1040 <mmult@plt>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    13c1 <hessenberg+0x271>
mov    -0x14(%rbp),%eax
imul   -0xc(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x10(%rbp),%eax
mov    %eax,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    13ae <hessenberg+0x25e>
mov    -0x38(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    -0x10(%rbp),%ecx
sub    $0x1,%ecx
mov    -0xc(%rbp),%edx
sub    -0x10(%rbp),%edx
imul   %edx,%ecx
add    -0x18(%rbp),%ecx
sub    -0x10(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    -0x18(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1355 <hessenberg+0x205>
jmp    13b3 <hessenberg+0x263>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1339 <hessenberg+0x1e9>
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1450 <hessenberg+0x300>
mov    -0x14(%rbp),%eax
imul   -0xc(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x14(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    -0x10(%rbp),%ecx
sub    $0x1,%ecx
imul   %ecx,%eax
mov    %eax,-0x20(%rbp)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    143d <hessenberg+0x2ed>
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    -0x18(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x38(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    -0x18(%rbp),%ecx
sub    -0x10(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    13f9 <hessenberg+0x2a9>
jmp    1442 <hessenberg+0x2f2>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    13c8 <hessenberg+0x278>
mov    -0x38(%rbp),%rdi
mov    -0x30(%rbp),%rsi
mov    -0x40(%rbp),%rdx
mov    -0xc(%rbp),%ecx
mov    -0xc(%rbp),%r8d
sub    -0x10(%rbp),%r8d
sub    $0x1,%r8d
mov    $0x1,%r9d
call   1040 <mmult@plt>
mov    -0x40(%rbp),%rdi
mov    -0xc(%rbp),%esi
movsd  -0x48(%rbp),%xmm0
mov    $0x1,%edx
call   1060 <scale@plt>
mov    -0x40(%rbp),%rdi
mov    -0x30(%rbp),%rsi
mov    -0x38(%rbp),%rdx
mov    -0xc(%rbp),%ecx
mov    -0xc(%rbp),%r9d
sub    -0x10(%rbp),%r9d
sub    $0x1,%r9d
mov    $0x1,%r8d
call   1040 <mmult@plt>
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    154a <hessenberg+0x3fa>
mov    -0x14(%rbp),%eax
imul   -0xc(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x14(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    -0x10(%rbp),%ecx
sub    $0x1,%ecx
imul   %ecx,%eax
mov    %eax,-0x20(%rbp)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1537 <hessenberg+0x3e7>
mov    -0x38(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    -0x18(%rbp),%ecx
sub    -0x10(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    -0x18(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    14ea <hessenberg+0x39a>
jmp    153c <hessenberg+0x3ec>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    14b9 <hessenberg+0x369>
jmp    154f <hessenberg+0x3ff>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    11ba <hessenberg+0x6a>
mov    -0x28(%rbp),%rdi
call   1030 <free@plt>
mov    -0x30(%rbp),%rdi
call   1030 <free@plt>
mov    -0x38(%rbp),%rdi
call   1030 <free@plt>
mov    -0x40(%rbp),%rdi
call   1030 <free@plt>
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<hessenberg>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x178,%rsp
mov    %esi,%r12d
mov    %rdi,0x18(%rsp)
lea    0x0(,%r12,8),%r14d
mov    %r14d,%edi
call   1050 <malloc@plt>
mov    %rax,%rbp
mov    %r14d,%edi
call   1050 <malloc@plt>
mov    %rax,%r15
mov    %r14d,%edi
imul   %r12d,%edi
call   1050 <malloc@plt>
mov    %rax,%rbx
mov    %r14d,%edi
call   1050 <malloc@plt>
mov    %r12,%r11
mov    %rax,0x60(%rsp)
cmp    $0x3,%r11d
jl     1e06 <hessenberg+0xcb6>
movabs $0xffffffff00000000,%r8
ff ff ff
movabs $0x200000000,%r14
lea    -0x2(%r11),%eax
mov    %rax,0x148(%rsp)
mov    %r11d,%r12d
lea    -0x2(%r12),%rax
mov    %rax,0x140(%rsp)
lea    -0x1(%r11),%eax
mov    %rax,0x138(%rsp)
lea    0x0(,%r12,8),%rdx
lea    0x8(,%r12,8),%rax
mov    %rax,0xd8(%rsp)
lea    -0x1(%r12),%r9
mov    %r12,%rax
neg    %rax
mov    %rax,0x130(%rsp)
mov    0x18(%rsp),%rax
lea    (%rax,%r12,8),%rdi
mov    %r12,%rcx
shl    $0x20,%rcx
lea    (%rax,%r12,8),%rsi
add    $0x10,%rsi
mov    %rsi,0xd0(%rsp)
lea    0x8(%rax,%r12,8),%rsi
mov    %rsi,0xc8(%rsp)
mov    %rcx,0x128(%rsp)
lea    (%rcx,%r8,1),%rcx
mov    %rcx,0xb8(%rsp)
lea    0x38(%rax),%rcx
mov    %rcx,0xf8(%rsp)
lea    -0x5(%r12),%rcx
lea    0x8(%rax),%rsi
mov    %rsi,0xf0(%rsp)
add    $0x18,%rax
mov    %rax,0xe8(%rsp)
mov    $0x1,%eax
mov    %rax,0xc0(%rsp)
movq   $0xffffffffffffffff,0xa0(%rsp)
xor    %eax,%eax
mov    %rax,0x68(%rsp)
mov    %r12,%rsi
mov    %r9,0x120(%rsp)
xor    %r8d,%r8d
mov    $0x1,%eax
mov    %rax,0x8(%rsp)
mov    %r11,0x80(%rsp)
mov    %rbp,0x150(%rsp)
mov    %r15,0x30(%rsp)
mov    %rdx,0x58(%rsp)
jmp    139e <hessenberg+0x24e>
nopl   0x0(%rax,%rax,1)
addq   $0x1,0x8(%rsp)
mov    0xd8(%rsp),%rcx
mov    0x100(%rsp),%rdi
add    %rcx,%rdi
mov    0x50(%rsp),%r9
add    $0xffffffffffffffff,%r9
mov    0xa8(%rsp),%rsi
add    $0xffffffffffffffff,%rsi
movabs $0xffffffff00000000,%rax
ff ff ff
add    %rax,0x128(%rsp)
add    %rcx,0xd0(%rsp)
addq   $0xffffffffffffffff,0x68(%rsp)
addq   $0xffffffffffffffff,0xc0(%rsp)
mov    0x58(%rsp),%rdx
add    %rdx,0xc8(%rsp)
add    %rax,0xb8(%rsp)
addq   $0x8,0xf8(%rsp)
mov    0x170(%rsp),%rcx
add    $0xffffffffffffffff,%rcx
addq   $0xffffffffffffffff,0xa0(%rsp)
addq   $0x8,0xe8(%rsp)
mov    0x48(%rsp),%rax
mov    %rax,%r8
cmp    0x148(%rsp),%rax
mov    0x80(%rsp),%r11
mov    0x150(%rsp),%rbp
mov    0x30(%rsp),%r15
je     1e06 <hessenberg+0xcb6>
mov    %rsi,0xa8(%rsp)
mov    %rcx,0x170(%rsp)
mov    %rcx,%rsi
shr    $0x2,%rsi
mov    0x120(%rsp),%rax
sub    %r8,%rax
add    $0xfffffffffffffffc,%rax
mov    %rax,0x118(%rsp)
mov    %r12,%r10
sub    %r8,%r10
lea    0x1(%r8),%rax
mov    %rax,0x48(%rsp)
cmp    %r12,%rax
jae    1406 <hessenberg+0x2b6>
xor    %eax,%eax
mov    %rdi,%rcx
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rcx),%xmm0
movsd  %xmm0,0x0(%rbp,%rax,8)
add    %rdx,%rcx
add    $0x1,%rax
cmp    %rax,%r9
jne    13f0 <hessenberg+0x2a0>
mov    %r9,0x50(%rsp)
mov    %rdi,0x100(%rsp)
add    $0x1,%rsi
mov    %rsi,0x108(%rsp)
mov    0x118(%rsp),%rax
shr    $0x2,%rax
mov    %rax,0x110(%rsp)
mov    %r8,%rcx
not    %rcx
mov    %rcx,0xb0(%rsp)
lea    (%rcx,%r12,1),%rax
mov    %rax,0x28(%rsp)
mov    0x140(%rsp),%rax
sub    %r8,%rax
mov    %rax,0x78(%rsp)
lea    -0x1(%r10),%rsi
mov    %r8d,%r13d
not    %r13d
add    %r11d,%r13d
mov    %rbp,%rdi
mov    %rsi,0x70(%rsp)
mov    %r15,%rdx
mov    %r10,%rbp
mov    %r8,0xe0(%rsp)
call   1070 <house@plt>
cmp    %r12,0x48(%rsp)
movsd  %xmm0,0x88(%rsp)
jae    1840 <hessenberg+0x6f0>
mov    %r13d,0x10(%rsp)
mov    0xe0(%rsp),%rdi
lea    0x2(%rdi),%rcx
mov    0xa8(%rsp),%rsi
mov    %rsi,%rax
and    $0xfffffffffffffffc,%rax
mov    %rax,0x168(%rsp)
mov    0x50(%rsp),%rax
shl    $0x20,%rax
mov    0xd8(%rsp),%rdx
imul   %rdi,%rdx
mov    %rdx,0x38(%rsp)
imul   %r12,%rcx
mov    %rcx,0x90(%rsp)
mov    %rbp,0x20(%rsp)
lea    (%rbx,%rbp,8),%rcx
mov    %rcx,0x98(%rsp)
mov    0x80(%rsp),%rcx
mov    %rdi,%r15
sub    %edi,%ecx
mov    %ecx,0x40(%rsp)
mov    0x100(%rsp),%rcx
xor    %r8d,%r8d
mov    0x8(%rsp),%r11
mov    0x58(%rsp),%r9
movabs $0x100000000,%r10
cs nopw 0x0(%rax,%rax,1)
nop
xor    %edi,%edi
mov    %r8,%rbp
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rcx,%rdi,8),%xmm0
mov    %rbp,%rdx
sar    $0x1d,%rdx
movsd  %xmm0,(%rbx,%rdx,1)
add    %rax,%rbp
add    $0x1,%rdi
cmp    %rdi,%rsi
jne    1540 <hessenberg+0x3f0>
add    $0x1,%r11
add    %r10,%r8
add    %r9,%rcx
cmp    %r12,%r11
jne    1530 <hessenberg+0x3e0>
mov    %rbx,%rdi
mov    0x30(%rsp),%rsi
mov    0x60(%rsp),%rbp
mov    %rbp,%rdx
mov    0x20(%rsp),%r13
mov    %r13d,%ecx
mov    0x10(%rsp),%r8d
mov    $0x1,%r9d
call   1040 <mmult@plt>
mov    %rbp,%rdi
mov    %r13d,%esi
mov    $0x1,%edx
movsd  0x88(%rsp),%xmm0
call   1060 <scale@plt>
mov    0x30(%rsp),%rdi
mov    %rbp,%rsi
mov    %rbx,%rdx
mov    0x10(%rsp),%ecx
mov    $0x1,%r8d
mov    %r13d,%r9d
call   1040 <mmult@plt>
cmp    %r12,0x48(%rsp)
mov    %r15,%r8
jae    189a <hessenberg+0x74a>
mov    0x20(%rsp),%rax
and    $0xfffffffffffffffc,%rax
mov    %rax,0x160(%rsp)
add    %r8,%rax
mov    %rax,0x158(%rsp)
mov    0xc8(%rsp),%rdx
mov    0xc0(%rsp),%r11
mov    0x68(%rsp),%r9
mov    0xd0(%rsp),%rbp
xor    %eax,%eax
mov    %rax,0x10(%rsp)
xor    %r13d,%r13d
mov    0x8(%rsp),%r10
jmp    166f <hessenberg+0x51f>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%r10
add    $0x1,%r13
mov    0x10(%rsp),%rax
add    0x128(%rsp),%rax
mov    %rax,0x10(%rsp)
mov    0x58(%rsp),%rax
add    %rax,%rbp
mov    0xa8(%rsp),%rcx
add    %rcx,%r9
add    %rcx,%r11
add    %rax,%rdx
cmp    %r12,%r10
mov    %r15,%r8
je     189a <hessenberg+0x74a>
cmpq   $0x4,0x20(%rsp)
jae    1740 <hessenberg+0x5f0>
mov    %r12d,%eax
sub    %r8d,%eax
mov    %r8,%rsi
test   $0x1,%al
je     16c5 <hessenberg+0x575>
mov    %r10,%rax
imul   %r12,%rax
mov    0xb0(%rsp),%rcx
lea    (%r10,%rcx,1),%rsi
imul   0x20(%rsp),%rsi
sub    %r15,%rsi
add    %r8d,%esi
movslq %esi,%rsi
add    %r8,%rax
mov    0x18(%rsp),%rcx
movsd  (%rcx,%rax,8),%xmm0
subsd  (%rbx,%rsi,8),%xmm0
movsd  %xmm0,(%rcx,%rax,8)
lea    0x1(%r8),%rsi
cmp    %r8,0x120(%rsp)
je     1630 <hessenberg+0x4e0>
lea    (%rsi,%r9,1),%ecx
shl    $0x20,%rcx
lea    (%r11,%rsi,1),%eax
shl    $0x20,%rax
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %rcx,%rdi
sar    $0x1d,%rdi
movsd  -0x8(%rdx,%rsi,8),%xmm0
movsd  (%rdx,%rsi,8),%xmm1
subsd  (%rbx,%rdi,1),%xmm0
movsd  %xmm0,-0x8(%rdx,%rsi,8)
mov    %rax,%rdi
sar    $0x1d,%rdi
subsd  (%rbx,%rdi,1),%xmm1
movsd  %xmm1,(%rdx,%rsi,8)
add    $0x2,%rsi
add    %r14,%rcx
add    %r14,%rax
cmp    %rsi,%r12
jne    16f0 <hessenberg+0x5a0>
jmp    1630 <hessenberg+0x4e0>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    0x40(%rsp),%eax
imul   %r13d,%eax
mov    0x28(%rsp),%rcx
lea    (%rax,%rcx,1),%esi
mov    %r15,%r8
cmp    %eax,%esi
jl     167b <hessenberg+0x52b>
mov    0x28(%rsp),%rax
shr    $0x20,%rax
mov    %r15,%r8
jne    167b <hessenberg+0x52b>
mov    %r13,%rax
imul   %r12,%rax
lea    (%rax,%r12,1),%rsi
mov    0x18(%rsp),%rcx
lea    (%rcx,%rsi,8),%rdi
add    0x38(%rsp),%rdi
mov    0x20(%rsp),%rcx
mov    %ecx,%esi
imul   %r13d,%esi
movslq %esi,%r8
mov    0x98(%rsp),%rcx
lea    (%rcx,%r8,8),%rsi
cmp    %rsi,%rdi
jae    17c6 <hessenberg+0x676>
add    0x90(%rsp),%rax
mov    0x18(%rsp),%rcx
lea    (%rcx,%rax,8),%rax
lea    (%rbx,%r8,8),%rsi
mov    %r15,%r8
cmp    %rax,%rsi
jb     167b <hessenberg+0x52b>
mov    0x10(%rsp),%rax
xor    %esi,%esi
movabs $0x400000000,%r8
mov    0x168(%rsp),%rcx
nop
mov    %rax,%rdi
sar    $0x1d,%rdi
movupd (%rbx,%rdi,1),%xmm0
movupd 0x10(%rbx,%rdi,1),%xmm1
movupd -0x10(%rbp,%rsi,8),%xmm2
subpd  %xmm0,%xmm2
movupd 0x0(%rbp,%rsi,8),%xmm0
subpd  %xmm1,%xmm0
movupd %xmm2,-0x10(%rbp,%rsi,8)
movupd %xmm0,0x0(%rbp,%rsi,8)
add    $0x4,%rsi
add    %r8,%rax
cmp    %rsi,%rcx
jne    17e0 <hessenberg+0x690>
mov    0x158(%rsp),%r8
mov    0x160(%rsp),%rax
cmp    %rax,0x20(%rsp)
je     1630 <hessenberg+0x4e0>
jmp    167b <hessenberg+0x52b>
xchg   %ax,%ax
mov    %rbx,%rdi
mov    %r15,%rsi
mov    0x60(%rsp),%r15
mov    %r15,%rdx
mov    %ebp,%ecx
mov    %r13d,%r8d
mov    $0x1,%r9d
call   1040 <mmult@plt>
mov    %r15,%rdi
mov    %ebp,%esi
mov    $0x1,%edx
movsd  0x88(%rsp),%xmm0
call   1060 <scale@plt>
mov    0x30(%rsp),%rdi
mov    %r15,%rsi
mov    %rbx,%rdx
mov    %r13d,%ecx
mov    $0x1,%r8d
mov    %ebp,%r9d
call   1040 <mmult@plt>
mov    0xe0(%rsp),%r8
mov    0x50(%rsp),%rax
and    $0xfffffffffffffffc,%rax
mov    %rax,0x98(%rsp)
andq   $0xfffffffffffffffe,0x108(%rsp)
addq   $0x1,0x110(%rsp)
mov    0x138(%rsp),%rax
sub    %r8,%rax
mov    %rax,0x10(%rsp)
mov    0x28(%rsp),%rax
lea    (%rbx,%rax,8),%rcx
mov    %rcx,0x20(%rsp)
mov    %rax,%rcx
and    $0xfffffffffffffffc,%rcx
mov    0x8(%rsp),%rax
mov    %rcx,0x40(%rsp)
add    %rcx,%rax
mov    %rax,0x90(%rsp)
mov    0x68(%rsp),%r13
mov    0xf0(%rsp),%rdx
mov    0xa0(%rsp),%r10
mov    0xf8(%rsp),%rdi
xor    %r11d,%r11d
xor    %r9d,%r9d
jmp    194b <hessenberg+0x7fb>
nopl   0x0(%rax)
add    $0x1,%r9
add    0xb8(%rsp),%r11
mov    0x58(%rsp),%rax
add    %rax,%rdi
mov    0x50(%rsp),%rcx
add    %rcx,%r10
add    %rax,%rdx
add    %rcx,%r13
cmp    %r12,%r9
je     1b70 <hessenberg+0xa20>
cmp    %r12,0x48(%rsp)
jae    1920 <hessenberg+0x7d0>
mov    %r9,%rbp
imul   %r12,%rbp
mov    %r9,%r15
imul   0x70(%rsp),%r15
add    0xb0(%rsp),%r15
mov    0x8(%rsp),%r8
cmpq   $0x4,0x28(%rsp)
jae    1a10 <hessenberg+0x8c0>
mov    %r12d,%ecx
sub    %r8d,%ecx
mov    %r8,%rax
test   $0x1,%cl
je     19a5 <hessenberg+0x855>
add    %r8,%rbp
mov    0x18(%rsp),%rax
movsd  (%rax,%rbp,8),%xmm0
add    %r8d,%r15d
movslq %r15d,%rax
movsd  %xmm0,(%rbx,%rax,8)
lea    0x1(%r8),%rax
not    %r8
cmp    0x130(%rsp),%r8
je     1920 <hessenberg+0x7d0>
lea    (%rax,%r10,1),%ebp
shl    $0x20,%rbp
lea    (%rax,%r13,1),%esi
shl    $0x20,%rsi
cs nopw 0x0(%rax,%rax,1)
movsd  -0x8(%rdx,%rax,8),%xmm0
mov    %rbp,%rcx
sar    $0x1d,%rcx
movsd  %xmm0,(%rbx,%rcx,1)
movsd  (%rdx,%rax,8),%xmm0
mov    %rsi,%rcx
sar    $0x1d,%rcx
movsd  %xmm0,(%rbx,%rcx,1)
add    $0x2,%rax
add    %r14,%rbp
add    %r14,%rsi
cmp    %rax,%r12
jne    19d0 <hessenberg+0x880>
jmp    1920 <hessenberg+0x7d0>
nopw   0x0(%rax,%rax,1)
mov    0x10(%rsp),%rax
imul   %r9d,%eax
mov    0x78(%rsp),%rcx
add    %eax,%ecx
mov    0x8(%rsp),%r8
cmp    %eax,%ecx
jl     197b <hessenberg+0x82b>
mov    0x78(%rsp),%rax
shr    $0x20,%rax
mov    0x8(%rsp),%r8
jne    197b <hessenberg+0x82b>
mov    0x8(%rsp),%rax
add    %rbp,%rax
mov    0x10(%rsp),%rcx
imul   %r9d,%ecx
mov    0x18(%rsp),%rsi
lea    (%rsi,%rax,8),%rsi
movslq %ecx,%rax
mov    0x20(%rsp),%rcx
lea    (%rcx,%rax,8),%rcx
cmp    %rcx,%rsi
jae    1a8b <hessenberg+0x93b>
lea    (%r12,%rbp,1),%rcx
mov    0x18(%rsp),%rsi
lea    (%rsi,%rcx,8),%rcx
lea    (%rbx,%rax,8),%rax
mov    0x8(%rsp),%r8
cmp    %rcx,%rax
jb     197b <hessenberg+0x82b>
cmpq   $0x4,0x118(%rsp)
mov    %rbp,0x38(%rsp)
jae    1a9f <hessenberg+0x94f>
xor    %eax,%eax
jmp    1b0e <hessenberg+0x9be>
mov    0x108(%rsp),%rcx
mov    %r11,%rbp
xor    %eax,%eax
movabs $0x800000000,%r8
cs nopw 0x0(%rax,%rax,1)
movups -0x30(%rdi,%rax,8),%xmm0
movups -0x20(%rdi,%rax,8),%xmm1
mov    %rbp,%rsi
sar    $0x1d,%rsi
movups %xmm0,(%rbx,%rsi,1)
movups %xmm1,0x10(%rbx,%rsi,1)
movupd -0x10(%rdi,%rax,8),%xmm0
movupd (%rdi,%rax,8),%xmm1
movabs $0x400000000,%rsi
add    %rbp,%rsi
sar    $0x1d,%rsi
movupd %xmm0,(%rbx,%rsi,1)
movupd %xmm1,0x10(%rbx,%rsi,1)
add    $0x8,%rax
add    %r8,%rbp
add    $0xfffffffffffffffe,%rcx
jne    1ac0 <hessenberg+0x970>
testb  $0x1,0x110(%rsp)
01
mov    0x38(%rsp),%rbp
je     1b46 <hessenberg+0x9f6>
add    0x8(%rsp),%rax
lea    (%rax,%rbp,1),%rcx
mov    0x18(%rsp),%rsi
movupd (%rsi,%rcx,8),%xmm0
movupd 0x10(%rsi,%rcx,8),%xmm1
add    %r15d,%eax
cltq
movupd %xmm0,(%rbx,%rax,8)
movupd %xmm1,0x10(%rbx,%rax,8)
mov    0x90(%rsp),%r8
mov    0x40(%rsp),%rax
cmp    %rax,0x28(%rsp)
je     1920 <hessenberg+0x7d0>
jmp    197b <hessenberg+0x82b>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %rbx,%rdi
mov    0x30(%rsp),%r15
mov    %r15,%rsi
mov    0x60(%rsp),%rbp
mov    %rbp,%rdx
mov    0x80(%rsp),%r13
mov    %r13d,%ecx
mov    0x70(%rsp),%r8
mov    $0x1,%r9d
call   1040 <mmult@plt>
mov    %rbp,%rdi
mov    %r13d,%esi
mov    $0x1,%edx
movsd  0x88(%rsp),%xmm0
call   1060 <scale@plt>
mov    %rbp,%rdi
mov    %r15,%rsi
mov    %rbx,%rdx
mov    %r13d,%ecx
mov    $0x1,%r8d
mov    0x70(%rsp),%r9
call   1040 <mmult@plt>
mov    0x8(%rsp),%rax
mov    0x40(%rsp),%rcx
add    %rcx,%rax
mov    %rax,0x38(%rsp)
mov    0xf0(%rsp),%rcx
mov    0x68(%rsp),%r11
mov    0xa0(%rsp),%r13
mov    0xe8(%rsp),%rdi
xor    %r10d,%r10d
xor    %r15d,%r15d
jmp    1c3b <hessenberg+0xaeb>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r15
add    0xb8(%rsp),%r10
mov    0x58(%rsp),%rax
add    %rax,%rdi
mov    0x50(%rsp),%rdx
add    %rdx,%r13
add    %rdx,%r11
add    %rax,%rcx
cmp    %r12,%r15
je     12e0 <hessenberg+0x190>
cmp    %r12,0x48(%rsp)
jae    1c10 <hessenberg+0xac0>
mov    %r15,%r9
imul   %r12,%r9
mov    0x8(%rsp),%rdx
cmpq   $0x4,0x28(%rsp)
jae    1d10 <hessenberg+0xbc0>
mov    %r12d,%esi
sub    %edx,%esi
mov    %rdx,%rbp
test   $0x1,%sil
je     1c9b <hessenberg+0xb4b>
mov    %r15,%rsi
imul   0x70(%rsp),%rsi
add    0xb0(%rsp),%rsi
add    %edx,%esi
movslq %esi,%rsi
add    %rdx,%r9
mov    0x18(%rsp),%rax
movsd  (%rax,%r9,8),%xmm0
subsd  (%rbx,%rsi,8),%xmm0
movsd  %xmm0,(%rax,%r9,8)
lea    0x1(%rdx),%rbp
not    %rdx
cmp    0x130(%rsp),%rdx
je     1c10 <hessenberg+0xac0>
lea    0x0(%rbp,%r13,1),%eax
shl    $0x20,%rax
lea    (%r11,%rbp,1),%esi
shl    $0x20,%rsi
nopl   (%rax)
mov    %rax,%rdx
sar    $0x1d,%rdx
movsd  -0x8(%rcx,%rbp,8),%xmm0
movsd  (%rcx,%rbp,8),%xmm1
subsd  (%rbx,%rdx,1),%xmm0
movsd  %xmm0,-0x8(%rcx,%rbp,8)
mov    %rsi,%rdx
sar    $0x1d,%rdx
subsd  (%rbx,%rdx,1),%xmm1
movsd  %xmm1,(%rcx,%rbp,8)
add    $0x2,%rbp
add    %r14,%rax
add    %r14,%rsi
cmp    %rbp,%r12
jne    1cc0 <hessenberg+0xb70>
jmp    1c10 <hessenberg+0xac0>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    0x10(%rsp),%rsi
imul   %r15d,%esi
mov    0x78(%rsp),%rax
lea    (%rsi,%rax,1),%ebp
mov    0x8(%rsp),%rdx
cmp    %esi,%ebp
jl     1c5a <hessenberg+0xb0a>
mov    0x78(%rsp),%rdx
shr    $0x20,%rdx
mov    0x8(%rsp),%rdx
jne    1c5a <hessenberg+0xb0a>
mov    0x8(%rsp),%rax
lea    (%rax,%r9,1),%rdx
mov    0x10(%rsp),%rax
mov    %eax,%esi
imul   %r15d,%esi
mov    0x18(%rsp),%rax
lea    (%rax,%rdx,8),%rbp
movslq %esi,%rdx
mov    0x20(%rsp),%rax
lea    (%rax,%rdx,8),%rsi
cmp    %rsi,%rbp
jae    1d8f <hessenberg+0xc3f>
lea    (%r9,%r12,1),%rsi
mov    0x18(%rsp),%rax
lea    (%rax,%rsi,8),%rsi
lea    (%rbx,%rdx,8),%rbp
mov    0x8(%rsp),%rdx
cmp    %rsi,%rbp
jb     1c5a <hessenberg+0xb0a>
mov    %r10,%rbp
xor    %esi,%esi
movabs $0x400000000,%rax
mov    0x98(%rsp),%r8
cs nopw 0x0(%rax,%rax,1)
mov    %rbp,%rdx
sar    $0x1d,%rdx
movupd (%rbx,%rdx,1),%xmm0
movupd 0x10(%rbx,%rdx,1),%xmm1
movupd -0x10(%rdi,%rsi,8),%xmm2
subpd  %xmm0,%xmm2
movupd (%rdi,%rsi,8),%xmm0
subpd  %xmm1,%xmm0
movupd %xmm2,-0x10(%rdi,%rsi,8)
movupd %xmm0,(%rdi,%rsi,8)
add    $0x4,%rsi
add    %rax,%rbp
cmp    %rsi,%r8
jne    1db0 <hessenberg+0xc60>
mov    0x38(%rsp),%rdx
mov    0x40(%rsp),%rax
cmp    %rax,0x28(%rsp)
je     1c10 <hessenberg+0xac0>
jmp    1c5a <hessenberg+0xb0a>
mov    %rbp,%rdi
call   1030 <free@plt>
mov    %r15,%rdi
call   1030 <free@plt>
mov    %rbx,%rdi
call   1030 <free@plt>
mov    0x60(%rsp),%rdi
add    $0x178,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

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

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