#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
double calculate_column_norm(float *, int);
int /*<<< orphan*/ fprintf(int /*<<< orphan*/, char *);
int /*<<< orphan*/ free(float *);
scalar_t__ malloc(int);
float *project_vector(float *, float *, int);
int /*<<< orphan*/ stderr;
float ***qr_decomposition(float **mat, int rows, int cols) {
  float **q, **r;
  int i, j, k;
  float ***qr_mat = (float ***)malloc(sizeof(float **) * 2);
  if (qr_mat == NULL) {
    fprintf(stderr,
            "ERROR: Allocating memory for decomposition matrices failed\n");
    return NULL;
  }
  q = (float **)malloc(sizeof(float *) * rows);
  r = (float **)malloc(sizeof(float *) * rows);
  if (q == NULL || r == NULL) {
    fprintf(stderr,
            "ERROR: Allocating memory for decompositing matrix failed\n");
    return NULL;
  }
  for (i = 0; i < rows; i++) {
    q[i] = (float *)malloc(sizeof(float) * cols);
    r[i] = (float *)malloc(sizeof(float) * rows);
    if (q[i] == NULL || r[i] == NULL) {
      fprintf(stderr,
              "ERROR: Allocating memory for decomposition matrix failed\n");
      return NULL;
    }
  }
  qr_mat[0] = q;
  qr_mat[1] = r;
  for (i = 0; i < cols; i++) {
    float *of, *on;
    of = (float *)malloc(sizeof(float) * rows);
    on = (float *)malloc(sizeof(float) * rows);
    if (of == NULL || on == NULL) {
      fprintf(stderr,
              "ERROR: Calculating the projection failed. Out of memory\n");
      return NULL;
    }
    if (i == 0) {
      for (j = 0; j < rows; j++)
        of[j] = mat[j][i];
      double norm = calculate_column_norm(of, rows);
      for (j = 0; j < rows; j++)
        q[j][i] = mat[j][i] / norm;
    } else {
      for (j = 0; j < rows; j++) {
        of[j] = mat[j][i];
        q[j][i] = mat[j][i];
      }
      for (k = 0; k < i; k++) {
        float *proj;
        for (j = 0; j < rows; j++)
          on[j] = q[j][k];
        proj = project_vector(on, of, rows);
        if (proj == NULL) {
          fprintf(stderr, "ERROR: Calculating the projection failed.\n");
          return NULL;
        }
        for (j = 0; j < rows; j++)
          q[j][i] -= proj[j];
        free(proj);
      }
      for (j = 0; j < rows; j++)
        of[j] = q[j][i];
      double norm = calculate_column_norm(of, rows);
      for (j = 0; j < rows; j++)
        q[j][i] /= norm;
    }
    free(on);
    free(of);
  }
  for (i = 0; i < rows; i++)
    for (j = 0; j < rows; j++) {
      r[i][j] = 0;
      for (k = 0; k < rows; k++)
        r[i][j] += (q[k][i] * mat[k][j]);
    }
  return qr_mat;
}


