#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 */
struct timeval {
  double tv_sec;
  double tv_usec;
};

/* Variables and functions */
int atoi(char *);
int /*<<< orphan*/ choleskyblk(double *, int, int, int *);
int /*<<< orphan*/ choleskytile(double *, int, int, int *);
int /*<<< orphan*/ convtile(double *, double *, int, int, int);
int /*<<< orphan*/ dgemv_(char *, int *, int *, double const *, double *, int *,
                          double *, int const *, double const *, double *,
                          int const *, int);
int /*<<< orphan*/ dpotrs_(char *, int *, int const *, double *, int *,
                           double *, int *, int *, int);
double drand48();
int /*<<< orphan*/ exit(int);
scalar_t__ fabs(double);
int /*<<< orphan*/ free(double *);
int /*<<< orphan*/ gettimeofday(struct timeval *, int /*<<< orphan*/ *);
scalar_t__ malloc(int);
int /*<<< orphan*/ printf(char *, ...);
int main(int argc, char **argv) {
  int i, j, n = 640, bs = 64, ft = 0, info;
  double *A, *B, *x, *b, err;
  const double done = 1.0, dzero = 0.0;
  const int one = 1;
  struct timeval t0, t1;
  double time;
  if (argc >= 2) {
    if ((n = atoi(argv[1])) < 0)
      n = 640;
  }
  if (argc >= 3) {
    if ((bs = atoi(argv[2])) < 0)
      bs = 64;
  }
  if (n < bs || n % bs) {
    printf("ERROR: el tamaño de la matriz ha de ser divisible por bs\n");
    exit(1);
  }
  if (argc >= 4)
    ft = atoi(argv[3]);
  A = (double *)malloc(n * n * sizeof(double));
  B = (double *)malloc(n * n * sizeof(double));
  x = (double *)malloc(n * sizeof(double));
  b = (double *)malloc(n * sizeof(double));
  for (i = 0; i < n; i++) {
    for (j = 0; j < n; j++) {
      A[i + n * j] = drand48();
      A[j + n * i] = A[i + n * j];
    }
  }
  for (i = 0; i < n; i++)
    A[i + i * n] += (double)n;
  for (i = 0; i < n; i++)
    x[i] = 1.0;
  dgemv_("N", &n, &n, &done, A, &n, x, &one, &dzero, b, &one, 1);
  if (ft)
    choleskyblk(A, n, bs, &info);
  else {
    convtile(A, B, n, bs, 0);
    gettimeofday(&t0, NULL);
    choleskytile(B, n, bs, &info);
    gettimeofday(&t1, NULL);
    convtile(B, A, n, bs, 1);
  }
  if (info != 0) {
    printf("ERROR: falló la factorización de Cholesky\n");
    exit(1);
  }
  dpotrs_("L", &n, &one, A, &n, b, &n, &info, 1);
  if (info != 0) {
    printf("ERROR: falló la resolución triangular\n");
    exit(1);
  }
  time = (t1.tv_sec - t0.tv_sec) + (t1.tv_usec - t0.tv_usec) / 1000000.0;
  printf("%d;%d;%f;", n, bs, time);
  err = 0.0;
  for (i = 0; i < n; i++)
    err += fabs(1.0 - b[i]);
  printf("%g\n", err);
  free(A);
  free(B);
  free(x);
  free(b);
  return 0;
}


