#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*/ exit(int);
int /*<<< orphan*/ free(double *);
double house(double *, int, double *);
scalar_t__ malloc(int);
int /*<<< orphan*/ mmult(double *, double *, double *, int, int, int);
int /*<<< orphan*/ printf(char *);
int /*<<< orphan*/ scale(double *, int, int, double);
void bidiag(double *A, int M, int N) {
  int j, i, k, u, t;
  double *x, *v, *AT, *w;
  double beta;
  if (M < N) {
    printf("M should be greater than or equal to N");
    exit(1);
  }
  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 (j = 0; j < N; ++j) {
    for (i = j; i < M; ++i) {
      x[i - j] = A[i * N + j];
    }
    beta = house(x, M - j, v);
    for (i = j; i < M; i++) {
      t = i * N;
      u = 0;
      for (k = j; k < N; k++) {
        AT[u + i - j] = A[k + t];
        u += (M - j);
      }
    }
    mmult(AT, v, w, N - j, M - j, 1);
    scale(w, N - j, 1, beta);
    mmult(v, w, AT, M - j, 1, N - j);
    for (i = j; i < M; i++) {
      t = i * N;
      for (k = j; k < N; k++) {
        A[t + k] -= AT[(i - j) * (N - j) + k - j];
      }
    }
    for (i = j + 1; i < M; ++i) {
      A[i * N + j] = v[i - j];
    }
    if (j < N - 2) {
      for (i = j + 1; i < N; ++i) {
        x[i - j - 1] = A[j * N + i];
      }
      beta = house(x, N - j - 1, v);
      for (i = j; i < M; i++) {
        t = i * N;
        u = (i - j) * (N - j - 1);
        for (k = j + 1; k < N; k++) {
          AT[u + k - j - 1] = A[k + t];
        }
      }
      mmult(AT, v, w, M - j, N - j - 1, 1);
      scale(w, M - j, 1, beta);
      mmult(w, v, AT, M - j, 1, N - j - 1);
      for (i = j; i < M; i++) {
        t = i * N;
        u = (i - j) * (N - j - 1);
        for (k = j + 1; k < N; k++) {
          A[k + t] -= AT[u + k - j - 1];
        }
      }
      u = 1;
      t = j * N;
      for (i = j + 2; i < N; ++i) {
        A[t + i] = v[u];
        u++;
      }
    }
  }
  free(x);
  free(v);
  free(AT);
}


