#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 *);
scalar_t__ malloc(int);
int /*<<< orphan*/ mmult(double *, double *, double *, int, int, int);
int /*<<< orphan*/ scale(double *, int, int, double);
void getQR_house(double *A, int M, int N, double *bvec, double *Q, double *R) {
  int i, j, k, t, u;
  double *x, *v, *AT, *w;
  x = (double *)malloc(sizeof(double) * M);
  v = (double *)malloc(sizeof(double) * M);
  AT = (double *)malloc(sizeof(double) * M * N);
  w = (double *)malloc(sizeof(double) * M * M);
  for (i = 0; i < N; ++i) {
    t = i * N;
    for (j = 0; j < N; ++j) {
      if (i > j) {
        R[t + j] = 0.;
      } else {
        R[t + j] = A[t + j];
      }
    }
  }
  for (i = 0; i < M; ++i) {
    t = i * N;
    for (j = 0; j < N; ++j) {
      if (i == j) {
        Q[t + j] = 1.;
      } else {
        Q[t + j] = 0.;
      }
    }
  }
  for (j = N - 1; j >= 0; --j) {
    v[0] = 1.;
    for (i = j + 1; i < M; ++i) {
      v[i - j] = A[i * N + j];
    }
    for (i = j; i < M; i++) {
      t = i * N;
      u = 0;
      for (k = j; k < N; k++) {
        AT[u + i - j] = Q[k + t];
        u += (M - j);
      }
    }
    mmult(AT, v, w, N - j, M - j, 1);
    scale(w, N - j, 1, bvec[j]);
    mmult(v, w, AT, M - j, 1, N - j);
    for (i = j; i < M; i++) {
      t = i * N;
      for (k = j; k < N; k++) {
        Q[t + k] -= AT[(i - j) * (N - j) + k - j];
      }
    }
  }
  free(x);
  free(v);
  free(AT);
  free(w);
}