x86 assembly (clang, -O0):
<qr_decomposition>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    %edx,-0x18(%rbp)
mov    $0x10,%edi
call   1070 <malloc@plt>
mov    %rax,-0x40(%rbp)
cmpq   $0x0,-0x40(%rbp)
jne    119d <qr_decomposition+0x4d>
mov    0x2e5e(%rip),%rax
mov    (%rax),%edi
lea    0xe75(%rip),%rsi
call   1060 <fprintf@plt>
movq   $0x0,-0x8(%rbp)
jmp    16ed <qr_decomposition+0x59d>
movslq -0x14(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   1070 <malloc@plt>
mov    %rax,-0x20(%rbp)
movslq -0x14(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   1070 <malloc@plt>
mov    %rax,-0x28(%rbp)
cmpq   $0x0,-0x20(%rbp)
je     11d9 <qr_decomposition+0x89>
cmpq   $0x0,-0x28(%rbp)
jne    11fb <qr_decomposition+0xab>
mov    0x2e00(%rip),%rax
mov    (%rax),%edi
lea    0xe53(%rip),%rsi
call   1060 <fprintf@plt>
movq   $0x0,-0x8(%rbp)
jmp    16ed <qr_decomposition+0x59d>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    12a5 <qr_decomposition+0x155>
movslq -0x18(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1070 <malloc@plt>
mov    %rax,%rdx
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
movslq -0x14(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1070 <malloc@plt>
mov    %rax,%rdx
mov    -0x28(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
cmpq   $0x0,(%rax,%rcx,8)
je     1270 <qr_decomposition+0x120>
mov    -0x28(%rbp),%rax
movslq -0x2c(%rbp),%rcx
cmpq   $0x0,(%rax,%rcx,8)
jne    1292 <qr_decomposition+0x142>
mov    0x2d69(%rip),%rax
mov    (%rax),%edi
lea    0xdf6(%rip),%rsi
call   1060 <fprintf@plt>
movq   $0x0,-0x8(%rbp)
jmp    16ed <qr_decomposition+0x59d>
jmp    1297 <qr_decomposition+0x147>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1202 <qr_decomposition+0xb2>
mov    -0x20(%rbp),%rcx
mov    -0x40(%rbp),%rax
mov    %rcx,(%rax)
mov    -0x28(%rbp),%rcx
mov    -0x40(%rbp),%rax
mov    %rcx,0x8(%rax)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    1614 <qr_decomposition+0x4c4>
movslq -0x14(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1070 <malloc@plt>
mov    %rax,-0x48(%rbp)
movslq -0x14(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1070 <malloc@plt>
mov    %rax,-0x50(%rbp)
cmpq   $0x0,-0x48(%rbp)
je     130b <qr_decomposition+0x1bb>
cmpq   $0x0,-0x50(%rbp)
jne    132d <qr_decomposition+0x1dd>
mov    0x2cce(%rip),%rax
mov    (%rax),%edi
lea    0xd95(%rip),%rsi
call   1060 <fprintf@plt>
movq   $0x0,-0x8(%rbp)
jmp    16ed <qr_decomposition+0x59d>
cmpl   $0x0,-0x2c(%rbp)
jne    13e8 <qr_decomposition+0x298>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    137a <qr_decomposition+0x22a>
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x48(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    133e <qr_decomposition+0x1ee>
mov    -0x48(%rbp),%rdi
mov    -0x14(%rbp),%esi
call   1040 <calculate_column_norm@plt>
movsd  %xmm0,-0x58(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    13e3 <qr_decomposition+0x293>
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
divsd  -0x58(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1392 <qr_decomposition+0x242>
jmp    15f4 <qr_decomposition+0x4a4>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1455 <qr_decomposition+0x305>
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x48(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    13ef <qr_decomposition+0x29f>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    154f <qr_decomposition+0x3ff>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    14ab <qr_decomposition+0x35b>
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x50(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    146f <qr_decomposition+0x31f>
mov    -0x50(%rbp),%rdi
mov    -0x48(%rbp),%rsi
mov    -0x14(%rbp),%edx
call   1050 <project_vector@plt>
mov    %rax,-0x60(%rbp)
cmpq   $0x0,-0x60(%rbp)
jne    14ec <qr_decomposition+0x39c>
mov    0x2b0f(%rip),%rax
mov    (%rax),%edi
lea    0xc0f(%rip),%rsi
call   1060 <fprintf@plt>
movq   $0x0,-0x8(%rbp)
jmp    16ed <qr_decomposition+0x59d>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1538 <qr_decomposition+0x3e8>
mov    -0x60(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    14f3 <qr_decomposition+0x3a3>
mov    -0x60(%rbp),%rdi
call   1030 <free@plt>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    145c <qr_decomposition+0x30c>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1592 <qr_decomposition+0x442>
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x48(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1556 <qr_decomposition+0x406>
mov    -0x48(%rbp),%rdi
mov    -0x14(%rbp),%esi
call   1040 <calculate_column_norm@plt>
movsd  %xmm0,-0x68(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    15ef <qr_decomposition+0x49f>
movsd  -0x68(%rbp),%xmm1
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    15aa <qr_decomposition+0x45a>
jmp    15f4 <qr_decomposition+0x4a4>
mov    -0x50(%rbp),%rdi
call   1030 <free@plt>
mov    -0x48(%rbp),%rdi
call   1030 <free@plt>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    12c3 <qr_decomposition+0x173>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    16e5 <qr_decomposition+0x595>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    16d2 <qr_decomposition+0x582>
mov    -0x28(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    16bf <qr_decomposition+0x56f>
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x28(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1659 <qr_decomposition+0x509>
jmp    16c4 <qr_decomposition+0x574>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    162e <qr_decomposition+0x4de>
jmp    16d7 <qr_decomposition+0x587>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    161b <qr_decomposition+0x4cb>
mov    -0x40(%rbp),%rax
mov    %rax,-0x8(%rbp)
mov    -0x8(%rbp),%rax
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<qr_decomposition>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x68,%rsp
mov    %edx,%r13d
mov    %esi,%ebx
mov    %rdi,%r12
mov    $0x10,%edi
call   1070 <malloc@plt>
test   %rax,%rax
je     181a <qr_decomposition+0x6ca>
mov    %rax,%r14
mov    %rbx,0x10(%rsp)
lea    0x0(,%rbx,8),%ebp
mov    %ebp,%edi
call   1070 <malloc@plt>
mov    %rax,%r15
mov    %ebp,%edi
call   1070 <malloc@plt>
lea    0xe9c(%rip),%rsi
test   %r15,%r15
je     1833 <qr_decomposition+0x6e3>
mov    %rax,%rcx
test   %rax,%rax
je     1833 <qr_decomposition+0x6e3>
mov    %r13,0x18(%rsp)
mov    %r14,0x38(%rsp)
cmpl   $0x0,0x10(%rsp)
jle    122c <qr_decomposition+0xdc>
mov    0x18(%rsp),%rax
lea    0x0(,%rax,4),%eax
mov    %eax,0x8(%rsp)
mov    0x10(%rsp),%rax
lea    0x0(,%rax,4),%r13d
mov    %eax,%r14d
xor    %ebp,%ebp
nopl   0x0(%rax,%rax,1)
mov    0x8(%rsp),%edi
mov    %rcx,%rbx
call   1070 <malloc@plt>
mov    %rax,(%r15,%rbp,8)
mov    %r13d,%edi
call   1070 <malloc@plt>
mov    %rax,(%rbx,%rbp,8)
cmpq   $0x0,(%r15,%rbp,8)
je     182c <qr_decomposition+0x6dc>
test   %rax,%rax
je     182c <qr_decomposition+0x6dc>
mov    %rbx,%rcx
add    $0x1,%rbp
cmp    %rbp,%r14
jne    11f0 <qr_decomposition+0xa0>
mov    0x38(%rsp),%rax
mov    %r15,(%rax)
mov    %rcx,0x58(%rsp)
mov    %rcx,0x8(%rax)
mov    0x18(%rsp),%rcx
test   %ecx,%ecx
jle    1852 <qr_decomposition+0x702>
mov    0x10(%rsp),%rax
lea    0x0(,%rax,4),%ebx
mov    %ecx,%ecx
mov    %rcx,0x60(%rsp)
mov    %eax,%eax
lea    -0x1(%rax),%rcx
mov    %rcx,0x8(%rsp)
mov    %eax,%ecx
and    $0xfffffffe,%ecx
mov    %rcx,0x28(%rsp)
mov    %eax,%ecx
and    $0x3,%ecx
mov    %rcx,0x30(%rsp)
mov    %rax,0x20(%rsp)
and    $0xfffffffc,%eax
mov    %rax,0x40(%rsp)
xor    %r13d,%r13d
mov    %ebx,0x4c(%rsp)
cmp    $0x2,%r13d
mov    $0x1,%ebp
cmovae %r13d,%ebp
mov    %ebx,%edi
call   1070 <malloc@plt>
mov    %rax,%r14
mov    %ebx,%edi
call   1070 <malloc@plt>
mov    %r14,%r8
test   %r14,%r14
je     1932 <qr_decomposition+0x7e2>
mov    %rax,%r14
test   %rax,%rax
je     1932 <qr_decomposition+0x7e2>
test   %r13,%r13
mov    %r8,0x18(%rsp)
je     153f <qr_decomposition+0x3ef>
mov    0x10(%rsp),%rdx
test   %edx,%edx
mov    0x30(%rsp),%rbx
jle    1373 <qr_decomposition+0x223>
cmpq   $0x0,0x8(%rsp)
je     155c <qr_decomposition+0x40c>
xor    %eax,%eax
mov    0x28(%rsp),%rsi
nopl   0x0(%rax)
mov    (%r12,%rax,8),%rcx
movss  (%rcx,%r13,4),%xmm0
movss  %xmm0,(%r8,%rax,4)
movss  (%rcx,%r13,4),%xmm0
mov    (%r15,%rax,8),%rcx
movss  %xmm0,(%rcx,%r13,4)
mov    0x8(%r12,%rax,8),%rcx
movss  (%rcx,%r13,4),%xmm0
movss  %xmm0,0x4(%r8,%rax,4)
movss  (%rcx,%r13,4),%xmm0
mov    0x8(%r15,%rax,8),%rcx
movss  %xmm0,(%rcx,%r13,4)
add    $0x2,%rax
cmp    %rax,%rsi
jne    1300 <qr_decomposition+0x1b0>
testb  $0x1,0x20(%rsp)
je     1373 <qr_decomposition+0x223>
mov    (%r12,%rax,8),%rcx
movss  (%rcx,%r13,4),%xmm0
movss  %xmm0,(%r8,%rax,4)
movss  (%rcx,%r13,4),%xmm0
mov    (%r15,%rax,8),%rax
movss  %xmm0,(%rax,%r13,4)
mov    %ebp,%eax
mov    %rax,0x50(%rsp)
xor    %ebp,%ebp
jmp    13a1 <qr_decomposition+0x251>
xchg   %ax,%ax
mov    %rax,%rdi
call   1030 <free@plt>
add    $0x1,%rbp
cmp    0x50(%rsp),%rbp
mov    0x10(%rsp),%rdx
mov    0x18(%rsp),%r8
je     156e <qr_decomposition+0x41e>
test   %edx,%edx
jle    13c0 <qr_decomposition+0x270>
cmpq   $0x3,0x8(%rsp)
jae    13e0 <qr_decomposition+0x290>
xor    %eax,%eax
jmp    143b <qr_decomposition+0x2eb>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    %r14,%rdi
mov    %r8,%rsi
call   1050 <project_vector@plt>
test   %rax,%rax
jne    14f5 <qr_decomposition+0x3a5>
jmp    1823 <qr_decomposition+0x6d3>
nopl   0x0(%rax)
xor    %eax,%eax
mov    0x40(%rsp),%rsi
nopw   0x0(%rax,%rax,1)
mov    (%r15,%rax,8),%rcx
movss  (%rcx,%rbp,4),%xmm0
movss  %xmm0,(%r14,%rax,4)
mov    0x8(%r15,%rax,8),%rcx
movss  (%rcx,%rbp,4),%xmm0
movss  %xmm0,0x4(%r14,%rax,4)
mov    0x10(%r15,%rax,8),%rcx
movss  (%rcx,%rbp,4),%xmm0
movss  %xmm0,0x8(%r14,%rax,4)
mov    0x18(%r15,%rax,8),%rcx
movss  (%rcx,%rbp,4),%xmm0
movss  %xmm0,0xc(%r14,%rax,4)
add    $0x4,%rax
cmp    %rax,%rsi
jne    13f0 <qr_decomposition+0x2a0>
test   %rbx,%rbx
je     1467 <qr_decomposition+0x317>
lea    (%r14,%rax,4),%rcx
lea    (%r15,%rax,8),%rax
xor    %edi,%edi
nopw   0x0(%rax,%rax,1)
mov    (%rax,%rdi,8),%rsi
movss  (%rsi,%rbp,4),%xmm0
movss  %xmm0,(%rcx,%rdi,4)
add    $0x1,%rdi
cmp    %rdi,%rbx
jne    1450 <qr_decomposition+0x300>
mov    %r14,%rdi
mov    %r8,%rsi
mov    %rdx,%rbx
call   1050 <project_vector@plt>
test   %rax,%rax
je     1823 <qr_decomposition+0x6d3>
test   %ebx,%ebx
jle    14f0 <qr_decomposition+0x3a0>
cmpq   $0x0,0x8(%rsp)
mov    0x28(%rsp),%rsi
mov    0x30(%rsp),%rbx
je     1518 <qr_decomposition+0x3c8>
xor    %ecx,%ecx
nopw   0x0(%rax,%rax,1)
mov    (%r15,%rcx,8),%rdx
movss  (%rdx,%r13,4),%xmm0
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rdx,%r13,4)
mov    0x8(%r15,%rcx,8),%rdx
movss  (%rdx,%r13,4),%xmm0
subss  0x4(%rax,%rcx,4),%xmm0
movss  %xmm0,(%rdx,%r13,4)
add    $0x2,%rcx
cmp    %rcx,%rsi
jne    14a0 <qr_decomposition+0x350>
testb  $0x1,0x20(%rsp)
je     1380 <qr_decomposition+0x230>
jmp    1525 <qr_decomposition+0x3d5>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    0x30(%rsp),%rbx
mov    %rax,%rdi
call   1030 <free@plt>
add    $0x1,%rbp
cmp    0x50(%rsp),%rbp
mov    0x10(%rsp),%rdx
mov    0x18(%rsp),%r8
jne    13a1 <qr_decomposition+0x251>
jmp    1585 <qr_decomposition+0x435>
xor    %ecx,%ecx
testb  $0x1,0x20(%rsp)
je     1380 <qr_decomposition+0x230>
mov    (%r15,%rcx,8),%rdx
movss  (%rdx,%r13,4),%xmm0
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rdx,%r13,4)
jmp    1380 <qr_decomposition+0x230>
mov    0x10(%rsp),%rdx
test   %edx,%edx
mov    0x30(%rsp),%rbp
jle    1585 <qr_decomposition+0x435>
cmpq   $0x3,0x8(%rsp)
jae    15ba <qr_decomposition+0x46a>
xor    %eax,%eax
jmp    1617 <qr_decomposition+0x4c7>
xor    %eax,%eax
testb  $0x1,0x20(%rsp)
jne    1353 <qr_decomposition+0x203>
jmp    1373 <qr_decomposition+0x223>
test   %edx,%edx
jle    1585 <qr_decomposition+0x435>
cmpq   $0x3,0x8(%rsp)
jae    16eb <qr_decomposition+0x59b>
xor    %eax,%eax
jmp    174f <qr_decomposition+0x5ff>
mov    %r8,%rdi
mov    %edx,%esi
mov    %r8,%rbp
call   1040 <calculate_column_norm@plt>
mov    %r14,%rdi
call   1030 <free@plt>
mov    %rbp,%rdi
call   1030 <free@plt>
add    $0x1,%r13
cmp    0x60(%rsp),%r13
mov    0x4c(%rsp),%ebx
jne    1290 <qr_decomposition+0x140>
jmp    1852 <qr_decomposition+0x702>
xor    %eax,%eax
mov    0x40(%rsp),%rsi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    (%r12,%rax,8),%rcx
movss  (%rcx),%xmm0
movss  %xmm0,(%r8,%rax,4)
mov    0x8(%r12,%rax,8),%rcx
movss  (%rcx),%xmm0
movss  %xmm0,0x4(%r8,%rax,4)
mov    0x10(%r12,%rax,8),%rcx
movss  (%rcx),%xmm0
movss  %xmm0,0x8(%r8,%rax,4)
mov    0x18(%r12,%rax,8),%rcx
movss  (%rcx),%xmm0
movss  %xmm0,0xc(%r8,%rax,4)
add    $0x4,%rax
cmp    %rax,%rsi
jne    15d0 <qr_decomposition+0x480>
testb  $0x3,0x20(%rsp)
je     1646 <qr_decomposition+0x4f6>
lea    (%r8,%rax,4),%rcx
lea    (%r12,%rax,8),%rax
xor    %edi,%edi
nopl   0x0(%rax,%rax,1)
mov    (%rax,%rdi,8),%rsi
movss  (%rsi),%xmm0
movss  %xmm0,(%rcx,%rdi,4)
add    $0x1,%rdi
cmp    %rdi,%rbp
jne    1630 <qr_decomposition+0x4e0>
mov    %r8,%rdi
mov    %edx,%esi
mov    %rdx,%rbx
mov    %r8,%rbp
call   1040 <calculate_column_norm@plt>
test   %ebx,%ebx
jle    1592 <qr_decomposition+0x442>
cmpq   $0x0,0x8(%rsp)
je     16b5 <qr_decomposition+0x565>
xor    %eax,%eax
mov    0x28(%rsp),%rdx
nopl   (%rax)
mov    (%r12,%rax,8),%rcx
movss  (%rcx),%xmm1
cvtss2sd %xmm1,%xmm1
divsd  %xmm0,%xmm1
cvtsd2ss %xmm1,%xmm1
mov    (%r15,%rax,8),%rcx
movss  %xmm1,(%rcx)
mov    0x8(%r12,%rax,8),%rcx
movss  (%rcx),%xmm1
cvtss2sd %xmm1,%xmm1
divsd  %xmm0,%xmm1
cvtsd2ss %xmm1,%xmm1
mov    0x8(%r15,%rax,8),%rcx
movss  %xmm1,(%rcx)
add    $0x2,%rax
cmp    %rax,%rdx
jne    1670 <qr_decomposition+0x520>
jmp    16b7 <qr_decomposition+0x567>
xor    %eax,%eax
testb  $0x1,0x20(%rsp)
mov    0x18(%rsp),%rbp
je     1592 <qr_decomposition+0x442>
mov    (%r12,%rax,8),%rcx
movss  (%rcx),%xmm1
cvtss2sd %xmm1,%xmm1
divsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
mov    (%r15,%rax,8),%rax
movss  %xmm0,(%rax)
jmp    1592 <qr_decomposition+0x442>
xor    %eax,%eax
mov    0x40(%rsp),%rsi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    (%r15,%rax,8),%rcx
movss  (%rcx,%r13,4),%xmm0
movss  %xmm0,(%r8,%rax,4)
mov    0x8(%r15,%rax,8),%rcx
movss  (%rcx,%r13,4),%xmm0
movss  %xmm0,0x4(%r8,%rax,4)
mov    0x10(%r15,%rax,8),%rcx
movss  (%rcx,%r13,4),%xmm0
movss  %xmm0,0x8(%r8,%rax,4)
mov    0x18(%r15,%rax,8),%rcx
movss  (%rcx,%r13,4),%xmm0
movss  %xmm0,0xc(%r8,%rax,4)
add    $0x4,%rax
cmp    %rax,%rsi
jne    1700 <qr_decomposition+0x5b0>
testb  $0x3,0x20(%rsp)
je     1778 <qr_decomposition+0x628>
lea    (%r8,%rax,4),%rcx
lea    (%r15,%rax,8),%rax
xor    %edi,%edi
mov    (%rax,%rdi,8),%rsi
movss  (%rsi,%r13,4),%xmm0
movss  %xmm0,(%rcx,%rdi,4)
add    $0x1,%rdi
cmp    %rdi,%rbx
jne    1760 <qr_decomposition+0x610>
mov    %r8,%rdi
mov    %edx,%esi
mov    %rdx,%rbx
mov    %r8,%rbp
call   1040 <calculate_column_norm@plt>
test   %ebx,%ebx
jle    1592 <qr_decomposition+0x442>
cmpq   $0x0,0x8(%rsp)
je     17e4 <qr_decomposition+0x694>
xor    %eax,%eax
mov    0x28(%rsp),%rdx
nop
mov    (%r15,%rax,8),%rcx
movss  (%rcx,%r13,4),%xmm1
cvtss2sd %xmm1,%xmm1
divsd  %xmm0,%xmm1
cvtsd2ss %xmm1,%xmm1
movss  %xmm1,(%rcx,%r13,4)
mov    0x8(%r15,%rax,8),%rcx
movss  (%rcx,%r13,4),%xmm1
cvtss2sd %xmm1,%xmm1
divsd  %xmm0,%xmm1
cvtsd2ss %xmm1,%xmm1
movss  %xmm1,(%rcx,%r13,4)
add    $0x2,%rax
cmp    %rax,%rdx
jne    17a0 <qr_decomposition+0x650>
jmp    17e6 <qr_decomposition+0x696>
xor    %eax,%eax
testb  $0x1,0x20(%rsp)
mov    0x18(%rsp),%rbp
je     1592 <qr_decomposition+0x442>
mov    (%r15,%rax,8),%rax
movss  (%rax,%r13,4),%xmm1
cvtss2sd %xmm1,%xmm1
divsd  %xmm0,%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
movss  %xmm0,(%rax,%r13,4)
jmp    1592 <qr_decomposition+0x442>
lea    0x7df(%rip),%rsi
jmp    1833 <qr_decomposition+0x6e3>
lea    0x8bf(%rip),%rsi
jmp    1833 <qr_decomposition+0x6e3>
lea    0x843(%rip),%rsi
mov    0x27a6(%rip),%rax
mov    (%rax),%edi
call   1060 <fprintf@plt>
xor    %eax,%eax
add    $0x68,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    0x10(%rsp),%r9
test   %r9d,%r9d
jle    1928 <qr_decomposition+0x7d8>
mov    %r9d,%r8d
mov    %r8d,%ecx
and    $0xfffffffe,%ecx
xor    %edx,%edx
mov    0x38(%rsp),%rax
mov    0x58(%rsp),%r10
jmp    1889 <qr_decomposition+0x739>
nopw   0x0(%rax,%rax,1)
add    $0x1,%rdx
cmp    %r8,%rdx
je     1843 <qr_decomposition+0x6f3>
mov    (%r10,%rdx,8),%rsi
xor    %edi,%edi
jmp    18a9 <qr_decomposition+0x759>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
add    $0x1,%rdi
cmp    %r8,%rdi
je     1880 <qr_decomposition+0x730>
movl   $0x0,(%rsi,%rdi,4)
xorps  %xmm0,%xmm0
xor    %ebp,%ebp
cmp    $0x1,%r9d
je     1901 <qr_decomposition+0x7b1>
nopl   0x0(%rax,%rax,1)
mov    (%r15,%rbp,8),%rbx
movss  (%rbx,%rdx,4),%xmm1
mov    (%r12,%rbp,8),%rbx
mulss  (%rbx,%rdi,4),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,(%rsi,%rdi,4)
mov    0x8(%r15,%rbp,8),%rbx
movss  (%rbx,%rdx,4),%xmm0
mov    0x8(%r12,%rbp,8),%rbx
mulss  (%rbx,%rdi,4),%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rsi,%rdi,4)
add    $0x2,%rbp
cmp    %rbp,%rcx
jne    18c0 <qr_decomposition+0x770>
test   $0x1,%r8b
je     18a0 <qr_decomposition+0x750>
mov    (%r15,%rbp,8),%rbx
movss  (%rbx,%rdx,4),%xmm1
mov    (%r12,%rbp,8),%rbp
mulss  0x0(%rbp,%rdi,4),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,(%rsi,%rdi,4)
jmp    18a0 <qr_decomposition+0x750>
mov    0x38(%rsp),%rax
jmp    1843 <qr_decomposition+0x6f3>
lea    0x777(%rip),%rsi
jmp    1833 <qr_decomposition+0x6e3>

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

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