x86 assembly (clang, -O0):
<main>:
push   %rbp
mov    %rsp,%rbp
push   %rbx
sub    $0xc8,%rsp
movl   $0x0,-0xc(%rbp)
mov    %edi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
movl   $0x280,-0x24(%rbp)
movl   $0x40,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
movsd  0xdf9(%rip),%xmm0
movsd  %xmm0,-0x60(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x68(%rbp)
movl   $0x1,-0x6c(%rbp)
cmpl   $0x2,-0x10(%rbp)
jl     124a <main+0x7a>
mov    -0x18(%rbp),%rax
mov    0x8(%rax),%rdi
call   10c0 <atoi@plt>
mov    %eax,-0x24(%rbp)
cmp    $0x0,%eax
jge    1245 <main+0x75>
movl   $0x280,-0x24(%rbp)
jmp    124a <main+0x7a>
cmpl   $0x3,-0x10(%rbp)
jl     1279 <main+0xa9>
mov    -0x18(%rbp),%rax
mov    0x10(%rax),%rdi
call   10c0 <atoi@plt>
mov    %eax,-0x28(%rbp)
cmp    $0x0,%eax
jge    1274 <main+0xa4>
movl   $0x40,-0x28(%rbp)
jmp    1279 <main+0xa9>
mov    -0x24(%rbp),%eax
cmp    -0x28(%rbp),%eax
jl     1295 <main+0xc5>
mov    -0x24(%rbp),%eax
cltd
idivl  -0x28(%rbp)
cmp    $0x0,%edx
je     12ad <main+0xdd>
lea    0xd74(%rip),%rdi
mov    $0x0,%al
call   1050 <printf@plt>
mov    $0x1,%edi
call   10d0 <exit@plt>
cmpl   $0x4,-0x10(%rbp)
jl     12c7 <main+0xf7>
mov    -0x18(%rbp),%rax
mov    0x18(%rax),%rdi
call   10c0 <atoi@plt>
mov    %eax,-0x2c(%rbp)
mov    -0x24(%rbp),%eax
imul   -0x24(%rbp),%eax
cltq
shl    $0x3,%rax
mov    %eax,%edi
call   10a0 <malloc@plt>
mov    %rax,-0x38(%rbp)
mov    -0x24(%rbp),%eax
imul   -0x24(%rbp),%eax
cltq
shl    $0x3,%rax
mov    %eax,%edi
call   10a0 <malloc@plt>
mov    %rax,-0x40(%rbp)
movslq -0x24(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   10a0 <malloc@plt>
mov    %rax,-0x48(%rbp)
movslq -0x24(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   10a0 <malloc@plt>
mov    %rax,-0x50(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    13b3 <main+0x1e3>
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    13a0 <main+0x1d0>
mov    $0x0,%al
call   10f0 <drand48@plt>
mov    -0x38(%rbp),%rax
mov    -0x1c(%rbp),%ecx
mov    -0x24(%rbp),%edx
imul   -0x20(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%rax
mov    -0x1c(%rbp),%ecx
mov    -0x24(%rbp),%edx
imul   -0x20(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x38(%rbp),%rax
mov    -0x20(%rbp),%ecx
mov    -0x24(%rbp),%edx
imul   -0x1c(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1337 <main+0x167>
jmp    13a5 <main+0x1d5>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1324 <main+0x154>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    13f6 <main+0x226>
cvtsi2sdl -0x24(%rbp),%xmm0
mov    -0x38(%rbp),%rax
mov    -0x1c(%rbp),%ecx
mov    -0x1c(%rbp),%edx
imul   -0x24(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    13ba <main+0x1ea>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    142c <main+0x25c>
mov    -0x48(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  0xbe7(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    13fd <main+0x22d>
mov    -0x38(%rbp),%r8
mov    -0x48(%rbp),%rbx
mov    -0x50(%rbp),%r10
lea    0xc0c(%rip),%rdi
lea    -0x24(%rbp),%r9
lea    -0x60(%rbp),%rcx
lea    -0x6c(%rbp),%rax
lea    -0x68(%rbp),%r11
mov    %r9,%rsi
mov    %r9,%rdx
mov    %rbx,(%rsp)
mov    %rax,0x8(%rsp)
mov    %r11,0x10(%rsp)
mov    %r10,0x18(%rsp)
mov    %rax,0x20(%rsp)
movl   $0x1,0x28(%rsp)
call   10b0 <dgemv_@plt>
cmpl   $0x0,-0x2c(%rbp)
je     149c <main+0x2cc>
mov    -0x38(%rbp),%rdi
mov    -0x24(%rbp),%esi
mov    -0x28(%rbp),%edx
lea    -0x30(%rbp),%rcx
call   1040 <choleskyblk@plt>
jmp    14fb <main+0x32b>
mov    -0x38(%rbp),%rdi
mov    -0x40(%rbp),%rsi
mov    -0x24(%rbp),%edx
mov    -0x28(%rbp),%ecx
xor    %r8d,%r8d
call   10e0 <convtile@plt>
lea    -0x80(%rbp),%rdi
xor    %eax,%eax
mov    %eax,%esi
call   1060 <gettimeofday@plt>
mov    -0x40(%rbp),%rdi
mov    -0x24(%rbp),%esi
mov    -0x28(%rbp),%edx
lea    -0x30(%rbp),%rcx
call   1090 <choleskytile@plt>
lea    -0x90(%rbp),%rdi
xor    %eax,%eax
mov    %eax,%esi
call   1060 <gettimeofday@plt>
mov    -0x40(%rbp),%rdi
mov    -0x38(%rbp),%rsi
mov    -0x24(%rbp),%edx
mov    -0x28(%rbp),%ecx
mov    $0x1,%r8d
call   10e0 <convtile@plt>
cmpl   $0x0,-0x30(%rbp)
je     151d <main+0x34d>
lea    0xb41(%rip),%rdi
mov    $0x0,%al
call   1050 <printf@plt>
mov    $0x1,%edi
call   10d0 <exit@plt>
mov    -0x38(%rbp),%rcx
mov    -0x50(%rbp),%r9
lea    0xb4e(%rip),%rdi
lea    -0x24(%rbp),%r10
lea    -0x6c(%rbp),%rdx
lea    -0x30(%rbp),%rax
mov    %r10,%rsi
mov    %r10,%r8
mov    %r10,(%rsp)
mov    %rax,0x8(%rsp)
movl   $0x1,0x10(%rsp)
call   1080 <dpotrs_@plt>
cmpl   $0x0,-0x30(%rbp)
je     1576 <main+0x3a6>
lea    0xb17(%rip),%rdi
mov    $0x0,%al
call   1050 <printf@plt>
mov    $0x1,%edi
call   10d0 <exit@plt>
movsd  -0x90(%rbp),%xmm0
ff
subsd  -0x80(%rbp),%xmm0
movsd  -0x88(%rbp),%xmm1
ff
subsd  -0x78(%rbp),%xmm1
movsd  0xa70(%rip),%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
mov    -0x24(%rbp),%esi
mov    -0x28(%rbp),%edx
movsd  -0x98(%rbp),%xmm0
ff
lea    0xae8(%rip),%rdi
mov    $0x1,%al
call   1050 <printf@plt>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1616 <main+0x446>
mov    -0x50(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  0xa11(%rip),%xmm0
subsd  (%rax,%rcx,8),%xmm0
call   1070 <fabs@plt>
cvtsi2sd %rax,%xmm0
addsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    15d3 <main+0x403>
movsd  -0x58(%rbp),%xmm0
lea    0xa8d(%rip),%rdi
mov    $0x1,%al
call   1050 <printf@plt>
mov    -0x38(%rbp),%rdi
call   1030 <free@plt>
mov    -0x40(%rbp),%rdi
call   1030 <free@plt>
mov    -0x48(%rbp),%rdi
call   1030 <free@plt>
mov    -0x50(%rbp),%rdi
call   1030 <free@plt>
xor    %eax,%eax
add    $0xc8,%rsp
pop    %rbx
pop    %rbp
ret

x86 assembly (clang, -O2):
<main>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x88,%rsp
mov    %rsi,%rbx
mov    %edi,%ebp
movl   $0x280,0x8(%rsp)
movabs $0x3ff0000000000000,%rax
mov    %rax,0x58(%rsp)
movq   $0x0,0x50(%rsp)
movl   $0x1,0x1c(%rsp)
cmp    $0x2,%edi
jl     126a <main+0x8a>
mov    0x8(%rbx),%rdi
call   10d0 <atoi@plt>
test   %eax,%eax
mov    $0x280,%esi
cmovns %eax,%esi
mov    %esi,0x8(%rsp)
mov    $0x40,%r14d
cmp    $0x2,%ebp
je     1260 <main+0x80>
mov    0x10(%rbx),%rdi
mov    %rsi,%r14
call   10d0 <atoi@plt>
mov    %r14,%rsi
test   %eax,%eax
mov    $0x40,%r14d
cmovns %eax,%r14d
cmp    %r14d,%esi
jae    1275 <main+0x95>
jmp    173d <main+0x55d>
mov    $0x280,%esi
mov    $0x40,%r14d
xor    %ecx,%ecx
mov    %esi,%eax
xor    %edx,%edx
div    %r14d
test   %edx,%edx
jne    173d <main+0x55d>
cmp    $0x4,%ebp
mov    %rsi,0x20(%rsp)
jl     12a0 <main+0xc0>
mov    0x18(%rbx),%rdi
call   10d0 <atoi@plt>
mov    0x20(%rsp),%rsi
mov    %eax,%ecx
mov    %ecx,0x2c(%rsp)
lea    0x0(,%rsi,8),%ebx
mov    %ebx,%ebp
imul   %esi,%ebp
mov    %ebp,%edi
call   10b0 <malloc@plt>
mov    %rax,%r12
mov    %ebp,%edi
call   10b0 <malloc@plt>
mov    %rax,%r15
mov    %ebx,%edi
call   10b0 <malloc@plt>
mov    %rax,%r13
mov    %ebx,%edi
call   10b0 <malloc@plt>
mov    0x20(%rsp),%rcx
mov    %rax,0x30(%rsp)
test   %ecx,%ecx
je     156e <main+0x38e>
mov    %r13,0x38(%rsp)
mov    %r15,0x40(%rsp)
mov    %r14d,0x28(%rsp)
mov    %ecx,%ebx
lea    0x0(,%rbx,8),%r15
mov    %r12,%r13
mov    %r12,0x48(%rsp)
mov    %r12,%rbp
xor    %r14d,%r14d
xchg   %ax,%ax
mov    %r13,0x10(%rsp)
xor    %r12d,%r12d
nopl   0x0(%rax,%rax,1)
xor    %eax,%eax
call   1100 <drand48@plt>
movsd  %xmm0,0x0(%r13)
movsd  %xmm0,0x0(%rbp,%r12,8)
add    $0x1,%r12
add    %r15,%r13
cmp    %r12,%rbx
jne    1320 <main+0x140>
add    $0x1,%r14
add    %r15,%rbp
mov    0x10(%rsp),%r13
add    $0x8,%r13
cmp    %rbx,%r14
jne    1310 <main+0x130>
mov    0x20(%rsp),%rbp
test   %ebp,%ebp
mov    0x28(%rsp),%r14d
mov    0x48(%rsp),%r12
mov    0x40(%rsp),%r15
mov    0x38(%rsp),%r13
movabs $0x3ff0000000000000,%r11
je     156e <main+0x38e>
xorps  %xmm0,%xmm0
cvtsi2sd %ebp,%xmm0
xor    %eax,%eax
cmp    $0x1,%ebp
jne    13a4 <main+0x1c4>
test   $0x1,%bl
jne    1418 <main+0x238>
test   %ebp,%ebp
jne    1437 <main+0x257>
jmp    156e <main+0x38e>
lea    0x0(%rbp,%rbp,1),%r9d
mov    %ebx,%r8d
and    $0xfffffffe,%r8d
neg    %r8
mov    $0x1,%edi
xor    %esi,%esi
xor    %edx,%edx
nopl   0x0(%rax)
movslq %ebp,%rbp
lea    (%rdi,%rbp,1),%r10
cltq
lea    (%rdx,%rax,1),%rcx
movsd  (%r12,%rcx,8),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%r12,%rcx,8)
movsd  (%r12,%r10,8),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%r12,%r10,8)
add    $0x2,%rdx
add    %r9d,%eax
add    $0x2,%rdi
add    %r9d,%ebp
add    $0xfffffffffffffffe,%rsi
cmp    %rsi,%r8
jne    13c0 <main+0x1e0>
neg    %rsi
mov    %rsi,%rax
mov    0x20(%rsp),%rbp
test   $0x1,%bl
je     1397 <main+0x1b7>
mov    %ebp,%ecx
imul   %eax,%ecx
movslq %ecx,%rcx
add    %rax,%rcx
addsd  (%r12,%rcx,8),%xmm0
movsd  %xmm0,(%r12,%rcx,8)
test   %ebp,%ebp
je     156e <main+0x38e>
cmp    $0x4,%ebp
jae    1443 <main+0x263>
xor    %eax,%eax
jmp    1560 <main+0x380>
mov    %ebx,%eax
and    $0xfffffffc,%eax
lea    -0x4(%rax),%rsi
mov    %rsi,%rdx
shr    $0x2,%rdx
add    $0x1,%rdx
mov    %edx,%ecx
and    $0x7,%ecx
cmp    $0x1c,%rsi
jae    1469 <main+0x289>
xor    %esi,%esi
jmp    1516 <main+0x336>
and    $0xfffffffffffffff8,%rdx
xor    %esi,%esi
movapd 0xb89(%rip),%xmm0
nopw   0x0(%rax,%rax,1)
movupd %xmm0,0x0(%r13,%rsi,8)
movupd %xmm0,0x10(%r13,%rsi,8)
movupd %xmm0,0x20(%r13,%rsi,8)
movupd %xmm0,0x30(%r13,%rsi,8)
movupd %xmm0,0x40(%r13,%rsi,8)
movupd %xmm0,0x50(%r13,%rsi,8)
movupd %xmm0,0x60(%r13,%rsi,8)
movupd %xmm0,0x70(%r13,%rsi,8)
movupd %xmm0,0x80(%r13,%rsi,8)
movupd %xmm0,0x90(%r13,%rsi,8)
movupd %xmm0,0xa0(%r13,%rsi,8)
movupd %xmm0,0xb0(%r13,%rsi,8)
movupd %xmm0,0xc0(%r13,%rsi,8)
movupd %xmm0,0xd0(%r13,%rsi,8)
movupd %xmm0,0xe0(%r13,%rsi,8)
movupd %xmm0,0xf0(%r13,%rsi,8)
add    $0x20,%rsi
add    $0xfffffffffffffff8,%rdx
jne    1480 <main+0x2a0>
test   %rcx,%rcx
je     1554 <main+0x374>
lea    0x10(,%rsi,8),%rdx
add    %r13,%rdx
shl    $0x5,%rcx
xor    %esi,%esi
movapd 0xacc(%rip),%xmm0
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movupd %xmm0,-0x10(%rdx,%rsi,1)
movupd %xmm0,(%rdx,%rsi,1)
add    $0x20,%rsi
cmp    %rsi,%rcx
jne    1540 <main+0x360>
cmp    %rbx,%rax
je     156e <main+0x38e>
nopl   0x0(%rax)
mov    %r11,0x0(%r13,%rax,8)
add    $0x1,%rax
cmp    %rax,%rbx
jne    1560 <main+0x380>
lea    0x1c(%rsp),%rax
lea    0x50(%rsp),%rbp
lea    0xaa1(%rip),%rdi
lea    0x58(%rsp),%rcx
lea    0x8(%rsp),%rbx
mov    %rbx,%rsi
mov    %rbx,%rdx
mov    %r12,%r8
mov    %rbx,%r9
push   $0x1
push   %rax
push   0x40(%rsp)
push   %rbp
push   %rax
push   %r13
call   10c0 <dgemv_@plt>
add    $0x30,%rsp
mov    0x8(%rsp),%edx
cmpl   $0x0,0x2c(%rsp)
je     15c8 <main+0x3e8>
lea    0xc(%rsp),%rcx
mov    %r12,%rdi
mov    %edx,%esi
mov    %r14d,%edx
call   1050 <choleskyblk@plt>
jmp    161d <main+0x43d>
mov    %r12,%rdi
mov    %r15,%rsi
mov    %r14d,%ecx
xor    %r8d,%r8d
call   10f0 <convtile@plt>
lea    0x70(%rsp),%rdi
xor    %esi,%esi
call   1070 <gettimeofday@plt>
mov    0x8(%rsp),%esi
lea    0xc(%rsp),%rcx
mov    %r15,%rdi
mov    %r14d,%edx
call   10a0 <choleskytile@plt>
lea    0x60(%rsp),%rdi
xor    %esi,%esi
call   1070 <gettimeofday@plt>
mov    0x8(%rsp),%edx
mov    %r15,%rdi
mov    %r12,%rsi
mov    %r14d,%ecx
mov    $0x1,%r8d
call   10f0 <convtile@plt>
cmpl   $0x0,0xc(%rsp)
mov    0x30(%rsp),%rbp
jne    1746 <main+0x566>
sub    $0x8,%rsp
lea    0x14(%rsp),%rax
lea    0x9e5(%rip),%rdi
lea    0x24(%rsp),%rdx
mov    %rbx,%rsi
mov    %r12,%rcx
mov    %rbx,%r8
mov    %rbp,%r9
push   $0x1
push   %rax
push   %rbx
call   1090 <dpotrs_@plt>
add    $0x20,%rsp
cmpl   $0x0,0xc(%rsp)
jne    174f <main+0x56f>
movapd 0x60(%rsp),%xmm1
subpd  0x70(%rsp),%xmm1
movapd %xmm1,%xmm0
unpckhpd %xmm1,%xmm0
divsd  0x98e(%rip),%xmm0
addsd  %xmm1,%xmm0
mov    0x8(%rsp),%esi
lea    0x993(%rip),%rdi
mov    %r14d,%edx
mov    $0x1,%al
call   1060 <printf@plt>
cmpl   $0x0,0x8(%rsp)
jle    16f7 <main+0x517>
xorpd  %xmm0,%xmm0
xor    %ebx,%ebx
nopl   0x0(%rax,%rax,1)
movsd  %xmm0,0x10(%rsp)
movsd  0x95a(%rip),%xmm0
subsd  0x0(%rbp,%rbx,8),%xmm0
call   1080 <fabs@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %rax,%xmm0
movsd  0x10(%rsp),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,0x10(%rsp)
movsd  0x10(%rsp),%xmm0
add    $0x1,%rbx
movslq 0x8(%rsp),%rax
cmp    %rax,%rbx
jl     16b0 <main+0x4d0>
jmp    16fb <main+0x51b>
xorpd  %xmm0,%xmm0
lea    0x92c(%rip),%rdi
mov    $0x1,%al
call   1060 <printf@plt>
mov    %r12,%rdi
call   1030 <free@plt>
mov    %r15,%rdi
call   1030 <free@plt>
mov    %r13,%rdi
call   1030 <free@plt>
mov    %rbp,%rdi
call   1030 <free@plt>
xor    %eax,%eax
add    $0x88,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
lea    0x942(%rip),%rdi
jmp    1756 <main+0x576>
lea    0x90d(%rip),%rdi
jmp    1756 <main+0x576>
lea    0x8dc(%rip),%rdi
call   1040 <puts@plt>
mov    $0x1,%edi
call   10e0 <exit@plt>

LLVM IR (-O0):
%struct.timeval = type { double, double }
@.str = [59 x i8] c"ERROR: el tama\C3\B1o de la matriz ha de ser divisible por bs\0A\00"
@.str.1 = [2 x i8] c"N\00"
@.str.2 = [45 x i8] c"ERROR: fall\C3\B3 la factorizaci\C3\B3n de Cholesky\0A\00"
@.str.3 = [2 x i8] c"L\00"
@.str.4 = [41 x i8] c"ERROR: fall\C3\B3 la resoluci\C3\B3n triangular\0A\00"
@.str.5 = [10 x i8] c"%d;%d;%f;\00"
@.str.6 = [4 x i8] c"%g\0A\00"
define i32 @main(i32 %0, i8** %1) {
%3 = alloca i32
%4 = alloca i32
%5 = alloca i8**
%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
%17 = alloca double
%18 = alloca double
%19 = alloca i32
%20 = alloca %struct.timeval
%21 = alloca %struct.timeval
%22 = alloca double
store i32 0, i32* %3
store i32 %0, i32* %4
store i8** %1, i8*** %5
store i32 640, i32* %8
store i32 64, i32* %9
store i32 0, i32* %10
store double 1.000000e+00, double* %17
store double 0.000000e+00, double* %18
store i32 1, i32* %19
%23 = load i32, i32* %4
%24 = icmp sge i32 %23, 2
br i1 %24, label %25, label %33
25: ; preds = %2
%26 = load i8**, i8*** %5
%27 = getelementptr inbounds i8*, i8** %26, i64 1
%28 = load i8*, i8** %27
%29 = call i32 @atoi(i8* %28)
store i32 %29, i32* %8
%30 = icmp slt i32 %29, 0
br i1 %30, label %31, label %32
31: ; preds = %25
store i32 640, i32* %8
br label %32
32: ; preds = %31, %25
br label %33
33: ; preds = %32, %2
%34 = load i32, i32* %4
%35 = icmp sge i32 %34, 3
br i1 %35, label %36, label %44
36: ; preds = %33
%37 = load i8**, i8*** %5
%38 = getelementptr inbounds i8*, i8** %37, i64 2
%39 = load i8*, i8** %38
%40 = call i32 @atoi(i8* %39)
store i32 %40, i32* %9
%41 = icmp slt i32 %40, 0
br i1 %41, label %42, label %43
42: ; preds = %36
store i32 64, i32* %9
br label %43
43: ; preds = %42, %36
br label %44
44: ; preds = %43, %33
%45 = load i32, i32* %8
%46 = load i32, i32* %9
%47 = icmp slt i32 %45, %46
br i1 %47, label %53, label %48
48: ; preds = %44
%49 = load i32, i32* %8
%50 = load i32, i32* %9
%51 = srem i32 %49, %50
%52 = icmp ne i32 %51, 0
br i1 %52, label %53, label %56
53: ; preds = %48, %44
%54 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([59 x i8], [59 x i8]* @.str, i64 0, i64 0))
%55 = call i32 @exit(i32 1)
unreachable
56: ; preds = %48
%57 = load i32, i32* %4
%58 = icmp sge i32 %57, 4
br i1 %58, label %59, label %64
59: ; preds = %56
%60 = load i8**, i8*** %5
%61 = getelementptr inbounds i8*, i8** %60, i64 3
%62 = load i8*, i8** %61
%63 = call i32 @atoi(i8* %62)
store i32 %63, i32* %10
br label %64
64: ; preds = %59, %56
%65 = load i32, i32* %8
%66 = load i32, i32* %8
%67 = mul nsw i32 %65, %66
%68 = sext i32 %67 to i64
%69 = mul i64 %68, 8
%70 = trunc i64 %69 to i32
%71 = call i64 @malloc(i32 %70)
%72 = inttoptr i64 %71 to double*
store double* %72, double** %12
%73 = load i32, i32* %8
%74 = load i32, i32* %8
%75 = mul nsw i32 %73, %74
%76 = sext i32 %75 to i64
%77 = mul i64 %76, 8
%78 = trunc i64 %77 to i32
%79 = call i64 @malloc(i32 %78)
%80 = inttoptr i64 %79 to double*
store double* %80, double** %13
%81 = load i32, i32* %8
%82 = sext i32 %81 to i64
%83 = mul i64 %82, 8
%84 = trunc i64 %83 to i32
%85 = call i64 @malloc(i32 %84)
%86 = inttoptr i64 %85 to double*
store double* %86, double** %14
%87 = load i32, i32* %8
%88 = sext i32 %87 to i64
%89 = mul i64 %88, 8
%90 = trunc i64 %89 to i32
%91 = call i64 @malloc(i32 %90)
%92 = inttoptr i64 %91 to double*
store double* %92, double** %15
store i32 0, i32* %6
br label %93
93: ; preds = %133, %64
%94 = load i32, i32* %6
%95 = load i32, i32* %8
%96 = icmp slt i32 %94, %95
br i1 %96, label %97, label %136
97: ; preds = %93
store i32 0, i32* %7
br label %98
98: ; preds = %129, %97
%99 = load i32, i32* %7
%100 = load i32, i32* %8
%101 = icmp slt i32 %99, %100
br i1 %101, label %102, label %132
102: ; preds = %98
%103 = call double (...) @drand48()
%104 = load double*, double** %12
%105 = load i32, i32* %6
%106 = load i32, i32* %8
%107 = load i32, i32* %7
%108 = mul nsw i32 %106, %107
%109 = add nsw i32 %105, %108
%110 = sext i32 %109 to i64
%111 = getelementptr inbounds double, double* %104, i64 %110
store double %103, double* %111
%112 = load double*, double** %12
%113 = load i32, i32* %6
%114 = load i32, i32* %8
%115 = load i32, i32* %7
%116 = mul nsw i32 %114, %115
%117 = add nsw i32 %113, %116
%118 = sext i32 %117 to i64
%119 = getelementptr inbounds double, double* %112, i64 %118
%120 = load double, double* %119
%121 = load double*, double** %12
%122 = load i32, i32* %7
%123 = load i32, i32* %8
%124 = load i32, i32* %6
%125 = mul nsw i32 %123, %124
%126 = add nsw i32 %122, %125
%127 = sext i32 %126 to i64
%128 = getelementptr inbounds double, double* %121, i64 %127
store double %120, double* %128
br label %129
129: ; preds = %102
%130 = load i32, i32* %7
%131 = add nsw i32 %130, 1
store i32 %131, i32* %7
br label %98
132: ; preds = %98
br label %133
133: ; preds = %132
%134 = load i32, i32* %6
%135 = add nsw i32 %134, 1
store i32 %135, i32* %6
br label %93
136: ; preds = %93
store i32 0, i32* %6
br label %137
137: ; preds = %154, %136
%138 = load i32, i32* %6
%139 = load i32, i32* %8
%140 = icmp slt i32 %138, %139
br i1 %140, label %141, label %157
141: ; preds = %137
%142 = load i32, i32* %8
%143 = sitofp i32 %142 to double
%144 = load double*, double** %12
%145 = load i32, i32* %6
%146 = load i32, i32* %6
%147 = load i32, i32* %8
%148 = mul nsw i32 %146, %147
%149 = add nsw i32 %145, %148
%150 = sext i32 %149 to i64
%151 = getelementptr inbounds double, double* %144, i64 %150
%152 = load double, double* %151
%153 = fadd double %152, %143
store double %153, double* %151
br label %154
154: ; preds = %141
%155 = load i32, i32* %6
%156 = add nsw i32 %155, 1
store i32 %156, i32* %6
br label %137
157: ; preds = %137
store i32 0, i32* %6
br label %158
158: ; preds = %167, %157
%159 = load i32, i32* %6
%160 = load i32, i32* %8
%161 = icmp slt i32 %159, %160
br i1 %161, label %162, label %170
162: ; preds = %158
%163 = load double*, double** %14
%164 = load i32, i32* %6
%165 = sext i32 %164 to i64
%166 = getelementptr inbounds double, double* %163, i64 %165
store double 1.000000e+00, double* %166
br label %167
167: ; preds = %162
%168 = load i32, i32* %6
%169 = add nsw i32 %168, 1
store i32 %169, i32* %6
br label %158
170: ; preds = %158
%171 = load double*, double** %12
%172 = load double*, double** %14
%173 = load double*, double** %15
%174 = call i32 @dgemv_(i8* getelementptr inbounds ([2 x i8], [2 x i8]* @.str.1, i64 0, i64 0), i32* %8, i32* %8, double* %17, double* %171, i32* %8, double* %172, i32* %19, double* %18, double* %173, i32* %19, i32 1)
%175 = load i32, i32* %10
%176 = icmp ne i32 %175, 0
br i1 %176, label %177, label %182
177: ; preds = %170
%178 = load double*, double** %12
%179 = load i32, i32* %8
%180 = load i32, i32* %9
%181 = call i32 @choleskyblk(double* %178, i32 %179, i32 %180, i32* %11)
br label %199
182: ; preds = %170
%183 = load double*, double** %12
%184 = load double*, double** %13
%185 = load i32, i32* %8
%186 = load i32, i32* %9
%187 = call i32 @convtile(double* %183, double* %184, i32 %185, i32 %186, i32 0)
%188 = call i32 @gettimeofday(%struct.timeval* %20, i32* null)
%189 = load double*, double** %13
%190 = load i32, i32* %8
%191 = load i32, i32* %9
%192 = call i32 @choleskytile(double* %189, i32 %190, i32 %191, i32* %11)
%193 = call i32 @gettimeofday(%struct.timeval* %21, i32* null)
%194 = load double*, double** %13
%195 = load double*, double** %12
%196 = load i32, i32* %8
%197 = load i32, i32* %9
%198 = call i32 @convtile(double* %194, double* %195, i32 %196, i32 %197, i32 1)
br label %199
199: ; preds = %182, %177
%200 = load i32, i32* %11
%201 = icmp ne i32 %200, 0
br i1 %201, label %202, label %205
202: ; preds = %199
%203 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([45 x i8], [45 x i8]* @.str.2, i64 0, i64 0))
%204 = call i32 @exit(i32 1)
unreachable
205: ; preds = %199
%206 = load double*, double** %12
%207 = load double*, double** %15
%208 = call i32 @dpotrs_(i8* getelementptr inbounds ([2 x i8], [2 x i8]* @.str.3, i64 0, i64 0), i32* %8, i32* %19, double* %206, i32* %8, double* %207, i32* %8, i32* %11, i32 1)
%209 = load i32, i32* %11
%210 = icmp ne i32 %209, 0
br i1 %210, label %211, label %214
211: ; preds = %205
%212 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([41 x i8], [41 x i8]* @.str.4, i64 0, i64 0))
%213 = call i32 @exit(i32 1)
unreachable
214: ; preds = %205
%215 = getelementptr inbounds %struct.timeval, %struct.timeval* %21, i32 0, i32 0
%216 = load double, double* %215
%217 = getelementptr inbounds %struct.timeval, %struct.timeval* %20, i32 0, i32 0
%218 = load double, double* %217
%219 = fsub double %216, %218
%220 = getelementptr inbounds %struct.timeval, %struct.timeval* %21, i32 0, i32 1
%221 = load double, double* %220
%222 = getelementptr inbounds %struct.timeval, %struct.timeval* %20, i32 0, i32 1
%223 = load double, double* %222
%224 = fsub double %221, %223
%225 = fdiv double %224, 1.000000e+06
%226 = fadd double %219, %225
store double %226, double* %22
%227 = load i32, i32* %8
%228 = load i32, i32* %9
%229 = load double, double* %22
%230 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([10 x i8], [10 x i8]* @.str.5, i64 0, i64 0), i32 %227, i32 %228, double %229)
store double 0.000000e+00, double* %16
store i32 0, i32* %6
br label %231
231: ; preds = %246, %214
%232 = load i32, i32* %6
%233 = load i32, i32* %8
%234 = icmp slt i32 %232, %233
br i1 %234, label %235, label %249
235: ; preds = %231
%236 = load double*, double** %15
%237 = load i32, i32* %6
%238 = sext i32 %237 to i64
%239 = getelementptr inbounds double, double* %236, i64 %238
%240 = load double, double* %239
%241 = fsub double 1.000000e+00, %240
%242 = call i64 @fabs(double %241)
%243 = sitofp i64 %242 to double
%244 = load double, double* %16
%245 = fadd double %244, %243
store double %245, double* %16
br label %246
246: ; preds = %235
%247 = load i32, i32* %6
%248 = add nsw i32 %247, 1
store i32 %248, i32* %6
br label %231
249: ; preds = %231
%250 = load double, double* %16
%251 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([4 x i8], [4 x i8]* @.str.6, i64 0, i64 0), double %250)
%252 = load double*, double** %12
%253 = call i32 @free(double* %252)
%254 = load double*, double** %13
%255 = call i32 @free(double* %254)
%256 = load double*, double** %14
%257 = call i32 @free(double* %256)
%258 = load double*, double** %15
%259 = call i32 @free(double* %258)
ret i32 0
}

LLVM IR (-O2):
%struct.timeval = type { double, double }
@.str.1 = [2 x i8] c"N\00"
@.str.3 = [2 x i8] c"L\00"
@.str.5 = [10 x i8] c"%d;%d;%f;\00"
@.str.6 = [4 x i8] c"%g\0A\00"
@str = [40 x i8] c"ERROR: fall\C3\B3 la resoluci\C3\B3n triangular\00"
@str.7 = [44 x i8] c"ERROR: fall\C3\B3 la factorizaci\C3\B3n de Cholesky\00"
@str.8 = [58 x i8] c"ERROR: el tama\C3\B1o de la matriz ha de ser divisible por bs\00"
define i32 @main(i32 %0, i8** nocapture %1) {
%3 = alloca i32
%4 = alloca i32
%5 = alloca double
%6 = alloca double
%7 = alloca i32
%8 = alloca <2 x double>
%9 = bitcast <2 x double>* %8 to %struct.timeval*
%10 = alloca <2 x double>
%11 = bitcast <2 x double>* %10 to %struct.timeval*
%12 = bitcast i32* %3 to i8*
store i32 640, i32* %3
%13 = bitcast i32* %4 to i8*
%14 = bitcast double* %5 to i8*
store double 1.000000e+00, double* %5
%15 = bitcast double* %6 to i8*
store double 0.000000e+00, double* %6
%16 = bitcast i32* %7 to i8*
store i32 1, i32* %7
%17 = bitcast <2 x double>* %8 to i8*
%18 = bitcast <2 x double>* %10 to i8*
%19 = icmp sgt i32 %0, 1
br i1 %19, label %20, label %36
20: ; preds = %2
%21 = getelementptr inbounds i8*, i8** %1, i64 1
%22 = load i8*, i8** %21
%23 = tail call i32 @atoi(i8* %22)
%24 = icmp slt i32 %23, 0
%25 = select i1 %24, i32 640, i32 %23
store i32 %25, i32* %3
%26 = icmp eq i32 %0, 2
br i1 %26, label %33, label %27
27: ; preds = %20
%28 = getelementptr inbounds i8*, i8** %1, i64 2
%29 = load i8*, i8** %28
%30 = tail call i32 @atoi(i8* %29)
%31 = icmp slt i32 %30, 0
%32 = select i1 %31, i32 64, i32 %30
br label %33
33: ; preds = %27, %20
%34 = phi i32 [ %32, %27 ], [ 64, %20 ]
%35 = icmp ult i32 %25, %34
br i1 %35, label %41, label %36
36: ; preds = %2, %33
%37 = phi i32 [ %34, %33 ], [ 64, %2 ]
%38 = phi i32 [ %25, %33 ], [ 640, %2 ]
%39 = urem i32 %38, %37
%40 = icmp eq i32 %39, 0
br i1 %40, label %44, label %41
41: ; preds = %36, %33
%42 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([58 x i8], [58 x i8]* @str.8, i64 0, i64 0))
%43 = tail call i32 @exit(i32 1)
unreachable
44: ; preds = %36
%45 = icmp sgt i32 %0, 3
br i1 %45, label %46, label %50
46: ; preds = %44
%47 = getelementptr inbounds i8*, i8** %1, i64 3
%48 = load i8*, i8** %47
%49 = tail call i32 @atoi(i8* %48)
br label %50
50: ; preds = %46, %44
%51 = phi i32 [ %49, %46 ], [ 0, %44 ]
%52 = shl i32 %38, 3
%53 = mul i32 %52, %38
%54 = tail call i64 @malloc(i32 %53)
%55 = inttoptr i64 %54 to double*
%56 = tail call i64 @malloc(i32 %53)
%57 = inttoptr i64 %56 to double*
%58 = tail call i64 @malloc(i32 %52)
%59 = inttoptr i64 %58 to double*
%60 = tail call i64 @malloc(i32 %52)
%61 = inttoptr i64 %60 to double*
%62 = icmp eq i32 %38, 0
br i1 %62, label %201, label %63
63: ; preds = %50
%64 = zext i32 %38 to i64
%65 = zext i32 %38 to i64
br label %66
66: ; preds = %86, %63
%67 = phi i64 [ 0, %63 ], [ %87, %86 ]
%68 = mul nsw i64 %67, %64
br label %76
69: ; preds = %86
br i1 %62, label %201, label %70
70: ; preds = %69
%71 = sitofp i32 %38 to double
%72 = and i64 %65, 1
%73 = icmp eq i32 %38, 1
br i1 %73, label %89, label %74
74: ; preds = %70
%75 = and i64 %65, 4294967294
br label %175
76: ; preds = %66, %76
%77 = phi i64 [ 0, %66 ], [ %84, %76 ]
%78 = tail call double (...) @drand48()
%79 = mul nsw i64 %77, %64
%80 = add nuw nsw i64 %79, %67
%81 = getelementptr inbounds double, double* %55, i64 %80
store double %78, double* %81
%82 = add nuw nsw i64 %68, %77
%83 = getelementptr inbounds double, double* %55, i64 %82
store double %78, double* %83
%84 = add nuw nsw i64 %77, 1
%85 = icmp eq i64 %84, %65
br i1 %85, label %86, label %76
86: ; preds = %76
%87 = add nuw nsw i64 %67, 1
%88 = icmp eq i64 %87, %65
br i1 %88, label %69, label %66
89: ; preds = %175, %70
%90 = phi i64 [ 0, %70 ], [ %193, %175 ]
%91 = icmp eq i64 %72, 0
br i1 %91, label %100, label %92
92: ; preds = %89
%93 = trunc i64 %90 to i32
%94 = mul nsw i32 %38, %93
%95 = sext i32 %94 to i64
%96 = add nsw i64 %90, %95
%97 = getelementptr inbounds double, double* %55, i64 %96
%98 = load double, double* %97
%99 = fadd double %98, %71
store double %99, double* %97
br label %100
100: ; preds = %89, %92
br i1 %62, label %201, label %101
101: ; preds = %100
%102 = zext i32 %38 to i64
%103 = icmp ult i32 %38, 4
br i1 %103, label %173, label %104
104: ; preds = %101
%105 = and i64 %65, 4294967292
%106 = add nsw i64 %105, -4
%107 = lshr exact i64 %106, 2
%108 = add nuw nsw i64 %107, 1
%109 = and i64 %108, 7
%110 = icmp ult i64 %106, 28
br i1 %110, label %158, label %111
111: ; preds = %104
%112 = and i64 %108, 9223372036854775800
br label %113
113: ; preds = %113, %111
%114 = phi i64 [ 0, %111 ], [ %155, %113 ]
%115 = phi i64 [ 0, %111 ], [ %156, %113 ]
%116 = getelementptr inbounds double, double* %59, i64 %114
%117 = bitcast double* %116 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %117
%118 = getelementptr inbounds double, double* %116, i64 2
%119 = bitcast double* %118 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %119
%120 = or i64 %114, 4
%121 = getelementptr inbounds double, double* %59, i64 %120
%122 = bitcast double* %121 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %122
%123 = getelementptr inbounds double, double* %121, i64 2
%124 = bitcast double* %123 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %124
%125 = or i64 %114, 8
%126 = getelementptr inbounds double, double* %59, i64 %125
%127 = bitcast double* %126 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %127
%128 = getelementptr inbounds double, double* %126, i64 2
%129 = bitcast double* %128 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %129
%130 = or i64 %114, 12
%131 = getelementptr inbounds double, double* %59, i64 %130
%132 = bitcast double* %131 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %132
%133 = getelementptr inbounds double, double* %131, i64 2
%134 = bitcast double* %133 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %134
%135 = or i64 %114, 16
%136 = getelementptr inbounds double, double* %59, i64 %135
%137 = bitcast double* %136 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %137
%138 = getelementptr inbounds double, double* %136, i64 2
%139 = bitcast double* %138 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %139
%140 = or i64 %114, 20
%141 = getelementptr inbounds double, double* %59, i64 %140
%142 = bitcast double* %141 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %142
%143 = getelementptr inbounds double, double* %141, i64 2
%144 = bitcast double* %143 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %144
%145 = or i64 %114, 24
%146 = getelementptr inbounds double, double* %59, i64 %145
%147 = bitcast double* %146 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %147
%148 = getelementptr inbounds double, double* %146, i64 2
%149 = bitcast double* %148 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %149
%150 = or i64 %114, 28
%151 = getelementptr inbounds double, double* %59, i64 %150
%152 = bitcast double* %151 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %152
%153 = getelementptr inbounds double, double* %151, i64 2
%154 = bitcast double* %153 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %154
%155 = add nuw i64 %114, 32
%156 = add i64 %115, 8
%157 = icmp eq i64 %156, %112
br i1 %157, label %158, label %113
158: ; preds = %113, %104
%159 = phi i64 [ 0, %104 ], [ %155, %113 ]
%160 = icmp eq i64 %109, 0
br i1 %160, label %171, label %161
161: ; preds = %158, %161
%162 = phi i64 [ %168, %161 ], [ %159, %158 ]
%163 = phi i64 [ %169, %161 ], [ 0, %158 ]
%164 = getelementptr inbounds double, double* %59, i64 %162
%165 = bitcast double* %164 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %165
%166 = getelementptr inbounds double, double* %164, i64 2
%167 = bitcast double* %166 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 1.000000e+00>, <2 x double>* %167
%168 = add nuw i64 %162, 4
%169 = add i64 %163, 1
%170 = icmp eq i64 %169, %109
br i1 %170, label %171, label %161
171: ; preds = %161, %158
%172 = icmp eq i64 %105, %65
br i1 %172, label %201, label %173
173: ; preds = %101, %171
%174 = phi i64 [ 0, %101 ], [ %105, %171 ]
br label %196
175: ; preds = %175, %74
%176 = phi i64 [ 0, %74 ], [ %193, %175 ]
%177 = phi i64 [ 0, %74 ], [ %194, %175 ]
%178 = trunc i64 %176 to i32
%179 = mul nsw i32 %38, %178
%180 = sext i32 %179 to i64
%181 = add nsw i64 %176, %180
%182 = getelementptr inbounds double, double* %55, i64 %181
%183 = load double, double* %182
%184 = fadd double %183, %71
store double %184, double* %182
%185 = or i64 %176, 1
%186 = trunc i64 %185 to i32
%187 = mul nsw i32 %38, %186
%188 = sext i32 %187 to i64
%189 = add nsw i64 %185, %188
%190 = getelementptr inbounds double, double* %55, i64 %189
%191 = load double, double* %190
%192 = fadd double %191, %71
store double %192, double* %190
%193 = add nuw nsw i64 %176, 2
%194 = add i64 %177, 2
%195 = icmp eq i64 %194, %75
br i1 %195, label %89, label %175
196: ; preds = %173, %196
%197 = phi i64 [ %199, %196 ], [ %174, %173 ]
%198 = getelementptr inbounds double, double* %59, i64 %197
store double 1.000000e+00, double* %198
%199 = add nuw nsw i64 %197, 1
%200 = icmp eq i64 %199, %102
br i1 %200, label %201, label %196
201: ; preds = %196, %171, %50, %69, %100
%202 = call i32 @dgemv_(i8* getelementptr inbounds ([2 x i8], [2 x i8]* @.str.1, i64 0, i64 0), i32* %3, i32* %3, double* %5, double* %55, i32* %3, double* %59, i32* %7, double* %6, double* %61, i32* %7, i32 1)
%203 = icmp eq i32 %51, 0
%204 = load i32, i32* %3
br i1 %203, label %207, label %205
205: ; preds = %201
%206 = call i32 @choleskyblk(double* %55, i32 %204, i32 %37, i32* %4)
br label %215
207: ; preds = %201
%208 = call i32 @convtile(double* %55, double* %57, i32 %204, i32 %37, i32 0)
%209 = call i32 @gettimeofday(%struct.timeval* %9, i32* null)
%210 = load i32, i32* %3
%211 = call i32 @choleskytile(double* %57, i32 %210, i32 %37, i32* %4)
%212 = call i32 @gettimeofday(%struct.timeval* %11, i32* null)
%213 = load i32, i32* %3
%214 = call i32 @convtile(double* %57, double* %55, i32 %213, i32 %37, i32 1)
br label %215
215: ; preds = %207, %205
%216 = load i32, i32* %4
%217 = icmp eq i32 %216, 0
br i1 %217, label %221, label %218
218: ; preds = %215
%219 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([44 x i8], [44 x i8]* @str.7, i64 0, i64 0))
%220 = call i32 @exit(i32 1)
unreachable
221: ; preds = %215
%222 = call i32 @dpotrs_(i8* getelementptr inbounds ([2 x i8], [2 x i8]* @.str.3, i64 0, i64 0), i32* %3, i32* %7, double* %55, i32* %3, double* %61, i32* %3, i32* %4, i32 1)
%223 = load i32, i32* %4
%224 = icmp eq i32 %223, 0
br i1 %224, label %228, label %225
225: ; preds = %221
%226 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([40 x i8], [40 x i8]* @str, i64 0, i64 0))
%227 = call i32 @exit(i32 1)
unreachable
228: ; preds = %221
%229 = load <2 x double>, <2 x double>* %10
%230 = load <2 x double>, <2 x double>* %8
%231 = fsub <2 x double> %229, %230
%232 = extractelement <2 x double> %231, i64 1
%233 = fdiv double %232, 1.000000e+06
%234 = extractelement <2 x double> %231, i64 0
%235 = fadd double %234, %233
%236 = load i32, i32* %3
%237 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([10 x i8], [10 x i8]* @.str.5, i64 0, i64 0), i32 %236, i32 %37, double %235)
%238 = load i32, i32* %3
%239 = icmp sgt i32 %238, 0
br i1 %239, label %240, label %253
240: ; preds = %228, %240
%241 = phi i64 [ %249, %240 ], [ 0, %228 ]
%242 = phi double [ %248, %240 ], [ 0.000000e+00, %228 ]
%243 = getelementptr inbounds double, double* %61, i64 %241
%244 = load double, double* %243
%245 = fsub double 1.000000e+00, %244
%246 = call i64 @fabs(double %245)
%247 = sitofp i64 %246 to double
%248 = fadd double %242, %247
%249 = add nuw nsw i64 %241, 1
%250 = load i32, i32* %3
%251 = sext i32 %250 to i64
%252 = icmp slt i64 %249, %251
br i1 %252, label %240, label %253
253: ; preds = %240, %228
%254 = phi double [ 0.000000e+00, %228 ], [ %248, %240 ]
%255 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([4 x i8], [4 x i8]* @.str.6, i64 0, i64 0), double %254)
%256 = call i32 @free(double* %55)
%257 = call i32 @free(double* %57)
%258 = call i32 @free(double* %59)
%259 = call i32 @free(double* %61)
ret i32 0
}