x86 assembly (clang, -O0):
<getQR_house>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
movslq -0xc(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   1050 <malloc@plt>
mov    %rax,-0x48(%rbp)
movslq -0xc(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   1050 <malloc@plt>
mov    %rax,-0x50(%rbp)
movslq -0xc(%rbp),%rax
shl    $0x3,%rax
movslq -0x10(%rbp),%rcx
imul   %rcx,%rax
mov    %eax,%edi
call   1050 <malloc@plt>
mov    %rax,-0x58(%rbp)
movslq -0xc(%rbp),%rax
shl    $0x3,%rax
movslq -0xc(%rbp),%rcx
imul   %rcx,%rax
mov    %eax,%edi
call   1050 <malloc@plt>
mov    %rax,-0x60(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    125a <getQR_house+0x11a>
mov    -0x2c(%rbp),%eax
imul   -0x10(%rbp),%eax
mov    %eax,-0x38(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1247 <getQR_house+0x107>
mov    -0x2c(%rbp),%eax
cmp    -0x30(%rbp),%eax
jle    1210 <getQR_house+0xd0>
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1234 <getQR_house+0xf4>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1239 <getQR_house+0xf9>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    11de <getQR_house+0x9e>
jmp    124c <getQR_house+0x10c>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11c1 <getQR_house+0x81>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12f0 <getQR_house+0x1b0>
mov    -0x2c(%rbp),%eax
imul   -0x10(%rbp),%eax
mov    %eax,-0x38(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    12dd <getQR_house+0x19d>
mov    -0x2c(%rbp),%eax
cmp    -0x30(%rbp),%eax
jne    12b5 <getQR_house+0x175>
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  0xd55(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    12ca <getQR_house+0x18a>
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    12cf <getQR_house+0x18f>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    127e <getQR_house+0x13e>
jmp    12e2 <getQR_house+0x1a2>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1261 <getQR_house+0x121>
mov    -0x10(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x30(%rbp)
cmpl   $0x0,-0x30(%rbp)
jl     14f0 <getQR_house+0x3b0>
mov    -0x50(%rbp),%rax
movsd  0xcf1(%rip),%xmm0
movsd  %xmm0,(%rax)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    135e <getQR_house+0x21e>
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x10(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x50(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    131c <getQR_house+0x1dc>
mov    -0x30(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    13e7 <getQR_house+0x2a7>
mov    -0x2c(%rbp),%eax
imul   -0x10(%rbp),%eax
mov    %eax,-0x38(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x30(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    13d4 <getQR_house+0x294>
mov    -0x20(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x58(%rbp),%rax
mov    -0x3c(%rbp),%ecx
add    -0x2c(%rbp),%ecx
sub    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xc(%rbp),%eax
sub    -0x30(%rbp),%eax
add    -0x3c(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1387 <getQR_house+0x247>
jmp    13d9 <getQR_house+0x299>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1364 <getQR_house+0x224>
mov    -0x58(%rbp),%rdi
mov    -0x50(%rbp),%rsi
mov    -0x60(%rbp),%rdx
mov    -0x10(%rbp),%ecx
sub    -0x30(%rbp),%ecx
mov    -0xc(%rbp),%r8d
sub    -0x30(%rbp),%r8d
mov    $0x1,%r9d
call   1040 <mmult@plt>
mov    -0x60(%rbp),%rdi
mov    -0x10(%rbp),%esi
sub    -0x30(%rbp),%esi
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    $0x1,%edx
call   1060 <scale@plt>
mov    -0x50(%rbp),%rdi
mov    -0x60(%rbp),%rsi
mov    -0x58(%rbp),%rdx
mov    -0xc(%rbp),%ecx
sub    -0x30(%rbp),%ecx
mov    -0x10(%rbp),%r9d
sub    -0x30(%rbp),%r9d
mov    $0x1,%r8d
call   1040 <mmult@plt>
mov    -0x30(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    14dd <getQR_house+0x39d>
mov    -0x2c(%rbp),%eax
imul   -0x10(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x30(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    14ca <getQR_house+0x38a>
mov    -0x58(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0x30(%rbp),%ecx
mov    -0x10(%rbp),%edx
sub    -0x30(%rbp),%edx
imul   %edx,%ecx
add    -0x34(%rbp),%ecx
sub    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1474 <getQR_house+0x334>
jmp    14cf <getQR_house+0x38f>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1458 <getQR_house+0x318>
jmp    14e2 <getQR_house+0x3a2>
mov    -0x30(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x30(%rbp)
jmp    12f9 <getQR_house+0x1b9>
mov    -0x48(%rbp),%rdi
call   1030 <free@plt>
mov    -0x50(%rbp),%rdi
call   1030 <free@plt>
mov    -0x58(%rbp),%rdi
call   1030 <free@plt>
mov    -0x60(%rbp),%rdi
call   1030 <free@plt>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<getQR_house>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x108,%rsp
mov    %r9,%rbp
mov    %r8,0x18(%rsp)
mov    %rcx,0xb8(%rsp)
mov    %edx,%r14d
mov    %esi,%ebx
mov    %rdi,%r15
lea    0x0(,%rbx,8),%r12d
mov    %r12d,%edi
call   1050 <malloc@plt>
mov    %rax,0xa8(%rsp)
mov    %r12d,%edi
call   1050 <malloc@plt>
mov    %rax,0x78(%rsp)
mov    %r12d,%edi
imul   %r14d,%edi
call   1050 <malloc@plt>
mov    %rax,%r13
mov    %rbx,0x30(%rsp)
imul   %ebx,%r12d
mov    %r12d,%edi
call   1050 <malloc@plt>
mov    %rax,0x70(%rsp)
mov    %r14d,%eax
mov    %rax,(%rsp)
test   %r14d,%r14d
jle    12a3 <getQR_house+0x163>
mov    %rbp,%r9
mov    (%rsp),%rdx
mov    %edx,%eax
and    $0xfffffffe,%eax
lea    0x8(%rbp),%rcx
lea    0x0(,%rdx,8),%r8
lea    0x8(%r15),%rsi
xor    %edi,%edi
jmp    120d <getQR_house+0xcd>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  %xmm0,(%r9,%rbx,8)
mov    (%rsp),%rdx
add    $0x1,%rdi
add    %r8,%rcx
add    %r8,%rsi
cmp    %rdx,%rdi
je     12a3 <getQR_house+0x163>
cmp    $0x1,%r14d
jne    1220 <getQR_house+0xe0>
xor    %ebp,%ebp
jmp    1274 <getQR_house+0x134>
nopw   0x0(%rax,%rax,1)
xor    %ebx,%ebx
jmp    123e <getQR_house+0xfe>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  %xmm0,(%rcx,%rbp,8)
add    $0x1,%rbx
cmp    %rax,%rbx
je     1270 <getQR_house+0x130>
mov    %rbx,%rbp
xorpd  %xmm0,%xmm0
xorpd  %xmm1,%xmm1
cmp    %rbx,%rdi
ja     1254 <getQR_house+0x114>
movsd  -0x8(%rsi,%rbp,8),%xmm1
movsd  %xmm1,-0x8(%rcx,%rbp,8)
lea    0x1(%rbp),%rbx
cmp    %rbx,%rdi
ja     1230 <getQR_house+0xf0>
movsd  (%rsi,%rbp,8),%xmm0
jmp    1230 <getQR_house+0xf0>
nopw   0x0(%rax,%rax,1)
add    $0x2,%rbp
mov    (%rsp),%rdx
test   $0x1,%dl
je     11fa <getQR_house+0xba>
mov    %rdi,%rbx
imul   %rdx,%rbx
add    %rbp,%rbx
xorpd  %xmm0,%xmm0
cmp    %rbp,%rdi
ja     11f0 <getQR_house+0xb0>
movsd  (%r15,%rbx,8),%xmm0
jmp    11f0 <getQR_house+0xb0>
movabs $0x3ff0000000000000,%r11
cmpl   $0x0,0x30(%rsp)
mov    0x18(%rsp),%r12
jle    1397 <getQR_house+0x257>
movslq %r14d,%r8
mov    0x30(%rsp),%r9d
mov    (%rsp),%rdx
and    $0xfffffffe,%edx
lea    0x8(%r12),%rsi
lea    0x0(,%r8,8),%r10
xor    %eax,%eax
jmp    1300 <getQR_house+0x1c0>
mov    %r11,(%r12,%rdi,8)
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
add    $0x1,%rax
add    %r10,%rsi
cmp    %r9,%rax
je     1397 <getQR_house+0x257>
test   %r14d,%r14d
jle    12f0 <getQR_house+0x1b0>
mov    %rax,%rcx
imul   %r8,%rcx
lea    (%rax,%rcx,1),%rdi
cmp    $0x1,%r14d
jne    1320 <getQR_house+0x1e0>
xor    %ebp,%ebp
testb  $0x1,(%rsp)
je     12f0 <getQR_house+0x1b0>
jmp    137e <getQR_house+0x23e>
xor    %ebx,%ebx
jmp    133d <getQR_house+0x1fd>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    %r11,(%r12,%rdi,8)
add    $0x1,%rbx
cmp    %rdx,%rbx
je     1370 <getQR_house+0x230>
mov    %rbx,%rbp
cmp    %rbx,%rax
jne    1350 <getQR_house+0x210>
mov    %r11,(%r12,%rdi,8)
jmp    1359 <getQR_house+0x219>
nopl   0x0(%rax,%rax,1)
movq   $0x0,-0x8(%rsi,%rbp,8)
lea    0x1(%rbp),%rbx
cmp    %rbx,%rax
je     1330 <getQR_house+0x1f0>
movq   $0x0,(%rsi,%rbp,8)
jmp    1334 <getQR_house+0x1f4>
nopl   0x0(%rax)
add    $0x2,%rbp
testb  $0x1,(%rsp)
je     12f0 <getQR_house+0x1b0>
cmp    %rbp,%rax
je     12dd <getQR_house+0x19d>
add    %rcx,%rbp
movq   $0x0,(%r12,%rbp,8)
jmp    12f0 <getQR_house+0x1b0>
test   %r14d,%r14d
jle    19cd <getQR_house+0x88d>
mov    0x30(%rsp),%rax
movslq %eax,%rbp
sub    %r14d,%eax
mov    (%rsp),%rdx
mov    %rdx,%rsi
imul   %rdx,%rsi
lea    -0x1(%rsi),%rcx
mov    %rcx,0xe8(%rsp)
lea    0x0(,%rdx,8),%rdi
mov    %rdi,%rcx
xor    $0xfffffffffffffff8,%rcx
mov    %rcx,0x98(%rsp)
mov    %rdx,%rcx
neg    %rcx
mov    %rcx,0x100(%rsp)
lea    0x1(%rdx),%rcx
imul   %rdx,%rcx
lea    (%r15,%rcx,8),%rbx
add    $0xfffffffffffffff8,%rbx
lea    (%rax,%rax,1),%r15d
add    $0x2,%r15d
add    $0x1,%eax
mov    %eax,0x2c(%rsp)
mov    $0x1,%eax
mov    %rax,0x38(%rsp)
mov    $0x1,%eax
sub    %rdx,%rax
mov    %rax,0x50(%rsp)
lea    (%r12,%rsi,8),%rax
mov    %rax,0x90(%rsp)
mov    %rsi,0xb0(%rsp)
lea    (%r12,%rsi,8),%rax
add    $0x8,%rax
mov    %rax,0x88(%rsp)
mov    $0x2,%eax
sub    %r14d,%eax
mov    %eax,0x28(%rsp)
lea    -0x1(%rdx),%rax
imul   %rdx,%rax
add    $0x1,%rax
mov    %rax,0x80(%rsp)
xor    %eax,%eax
mov    %rax,0x10(%rsp)
mov    %rdx,0x8(%rsp)
mov    %rbp,0x40(%rsp)
mov    %rdi,0x60(%rsp)
jmp    14e5 <getQR_house+0x3a5>
nopl   0x0(%rax)
addl   $0x1,0x2c(%rsp)
mov    0x98(%rsp),%rcx
mov    0xe0(%rsp),%rbx
add    %rcx,%rbx
add    $0x2,%r15d
addq   $0x1,0x50(%rsp)
add    %rcx,0x90(%rsp)
addq   $0x1,0x38(%rsp)
add    %rcx,0x88(%rsp)
addl   $0x1,0x28(%rsp)
sub    %rax,0x80(%rsp)
cmpq   $0x2,0x48(%rsp)
mov    0x18(%rsp),%r12
movabs $0x3ff0000000000000,%r11
mov    0x60(%rsp),%rdi
jl     19cd <getQR_house+0x88d>
mov    0x8(%rsp),%rsi
mov    0x10(%rsp),%rax
mov    %rax,0x20(%rsp)
mov    0x78(%rsp),%r14
mov    %r11,(%r14)
cmp    %rbp,%rsi
jge    152f <getQR_house+0x3ef>
mov    %rbx,%rax
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
movsd  (%rax),%xmm0
movsd  %xmm0,0x8(%r14,%rcx,8)
add    %rdi,%rax
lea    (%rsi,%rcx,1),%rdx
add    $0x1,%rdx
add    $0x1,%rcx
cmp    %rbp,%rdx
jl     1510 <getQR_house+0x3d0>
mov    %rbx,0xe0(%rsp)
mov    0x20(%rsp),%rax
add    $0x1,%rax
mov    %rax,0x10(%rsp)
lea    -0x1(%rsi),%rax
mov    %rax,0x8(%rsp)
mov    $0x1,%eax
sub    %esi,%eax
mov    %rax,0xa0(%rsp)
mov    %rsi,0x48(%rsp)
cmp    %rbp,%rsi
jle    18d0 <getQR_house+0x790>
mov    (%rsp),%rbx
mov    0x8(%rsp),%r12
sub    %r12,%rbx
mov    0x30(%rsp),%rax
mov    %eax,%r8d
sub    %r12d,%r8d
mov    %r8d,0x58(%rsp)
mov    %r13,%rdi
mov    %r14,%rsi
mov    0x70(%rsp),%rbp
mov    %rbp,%rdx
mov    %ebx,%ecx
mov    $0x1,%r9d
call   1040 <mmult@plt>
mov    0xb8(%rsp),%rax
movsd  (%rax,%r12,8),%xmm0
mov    %rbp,%rdi
mov    %ebx,%esi
mov    $0x1,%edx
call   1060 <scale@plt>
mov    %r14,%rdi
mov    %rbp,%rsi
mov    0x40(%rsp),%rbp
mov    %r13,%rdx
mov    0x58(%rsp),%ecx
mov    $0x1,%r8d
mov    %rbx,0x58(%rsp)
mov    %ebx,%r9d
call   1040 <mmult@plt>
mov    (%rsp),%rax
cmp    %rbp,0x48(%rsp)
jg     1480 <getQR_house+0x340>
mov    0x20(%rsp),%rdx
mov    %rdx,%rcx
imul   %rax,%rcx
mov    0x38(%rsp),%rsi
and    $0xfffffffffffffffc,%rsi
mov    %rsi,0xd0(%rsp)
mov    0x98(%rsp),%rsi
imul   %rdx,%rsi
mov    %rsi,0xf8(%rsp)
mov    0xb0(%rsp),%rsi
sub    %rcx,%rsi
mov    %rsi,0xd8(%rsp)
lea    0x8(,%rdx,8),%rcx
add    %r13,%rcx
mov    %rcx,0xf0(%rsp)
lea    0x1(%rdx),%ecx
mov    %ecx,0x6c(%rsp)
mov    0x10(%rsp),%rcx
and    $0xfffffffffffffffc,%rcx
mov    0x8(%rsp),%r11
mov    %rcx,0xc8(%rsp)
add    %r11,%rcx
mov    %rcx,0xc0(%rsp)
mov    0x80(%rsp),%rsi
mov    0x50(%rsp),%rdx
mov    0x28(%rsp),%ecx
mov    %ecx,%ebx
mov    0x88(%rsp),%r8
xor    %r9d,%r9d
xor    %r10d,%r10d
jmp    16ca <getQR_house+0x58a>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x1,%r11
add    $0x1,%r10
mov    0x38(%rsp),%rcx
add    %rcx,%r9
add    0x60(%rsp),%r8
add    %ecx,%ebx
add    %ecx,%edx
add    %rax,%rsi
mov    0x40(%rsp),%rbp
cmp    %rbp,%r11
je     1480 <getQR_house+0x340>
cmp    %rax,0x48(%rsp)
jg     16a0 <getQR_house+0x560>
mov    0x8(%rsp),%rcx
cmpq   $0x4,0x10(%rsp)
jae    17b0 <getQR_house+0x670>
sub    %ecx,%eax
mov    %rcx,%rbp
test   $0x1,%al
je     172e <getQR_house+0x5ee>
mov    %r11,%rax
imul   (%rsp),%rax
mov    %r11d,%edi
sub    0x8(%rsp),%edi
mov    0x58(%rsp),%rbp
imul   %edi,%ebp
add    0xa0(%rsp),%ebp
lea    (%rcx,%rbp,1),%edi
movslq %edi,%rdi
add    %rcx,%rax
mov    0x18(%rsp),%rbp
movsd  0x0(%rbp,%rax,8),%xmm0
subsd  0x0(%r13,%rdi,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rax,8)
lea    0x1(%rcx),%rbp
not    %rcx
cmp    0x100(%rsp),%rcx
mov    (%rsp),%rax
je     16a0 <getQR_house+0x560>
mov    %rax,%rcx
sub    %rbp,%rcx
lea    (%rbx,%rbp,1),%edi
lea    (%rdx,%rbp,1),%r14d
add    %rsi,%rbp
mov    0x18(%rsp),%rax
lea    (%rax,%rbp,8),%r12
mov    (%rsp),%rax
movslq %r14d,%r14
movsd  -0x8(%r12),%xmm0
movsd  (%r12),%xmm1
subsd  0x0(%r13,%r14,8),%xmm0
movsd  %xmm0,-0x8(%r12)
movslq %edi,%rdi
subsd  0x0(%r13,%rdi,8),%xmm1
movsd  %xmm1,(%r12)
add    $0x2,%edi
add    $0x2,%r14d
add    $0x10,%r12
add    $0xfffffffffffffffe,%rcx
jne    1760 <getQR_house+0x620>
jmp    16a0 <getQR_house+0x560>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    0x6c(%rsp),%ecx
mov    %ecx,%edi
imul   %r10d,%edi
mov    0x20(%rsp),%rcx
lea    (%rdi,%rcx,1),%ebp
mov    0x8(%rsp),%rcx
cmp    %edi,%ebp
jl     16e2 <getQR_house+0x5a2>
mov    0x20(%rsp),%rcx
shr    $0x20,%rcx
mov    0x8(%rsp),%rcx
jne    16e2 <getQR_house+0x5a2>
mov    %r10,%rcx
imul   %rax,%rcx
mov    0xe8(%rsp),%rax
lea    (%rax,%rcx,1),%rdi
mov    0x18(%rsp),%rax
lea    (%rax,%rdi,8),%rbp
add    0xf8(%rsp),%rbp
mov    0x10(%rsp),%rax
mov    %eax,%edi
imul   %r10d,%edi
movslq %edi,%rdi
mov    0xf0(%rsp),%rax
lea    (%rax,%rdi,8),%rax
cmp    %rax,%rbp
jae    1854 <getQR_house+0x714>
add    0xd8(%rsp),%rcx
mov    0x18(%rsp),%rax
lea    (%rax,%rcx,8),%rax
lea    0x0(,%rdi,8),%rdi
add    %r13,%rdi
mov    0x8(%rsp),%rcx
cmp    %rax,%rdi
mov    (%rsp),%rax
jb     16e2 <getQR_house+0x5a2>
xor    %ecx,%ecx
mov    0xd0(%rsp),%rdi
xchg   %ax,%ax
lea    (%r9,%rcx,1),%eax
cltq
movupd 0x0(%r13,%rax,8),%xmm0
movupd 0x10(%r13,%rax,8),%xmm1
movupd -0x10(%r8,%rcx,8),%xmm2
subpd  %xmm0,%xmm2
movupd (%r8,%rcx,8),%xmm0
subpd  %xmm1,%xmm0
movupd %xmm2,-0x10(%r8,%rcx,8)
movupd %xmm0,(%r8,%rcx,8)
add    $0x4,%rcx
cmp    %rcx,%rdi
jne    1860 <getQR_house+0x720>
mov    0xc0(%rsp),%rcx
mov    0xc8(%rsp),%rax
cmp    %rax,0x10(%rsp)
mov    (%rsp),%rax
je     16a0 <getQR_house+0x560>
jmp    16e2 <getQR_house+0x5a2>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    0x38(%rsp),%rax
and    $0xfffffffffffffffe,%rax
mov    0x2c(%rsp),%r10d
mov    %r10d,%ecx
mov    0x50(%rsp),%r8
neg    %r8
xor    %r9d,%r9d
mov    0x90(%rsp),%rdi
mov    0x8(%rsp),%r11
jmp    191a <getQR_house+0x7da>
nopl   0x0(%rax,%rax,1)
add    $0x1,%r11
add    $0x1,%r9d
add    $0x1,%r10d
add    0x60(%rsp),%rdi
cmp    %rbp,%r11
je     156b <getQR_house+0x42b>
mov    0x48(%rsp),%rdx
cmp    (%rsp),%rdx
jg     1900 <getQR_house+0x7c0>
cmpq   $0x0,0x20(%rsp)
je     1989 <getQR_house+0x849>
mov    %r10d,%ebp
mov    %r9d,%edx
xor    %esi,%esi
xor    %ebx,%ebx
nopw   0x0(%rax,%rax,1)
movsd  -0x8(%rdi,%rbx,8),%xmm0
movslq %edx,%rdx
movsd  %xmm0,0x0(%r13,%rdx,8)
add    %rcx,%rsi
movsd  (%rdi,%rbx,8),%xmm0
movslq %ebp,%rbp
movsd  %xmm0,0x0(%r13,%rbp,8)
add    %rcx,%rsi
add    $0x2,%rbx
add    %r15d,%edx
add    %r15d,%ebp
cmp    %rbx,%rax
jne    1940 <getQR_house+0x800>
add    %r8,%rbx
mov    0x40(%rsp),%rbp
testb  $0x1,0x10(%rsp)
je     1900 <getQR_house+0x7c0>
jmp    199b <getQR_house+0x85b>
xor    %esi,%esi
mov    0x8(%rsp),%rbx
testb  $0x1,0x10(%rsp)
je     1900 <getQR_house+0x7c0>
mov    %r11,%rdx
imul   (%rsp),%rdx
mov    0xa0(%rsp),%rbp
add    %r11d,%ebp
add    %rdx,%rbx
movsd  (%r12,%rbx,8),%xmm0
add    %ebp,%esi
mov    0x40(%rsp),%rbp
movslq %esi,%rdx
movsd  %xmm0,0x0(%r13,%rdx,8)
jmp    1900 <getQR_house+0x7c0>
mov    0xa8(%rsp),%rdi
call   1030 <free@plt>
mov    0x78(%rsp),%rdi
call   1030 <free@plt>
mov    %r13,%rdi
call   1030 <free@plt>
mov    0x70(%rsp),%rdi
add    $0x108,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

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

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