x86 assembly (clang, -O0):
<bidiag>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    -0xc(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    11a4 <bidiag+0x34>
lea    0xe6b(%rip),%rdi
call   1050 <printf@plt>
mov    $0x1,%edi
call   1090 <exit@plt>
movslq -0xc(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   1060 <malloc@plt>
mov    %rax,-0x30(%rbp)
movslq -0xc(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   1060 <malloc@plt>
mov    %rax,-0x38(%rbp)
movslq -0xc(%rbp),%rax
shl    $0x3,%rax
movslq -0x10(%rbp),%rcx
imul   %rcx,%rax
mov    %eax,%edi
call   1060 <malloc@plt>
mov    %rax,-0x40(%rbp)
movslq -0xc(%rbp),%rax
shl    $0x3,%rax
movslq -0xc(%rbp),%rcx
imul   %rcx,%rax
mov    %eax,%edi
call   1060 <malloc@plt>
mov    %rax,-0x48(%rbp)
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    16b9 <bidiag+0x549>
mov    -0x14(%rbp),%eax
mov    %eax,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    125b <bidiag+0xeb>
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%ecx
imul   -0x10(%rbp),%ecx
add    -0x14(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
mov    -0x18(%rbp),%ecx
sub    -0x14(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1219 <bidiag+0xa9>
mov    -0x30(%rbp),%rdi
mov    -0xc(%rbp),%esi
sub    -0x14(%rbp),%esi
mov    -0x38(%rbp),%rdx
call   1080 <house@plt>
movsd  %xmm0,-0x50(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12fc <bidiag+0x18c>
mov    -0x18(%rbp),%eax
imul   -0x10(%rbp),%eax
mov    %eax,-0x24(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    12e9 <bidiag+0x179>
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x40(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    -0x18(%rbp),%ecx
sub    -0x14(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xc(%rbp),%eax
sub    -0x14(%rbp),%eax
add    -0x20(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    129c <bidiag+0x12c>
jmp    12ee <bidiag+0x17e>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1279 <bidiag+0x109>
mov    -0x40(%rbp),%rdi
mov    -0x38(%rbp),%rsi
mov    -0x48(%rbp),%rdx
mov    -0x10(%rbp),%ecx
sub    -0x14(%rbp),%ecx
mov    -0xc(%rbp),%r8d
sub    -0x14(%rbp),%r8d
mov    $0x1,%r9d
call   1040 <mmult@plt>
mov    -0x48(%rbp),%rdi
mov    -0x10(%rbp),%esi
sub    -0x14(%rbp),%esi
movsd  -0x50(%rbp),%xmm0
mov    $0x1,%edx
call   1070 <scale@plt>
mov    -0x38(%rbp),%rdi
mov    -0x48(%rbp),%rsi
mov    -0x40(%rbp),%rdx
mov    -0xc(%rbp),%ecx
sub    -0x14(%rbp),%ecx
mov    -0x10(%rbp),%r9d
sub    -0x14(%rbp),%r9d
mov    $0x1,%r8d
call   1040 <mmult@plt>
mov    -0x14(%rbp),%eax
mov    %eax,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    13ea <bidiag+0x27a>
mov    -0x18(%rbp),%eax
imul   -0x10(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    13d7 <bidiag+0x267>
mov    -0x40(%rbp),%rax
mov    -0x18(%rbp),%ecx
sub    -0x14(%rbp),%ecx
mov    -0x10(%rbp),%edx
sub    -0x14(%rbp),%edx
imul   %edx,%ecx
add    -0x1c(%rbp),%ecx
sub    -0x14(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1381 <bidiag+0x211>
jmp    13dc <bidiag+0x26c>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1365 <bidiag+0x1f5>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1435 <bidiag+0x2c5>
mov    -0x38(%rbp),%rax
mov    -0x18(%rbp),%ecx
sub    -0x14(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%ecx
imul   -0x10(%rbp),%ecx
add    -0x14(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    13f3 <bidiag+0x283>
mov    -0x14(%rbp),%eax
mov    -0x10(%rbp),%ecx
sub    $0x2,%ecx
cmp    %ecx,%eax
jge    16a6 <bidiag+0x536>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1494 <bidiag+0x324>
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   -0x10(%rbp),%ecx
add    -0x18(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
mov    -0x18(%rbp),%ecx
sub    -0x14(%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    144f <bidiag+0x2df>
mov    -0x30(%rbp),%rdi
mov    -0x10(%rbp),%esi
sub    -0x14(%rbp),%esi
sub    $0x1,%esi
mov    -0x38(%rbp),%rdx
call   1080 <house@plt>
movsd  %xmm0,-0x50(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1540 <bidiag+0x3d0>
mov    -0x18(%rbp),%eax
imul   -0x10(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x18(%rbp),%eax
sub    -0x14(%rbp),%eax
mov    -0x10(%rbp),%ecx
sub    -0x14(%rbp),%ecx
sub    $0x1,%ecx
imul   %ecx,%eax
mov    %eax,-0x20(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    152d <bidiag+0x3bd>
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x40(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    -0x1c(%rbp),%ecx
sub    -0x14(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    14e9 <bidiag+0x379>
jmp    1532 <bidiag+0x3c2>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    14b5 <bidiag+0x345>
mov    -0x40(%rbp),%rdi
mov    -0x38(%rbp),%rsi
mov    -0x48(%rbp),%rdx
mov    -0xc(%rbp),%ecx
sub    -0x14(%rbp),%ecx
mov    -0x10(%rbp),%r8d
sub    -0x14(%rbp),%r8d
sub    $0x1,%r8d
mov    $0x1,%r9d
call   1040 <mmult@plt>
mov    -0x48(%rbp),%rdi
mov    -0xc(%rbp),%esi
sub    -0x14(%rbp),%esi
movsd  -0x50(%rbp),%xmm0
mov    $0x1,%edx
call   1070 <scale@plt>
mov    -0x48(%rbp),%rdi
mov    -0x38(%rbp),%rsi
mov    -0x40(%rbp),%rdx
mov    -0xc(%rbp),%ecx
sub    -0x14(%rbp),%ecx
mov    -0x10(%rbp),%r9d
sub    -0x14(%rbp),%r9d
sub    $0x1,%r9d
mov    $0x1,%r8d
call   1040 <mmult@plt>
mov    -0x14(%rbp),%eax
mov    %eax,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1645 <bidiag+0x4d5>
mov    -0x18(%rbp),%eax
imul   -0x10(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x18(%rbp),%eax
sub    -0x14(%rbp),%eax
mov    -0x10(%rbp),%ecx
sub    -0x14(%rbp),%ecx
sub    $0x1,%ecx
imul   %ecx,%eax
mov    %eax,-0x20(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1632 <bidiag+0x4c2>
mov    -0x40(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    -0x1c(%rbp),%ecx
sub    -0x14(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    15e5 <bidiag+0x475>
jmp    1637 <bidiag+0x4c7>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    15b1 <bidiag+0x441>
movl   $0x1,-0x20(%rbp)
mov    -0x14(%rbp),%eax
imul   -0x10(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x14(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    16a1 <bidiag+0x531>
mov    -0x38(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    -0x18(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    165f <bidiag+0x4ef>
jmp    16a6 <bidiag+0x536>
jmp    16ab <bidiag+0x53b>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1207 <bidiag+0x97>
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):
<bidiag>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x218,%rsp
mov    %rdi,0x18(%rsp)
cmp    %edx,%esi
jl     2535 <bidiag+0x13c5>
mov    %edx,%ebx
mov    %esi,%ebp
lea    0x0(,%rbp,8),%r14d
mov    %r14d,%edi
call   1060 <malloc@plt>
mov    %rax,%r13
mov    %r14d,%edi
call   1060 <malloc@plt>
mov    %rax,%r15
mov    %r14d,%edi
imul   %ebx,%edi
call   1060 <malloc@plt>
mov    %rax,%r12
mov    %rbp,0xa0(%rsp)
imul   %ebp,%r14d
mov    %r14d,%edi
mov    %r15,%r14
call   1060 <malloc@plt>
mov    %rax,0xf8(%rsp)
mov    %rbx,0x130(%rsp)
test   %ebx,%ebx
jle    250c <bidiag+0x139c>
mov    0x130(%rsp),%rcx
lea    -0x2(%rcx),%eax
mov    %ecx,%r15d
mov    0xa0(%rsp),%rsi
movslq %esi,%rdx
mov    %rdx,0x50(%rsp)
cltq
mov    %rax,0x208(%rsp)
mov    %esi,%eax
mov    %rax,0x88(%rsp)
mov    0x18(%rsp),%rax
lea    -0x8(%rax),%rbx
mov    %rbx,0x1d0(%rsp)
lea    -0x1(%rcx),%r11d
lea    0x8(%r14),%rsi
mov    %rsi,0x1c8(%rsp)
lea    -0x2(%r15),%rsi
mov    %rsi,0x1e0(%rsp)
lea    0x0(,%r15,8),%rdx
lea    0x8(,%r15,8),%rsi
mov    %rsi,0x188(%rsp)
mov    %r15,%rsi
neg    %rsi
mov    %rsi,0x1b0(%rsp)
mov    %r15,%rsi
shl    $0x20,%rsi
lea    0x8(%rax),%r8
lea    -0x1(%r15),%rbx
lea    0x10(%rax),%r10
lea    (%rax,%r15,8),%r9
mov    %r9,0x170(%rsp)
lea    -0x5(%r15),%rdi
mov    %rdi,0x168(%rsp)
lea    -0x8(%r13),%rdi
lea    0x8(%r13),%r9
mov    %rsi,0x1b8(%rsp)
movabs $0xffffffff00000000,%rbp
ff ff ff
lea    (%rsi,%rbp,1),%rsi
mov    %rsi,0x108(%rsp)
mov    %rdi,%rsi
lea    0x38(%rax),%rdi
mov    %rdi,0x148(%rsp)
lea    0x18(%rax),%rdi
mov    %rdi,0x140(%rsp)
lea    0x18(%r14),%rdi
mov    %rdi,0x1d8(%rsp)
movl   $0x2,0x94(%rsp)
mov    $0x1,%edi
mov    %rdi,0x160(%rsp)
movq   $0xffffffffffffffff,0xf0(%rsp)
mov    $0x2,%edi
mov    %rdi,0x138(%rsp)
xor    %edi,%edi
mov    %rdi,0x60(%rsp)
mov    %r11,0x118(%rsp)
mov    %r11d,%edi
mov    %r11d,0x3c(%rsp)
mov    %rbx,0xd8(%rsp)
mov    %rbx,0xd0(%rsp)
mov    %rax,%r11
mov    %r8,0x1c0(%rsp)
mov    %r8,0xb0(%rsp)
mov    %r10,0x128(%rsp)
mov    %r10,0x158(%rsp)
mov    %r15,0xe0(%rsp)
mov    0x88(%rsp),%r10
mov    %rax,0xa8(%rsp)
mov    %ecx,%eax
mov    %ecx,0x9c(%rsp)
mov    $0x1,%eax
mov    %rax,0x10(%rsp)
mov    0xa0(%rsp),%rbx
xor    %ebp,%ebp
mov    %r13,0x150(%rsp)
mov    %rdx,0x68(%rsp)
mov    %r14,0x100(%rsp)
mov    %r15,0x30(%rsp)
jmp    14a2 <bidiag+0x332>
nopl   0x0(%rax)
addl   $0xffffffff,0x9c(%rsp)
ff
addq   $0xffffffffffffffff,0x60(%rsp)
mov    0x210(%rsp),%rbx
add    $0xffffffff,%ebx
addq   $0x1,0x10(%rsp)
addl   $0x1,0x94(%rsp)
01
mov    0x188(%rsp),%rcx
add    %rcx,0xa8(%rsp)
add    $0xffffffffffffffff,%r10
addq   $0xffffffffffffffff,0xe0(%rsp)
movabs $0xffffffff00000000,%rax
ff ff ff
add    %rax,0x1b8(%rsp)
add    %rcx,0x158(%rsp)
addq   $0xffffffffffffffff,0x160(%rsp)
add    %rdx,0xb0(%rsp)
add    %rcx,0x170(%rsp)
addq   $0xffffffffffffffff,0x168(%rsp)
mov    0x190(%rsp),%rsi
add    $0xfffffffffffffff8,%rsi
add    %rdx,%r11
add    $0xfffffffffffffff8,%r9
add    %rax,0x108(%rsp)
add    %rcx,0x148(%rsp)
addq   $0xffffffffffffffff,0xf0(%rsp)
addq   $0xffffffffffffffff,0xd0(%rsp)
add    %rcx,0x140(%rsp)
add    $0xffffffff,%r8d
add    %rdx,0x128(%rsp)
add    %r15,0x138(%rsp)
mov    %r8d,%eax
mov    %r8d,0x3c(%rsp)
cmp    %r15,%rbp
je     250c <bidiag+0x139c>
mov    %rbp,%rdi
mov    %rbp,%rax
imul   %r15,%rax
mov    %rax,0xb8(%rsp)
cmp    0x50(%rsp),%rbp
mov    %rbp,0x8(%rsp)
mov    %rsi,0x190(%rsp)
mov    %r9,0x1a8(%rsp)
mov    %r11,0x1a0(%rsp)
mov    %r10,0x198(%rsp)
mov    %rbx,0x210(%rsp)
jge    15c0 <bidiag+0x450>
shl    $0x20,%rbx
mov    0xa8(%rsp),%rax
xor    %ecx,%ecx
nopw   0x0(%rax,%rax,1)
movsd  (%rax),%xmm0
movsd  %xmm0,0x0(%r13,%rcx,8)
add    $0x1,%rcx
add    %rdx,%rax
cmp    %rcx,%r10
jne    1500 <bidiag+0x390>
mov    0xa0(%rsp),%rax
sub    %edi,%eax
mov    %r13,%rdi
mov    %eax,%r13d
mov    %eax,%esi
mov    %r14,%rdx
call   1080 <house@plt>
mov    0x8(%rsp),%rdx
cmp    0x50(%rsp),%rdx
movsd  %xmm0,0x20(%rsp)
jge    15e5 <bidiag+0x475>
mov    0xa8(%rsp),%rax
xor    %r10d,%r10d
mov    0x68(%rsp),%r8
movabs $0x100000000,%r9
mov    0xe0(%rsp),%rbp
nopl   0x0(%rax)
xor    %esi,%esi
mov    %r10,%rcx
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rax,%rsi,8),%xmm1
mov    %rcx,%rdi
sar    $0x1d,%rdi
movsd  %xmm1,(%r12,%rdi,1)
add    %rbx,%rcx
add    $0x1,%rsi
cmp    %rsi,%rbp
jne    1580 <bidiag+0x410>
add    $0x1,%rdx
add    %r9,%r10
add    %r8,%rax
cmp    0x88(%rsp),%rdx
jne    1570 <bidiag+0x400>
mov    0x8(%rsp),%rdx
jmp    15e5 <bidiag+0x475>
nopl   0x0(%rax)
mov    0xa0(%rsp),%rax
sub    %edi,%eax
mov    %r13,%rdi
mov    %eax,%r13d
mov    %eax,%esi
mov    %r14,%rdx
call   1080 <house@plt>
mov    0x8(%rsp),%rdx
movsd  %xmm0,0x20(%rsp)
mov    0x188(%rsp),%rax
imul   %rdx,%rax
mov    %rax,0xc8(%rsp)
mov    0xb8(%rsp),%rax
add    %r15,%rax
mov    %rax,0xc0(%rsp)
sub    %rdx,%r15
mov    0x130(%rsp),%rax
mov    %eax,%ebp
sub    %edx,%ebp
mov    %r12,%rdi
mov    %r14,%rsi
mov    0xf8(%rsp),%rbx
mov    %rbx,%rdx
mov    %ebp,%ecx
mov    %r13d,%r8d
mov    $0x1,%r9d
call   1040 <mmult@plt>
mov    %rbx,%rdi
mov    %r15,0x80(%rsp)
mov    %r15d,%esi
mov    $0x1,%edx
movsd  0x20(%rsp),%xmm0
call   1070 <scale@plt>
mov    %r14,%rdi
mov    %rbx,%rsi
mov    %r12,%rdx
mov    %r13d,0x98(%rsp)
mov    %r13d,%ecx
mov    $0x1,%r8d
mov    %ebp,0x38(%rsp)
mov    %ebp,%r9d
call   1040 <mmult@plt>
mov    0x8(%rsp),%rax
cmp    0x50(%rsp),%rax
movabs $0x200000000,%rdi
jge    1930 <bidiag+0x7c0>
mov    0x30(%rsp),%rax
mov    0x60(%rsp),%r13
lea    (%rax,%r13,1),%rcx
mov    0xe0(%rsp),%rax
and    $0xfffffffffffffffc,%rax
mov    %rax,0x78(%rsp)
lea    (%r12,%rcx,8),%rax
mov    %rax,0x58(%rsp)
mov    0xd8(%rsp),%rax
add    %r13,%rax
mov    %rax,0x28(%rsp)
mov    %rcx,0x20(%rsp)
and    $0xfffffffffffffffc,%rcx
mov    0x8(%rsp),%r10
mov    %rcx,0x48(%rsp)
lea    (%r10,%rcx,1),%rax
mov    %rax,0x40(%rsp)
mov    0xb0(%rsp),%rbx
mov    0x160(%rsp),%r15
mov    0x158(%rsp),%rsi
xor    %r8d,%r8d
xor    %r9d,%r9d
jmp    1747 <bidiag+0x5d7>
xchg   %ax,%ax
add    $0x1,%r10
add    $0x1,%r9
add    0x1b8(%rsp),%r8
mov    0x68(%rsp),%rax
add    %rax,%rsi
mov    0xe0(%rsp),%rcx
add    %rcx,%r13
add    %rcx,%r15
add    %rax,%rbx
cmp    0x88(%rsp),%r10
je     1930 <bidiag+0x7c0>
mov    0x8(%rsp),%r14
cmpq   $0x4,0x20(%rsp)
mov    0x30(%rsp),%rdx
jae    1820 <bidiag+0x6b0>
mov    %edx,%eax
sub    %r14d,%eax
mov    %r14,%r11
test   $0x1,%al
je     17ac <bidiag+0x63c>
mov    %r10,%rax
imul   0x30(%rsp),%rax
mov    %r10,%rcx
sub    0x8(%rsp),%rcx
imul   0x80(%rsp),%rcx
add    0x60(%rsp),%rcx
add    %r14d,%ecx
movslq %ecx,%rcx
add    %r14,%rax
mov    0x18(%rsp),%rbp
movsd  0x0(%rbp,%rax,8),%xmm0
subsd  (%r12,%rcx,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rax,8)
lea    0x1(%r14),%r11
cmp    %r14,0xd8(%rsp)
mov    0x30(%rsp),%rdx
je     1710 <bidiag+0x5a0>
lea    (%r11,%r13,1),%ebp
shl    $0x20,%rbp
lea    (%r15,%r11,1),%ecx
shl    $0x20,%rcx
nop
mov    %rbp,%rax
sar    $0x1d,%rax
movsd  -0x8(%rbx,%r11,8),%xmm0
movsd  (%rbx,%r11,8),%xmm1
subsd  (%r12,%rax,1),%xmm0
movsd  %xmm0,-0x8(%rbx,%r11,8)
mov    %rcx,%rax
sar    $0x1d,%rax
subsd  (%r12,%rax,1),%xmm1
movsd  %xmm1,(%rbx,%r11,8)
add    $0x2,%r11
add    %rdi,%rbp
add    %rdi,%rcx
cmp    %r11,%rdx
jne    17d0 <bidiag+0x660>
jmp    1710 <bidiag+0x5a0>
nopl   0x0(%rax,%rax,1)
mov    0x9c(%rsp),%eax
imul   %r9d,%eax
mov    0x28(%rsp),%rcx
add    %eax,%ecx
mov    0x8(%rsp),%r14
cmp    %eax,%ecx
jl     175d <bidiag+0x5ed>
mov    0x28(%rsp),%rax
shr    $0x20,%rax
mov    0x8(%rsp),%r14
jne    175d <bidiag+0x5ed>
mov    %r9,%rax
imul   %rdx,%rax
mov    0x18(%rsp),%rcx
lea    (%rcx,%rax,8),%rbp
add    0xc8(%rsp),%rbp
mov    0x20(%rsp),%rcx
imul   %r9d,%ecx
movslq %ecx,%r11
mov    0x58(%rsp),%rcx
lea    (%rcx,%r11,8),%rcx
cmp    %rcx,%rbp
jae    18a8 <bidiag+0x738>
add    0xc0(%rsp),%rax
mov    0x18(%rsp),%rcx
lea    (%rcx,%rax,8),%rax
lea    (%r12,%r11,8),%rcx
mov    0x8(%rsp),%r14
cmp    %rax,%rcx
jb     175d <bidiag+0x5ed>
mov    %r8,%rbp
xor    %ecx,%ecx
movabs $0x400000000,%rdi
mov    0x78(%rsp),%rdx
nopl   0x0(%rax)
mov    %rbp,%rax
sar    $0x1d,%rax
movupd (%r12,%rax,1),%xmm0
movupd 0x10(%r12,%rax,1),%xmm1
movupd -0x10(%rsi,%rcx,8),%xmm2
subpd  %xmm0,%xmm2
movupd (%rsi,%rcx,8),%xmm0
subpd  %xmm1,%xmm0
movupd %xmm2,-0x10(%rsi,%rcx,8)
movupd %xmm0,(%rsi,%rcx,8)
add    $0x4,%rcx
add    %rdi,%rbp
cmp    %rcx,%rdx
jne    18c0 <bidiag+0x750>
mov    0x40(%rsp),%r14
mov    0x48(%rsp),%rax
cmp    %rax,0x20(%rsp)
movabs $0x200000000,%rdi
mov    0x30(%rsp),%rdx
je     1710 <bidiag+0x5a0>
jmp    175d <bidiag+0x5ed>
nopw   0x0(%rax,%rax,1)
mov    0x8(%rsp),%rsi
lea    0x1(%rsi),%rbp
cmp    0x50(%rsp),%rbp
mov    0x100(%rsp),%r14
mov    0x68(%rsp),%rdx
mov    0x198(%rsp),%r10
jge    1986 <bidiag+0x816>
mov    $0x1,%eax
mov    0x170(%rsp),%rcx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  (%r14,%rax,8),%xmm0
movsd  %xmm0,(%rcx)
add    %rdx,%rcx
add    $0x1,%rax
cmp    %rax,%r10
jne    1970 <bidiag+0x800>
mov    0x3c(%rsp),%r8d
cmp    0x208(%rsp),%rsi
mov    0x150(%rsp),%r13
mov    0x30(%rsp),%r15
mov    0x1a8(%rsp),%r9
mov    0x1a0(%rsp),%r11
jge    13c0 <bidiag+0x250>
mov    0x168(%rsp),%rax
shr    $0x2,%rax
add    $0x1,%rax
and    $0xfffffffffffffffe,%rax
mov    %rax,0x110(%rsp)
not    %rsi
mov    %rsi,0x78(%rsp)
lea    (%rsi,%r15,1),%rax
mov    %rax,0x28(%rsp)
cmp    %r15,%rbp
mov    0x118(%rsp),%rdi
mov    0x38(%rsp),%ebx
jae    1b9f <bidiag+0xa2f>
mov    0x10(%rsp),%rax
cmpq   $0x4,0x28(%rsp)
jb     1b20 <bidiag+0x9b0>
mov    0x18(%rsp),%rax
mov    0xc0(%rsp),%rcx
lea    (%rax,%rcx,8),%rax
cmp    %rax,%r13
jae    1a50 <bidiag+0x8e0>
mov    0x28(%rsp),%rax
lea    0x0(,%rax,8),%rcx
add    %r13,%rcx
mov    0x1c0(%rsp),%rax
mov    0xc8(%rsp),%rdx
add    %rax,%rdx
mov    0x10(%rsp),%rax
cmp    %rcx,%rdx
jb     1b20 <bidiag+0x9b0>
mov    0xd8(%rsp),%rax
mov    0x60(%rsp),%rcx
add    %rcx,%rax
add    $0xfffffffffffffffc,%rax
mov    %rax,%rcx
shr    $0x2,%rcx
add    $0x1,%rcx
cmp    $0x4,%rax
mov    0xa8(%rsp),%rsi
jae    1a81 <bidiag+0x911>
xor    %edx,%edx
jmp    1aca <bidiag+0x95a>
mov    0x110(%rsp),%rax
xor    %edx,%edx
nopl   0x0(%rax,%rax,1)
movups 0x8(%rsi,%rdx,8),%xmm0
movups 0x18(%rsi,%rdx,8),%xmm1
movups %xmm0,0x0(%r13,%rdx,8)
movups %xmm1,0x10(%r13,%rdx,8)
movupd 0x28(%rsi,%rdx,8),%xmm0
movupd 0x38(%rsi,%rdx,8),%xmm1
movupd %xmm0,0x20(%r13,%rdx,8)
movupd %xmm1,0x30(%r13,%rdx,8)
add    $0x8,%rdx
add    $0xfffffffffffffffe,%rax
jne    1a90 <bidiag+0x920>
mov    0x28(%rsp),%rax
and    $0xfffffffffffffffc,%rax
test   $0x1,%cl
je     1b0b <bidiag+0x99b>
add    0x10(%rsp),%rdx
mov    0xb8(%rsp),%rcx
add    %rdx,%rcx
mov    0x18(%rsp),%rsi
movupd (%rsi,%rcx,8),%xmm0
movupd 0x10(%rsi,%rcx,8),%xmm1
add    0x78(%rsp),%rdx
movupd %xmm0,0x0(%r13,%rdx,8)
movupd %xmm1,0x10(%r13,%rdx,8)
cmp    %rax,0x28(%rsp)
je     1b9f <bidiag+0xa2f>
add    0x10(%rsp),%rax
nopl   0x0(%rax,%rax,1)
mov    %r15d,%edx
sub    %eax,%edx
mov    %rax,%rcx
not    %rcx
add    %r15,%rcx
and    $0x3,%rdx
mov    0x190(%rsp),%rsi
je     1b55 <bidiag+0x9e5>
nopl   0x0(%rax)
movsd  (%r11,%rax,8),%xmm0
movsd  %xmm0,(%rsi,%rax,8)
add    $0x1,%rax
add    $0xffffffffffffffff,%rdx
jne    1b40 <bidiag+0x9d0>
cmp    $0x3,%rcx
jb     1b9f <bidiag+0xa2f>
nopl   0x0(%rax,%rax,1)
movsd  (%r11,%rax,8),%xmm0
movsd  %xmm0,-0x10(%r9,%rax,8)
movsd  0x8(%r11,%rax,8),%xmm0
movsd  %xmm0,-0x8(%r9,%rax,8)
movsd  0x10(%r11,%rax,8),%xmm0
movsd  %xmm0,(%r9,%rax,8)
movsd  0x18(%r11,%rax,8),%xmm0
movsd  %xmm0,0x8(%r9,%rax,8)
add    $0x4,%rax
cmp    %rax,%r15
jne    1b60 <bidiag+0x9f0>
mov    %r8,0x70(%rsp)
mov    %rbp,0x20(%rsp)
mov    %edi,%eax
mov    0x8(%rsp),%rbp
sub    %ebp,%eax
mov    %rax,0x48(%rsp)
mov    0x80(%rsp),%rsi
add    $0xffffffffffffffff,%rsi
add    $0xffffffff,%ebx
mov    %r13,%rdi
mov    %rsi,0x80(%rsp)
mov    %r14,%rdx
call   1080 <house@plt>
cmp    0x50(%rsp),%rbp
movsd  %xmm0,0x120(%rsp)
jge    1f40 <bidiag+0xdd0>
mov    %rbp,%rcx
mov    %ebx,0x38(%rsp)
mov    0xd8(%rsp),%rax
sub    %rbp,%rax
add    $0xfffffffffffffffc,%rax
mov    %rax,0x200(%rsp)
mov    %rax,%rdx
shr    $0x2,%rdx
mov    0xd0(%rsp),%rax
and    $0xfffffffffffffffc,%rax
mov    %rax,0x1e8(%rsp)
add    $0x1,%rdx
mov    %rdx,0x1f8(%rsp)
mov    0x118(%rsp),%rbx
sub    %rbp,%rbx
mov    %rbx,0x180(%rsp)
mov    0x28(%rsp),%rax
lea    (%r12,%rax,8),%rdx
mov    %rdx,0x178(%rsp)
mov    0x1e0(%rsp),%rdx
sub    %rbp,%rdx
mov    %rdx,0x58(%rsp)
mov    %rax,%rdx
and    $0xfffffffffffffffc,%rdx
mov    0x10(%rsp),%rax
mov    %rdx,0xe8(%rsp)
add    %rdx,%rax
mov    %rax,0x1f0(%rsp)
mov    0x60(%rsp),%r14
mov    %rbp,%rax
mov    0xb0(%rsp),%rdi
mov    0xf0(%rsp),%r13
mov    0x148(%rsp),%rcx
xor    %r11d,%r11d
xor    %r9d,%r9d
jmp    1ce7 <bidiag+0xb77>
nopl   (%rax)
add    $0x1,%rax
add    $0x1,%r9
add    0x108(%rsp),%r11
mov    0x68(%rsp),%rdx
add    %rdx,%rcx
mov    0xd0(%rsp),%rsi
add    %rsi,%r13
add    %rdx,%rdi
add    %rsi,%r14
cmp    0x88(%rsp),%rax
je     1fb0 <bidiag+0xe40>
cmp    %r15,0x20(%rsp)
jae    1cb0 <bidiag+0xb40>
mov    %rax,%rbx
imul   %r15,%rbx
mov    %rax,%r8
sub    0x8(%rsp),%r8
imul   0x80(%rsp),%r8
add    0x78(%rsp),%r8
mov    0x10(%rsp),%rsi
cmpq   $0x4,0x28(%rsp)
jb     1ea0 <bidiag+0xd30>
mov    0x48(%rsp),%rdx
imul   %r9d,%edx
mov    0x58(%rsp),%rsi
lea    (%rdx,%rsi,1),%ebp
mov    0x10(%rsp),%rsi
cmp    %edx,%ebp
jl     1ea0 <bidiag+0xd30>
mov    0x58(%rsp),%rdx
shr    $0x20,%rdx
mov    0x10(%rsp),%rsi
jne    1ea0 <bidiag+0xd30>
mov    %r9,%rdx
imul   %r15,%rdx
mov    0x180(%rsp),%rsi
imul   %r9d,%esi
mov    0x18(%rsp),%rbp
lea    0x0(%rbp,%rdx,8),%rbp
mov    %rbx,0x40(%rsp)
mov    0xc8(%rsp),%rbx
lea    0x8(%rbx,%rbp,1),%rbp
movslq %esi,%rsi
mov    0x178(%rsp),%rbx
lea    (%rbx,%rsi,8),%rbx
cmp    %rbx,%rbp
mov    0x40(%rsp),%rbx
jae    1dc0 <bidiag+0xc50>
add    0xc0(%rsp),%rdx
mov    0x18(%rsp),%rbp
lea    0x0(,%rdx,8),%rdx
add    %rbp,%rdx
lea    (%r12,%rsi,8),%rbp
mov    0x10(%rsp),%rsi
cmp    %rdx,%rbp
jb     1ea0 <bidiag+0xd30>
cmpq   $0x4,0x200(%rsp)
jae    1dcf <bidiag+0xc5f>
xor    %esi,%esi
jmp    1e42 <bidiag+0xcd2>
mov    0x110(%rsp),%rdx
mov    %r11,%rbp
xor    %esi,%esi
movabs $0x800000000,%r10
cs nopw 0x0(%rax,%rax,1)
movups -0x30(%rcx,%rsi,8),%xmm0
movups -0x20(%rcx,%rsi,8),%xmm1
mov    %rbp,%rbx
sar    $0x1d,%rbx
movups %xmm0,(%r12,%rbx,1)
movups %xmm1,0x10(%r12,%rbx,1)
movupd -0x10(%rcx,%rsi,8),%xmm0
movupd (%rcx,%rsi,8),%xmm1
movabs $0x400000000,%rbx
add    %rbp,%rbx
sar    $0x1d,%rbx
movupd %xmm0,(%r12,%rbx,1)
movupd %xmm1,0x10(%r12,%rbx,1)
add    $0x8,%rsi
add    %r10,%rbp
add    $0xfffffffffffffffe,%rdx
jne    1df0 <bidiag+0xc80>
testb  $0x1,0x1f8(%rsp)
01
mov    0x40(%rsp),%rbx
je     1e7e <bidiag+0xd0e>
add    0x10(%rsp),%rsi
lea    (%rsi,%rbx,1),%rdx
mov    0x18(%rsp),%rbp
movupd 0x0(%rbp,%rdx,8),%xmm0
movupd 0x10(%rbp,%rdx,8),%xmm1
add    %r8d,%esi
movslq %esi,%rdx
movupd %xmm0,(%r12,%rdx,8)
movupd %xmm1,0x10(%r12,%rdx,8)
mov    0x1f0(%rsp),%rsi
mov    0xe8(%rsp),%rdx
cmp    %rdx,0x28(%rsp)
je     1cb0 <bidiag+0xb40>
nopl   0x0(%rax)
mov    %r15d,%edx
sub    %esi,%edx
mov    %rsi,%r10
test   $0x1,%dl
je     1eca <bidiag+0xd5a>
add    %rsi,%rbx
mov    0x18(%rsp),%rdx
movsd  (%rdx,%rbx,8),%xmm0
add    %esi,%r8d
movslq %r8d,%rdx
movsd  %xmm0,(%r12,%rdx,8)
lea    0x1(%rsi),%r10
not    %rsi
cmp    0x1b0(%rsp),%rsi
movabs $0x200000000,%rbx
je     1cb0 <bidiag+0xb40>
lea    (%r10,%r13,1),%ebp
shl    $0x20,%rbp
lea    (%r14,%r10,1),%edx
shl    $0x20,%rdx
cs nopw 0x0(%rax,%rax,1)
nop
movsd  -0x8(%rdi,%r10,8),%xmm0
mov    %rbp,%rsi
sar    $0x1d,%rsi
movsd  %xmm0,(%r12,%rsi,1)
movsd  (%rdi,%r10,8),%xmm0
mov    %rdx,%rsi
sar    $0x1d,%rsi
movsd  %xmm0,(%r12,%rsi,1)
add    $0x2,%r10
add    %rbx,%rbp
add    %rbx,%rdx
cmp    %r10,%r15
jne    1f00 <bidiag+0xd90>
jmp    1cb0 <bidiag+0xb40>
nopl   0x0(%rax,%rax,1)
mov    %r12,%rdi
mov    %r14,%rsi
mov    0xf8(%rsp),%rbp
mov    %rbp,%rdx
mov    0x98(%rsp),%r15d
mov    %r15d,%ecx
mov    %ebx,%r8d
mov    $0x1,%r9d
call   1040 <mmult@plt>
mov    %rbp,%rdi
mov    %r15d,%esi
mov    $0x1,%edx
movsd  0x120(%rsp),%xmm0
call   1070 <scale@plt>
mov    %rbp,%rdi
mov    %r14,%rsi
mov    %r12,%rdx
mov    %r15d,%ecx
mov    $0x1,%r8d
mov    %ebx,%r9d
call   1040 <mmult@plt>
mov    0x20(%rsp),%rbp
jmp    22b0 <bidiag+0x1140>
nopw   0x0(%rax,%rax,1)
mov    %r12,%rdi
mov    0x100(%rsp),%r14
mov    %r14,%rsi
mov    0xf8(%rsp),%rbp
mov    %rbp,%rdx
mov    0x98(%rsp),%ebx
mov    %ebx,%ecx
mov    0x38(%rsp),%r15d
mov    %r15d,%r8d
mov    $0x1,%r9d
call   1040 <mmult@plt>
mov    %rbp,%rdi
mov    %ebx,%esi
mov    $0x1,%edx
movsd  0x120(%rsp),%xmm0
call   1070 <scale@plt>
mov    %rbp,%rdi
mov    %r14,%rsi
mov    %r12,%rdx
mov    %ebx,%ecx
mov    $0x1,%r8d
mov    %r15d,%r9d
call   1040 <mmult@plt>
mov    0x8(%rsp),%rax
cmp    0x50(%rsp),%rax
mov    0x150(%rsp),%r13
mov    0x20(%rsp),%rbp
jge    22b0 <bidiag+0x1140>
mov    0x10(%rsp),%rax
mov    0xe8(%rsp),%rcx
add    %rcx,%rax
mov    %rax,0x40(%rsp)
mov    0xb0(%rsp),%rcx
mov    0x60(%rsp),%r10
mov    0xf0(%rsp),%r11
mov    0x140(%rsp),%rdi
xor    %r9d,%r9d
xor    %r14d,%r14d
mov    0x8(%rsp),%r15
jmp    20bc <bidiag+0xf4c>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%r15
add    $0x1,%r14
add    0x108(%rsp),%r9
mov    0x68(%rsp),%rax
add    %rax,%rdi
mov    0xd0(%rsp),%rdx
add    %rdx,%r11
add    %rdx,%r10
add    %rax,%rcx
cmp    0x88(%rsp),%r15
mov    0x20(%rsp),%rbp
je     22b0 <bidiag+0x1140>
mov    0x30(%rsp),%rbx
cmp    %rbx,%rbp
jae    2080 <bidiag+0xf10>
mov    0x10(%rsp),%rsi
cmpq   $0x4,0x28(%rsp)
jb     21e0 <bidiag+0x1070>
mov    0x48(%rsp),%rax
imul   %r14d,%eax
mov    0x58(%rsp),%rdx
add    %eax,%edx
mov    0x10(%rsp),%rsi
cmp    %eax,%edx
jl     21e0 <bidiag+0x1070>
mov    0x58(%rsp),%rax
shr    $0x20,%rax
mov    0x10(%rsp),%rsi
jne    21e0 <bidiag+0x1070>
mov    %r14,%rax
imul   %rbx,%rax
mov    0x180(%rsp),%rdx
imul   %r14d,%edx
mov    0x18(%rsp),%rsi
lea    (%rsi,%rax,8),%rsi
mov    0xc8(%rsp),%rbp
add    %rbp,%rsi
add    $0x8,%rsi
movslq %edx,%rdx
mov    0x178(%rsp),%rbp
lea    0x0(%rbp,%rdx,8),%rbp
cmp    %rbp,%rsi
jae    2167 <bidiag+0xff7>
add    0xc0(%rsp),%rax
mov    0x18(%rsp),%rsi
lea    (%rsi,%rax,8),%rax
lea    (%r12,%rdx,8),%rdx
mov    0x10(%rsp),%rsi
cmp    %rax,%rdx
jb     21e0 <bidiag+0x1070>
mov    %r9,%rax
xor    %edx,%edx
movabs $0x400000000,%rbp
mov    0x1e8(%rsp),%rbx
xchg   %ax,%ax
mov    %rax,%rsi
sar    $0x1d,%rsi
movupd (%r12,%rsi,1),%xmm0
movupd 0x10(%r12,%rsi,1),%xmm1
movupd -0x10(%rdi,%rdx,8),%xmm2
subpd  %xmm0,%xmm2
movupd (%rdi,%rdx,8),%xmm0
subpd  %xmm1,%xmm0
movupd %xmm2,-0x10(%rdi,%rdx,8)
movupd %xmm0,(%rdi,%rdx,8)
add    $0x4,%rdx
add    %rbp,%rax
cmp    %rdx,%rbx
jne    2180 <bidiag+0x1010>
mov    0x40(%rsp),%rsi
mov    0xe8(%rsp),%rax
cmp    %rax,0x28(%rsp)
mov    0x30(%rsp),%rbx
je     2080 <bidiag+0xf10>
nopl   0x0(%rax,%rax,1)
mov    %ebx,%eax
sub    %esi,%eax
mov    %rsi,%rdx
test   $0x1,%al
je     222d <bidiag+0x10bd>
mov    %r15,%rax
imul   0x30(%rsp),%rax
mov    %r15,%rdx
sub    0x8(%rsp),%rdx
imul   0x80(%rsp),%rdx
add    0x78(%rsp),%rdx
add    %esi,%edx
movslq %edx,%rdx
add    %rsi,%rax
mov    0x18(%rsp),%rbp
movsd  0x0(%rbp,%rax,8),%xmm0
subsd  (%r12,%rdx,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rax,8)
lea    0x1(%rsi),%rdx
not    %rsi
cmp    0x1b0(%rsp),%rsi
movabs $0x200000000,%rbx
mov    0x30(%rsp),%r8
je     2080 <bidiag+0xf10>
lea    (%rdx,%r11,1),%ebp
shl    $0x20,%rbp
lea    (%r10,%rdx,1),%eax
shl    $0x20,%rax
nopl   (%rax)
mov    %rbp,%rsi
sar    $0x1d,%rsi
movsd  -0x8(%rcx,%rdx,8),%xmm0
movsd  (%rcx,%rdx,8),%xmm1
subsd  (%r12,%rsi,1),%xmm0
movsd  %xmm0,-0x8(%rcx,%rdx,8)
mov    %rax,%rsi
sar    $0x1d,%rsi
subsd  (%r12,%rsi,1),%xmm1
movsd  %xmm1,(%rcx,%rdx,8)
add    $0x2,%rdx
add    %rbx,%rbp
add    %rbx,%rax
cmp    %rdx,%r8
jne    2260 <bidiag+0x10f0>
jmp    2080 <bidiag+0xf10>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    0x8(%rsp),%rax
add    $0x2,%rax
mov    0x30(%rsp),%r15
cmp    %r15,%rax
mov    0x100(%rsp),%r14
mov    0x68(%rsp),%rdx
mov    0x1a8(%rsp),%r9
mov    0x1a0(%rsp),%r11
mov    0x198(%rsp),%r10
mov    0x70(%rsp),%r8
jae    13c0 <bidiag+0x250>
movslq 0x94(%rsp),%rax
mov    0x48(%rsp),%rbx
lea    -0x1(%rbx),%r8
mov    $0x1,%ecx
cmp    $0x4,%r8
jb     2440 <bidiag+0x12d0>
mov    0xb8(%rsp),%rsi
add    %rax,%rsi
mov    0x18(%rsp),%rdi
lea    (%rdi,%rsi,8),%rdi
lea    (%r14,%rbx,8),%rbp
cmp    %rbp,%rdi
jae    2351 <bidiag+0x11e1>
add    %rbx,%rsi
mov    0x1d0(%rsp),%rdi
lea    (%rdi,%rsi,8),%rsi
cmp    %rsi,0x1c8(%rsp)
jae    2351 <bidiag+0x11e1>
mov    0x20(%rsp),%rbp
jmp    2440 <bidiag+0x12d0>
lea    -0x5(%rbx),%rsi
mov    %rsi,%rcx
shr    $0x2,%rcx
add    $0x1,%rcx
cmp    $0x4,%rsi
jae    236f <bidiag+0x11ff>
xor    %esi,%esi
mov    0x20(%rsp),%rbp
jmp    23d9 <bidiag+0x1269>
mov    0x70(%rsp),%rsi
lea    -0x5(%rsi),%rdi
shr    $0x2,%rdi
add    $0x1,%rdi
and    $0xfffffffffffffffe,%rdi
mov    0x128(%rsp),%rsi
lea    (%rsi,%rax,8),%rbx
xor    %esi,%esi
mov    0x20(%rsp),%rbp
nopw   0x0(%rax,%rax,1)
movups 0x8(%r14,%rsi,8),%xmm0
movups 0x18(%r14,%rsi,8),%xmm1
movups %xmm0,-0x10(%rbx,%rsi,8)
movups %xmm1,(%rbx,%rsi,8)
movupd 0x28(%r14,%rsi,8),%xmm0
movupd 0x38(%r14,%rsi,8),%xmm1
movupd %xmm0,0x10(%rbx,%rsi,8)
movupd %xmm1,0x20(%rbx,%rsi,8)
add    $0x8,%rsi
add    $0xfffffffffffffffe,%rdi
jne    23a0 <bidiag+0x1230>
mov    %r8,%rdi
and    $0xfffffffffffffffc,%rdi
test   $0x1,%cl
je     2419 <bidiag+0x12a9>
lea    0x0(,%rsi,8),%rcx
add    %rax,%rsi
or     $0x8,%rcx
movupd (%r14,%rcx,1),%xmm0
movupd 0x10(%r14,%rcx,1),%xmm1
add    0xb8(%rsp),%rsi
mov    0x18(%rsp),%rcx
movupd %xmm0,(%rcx,%rsi,8)
movupd %xmm1,0x10(%rcx,%rsi,8)
cmp    %rdi,%r8
mov    0x70(%rsp),%r8
je     13c0 <bidiag+0x250>
lea    0x1(%rdi),%rcx
add    %rdi,%rax
mov    0x48(%rsp),%rbx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %ebx,%esi
sub    %ecx,%esi
mov    %rcx,%rdi
not    %rdi
add    %rbx,%rdi
test   $0x3,%sil
je     2496 <bidiag+0x1326>
mov    0x3c(%rsp),%esi
sub    %cl,%sil
movzbl %sil,%esi
and    $0x3,%esi
mov    0x70(%rsp),%r8
cs nopw 0x0(%rax,%rax,1)
movsd  (%r14,%rcx,8),%xmm0
movsd  %xmm0,(%r11,%rax,8)
add    $0x1,%rcx
add    $0x1,%rax
add    $0xffffffffffffffff,%rsi
jne    2470 <bidiag+0x1300>
cmp    $0x3,%rdi
jb     13c0 <bidiag+0x250>
jmp    24a5 <bidiag+0x1335>
mov    0x70(%rsp),%r8
cmp    $0x3,%rdi
jb     13c0 <bidiag+0x250>
mov    %r8,%rdi
sub    %rcx,%rdi
mov    0x1d8(%rsp),%rsi
lea    (%rsi,%rcx,8),%rcx
add    0x138(%rsp),%rax
mov    0x18(%rsp),%rsi
lea    (%rsi,%rax,8),%rax
xor    %esi,%esi
nopw   0x0(%rax,%rax,1)
movsd  -0x18(%rcx,%rsi,8),%xmm0
movsd  %xmm0,-0x10(%rax,%rsi,8)
movsd  -0x10(%rcx,%rsi,8),%xmm0
movsd  %xmm0,-0x8(%rax,%rsi,8)
movsd  -0x8(%rcx,%rsi,8),%xmm0
movsd  %xmm0,(%rax,%rsi,8)
movsd  (%rcx,%rsi,8),%xmm0
movsd  %xmm0,0x8(%rax,%rsi,8)
add    $0x4,%rsi
cmp    %rsi,%rdi
jne    24d0 <bidiag+0x1360>
jmp    13c0 <bidiag+0x250>
mov    %r13,%rdi
call   1030 <free@plt>
mov    %r14,%rdi
call   1030 <free@plt>
mov    %r12,%rdi
add    $0x218,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>
lea    0xac4(%rip),%rdi
call   1050 <printf@plt>
mov    $0x1,%edi
call   1090 <exit@plt>

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

LLVM IR (-O2):
@.str = [39 x i8] c"M should be greater than or equal to N\00"
define void @bidiag(double* nocapture %0, i32 %1, i32 %2) {
%4 = icmp slt i32 %1, %2
br i1 %4, label %5, label %8
5: ; preds = %3
%6 = tail call i32 @printf(i8* dereferenceable(1) getelementptr inbounds ([39 x i8], [39 x i8]* @.str, i64 0, i64 0))
%7 = tail call i32 @exit(i32 1)
unreachable
8: ; preds = %3
%9 = shl i32 %1, 3
%10 = tail call i64 @malloc(i32 %9)
%11 = inttoptr i64 %10 to double*
%12 = inttoptr i64 %10 to i8*
%13 = tail call i64 @malloc(i32 %9)
%14 = inttoptr i64 %13 to double*
%15 = mul i32 %9, %2
%16 = tail call i64 @malloc(i32 %15)
%17 = inttoptr i64 %16 to double*
%18 = mul i32 %9, %1
%19 = tail call i64 @malloc(i32 %18)
%20 = inttoptr i64 %19 to double*
%21 = icmp sgt i32 %2, 0
br i1 %21, label %22, label %828
22: ; preds = %8
%23 = add nsw i32 %2, -2
%24 = zext i32 %2 to i64
%25 = sext i32 %1 to i64
%26 = sext i32 %23 to i64
%27 = zext i32 %2 to i64
%28 = zext i32 %1 to i64
%29 = zext i32 %1 to i64
%30 = zext i32 %1 to i64
%31 = zext i32 %1 to i64
%32 = zext i32 %1 to i64
%33 = zext i32 %1 to i64
%34 = getelementptr double, double* %0, i64 -1
%35 = add i32 %2, -1
%36 = getelementptr double, double* %14, i64 1
%37 = add nsw i64 %27, -2
%38 = shl nuw nsw i64 %27, 3
%39 = add nuw nsw i64 %38, 8
%40 = add i32 %2, -1
%41 = zext i32 %40 to i64
%42 = add nsw i64 %27, -2
%43 = shl nuw nsw i64 %27, 3
%44 = add nuw nsw i64 %43, 8
%45 = add i32 %2, -1
%46 = zext i32 %45 to i64
%47 = getelementptr double, double* %0, i64 1
%48 = bitcast double* %47 to i8*
%49 = shl nuw nsw i64 %27, 3
%50 = add nuw nsw i64 %49, 8
%51 = add nsw i64 %27, -1
%52 = shl nuw nsw i64 %27, 3
%53 = add nuw nsw i64 %52, 8
%54 = add nsw i64 %27, -1
%55 = sub nsw i64 0, %27
%56 = sub nsw i64 0, %27
br label %57
57: ; preds = %22, %822
%58 = phi i64 [ 0, %22 ], [ %823, %822 ]
%59 = phi i32 [ %2, %22 ], [ %138, %822 ]
%60 = phi i32 [ 2, %22 ], [ %826, %822 ]
%61 = phi i64 [ 1, %22 ], [ %825, %822 ]
%62 = phi i32 [ %1, %22 ], [ %824, %822 ]
%63 = phi i64 [ 0, %22 ], [ %300, %822 ]
%64 = trunc i64 %63 to i32
%65 = sub i32 %40, %64
%66 = zext i32 %65 to i64
%67 = trunc i64 %63 to i32
%68 = sub i32 %40, %67
%69 = zext i32 %68 to i64
%70 = add nsw i64 %69, -5
%71 = lshr i64 %70, 2
%72 = add nuw nsw i64 %71, 1
%73 = sub i64 %54, %63
%74 = add i64 %73, -4
%75 = lshr i64 %74, 2
%76 = add nuw nsw i64 %75, 1
%77 = add i64 %54, %58
%78 = add i64 %77, -4
%79 = lshr i64 %78, 2
%80 = add nuw nsw i64 %79, 1
%81 = add i64 %58, %27
%82 = mul i64 %53, %63
%83 = mul i64 %63, %27
%84 = add i64 %83, %27
%85 = add i64 %58, %27
%86 = getelementptr double, double* %17, i64 %85
%87 = add i64 %51, %58
%88 = xor i64 %63, -1
%89 = add nsw i64 %88, %27
%90 = xor i64 %63, -1
%91 = add nsw i64 %90, %27
%92 = getelementptr double, double* %11, i64 %91
%93 = bitcast double* %92 to i8*
%94 = mul i64 %50, %63
%95 = getelementptr i8, i8* %48, i64 %94
%96 = mul i64 %63, %27
%97 = add i64 %96, %27
%98 = getelementptr double, double* %0, i64 %97
%99 = bitcast double* %98 to i8*
%100 = xor i64 %63, -1
%101 = add nsw i64 %100, %27
%102 = mul i64 %44, %63
%103 = mul i64 %63, %27
%104 = add i64 %103, %27
%105 = sub nsw i64 %46, %63
%106 = xor i64 %63, -1
%107 = add nsw i64 %106, %27
%108 = getelementptr double, double* %17, i64 %107
%109 = sub i64 %42, %63
%110 = trunc i64 %63 to i32
%111 = sub i32 %40, %110
%112 = xor i64 %63, -1
%113 = add nsw i64 %112, %27
%114 = mul i64 %39, %63
%115 = mul i64 %63, %27
%116 = add i64 %115, %27
%117 = sub nsw i64 %41, %63
%118 = xor i64 %63, -1
%119 = add nsw i64 %118, %27
%120 = getelementptr double, double* %17, i64 %119
%121 = sub i64 %37, %63
%122 = trunc i64 %63 to i32
%123 = sub i32 %35, %122
%124 = trunc i64 %63 to i32
%125 = sub i32 %35, %124
%126 = zext i32 %125 to i64
%127 = add nsw i64 %126, -1
%128 = mul i64 %63, %27
%129 = sext i32 %60 to i64
%130 = add i64 %128, %129
%131 = getelementptr double, double* %0, i64 %130
%132 = trunc i64 %63 to i32
%133 = sub i32 %35, %132
%134 = zext i32 %133 to i64
%135 = add i64 %130, %134
%136 = getelementptr double, double* %34, i64 %135
%137 = getelementptr double, double* %14, i64 %134
%138 = add i32 %59, -1
%139 = trunc i64 %63 to i32
%140 = sext i32 %60 to i64
%141 = zext i32 %62 to i64
%142 = icmp slt i64 %63, %25
br i1 %142, label %147, label %143
143: ; preds = %57
%144 = trunc i64 %63 to i32
%145 = sub i32 %1, %144
%146 = tail call double @house(double* %11, i32 %145, double* %14)
br label %181
147: ; preds = %57, %147
%148 = phi i64 [ %155, %147 ], [ %63, %57 ]
%149 = mul nsw i64 %148, %24
%150 = add nuw nsw i64 %149, %63
%151 = getelementptr inbounds double, double* %0, i64 %150
%152 = load double, double* %151
%153 = sub nuw nsw i64 %148, %63
%154 = getelementptr inbounds double, double* %11, i64 %153
store double %152, double* %154
%155 = add nuw nsw i64 %148, 1
%156 = icmp eq i64 %155, %28
br i1 %156, label %157, label %147
157: ; preds = %147
%158 = trunc i64 %63 to i32
%159 = sub i32 %1, %158
%160 = tail call double @house(double* %11, i32 %159, double* %14)
br i1 %142, label %161, label %181
161: ; preds = %157, %178
%162 = phi i64 [ %179, %178 ], [ %63, %157 ]
%163 = mul nsw i64 %162, %24
%164 = add nsw i64 %162, %58
br label %165
165: ; preds = %161, %165
%166 = phi i64 [ 0, %161 ], [ %175, %165 ]
%167 = phi i64 [ %63, %161 ], [ %176, %165 ]
%168 = add nuw nsw i64 %167, %163
%169 = getelementptr inbounds double, double* %0, i64 %168
%170 = load double, double* %169
%171 = add i64 %164, %166
%172 = shl i64 %171, 32
%173 = ashr exact i64 %172, 32
%174 = getelementptr inbounds double, double* %17, i64 %173
store double %170, double* %174
%175 = add i64 %166, %141
%176 = add nuw nsw i64 %167, 1
%177 = icmp eq i64 %176, %27
br i1 %177, label %178, label %165
178: ; preds = %165
%179 = add nuw nsw i64 %162, 1
%180 = icmp eq i64 %179, %29
br i1 %180, label %181, label %161
181: ; preds = %178, %143, %157
%182 = phi double [ %146, %143 ], [ %160, %157 ], [ %160, %178 ]
%183 = phi i32 [ %145, %143 ], [ %159, %157 ], [ %159, %178 ]
%184 = sub nsw i64 %24, %63
%185 = sub nsw i32 %2, %139
%186 = tail call i32 @mmult(double* %17, double* %14, double* %20, i32 %185, i32 %183, i32 1)
%187 = trunc i64 %184 to i32
%188 = tail call i32 @scale(double* %20, i32 %187, i32 1, double %182)
%189 = tail call i32 @mmult(double* %14, double* %20, double* %17, i32 %183, i32 1, i32 %185)
br i1 %142, label %190, label %299
190: ; preds = %181
%191 = icmp ult i64 %81, 4
%192 = trunc i64 %87 to i32
%193 = icmp ugt i64 %87, 4294967295
%194 = and i64 %81, -4
%195 = add i64 %63, %194
%196 = icmp eq i64 %81, %194
br label %197
197: ; preds = %190, %294
%198 = phi i64 [ %298, %294 ], [ 0, %190 ]
%199 = phi i32 [ %297, %294 ], [ 0, %190 ]
%200 = phi i64 [ %295, %294 ], [ %63, %190 ]
%201 = mul i64 %198, %27
%202 = getelementptr double, double* %0, i64 %201
%203 = bitcast double* %202 to i8*
%204 = getelementptr i8, i8* %203, i64 %82
%205 = add i64 %84, %201
%206 = getelementptr double, double* %0, i64 %205
%207 = mul i64 %85, %198
%208 = shl i64 %207, 32
%209 = ashr exact i64 %208, 32
%210 = getelementptr double, double* %17, i64 %209
%211 = getelementptr double, double* %86, i64 %209
%212 = bitcast double* %211 to i8*
%213 = mul i32 %59, %199
%214 = mul nsw i64 %200, %24
%215 = sub nsw i64 %200, %63
%216 = mul i64 %215, %184
%217 = add i64 %216, %58
br i1 %191, label %252, label %218
218: ; preds = %197
%219 = add i32 %213, %192
%220 = icmp slt i32 %219, %213
%221 = or i1 %220, %193
br i1 %221, label %252, label %222
222: ; preds = %218
%223 = icmp ult i8* %204, %212
%224 = icmp ult double* %210, %206
%225 = and i1 %223, %224
br i1 %225, label %252, label %226
226: ; preds = %222, %226
%227 = phi i64 [ %249, %226 ], [ 0, %222 ]
%228 = add i64 %63, %227
%229 = add i64 %217, %228
%230 = shl i64 %229, 32
%231 = ashr exact i64 %230, 32
%232 = getelementptr inbounds double, double* %17, i64 %231
%233 = bitcast double* %232 to <2 x double>*
%234 = load <2 x double>, <2 x double>* %233
%235 = getelementptr inbounds double, double* %232, i64 2
%236 = bitcast double* %235 to <2 x double>*
%237 = load <2 x double>, <2 x double>* %236
%238 = add nuw nsw i64 %228, %214
%239 = getelementptr inbounds double, double* %0, i64 %238
%240 = bitcast double* %239 to <2 x double>*
%241 = load <2 x double>, <2 x double>* %240
%242 = getelementptr inbounds double, double* %239, i64 2
%243 = bitcast double* %242 to <2 x double>*
%244 = load <2 x double>, <2 x double>* %243
%245 = fsub <2 x double> %241, %234
%246 = fsub <2 x double> %244, %237
%247 = bitcast double* %239 to <2 x double>*
store <2 x double> %245, <2 x double>* %247
%248 = bitcast double* %242 to <2 x double>*
store <2 x double> %246, <2 x double>* %248
%249 = add nuw i64 %227, 4
%250 = icmp eq i64 %249, %194
br i1 %250, label %251, label %226
251: ; preds = %226
br i1 %196, label %294, label %252
252: ; preds = %222, %218, %197, %251
%253 = phi i64 [ %63, %222 ], [ %63, %218 ], [ %63, %197 ], [ %195, %251 ]
%254 = sub i64 %27, %253
%255 = and i64 %254, 1
%256 = icmp eq i64 %255, 0
br i1 %256, label %268, label %257
257: ; preds = %252
%258 = add i64 %217, %253
%259 = shl i64 %258, 32
%260 = ashr exact i64 %259, 32
%261 = getelementptr inbounds double, double* %17, i64 %260
%262 = load double, double* %261
%263 = add nuw nsw i64 %253, %214
%264 = getelementptr inbounds double, double* %0, i64 %263
%265 = load double, double* %264
%266 = fsub double %265, %262
store double %266, double* %264
%267 = add nuw nsw i64 %253, 1
br label %268
268: ; preds = %257, %252
%269 = phi i64 [ %253, %252 ], [ %267, %257 ]
%270 = icmp eq i64 %54, %253
br i1 %270, label %294, label %271
271: ; preds = %268, %271
%272 = phi i64 [ %292, %271 ], [ %269, %268 ]
%273 = add i64 %217, %272
%274 = shl i64 %273, 32
%275 = ashr exact i64 %274, 32
%276 = getelementptr inbounds double, double* %17, i64 %275
%277 = load double, double* %276
%278 = add nuw nsw i64 %272, %214
%279 = getelementptr inbounds double, double* %0, i64 %278
%280 = load double, double* %279
%281 = fsub double %280, %277
store double %281, double* %279
%282 = add nuw nsw i64 %272, 1
%283 = add i64 %217, %282
%284 = shl i64 %283, 32
%285 = ashr exact i64 %284, 32
%286 = getelementptr inbounds double, double* %17, i64 %285
%287 = load double, double* %286
%288 = add nuw nsw i64 %282, %214
%289 = getelementptr inbounds double, double* %0, i64 %288
%290 = load double, double* %289
%291 = fsub double %290, %287
store double %291, double* %289
%292 = add nuw nsw i64 %272, 2
%293 = icmp eq i64 %292, %27
br i1 %293, label %294, label %271
294: ; preds = %268, %271, %251
%295 = add nuw nsw i64 %200, 1
%296 = icmp eq i64 %295, %30
%297 = add i32 %199, 1
%298 = add i64 %198, 1
br i1 %296, label %299, label %197
299: ; preds = %294, %181
%300 = add nuw nsw i64 %63, 1
%301 = icmp slt i64 %300, %25
br i1 %301, label %302, label %312
302: ; preds = %299, %302
%303 = phi i64 [ %310, %302 ], [ %61, %299 ]
%304 = sub nuw nsw i64 %303, %63
%305 = getelementptr inbounds double, double* %14, i64 %304
%306 = load double, double* %305
%307 = mul nsw i64 %303, %24
%308 = add nuw nsw i64 %307, %63
%309 = getelementptr inbounds double, double* %0, i64 %308
store double %306, double* %309
%310 = add nuw nsw i64 %303, 1
%311 = icmp eq i64 %310, %31
br i1 %311, label %312, label %302
312: ; preds = %302, %299
%313 = icmp slt i64 %63, %26
br i1 %313, label %314, label %822
314: ; preds = %312
%315 = icmp ult i64 %300, %24
br i1 %315, label %316, label %431
316: ; preds = %314
%317 = mul nsw i64 %63, %24
%318 = xor i64 %63, -1
%319 = icmp ult i64 %89, 4
br i1 %319, label %383, label %320
320: ; preds = %316
%321 = icmp ult i8* %12, %99
%322 = icmp ult i8* %95, %93
%323 = and i1 %321, %322
br i1 %323, label %383, label %324
324: ; preds = %320
%325 = and i64 %89, -4
%326 = add i64 %61, %325
%327 = and i64 %80, 1
%328 = icmp ult i64 %78, 4
br i1 %328, label %364, label %329
329: ; preds = %324
%330 = and i64 %80, 9223372036854775806
br label %331
331: ; preds = %331, %329
%332 = phi i64 [ 0, %329 ], [ %361, %331 ]
%333 = phi i64 [ 0, %329 ], [ %362, %331 ]
%334 = add i64 %61, %332
%335 = add nuw nsw i64 %334, %317
%336 = getelementptr inbounds double, double* %0, i64 %335
%337 = bitcast double* %336 to <2 x double>*
%338 = load <2 x double>, <2 x double>* %337
%339 = getelementptr inbounds double, double* %336, i64 2
%340 = bitcast double* %339 to <2 x double>*
%341 = load <2 x double>, <2 x double>* %340
%342 = add nsw i64 %334, %318
%343 = getelementptr inbounds double, double* %11, i64 %342
%344 = bitcast double* %343 to <2 x double>*
store <2 x double> %338, <2 x double>* %344
%345 = getelementptr inbounds double, double* %343, i64 2
%346 = bitcast double* %345 to <2 x double>*
store <2 x double> %341, <2 x double>* %346
%347 = or i64 %332, 4
%348 = add i64 %61, %347
%349 = add nuw nsw i64 %348, %317
%350 = getelementptr inbounds double, double* %0, i64 %349
%351 = bitcast double* %350 to <2 x double>*
%352 = load <2 x double>, <2 x double>* %351
%353 = getelementptr inbounds double, double* %350, i64 2
%354 = bitcast double* %353 to <2 x double>*
%355 = load <2 x double>, <2 x double>* %354
%356 = add nsw i64 %348, %318
%357 = getelementptr inbounds double, double* %11, i64 %356
%358 = bitcast double* %357 to <2 x double>*
store <2 x double> %352, <2 x double>* %358
%359 = getelementptr inbounds double, double* %357, i64 2
%360 = bitcast double* %359 to <2 x double>*
store <2 x double> %355, <2 x double>* %360
%361 = add nuw i64 %332, 8
%362 = add i64 %333, 2
%363 = icmp eq i64 %362, %330
br i1 %363, label %364, label %331
364: ; preds = %331, %324
%365 = phi i64 [ 0, %324 ], [ %361, %331 ]
%366 = icmp eq i64 %327, 0
br i1 %366, label %381, label %367
367: ; preds = %364
%368 = add i64 %61, %365
%369 = add nuw nsw i64 %368, %317
%370 = getelementptr inbounds double, double* %0, i64 %369
%371 = bitcast double* %370 to <2 x double>*
%372 = load <2 x double>, <2 x double>* %371
%373 = getelementptr inbounds double, double* %370, i64 2
%374 = bitcast double* %373 to <2 x double>*
%375 = load <2 x double>, <2 x double>* %374
%376 = add nsw i64 %368, %318
%377 = getelementptr inbounds double, double* %11, i64 %376
%378 = bitcast double* %377 to <2 x double>*
store <2 x double> %372, <2 x double>* %378
%379 = getelementptr inbounds double, double* %377, i64 2
%380 = bitcast double* %379 to <2 x double>*
store <2 x double> %375, <2 x double>* %380
br label %381
381: ; preds = %364, %367
%382 = icmp eq i64 %89, %325
br i1 %382, label %431, label %383
383: ; preds = %320, %316, %381
%384 = phi i64 [ %61, %320 ], [ %61, %316 ], [ %326, %381 ]
%385 = sub i64 %27, %384
%386 = xor i64 %384, -1
%387 = add i64 %386, %27
%388 = and i64 %385, 3
%389 = icmp eq i64 %388, 0
br i1 %389, label %401, label %390
390: ; preds = %383, %390
%391 = phi i64 [ %398, %390 ], [ %384, %383 ]
%392 = phi i64 [ %399, %390 ], [ 0, %383 ]
%393 = add nuw nsw i64 %391, %317
%394 = getelementptr inbounds double, double* %0, i64 %393
%395 = load double, double* %394
%396 = add nsw i64 %391, %318
%397 = getelementptr inbounds double, double* %11, i64 %396
store double %395, double* %397
%398 = add nuw nsw i64 %391, 1
%399 = add i64 %392, 1
%400 = icmp eq i64 %399, %388
br i1 %400, label %401, label %390
401: ; preds = %390, %383
%402 = phi i64 [ %384, %383 ], [ %398, %390 ]
%403 = icmp ult i64 %387, 3
br i1 %403, label %431, label %404
404: ; preds = %401, %404
%405 = phi i64 [ %429, %404 ], [ %402, %401 ]
%406 = add nuw nsw i64 %405, %317
%407 = getelementptr inbounds double, double* %0, i64 %406
%408 = load double, double* %407
%409 = add nsw i64 %405, %318
%410 = getelementptr inbounds double, double* %11, i64 %409
store double %408, double* %410
%411 = add nuw nsw i64 %405, 1
%412 = add nuw nsw i64 %411, %317
%413 = getelementptr inbounds double, double* %0, i64 %412
%414 = load double, double* %413
%415 = sub i64 %405, %63
%416 = getelementptr inbounds double, double* %11, i64 %415
store double %414, double* %416
%417 = add nuw nsw i64 %405, 2
%418 = add nuw nsw i64 %417, %317
%419 = getelementptr inbounds double, double* %0, i64 %418
%420 = load double, double* %419
%421 = add nsw i64 %417, %318
%422 = getelementptr inbounds double, double* %11, i64 %421
store double %420, double* %422
%423 = add nuw nsw i64 %405, 3
%424 = add nuw nsw i64 %423, %317
%425 = getelementptr inbounds double, double* %0, i64 %424
%426 = load double, double* %425
%427 = add nsw i64 %423, %318
%428 = getelementptr inbounds double, double* %11, i64 %427
store double %426, double* %428
%429 = add nuw nsw i64 %405, 4
%430 = icmp eq i64 %429, %27
br i1 %430, label %431, label %404
431: ; preds = %401, %404, %381, %314
%432 = add nsw i64 %184, -1
%433 = add nsw i32 %185, -1
%434 = trunc i64 %432 to i32
%435 = tail call double @house(double* %11, i32 %434, double* %14)
br i1 %142, label %436, label %695
436: ; preds = %431
%437 = xor i64 %63, -1
%438 = icmp ult i64 %101, 4
%439 = trunc i64 %109 to i32
%440 = icmp ugt i64 %109, 4294967295
%441 = and i64 %101, -4
%442 = add i64 %61, %441
%443 = and i64 %76, 1
%444 = icmp ult i64 %74, 4
%445 = and i64 %76, 9223372036854775806
%446 = icmp eq i64 %443, 0
%447 = icmp eq i64 %101, %441
br label %448
448: ; preds = %436, %573
%449 = phi i64 [ 0, %436 ], [ %577, %573 ]
%450 = phi i32 [ 0, %436 ], [ %576, %573 ]
%451 = phi i64 [ %63, %436 ], [ %574, %573 ]
%452 = mul i64 %449, %27
%453 = add i64 %452, 1
%454 = getelementptr double, double* %0, i64 %453
%455 = bitcast double* %454 to i8*
%456 = getelementptr i8, i8* %455, i64 %102
%457 = add i64 %104, %452
%458 = getelementptr double, double* %0, i64 %457
%459 = mul i64 %105, %449
%460 = shl i64 %459, 32
%461 = ashr exact i64 %460, 32
%462 = getelementptr double, double* %17, i64 %461
%463 = getelementptr double, double* %108, i64 %461
%464 = bitcast double* %463 to i8*
%465 = mul i32 %111, %450
%466 = mul nsw i64 %451, %24
br i1 %315, label %467, label %573
467: ; preds = %448
%468 = sub nuw nsw i64 %451, %63
%469 = mul nsw i64 %468, %432
%470 = add i64 %469, %437
br i1 %438, label %536, label %471
471: ; preds = %467
%472 = add i32 %465, %439
%473 = icmp slt i32 %472, %465
%474 = or i1 %473, %440
br i1 %474, label %536, label %475
475: ; preds = %471
%476 = icmp ult i8* %456, %464
%477 = icmp ult double* %462, %458
%478 = and i1 %476, %477
br i1 %478, label %536, label %479
479: ; preds = %475
br i1 %444, label %517, label %480
480: ; preds = %479, %480
%481 = phi i64 [ %514, %480 ], [ 0, %479 ]
%482 = phi i64 [ %515, %480 ], [ 0, %479 ]
%483 = add i64 %61, %481
%484 = add nuw nsw i64 %483, %466
%485 = getelementptr inbounds double, double* %0, i64 %484
%486 = bitcast double* %485 to <2 x double>*
%487 = load <2 x double>, <2 x double>* %486
%488 = getelementptr inbounds double, double* %485, i64 2
%489 = bitcast double* %488 to <2 x double>*
%490 = load <2 x double>, <2 x double>* %489
%491 = add i64 %470, %483
%492 = shl i64 %491, 32
%493 = ashr exact i64 %492, 32
%494 = getelementptr inbounds double, double* %17, i64 %493
%495 = bitcast double* %494 to <2 x double>*
store <2 x double> %487, <2 x double>* %495
%496 = getelementptr inbounds double, double* %494, i64 2
%497 = bitcast double* %496 to <2 x double>*
store <2 x double> %490, <2 x double>* %497
%498 = or i64 %481, 4
%499 = add i64 %61, %498
%500 = add nuw nsw i64 %499, %466
%501 = getelementptr inbounds double, double* %0, i64 %500
%502 = bitcast double* %501 to <2 x double>*
%503 = load <2 x double>, <2 x double>* %502
%504 = getelementptr inbounds double, double* %501, i64 2
%505 = bitcast double* %504 to <2 x double>*
%506 = load <2 x double>, <2 x double>* %505
%507 = add i64 %470, %499
%508 = shl i64 %507, 32
%509 = ashr exact i64 %508, 32
%510 = getelementptr inbounds double, double* %17, i64 %509
%511 = bitcast double* %510 to <2 x double>*
store <2 x double> %503, <2 x double>* %511
%512 = getelementptr inbounds double, double* %510, i64 2
%513 = bitcast double* %512 to <2 x double>*
store <2 x double> %506, <2 x double>* %513
%514 = add nuw i64 %481, 8
%515 = add i64 %482, 2
%516 = icmp eq i64 %515, %445
br i1 %516, label %517, label %480
517: ; preds = %480, %479
%518 = phi i64 [ 0, %479 ], [ %514, %480 ]
br i1 %446, label %535, label %519
519: ; preds = %517
%520 = add i64 %61, %518
%521 = add nuw nsw i64 %520, %466
%522 = getelementptr inbounds double, double* %0, i64 %521
%523 = bitcast double* %522 to <2 x double>*
%524 = load <2 x double>, <2 x double>* %523
%525 = getelementptr inbounds double, double* %522, i64 2
%526 = bitcast double* %525 to <2 x double>*
%527 = load <2 x double>, <2 x double>* %526
%528 = add i64 %470, %520
%529 = shl i64 %528, 32
%530 = ashr exact i64 %529, 32
%531 = getelementptr inbounds double, double* %17, i64 %530
%532 = bitcast double* %531 to <2 x double>*
store <2 x double> %524, <2 x double>* %532
%533 = getelementptr inbounds double, double* %531, i64 2
%534 = bitcast double* %533 to <2 x double>*
store <2 x double> %527, <2 x double>* %534
br label %535
535: ; preds = %517, %519
br i1 %447, label %573, label %536
536: ; preds = %475, %471, %467, %535
%537 = phi i64 [ %61, %475 ], [ %61, %471 ], [ %61, %467 ], [ %442, %535 ]
%538 = sub i64 %27, %537
%539 = xor i64 %537, -1
%540 = and i64 %538, 1
%541 = icmp eq i64 %540, 0
br i1 %541, label %551, label %542
542: ; preds = %536
%543 = add nuw nsw i64 %537, %466
%544 = getelementptr inbounds double, double* %0, i64 %543
%545 = load double, double* %544
%546 = add i64 %470, %537
%547 = shl i64 %546, 32
%548 = ashr exact i64 %547, 32
%549 = getelementptr inbounds double, double* %17, i64 %548
store double %545, double* %549
%550 = add nuw nsw i64 %537, 1
br label %551
551: ; preds = %542, %536
%552 = phi i64 [ %537, %536 ], [ %550, %542 ]
%553 = icmp eq i64 %539, %55
br i1 %553, label %573, label %554
554: ; preds = %551, %554
%555 = phi i64 [ %571, %554 ], [ %552, %551 ]
%556 = add nuw nsw i64 %555, %466
%557 = getelementptr inbounds double, double* %0, i64 %556
%558 = load double, double* %557
%559 = add i64 %470, %555
%560 = shl i64 %559, 32
%561 = ashr exact i64 %560, 32
%562 = getelementptr inbounds double, double* %17, i64 %561
store double %558, double* %562
%563 = add nuw nsw i64 %555, 1
%564 = add nuw nsw i64 %563, %466
%565 = getelementptr inbounds double, double* %0, i64 %564
%566 = load double, double* %565
%567 = add i64 %470, %563
%568 = shl i64 %567, 32
%569 = ashr exact i64 %568, 32
%570 = getelementptr inbounds double, double* %17, i64 %569
store double %566, double* %570
%571 = add nuw nsw i64 %555, 2
%572 = icmp eq i64 %571, %27
br i1 %572, label %573, label %554
573: ; preds = %551, %554, %535, %448
%574 = add nuw nsw i64 %451, 1
%575 = icmp eq i64 %574, %32
%576 = add i32 %450, 1
%577 = add i64 %449, 1
br i1 %575, label %578, label %448
578: ; preds = %573
%579 = tail call i32 @mmult(double* %17, double* %14, double* %20, i32 %183, i32 %433, i32 1)
%580 = tail call i32 @scale(double* %20, i32 %183, i32 1, double %435)
%581 = tail call i32 @mmult(double* %20, double* %14, double* %17, i32 %183, i32 1, i32 %433)
br i1 %142, label %582, label %699
582: ; preds = %578
%583 = xor i64 %63, -1
%584 = icmp ult i64 %113, 4
%585 = trunc i64 %121 to i32
%586 = icmp ugt i64 %121, 4294967295
%587 = and i64 %113, -4
%588 = add i64 %61, %587
%589 = icmp eq i64 %113, %587
br label %590
590: ; preds = %582, %690
%591 = phi i64 [ 0, %582 ], [ %694, %690 ]
%592 = phi i32 [ 0, %582 ], [ %693, %690 ]
%593 = phi i64 [ %63, %582 ], [ %691, %690 ]
%594 = mul i64 %591, %27
%595 = add i64 %594, 1
%596 = getelementptr double, double* %0, i64 %595
%597 = bitcast double* %596 to i8*
%598 = getelementptr i8, i8* %597, i64 %114
%599 = add i64 %116, %594
%600 = getelementptr double, double* %0, i64 %599
%601 = mul i64 %117, %591
%602 = shl i64 %601, 32
%603 = ashr exact i64 %602, 32
%604 = getelementptr double, double* %17, i64 %603
%605 = getelementptr double, double* %120, i64 %603
%606 = bitcast double* %605 to i8*
%607 = mul i32 %123, %592
%608 = mul nsw i64 %593, %24
br i1 %315, label %609, label %690
609: ; preds = %590
%610 = sub nuw nsw i64 %593, %63
%611 = mul nsw i64 %610, %432
%612 = add i64 %611, %583
br i1 %584, label %647, label %613
613: ; preds = %609
%614 = add i32 %607, %585
%615 = icmp slt i32 %614, %607
%616 = or i1 %615, %586
br i1 %616, label %647, label %617
617: ; preds = %613
%618 = icmp ult i8* %598, %606
%619 = icmp ult double* %604, %600
%620 = and i1 %618, %619
br i1 %620, label %647, label %621
621: ; preds = %617, %621
%622 = phi i64 [ %644, %621 ], [ 0, %617 ]
%623 = add i64 %61, %622
%624 = add i64 %612, %623
%625 = shl i64 %624, 32
%626 = ashr exact i64 %625, 32
%627 = getelementptr inbounds double, double* %17, i64 %626
%628 = bitcast double* %627 to <2 x double>*
%629 = load <2 x double>, <2 x double>* %628
%630 = getelementptr inbounds double, double* %627, i64 2
%631 = bitcast double* %630 to <2 x double>*
%632 = load <2 x double>, <2 x double>* %631
%633 = add nuw nsw i64 %623, %608
%634 = getelementptr inbounds double, double* %0, i64 %633
%635 = bitcast double* %634 to <2 x double>*
%636 = load <2 x double>, <2 x double>* %635
%637 = getelementptr inbounds double, double* %634, i64 2
%638 = bitcast double* %637 to <2 x double>*
%639 = load <2 x double>, <2 x double>* %638
%640 = fsub <2 x double> %636, %629
%641 = fsub <2 x double> %639, %632
%642 = bitcast double* %634 to <2 x double>*
store <2 x double> %640, <2 x double>* %642
%643 = bitcast double* %637 to <2 x double>*
store <2 x double> %641, <2 x double>* %643
%644 = add nuw i64 %622, 4
%645 = icmp eq i64 %644, %587
br i1 %645, label %646, label %621
646: ; preds = %621
br i1 %589, label %690, label %647
647: ; preds = %617, %613, %609, %646
%648 = phi i64 [ %61, %617 ], [ %61, %613 ], [ %61, %609 ], [ %588, %646 ]
%649 = sub i64 %27, %648
%650 = xor i64 %648, -1
%651 = and i64 %649, 1
%652 = icmp eq i64 %651, 0
br i1 %652, label %664, label %653
653: ; preds = %647
%654 = add i64 %612, %648
%655 = shl i64 %654, 32
%656 = ashr exact i64 %655, 32
%657 = getelementptr inbounds double, double* %17, i64 %656
%658 = load double, double* %657
%659 = add nuw nsw i64 %648, %608
%660 = getelementptr inbounds double, double* %0, i64 %659
%661 = load double, double* %660
%662 = fsub double %661, %658
store double %662, double* %660
%663 = add nuw nsw i64 %648, 1
br label %664
664: ; preds = %653, %647
%665 = phi i64 [ %648, %647 ], [ %663, %653 ]
%666 = icmp eq i64 %650, %56
br i1 %666, label %690, label %667
667: ; preds = %664, %667
%668 = phi i64 [ %688, %667 ], [ %665, %664 ]
%669 = add i64 %612, %668
%670 = shl i64 %669, 32
%671 = ashr exact i64 %670, 32
%672 = getelementptr inbounds double, double* %17, i64 %671
%673 = load double, double* %672
%674 = add nuw nsw i64 %668, %608
%675 = getelementptr inbounds double, double* %0, i64 %674
%676 = load double, double* %675
%677 = fsub double %676, %673
store double %677, double* %675
%678 = add nuw nsw i64 %668, 1
%679 = add i64 %612, %678
%680 = shl i64 %679, 32
%681 = ashr exact i64 %680, 32
%682 = getelementptr inbounds double, double* %17, i64 %681
%683 = load double, double* %682
%684 = add nuw nsw i64 %678, %608
%685 = getelementptr inbounds double, double* %0, i64 %684
%686 = load double, double* %685
%687 = fsub double %686, %683
store double %687, double* %685
%688 = add nuw nsw i64 %668, 2
%689 = icmp eq i64 %688, %27
br i1 %689, label %690, label %667
690: ; preds = %664, %667, %646, %590
%691 = add nuw nsw i64 %593, 1
%692 = icmp eq i64 %691, %33
%693 = add i32 %592, 1
%694 = add i64 %591, 1
br i1 %692, label %699, label %590
695: ; preds = %431
%696 = tail call i32 @mmult(double* %17, double* %14, double* %20, i32 %183, i32 %433, i32 1)
%697 = tail call i32 @scale(double* %20, i32 %183, i32 1, double %435)
%698 = tail call i32 @mmult(double* %20, double* %14, double* %17, i32 %183, i32 1, i32 %433)
br label %699
699: ; preds = %690, %695, %578
%700 = mul nsw i64 %63, %24
%701 = add nuw nsw i64 %63, 2
%702 = icmp ult i64 %701, %24
br i1 %702, label %703, label %822
703: ; preds = %699
%704 = zext i32 %138 to i64
%705 = icmp ult i64 %127, 4
br i1 %705, label %770, label %706
706: ; preds = %703
%707 = icmp ult double* %131, %137
%708 = icmp ult double* %36, %136
%709 = and i1 %707, %708
br i1 %709, label %770, label %710
710: ; preds = %706
%711 = and i64 %127, -4
%712 = or i64 %711, 1
%713 = add nsw i64 %711, %140
%714 = and i64 %72, 1
%715 = icmp ult i64 %70, 4
br i1 %715, label %751, label %716
716: ; preds = %710
%717 = and i64 %72, 9223372036854775806
br label %718
718: ; preds = %718, %716
%719 = phi i64 [ 0, %716 ], [ %748, %718 ]
%720 = phi i64 [ 0, %716 ], [ %749, %718 ]
%721 = or i64 %719, 1
%722 = add i64 %719, %140
%723 = getelementptr inbounds double, double* %14, i64 %721
%724 = bitcast double* %723 to <2 x double>*
%725 = load <2 x double>, <2 x double>* %724
%726 = getelementptr inbounds double, double* %723, i64 2
%727 = bitcast double* %726 to <2 x double>*
%728 = load <2 x double>, <2 x double>* %727
%729 = add nsw i64 %722, %700
%730 = getelementptr inbounds double, double* %0, i64 %729
%731 = bitcast double* %730 to <2 x double>*
store <2 x double> %725, <2 x double>* %731
%732 = getelementptr inbounds double, double* %730, i64 2
%733 = bitcast double* %732 to <2 x double>*
store <2 x double> %728, <2 x double>* %733
%734 = or i64 %719, 4
%735 = or i64 %719, 5
%736 = add i64 %734, %140
%737 = getelementptr inbounds double, double* %14, i64 %735
%738 = bitcast double* %737 to <2 x double>*
%739 = load <2 x double>, <2 x double>* %738
%740 = getelementptr inbounds double, double* %737, i64 2
%741 = bitcast double* %740 to <2 x double>*
%742 = load <2 x double>, <2 x double>* %741
%743 = add nsw i64 %736, %700
%744 = getelementptr inbounds double, double* %0, i64 %743
%745 = bitcast double* %744 to <2 x double>*
store <2 x double> %739, <2 x double>* %745
%746 = getelementptr inbounds double, double* %744, i64 2
%747 = bitcast double* %746 to <2 x double>*
store <2 x double> %742, <2 x double>* %747
%748 = add nuw i64 %719, 8
%749 = add i64 %720, 2
%750 = icmp eq i64 %749, %717
br i1 %750, label %751, label %718
751: ; preds = %718, %710
%752 = phi i64 [ 0, %710 ], [ %748, %718 ]
%753 = icmp eq i64 %714, 0
br i1 %753, label %768, label %754
754: ; preds = %751
%755 = or i64 %752, 1
%756 = add i64 %752, %140
%757 = getelementptr inbounds double, double* %14, i64 %755
%758 = bitcast double* %757 to <2 x double>*
%759 = load <2 x double>, <2 x double>* %758
%760 = getelementptr inbounds double, double* %757, i64 2
%761 = bitcast double* %760 to <2 x double>*
%762 = load <2 x double>, <2 x double>* %761
%763 = add nsw i64 %756, %700
%764 = getelementptr inbounds double, double* %0, i64 %763
%765 = bitcast double* %764 to <2 x double>*
store <2 x double> %759, <2 x double>* %765
%766 = getelementptr inbounds double, double* %764, i64 2
%767 = bitcast double* %766 to <2 x double>*
store <2 x double> %762, <2 x double>* %767
br label %768
768: ; preds = %751, %754
%769 = icmp eq i64 %127, %711
br i1 %769, label %822, label %770
770: ; preds = %706, %703, %768
%771 = phi i64 [ 1, %706 ], [ 1, %703 ], [ %712, %768 ]
%772 = phi i64 [ %140, %706 ], [ %140, %703 ], [ %713, %768 ]
%773 = sub nsw i64 %66, %771
%774 = xor i64 %771, -1
%775 = add nsw i64 %774, %66
%776 = and i64 %773, 3
%777 = icmp eq i64 %776, 0
br i1 %777, label %790, label %778
778: ; preds = %770, %778
%779 = phi i64 [ %786, %778 ], [ %771, %770 ]
%780 = phi i64 [ %787, %778 ], [ %772, %770 ]
%781 = phi i64 [ %788, %778 ], [ 0, %770 ]
%782 = getelementptr inbounds double, double* %14, i64 %779
%783 = load double, double* %782
%784 = add nsw i64 %780, %700
%785 = getelementptr inbounds double, double* %0, i64 %784
store double %783, double* %785
%786 = add nuw nsw i64 %779, 1
%787 = add nsw i64 %780, 1
%788 = add i64 %781, 1
%789 = icmp eq i64 %788, %776
br i1 %789, label %790, label %778
790: ; preds = %778, %770
%791 = phi i64 [ %771, %770 ], [ %786, %778 ]
%792 = phi i64 [ %772, %770 ], [ %787, %778 ]
%793 = icmp ult i64 %775, 3
br i1 %793, label %822, label %794
794: ; preds = %790, %794
%795 = phi i64 [ %819, %794 ], [ %791, %790 ]
%796 = phi i64 [ %820, %794 ], [ %792, %790 ]
%797 = getelementptr inbounds double, double* %14, i64 %795
%798 = load double, double* %797
%799 = add nsw i64 %796, %700
%800 = getelementptr inbounds double, double* %0, i64 %799
store double %798, double* %800
%801 = add nuw nsw i64 %795, 1
%802 = add nsw i64 %796, 1
%803 = getelementptr inbounds double, double* %14, i64 %801
%804 = load double, double* %803
%805 = add nsw i64 %802, %700
%806 = getelementptr inbounds double, double* %0, i64 %805
store double %804, double* %806
%807 = add nuw nsw i64 %795, 2
%808 = add nsw i64 %796, 2
%809 = getelementptr inbounds double, double* %14, i64 %807
%810 = load double, double* %809
%811 = add nsw i64 %808, %700
%812 = getelementptr inbounds double, double* %0, i64 %811
store double %810, double* %812
%813 = add nuw nsw i64 %795, 3
%814 = add nsw i64 %796, 3
%815 = getelementptr inbounds double, double* %14, i64 %813
%816 = load double, double* %815
%817 = add nsw i64 %814, %700
%818 = getelementptr inbounds double, double* %0, i64 %817
store double %816, double* %818
%819 = add nuw nsw i64 %795, 4
%820 = add nsw i64 %796, 4
%821 = icmp eq i64 %819, %704
br i1 %821, label %822, label %794
822: ; preds = %790, %794, %768, %699, %312
%823 = add nsw i64 %58, -1
%824 = add i32 %62, -1
%825 = add nuw nsw i64 %61, 1
%826 = add nuw i32 %60, 1
%827 = icmp eq i64 %300, %27
br i1 %827, label %828, label %57
828: ; preds = %822, %8
%829 = tail call i32 @free(double* %11)
%830 = tail call i32 @free(double* %14)
%831 = tail call i32 @free(double* %17)
ret void
}