double *vectores(double *arr, int row, int col, int c_actual) {
  int i, j, k;
  int sum_prod_c = 0;
  int sum_prod_p = 0;
  double beta = 0;
  double mu, sigma, v_1, acumula, valor_actualizado;
  int vectorc[4] = {0, 0, 0, 0};
  double vector[4] = {0, 0, 0, 0};
  for (i = 0; i < row - c_actual; i++) {
    vectorc[i] = *(arr + ((i + c_actual) * col) + c_actual);
    sum_prod_c += vectorc[i] * vectorc[i];
    if (i > 0)
      sum_prod_p += vectorc[i] * vectorc[i];
    if (i == 0)
      vector[i] = 1;
    else
      vector[i] = vectorc[i];
    printf(" Vector completo: %d \n", vectorc[i]);
    printf(" Vector V : %f \n", vector[i]);
  }
  mu = sqrt(sum_prod_c);
  sigma = sum_prod_p;
  printf(" MU: %f \n", mu);
  printf(" sigma: %f \n\n", sigma);
  if (sigma == 0 && vectorc[0] >= 0)
    beta = 0;
  else if (sigma == 0 && vectorc[0] < 0)
    beta = -2;
  else if (vectorc[0] <= 0)
    v_1 = vectorc[0] - mu;
  else
    v_1 = -sigma / (vectorc[0] + mu);
  beta = 2 * pow(v_1, 2) / (sigma + pow(v_1, 2));
  vector[0] = v_1;
  for (i = 0; i < row - c_actual; i++) {
    vector[i] = vector[i] / v_1;
    printf(" Vector V de householder : %f \n", vector[i]);
  }
  printf(" \n Beta: %f\n", beta);
  printf(" \n V_1: %f \n\n", v_1);
  double I_Bvvt[4][4] = {{0, 0, 0}, {0, 0, 0}, {0, 0, 0}, {0, 0, 0}};
  for (i = 0; i < row - c_actual; i++) {
    for (j = 0; j < row - c_actual; j++) {
      if (i != j)
        I_Bvvt[i][j] = -1 * (beta * vector[i] * vector[j]);
      else
        I_Bvvt[i][j] = 1 - (beta * vector[i] * vector[j]);
    }
  }
  printf("I_Bvvt \n");
  show((double *)I_Bvvt, 4, 4);
  for (i = 0; i <= 4; i++) {
    for (j = 0; j < col; j++) {
      for (k = 0; k < 4; k++) {
        acumula += ((I_Bvvt[j][k]) * (*(arr + k * col + i)));
        printf("arreglo a I_Bvvt: %f\n", I_Bvvt[j][k]);
        printf("arreglo a arreglo : %f\n", (*(arr + k * col + i)));
        printf("acumula uno : %f\n", acumula);
      }
      valor_actualizado = acumula;
      acumula = 0;
    }
    *(arr + (j * col + i)) = valor_actualizado;
  }
  printf("acaboooooo \n");
  printf("Nuevo Arreglo \n");
  show((double *)arr, 4, 3);
}


x86 assembly (clang, -O0):
<vectores>:
push   %rbp
mov    %rsp,%rbp
sub    $0x120,%rsp
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    %edx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
movl   $0x0,-0x2c(%rbp)
movl   $0x0,-0x30(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
lea    -0x70(%rbp),%rdi
xor    %esi,%esi
mov    $0x10,%edx
call   1040 <memset@plt>
lea    -0x90(%rbp),%rdi
xor    %esi,%esi
mov    $0x20,%edx
call   1040 <memset@plt>
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
cmp    %ecx,%eax
jge    1294 <vectores+0x144>
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    -0x1c(%rbp),%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x1c(%rbp),%rcx
cvttsd2si (%rax,%rcx,8),%ecx
movslq -0x20(%rbp),%rax
mov    %ecx,-0x70(%rbp,%rax,4)
movslq -0x20(%rbp),%rax
mov    -0x70(%rbp,%rax,4),%eax
movslq -0x20(%rbp),%rcx
imul   -0x70(%rbp,%rcx,4),%eax
add    -0x2c(%rbp),%eax
mov    %eax,-0x2c(%rbp)
cmpl   $0x0,-0x20(%rbp)
jle    121a <vectores+0xca>
movslq -0x20(%rbp),%rax
mov    -0x70(%rbp,%rax,4),%eax
movslq -0x20(%rbp),%rcx
imul   -0x70(%rbp,%rcx,4),%eax
add    -0x30(%rbp),%eax
mov    %eax,-0x30(%rbp)
cmpl   $0x0,-0x20(%rbp)
jne    123e <vectores+0xee>
movslq -0x20(%rbp),%rax
movsd  0xde0(%rip),%xmm0
movsd  %xmm0,-0x90(%rbp,%rax,8)
ff ff
jmp    1255 <vectores+0x105>
movslq -0x20(%rbp),%rax
cvtsi2sdl -0x70(%rbp,%rax,4),%xmm0
movslq -0x20(%rbp),%rax
movsd  %xmm0,-0x90(%rbp,%rax,8)
ff ff
movslq -0x20(%rbp),%rax
mov    -0x70(%rbp,%rax,4),%esi
lea    0xdbc(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
movslq -0x20(%rbp),%rax
movsd  -0x90(%rbp,%rax,8),%xmm0
ff ff
lea    0xdb8(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    11a8 <vectores+0x58>
cvtsi2sdl -0x2c(%rbp),%xmm0
call   1070 <sqrt@plt>
movsd  %xmm0,-0x40(%rbp)
cvtsi2sdl -0x30(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x40(%rbp),%xmm0
lea    0xd8f(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
movsd  -0x48(%rbp),%xmm0
lea    0xd86(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
movsd  -0x48(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1302 <vectores+0x1b2>
jp     1302 <vectores+0x1b2>
cmpl   $0x0,-0x70(%rbp)
jl     1302 <vectores+0x1b2>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
jmp    138d <vectores+0x23d>
movsd  -0x48(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1336 <vectores+0x1e6>
jp     1336 <vectores+0x1e6>
cmpl   $0x0,-0x70(%rbp)
jge    1336 <vectores+0x1e6>
movsd  0xcd4(%rip),%xmm0
movsd  %xmm0,-0x38(%rbp)
jmp    1388 <vectores+0x238>
cmpl   $0x0,-0x70(%rbp)
jg     1354 <vectores+0x204>
cvtsi2sdl -0x70(%rbp),%xmm0
subsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
jmp    1383 <vectores+0x233>
movsd  -0x48(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
cvtsi2sdl -0x70(%rbp),%xmm1
addsd  -0x40(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
jmp    1388 <vectores+0x238>
jmp    138d <vectores+0x23d>
movsd  -0x50(%rbp),%xmm0
movsd  0xc6e(%rip),%xmm1
call   1060 <pow@plt>
movaps %xmm0,%xmm1
movsd  0xc5e(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x118(%rbp)
ff
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x120(%rbp)
ff
movsd  -0x50(%rbp),%xmm0
movsd  0xc38(%rip),%xmm1
call   1060 <pow@plt>
movsd  -0x120(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0x118(%rbp),%xmm0
ff
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
cmp    %ecx,%eax
jge    1462 <vectores+0x312>
movslq -0x20(%rbp),%rax
movsd  -0x90(%rbp,%rax,8),%xmm0
ff ff
divsd  -0x50(%rbp),%xmm0
movslq -0x20(%rbp),%rax
movsd  %xmm0,-0x90(%rbp,%rax,8)
ff ff
movslq -0x20(%rbp),%rax
movsd  -0x90(%rbp,%rax,8),%xmm0
ff ff
lea    0xc13(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1409 <vectores+0x2b9>
movsd  -0x38(%rbp),%xmm0
lea    0xc12(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
movsd  -0x50(%rbp),%xmm0
lea    0xc0c(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
lea    -0x110(%rbp),%rdi
xor    %esi,%esi
mov    $0x80,%edx
call   1040 <memset@plt>
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
cmp    %ecx,%eax
jge    15a9 <vectores+0x459>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
cmp    %ecx,%eax
jge    1596 <vectores+0x446>
mov    -0x20(%rbp),%eax
cmp    -0x24(%rbp),%eax
je     1522 <vectores+0x3d2>
movsd  -0x38(%rbp),%xmm1
movslq -0x20(%rbp),%rax
mulsd  -0x90(%rbp,%rax,8),%xmm1
ff ff
movslq -0x24(%rbp),%rax
mulsd  -0x90(%rbp,%rax,8),%xmm1
ff ff
movsd  0xb1a(%rip),%xmm0
mulsd  %xmm1,%xmm0
movslq -0x20(%rbp),%rcx
lea    -0x110(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1583 <vectores+0x433>
movsd  -0x38(%rbp),%xmm0
movslq -0x20(%rbp),%rax
mulsd  -0x90(%rbp,%rax,8),%xmm0
ff ff
movslq -0x24(%rbp),%rax
movsd  -0x90(%rbp,%rax,8),%xmm1
ff ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm1,%xmm0
movsd  0xaac(%rip),%xmm1
addsd  %xmm1,%xmm0
movslq -0x20(%rbp),%rcx
lea    -0x110(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1588 <vectores+0x438>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    14ba <vectores+0x36a>
jmp    159b <vectores+0x44b>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    14a2 <vectores+0x352>
lea    0xaeb(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
lea    -0x110(%rbp),%rdi
mov    $0x4,%edx
mov    %edx,%esi
mov    $0x0,%al
call   1050 <show@plt>
movl   $0x0,-0x20(%rbp)
cmpl   $0x4,-0x20(%rbp)
jg     170b <vectores+0x5bb>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    16e2 <vectores+0x592>
movl   $0x0,-0x28(%rbp)
cmpl   $0x4,-0x28(%rbp)
jge    16c2 <vectores+0x572>
movslq -0x24(%rbp),%rcx
lea    -0x110(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x58(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movslq -0x24(%rbp),%rcx
lea    -0x110(%rbp),%rax
shl    $0x5,%rcx
add    %rcx,%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
lea    0xa36(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%ecx
imul   -0x18(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
lea    0xa20(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
movsd  -0x58(%rbp),%xmm0
lea    0xa25(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    15f7 <vectores+0x4a7>
movsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    15e4 <vectores+0x494>
movsd  -0x60(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x18(%rbp),%ecx
add    -0x20(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    15d3 <vectores+0x483>
lea    0x9d2(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
lea    0x9d1(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
mov    -0x10(%rbp),%rdi
mov    $0x4,%esi
mov    $0x3,%edx
mov    $0x0,%al
call   1050 <show@plt>
mov    -0x8(%rbp),%rax
add    $0x120,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<vectores>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xf8,%rsp
mov    %esi,%r14d
mov    %rdi,%r15
xorps  %xmm0,%xmm0
movaps %xmm0,0x50(%rsp)
movaps %xmm0,0x40(%rsp)
sub    %ecx,%r14d
mov    %edx,0x24(%rsp)
movslq %edx,%r12
mov    %r14d,%r13d
test   %r14d,%r14d
mov    %rdi,0x38(%rsp)
jle    1261 <vectores+0x121>
mov    %r14d,0x10(%rsp)
movslq %ecx,%r14
mov    %r14,%rax
imul   %r12,%rax
lea    (%r15,%rax,8),%rax
cvttsd2si (%rax,%r14,8),%esi
mov    %esi,%ebp
imul   %esi,%ebp
lea    0xe7c(%rip),%rdi
xor    %ebx,%ebx
mov    %esi,0x28(%rsp)
xor    %eax,%eax
call   1040 <printf@plt>
lea    0xe7f(%rip),%rdi
movsd  0xe40(%rip),%xmm0
mov    $0x1,%al
call   1040 <printf@plt>
cmpl   $0x1,0x10(%rsp)
je     125a <vectores+0x11a>
lea    0x1(%r14),%rax
imul   %r12,%rax
add    %r14,%rax
lea    (%r15,%rax,8),%r14
mov    %r12,0x30(%rsp)
lea    0x0(,%r12,8),%r15
mov    $0x1,%r12d
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
movsd  (%r14),%xmm0
cvttsd2si %xmm0,%esi
mov    %esi,%eax
imul   %esi,%eax
add    %eax,%ebp
cvttpd2dq %xmm0,%xmm0
add    %eax,%ebx
cvtdq2pd %xmm0,%xmm0
movaps %xmm0,(%rsp)
movlps %xmm0,0x40(%rsp,%r12,8)
lea    0xdf5(%rip),%rdi
xor    %eax,%eax
call   1040 <printf@plt>
lea    0xdfe(%rip),%rdi
movaps (%rsp),%xmm0
mov    $0x1,%al
call   1040 <printf@plt>
add    $0x1,%r12
add    %r15,%r14
cmp    %r12,%r13
jne    1200 <vectores+0xc0>
mov    0x38(%rsp),%r15
mov    0x30(%rsp),%r12
mov    0x10(%rsp),%r14d
jmp    126d <vectores+0x12d>
movl   $0x0,0x28(%rsp)
xor    %ebp,%ebp
xor    %ebx,%ebx
xorps  %xmm0,%xmm0
cvtsi2sd %ebp,%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     1284 <vectores+0x144>
sqrtsd %xmm0,%xmm0
jmp    1289 <vectores+0x149>
call   1060 <sqrt@plt>
xorps  %xmm1,%xmm1
cvtsi2sd %ebx,%xmm1
movapd %xmm1,(%rsp)
lea    0xdac(%rip),%rdi
movsd  %xmm0,0x10(%rsp)
mov    $0x1,%al
call   1040 <printf@plt>
lea    0xda2(%rip),%rdi
movapd (%rsp),%xmm0
mov    $0x1,%al
call   1040 <printf@plt>
test   %ebx,%ebx
je     12f0 <vectores+0x1b0>
mov    0x28(%rsp),%eax
test   %eax,%eax
jle    12f2 <vectores+0x1b2>
movapd 0xd40(%rip),%xmm2
movapd (%rsp),%xmm3
xorpd  %xmm3,%xmm2
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  0x10(%rsp),%xmm1
addsd  %xmm0,%xmm1
divsd  %xmm1,%xmm2
jmp    1301 <vectores+0x1c1>
jmp    12fc <vectores+0x1bc>
cvtsi2sd %eax,%xmm2
subsd  0x10(%rsp),%xmm2
movapd (%rsp),%xmm3
movapd %xmm2,(%rsp)
movapd %xmm2,%xmm0
mulsd  %xmm2,%xmm0
movapd %xmm0,%xmm1
addsd  %xmm0,%xmm1
addsd  %xmm3,%xmm0
divsd  %xmm0,%xmm1
movapd %xmm1,0x10(%rsp)
movsd  %xmm2,0x40(%rsp)
test   %r14d,%r14d
jle    1384 <vectores+0x244>
movapd (%rsp),%xmm0
divsd  %xmm0,%xmm0
movsd  %xmm0,0x40(%rsp)
lea    0xd1b(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
cmp    $0x1,%r14d
je     1384 <vectores+0x244>
mov    $0x1,%ebp
lea    0xd02(%rip),%rbx
xchg   %ax,%ax
movsd  0x40(%rsp,%rbp,8),%xmm0
divsd  (%rsp),%xmm0
movsd  %xmm0,0x40(%rsp,%rbp,8)
mov    %rbx,%rdi
mov    $0x1,%al
call   1040 <printf@plt>
add    $0x1,%rbp
cmp    %rbp,%r13
jne    1360 <vectores+0x220>
lea    0xcf5(%rip),%rdi
movaps 0x10(%rsp),%xmm0
mov    $0x1,%al
call   1040 <printf@plt>
lea    0xcef(%rip),%rdi
movapd (%rsp),%xmm0
mov    $0x1,%al
call   1040 <printf@plt>
movapd 0x10(%rsp),%xmm4
xorpd  %xmm0,%xmm0
movapd %xmm0,0xe0(%rsp)
movapd %xmm0,0xd0(%rsp)
movapd %xmm0,0xc0(%rsp)
movapd %xmm0,0xb0(%rsp)
movapd %xmm0,0xa0(%rsp)
movapd %xmm0,0x90(%rsp)
movapd %xmm0,0x80(%rsp)
movapd %xmm0,0x70(%rsp)
test   %r14d,%r14d
jle    1473 <vectores+0x333>
xorpd  0xc0a(%rip),%xmm4
lea    0x70(%rsp),%rax
xor    %ecx,%ecx
movsd  0xbeb(%rip),%xmm0
jmp    142d <vectores+0x2ed>
nopw   0x0(%rax,%rax,1)
add    $0x1,%rcx
add    $0x20,%rax
cmp    %r13,%rcx
je     1473 <vectores+0x333>
movsd  0x40(%rsp,%rcx,8),%xmm1
movapd %xmm1,%xmm2
mulsd  %xmm4,%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
xor    %edx,%edx
jmp    145e <vectores+0x31e>
nopw   0x0(%rax,%rax,1)
movsd  %xmm3,(%rax,%rdx,8)
add    $0x1,%rdx
cmp    %rdx,%r13
je     1420 <vectores+0x2e0>
movapd %xmm1,%xmm3
cmp    %rdx,%rcx
je     1450 <vectores+0x310>
movsd  0x40(%rsp,%rdx,8),%xmm3
mulsd  %xmm2,%xmm3
jmp    1450 <vectores+0x310>
lea    0xc61(%rip),%rdi
call   1030 <puts@plt>
xor    %ebp,%ebp
lea    0x70(%rsp),%rdi
mov    $0x4,%esi
mov    $0x4,%edx
xor    %eax,%eax
call   1050 <show@plt>
mov    %r12d,%edi
lea    (%r15,%r12,8),%r8
mov    %r12,%rdx
shl    $0x4,%rdx
add    %r15,%rdx
lea    (%r12,%r12,2),%rax
lea    (%r15,%rax,8),%rsi
shl    $0x5,%rdi
xorpd  %xmm1,%xmm1
mov    0x24(%rsp),%ecx
mov    %rdx,0x10(%rsp)
mov    %r8,0x28(%rsp)
mov    %rsi,0x30(%rsp)
mov    %rdi,0x68(%rsp)
mov    0x28(%rsp),%r12
jmp    14fb <vectores+0x3bb>
nopw   0x0(%rax,%rax,1)
imul   %ecx,%eax
add    %ebp,%eax
cltq
movsd  %xmm0,(%r15,%rax,8)
add    $0x1,%rbp
cmp    $0x5,%rbp
je     167b <vectores+0x53b>
mov    $0x0,%eax
test   %ecx,%ecx
jle    14e0 <vectores+0x3a0>
xor    %r15d,%r15d
nopw   0x0(%rax,%rax,1)
movsd  0x70(%rsp,%r15,1),%xmm0
mov    0x38(%rsp),%r13
movsd  0x0(%r13,%rbp,8),%xmm2
mulsd  %xmm0,%xmm2
addsd  %xmm1,%xmm2
movsd  %xmm2,(%rsp)
lea    0xb64(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
movsd  0x0(%r13,%rbp,8),%xmm0
lea    0xb65(%rip),%r13
mov    %r13,%rdi
mov    $0x1,%al
call   1040 <printf@plt>
lea    0xb6c(%rip),%r14
mov    %r14,%rdi
movsd  (%rsp),%xmm0
mov    $0x1,%al
call   1040 <printf@plt>
movsd  0x78(%rsp,%r15,1),%xmm0
movsd  (%r12,%rbp,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rsp),%xmm1
movsd  %xmm1,(%rsp)
lea    0xb0d(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
movsd  (%r12,%rbp,8),%xmm0
mov    %r13,%rdi
mov    $0x1,%al
call   1040 <printf@plt>
mov    %r14,%rdi
movsd  (%rsp),%xmm0
mov    $0x1,%al
call   1040 <printf@plt>
movsd  0x80(%rsp,%r15,1),%xmm0
mov    0x10(%rsp),%rbx
movsd  (%rbx,%rbp,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rsp),%xmm1
movsd  %xmm1,(%rsp)
lea    0xabe(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
movsd  (%rbx,%rbp,8),%xmm0
mov    %r13,%rdi
mov    $0x1,%al
call   1040 <printf@plt>
mov    %r14,%rdi
movsd  (%rsp),%xmm0
mov    $0x1,%al
call   1040 <printf@plt>
movsd  0x88(%rsp,%r15,1),%xmm0
mov    0x30(%rsp),%rbx
movsd  (%rbx,%rbp,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rsp),%xmm1
movsd  %xmm1,(%rsp)
lea    0xa70(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
movsd  (%rbx,%rbp,8),%xmm0
mov    %r13,%rdi
mov    $0x1,%al
call   1040 <printf@plt>
mov    %r14,%rdi
movsd  (%rsp),%xmm0
mov    $0x1,%al
call   1040 <printf@plt>
xorpd  %xmm1,%xmm1
mov    0x68(%rsp),%rdi
add    $0x20,%r15
cmp    %r15,%rdi
jne    1510 <vectores+0x3d0>
mov    0x24(%rsp),%ecx
mov    %ecx,%eax
mov    0x38(%rsp),%r15
movsd  (%rsp),%xmm0
jmp    14e0 <vectores+0x3a0>
lea    0xa61(%rip),%rdi
call   1030 <puts@plt>
lea    0xa61(%rip),%rdi
call   1030 <puts@plt>
mov    %r15,%rdi
mov    $0x4,%esi
mov    $0x3,%edx
xor    %eax,%eax
call   1050 <show@plt>
add    $0xf8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
@.str = [23 x i8] c" Vector completo: %d \0A\00"
@.str.1 = [17 x i8] c" Vector V : %f \0A\00"
@.str.2 = [10 x i8] c" MU: %f \0A\00"
@.str.3 = [14 x i8] c" sigma: %f \0A\0A\00"
@.str.4 = [32 x i8] c" Vector V de householder : %f \0A\00"
@.str.5 = [13 x i8] c" \0A Beta: %f\0A\00"
@.str.6 = [14 x i8] c" \0A V_1: %f \0A\0A\00"
@.str.8 = [22 x i8] c"arreglo a I_Bvvt: %f\0A\00"
@.str.9 = [24 x i8] c"arreglo a arreglo : %f\0A\00"
@.str.10 = [18 x i8] c"acumula uno : %f\0A\00"
@str = [8 x i8] c"I_Bvvt \00"
@str.13 = [12 x i8] c"acaboooooo \00"
@str.14 = [15 x i8] c"Nuevo Arreglo \00"
define double* @vectores(double* %0, i32 %1, i32 %2, i32 %3) {
%5 = alloca [4 x double]
%6 = alloca [4 x [4 x double]]
%7 = bitcast [4 x double]* %5 to i8*
%8 = sub i32 %1, %3
%9 = icmp sgt i32 %8, 0
br i1 %9, label %10, label %43
10: ; preds = %4
%11 = sext i32 %3 to i64
%12 = sext i32 %2 to i64
%13 = zext i32 %8 to i64
%14 = mul nsw i64 %11, %12
%15 = getelementptr inbounds double, double* %0, i64 %14
%16 = getelementptr inbounds double, double* %15, i64 %11
%17 = load double, double* %16
%18 = fptosi double %17 to i32
%19 = mul nsw i32 %18, %18
%20 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([23 x i8], [23 x i8]* @.str, i64 0, i64 0), i32 %18)
%21 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([17 x i8], [17 x i8]* @.str.1, i64 0, i64 0), double 1.000000e+00)
%22 = icmp eq i32 %8, 1
br i1 %22, label %43, label %23
23: ; preds = %10
%24 = getelementptr inbounds double, double* %0, i64 %11
br label %25
25: ; preds = %23, %25
%26 = phi i64 [ 1, %23 ], [ %41, %25 ]
%27 = phi i32 [ 0, %23 ], [ %36, %25 ]
%28 = phi i32 [ %19, %23 ], [ %35, %25 ]
%29 = add nsw i64 %26, %11
%30 = mul nsw i64 %29, %12
%31 = getelementptr inbounds double, double* %24, i64 %30
%32 = load double, double* %31
%33 = fptosi double %32 to i32
%34 = mul nsw i32 %33, %33
%35 = add nuw nsw i32 %34, %28
%36 = add nuw nsw i32 %34, %27
%37 = sitofp i32 %33 to double
%38 = getelementptr inbounds [4 x double], [4 x double]* %5, i64 0, i64 %26
store double %37, double* %38
%39 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([23 x i8], [23 x i8]* @.str, i64 0, i64 0), i32 %33)
%40 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([17 x i8], [17 x i8]* @.str.1, i64 0, i64 0), double %37)
%41 = add nuw nsw i64 %26, 1
%42 = icmp eq i64 %41, %13
br i1 %42, label %43, label %25
43: ; preds = %25, %10, %4
%44 = phi i32 [ 0, %4 ], [ %18, %10 ], [ %18, %25 ]
%45 = phi i32 [ 0, %4 ], [ %19, %10 ], [ %35, %25 ]
%46 = phi i32 [ 0, %4 ], [ 0, %10 ], [ %36, %25 ]
%47 = sitofp i32 %45 to double
%48 = tail call double @sqrt(double %47)
%49 = sitofp i32 %46 to double
%50 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([10 x i8], [10 x i8]* @.str.2, i64 0, i64 0), double %48)
%51 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([14 x i8], [14 x i8]* @.str.3, i64 0, i64 0), double %49)
%52 = icmp eq i32 %46, 0
br i1 %52, label %63, label %53
53: ; preds = %43
%54 = icmp slt i32 %44, 1
br i1 %54, label %55, label %58
55: ; preds = %53
%56 = sitofp i32 %44 to double
%57 = fsub double %56, %48
br label %63
58: ; preds = %53
%59 = fneg double %49
%60 = sitofp i32 %44 to double
%61 = fadd double %48, %60
%62 = fdiv double %59, %61
br label %63
63: ; preds = %43, %58, %55
%64 = phi double [ %57, %55 ], [ %62, %58 ], [ undef, %43 ]
%65 = fmul double %64, %64
%66 = fmul double %65, 2.000000e+00
%67 = fadd double %65, %49
%68 = fdiv double %66, %67
%69 = getelementptr inbounds [4 x double], [4 x double]* %5, i64 0, i64 0
store double %64, double* %69
br i1 %9, label %70, label %85
70: ; preds = %63
%71 = zext i32 %8 to i64
%72 = getelementptr inbounds [4 x double], [4 x double]* %5, i64 0, i64 0
%73 = fdiv double %64, %64
store double %73, double* %72
%74 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([32 x i8], [32 x i8]* @.str.4, i64 0, i64 0), double %73)
%75 = icmp eq i32 %8, 1
br i1 %75, label %85, label %76
76: ; preds = %70, %76
%77 = phi i64 [ %83, %76 ], [ 1, %70 ]
%78 = getelementptr inbounds [4 x double], [4 x double]* %5, i64 0, i64 %77
%79 = load double, double* %78
%80 = getelementptr inbounds [4 x double], [4 x double]* %5, i64 0, i64 %77
%81 = fdiv double %79, %64
store double %81, double* %80
%82 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([32 x i8], [32 x i8]* @.str.4, i64 0, i64 0), double %81)
%83 = add nuw nsw i64 %77, 1
%84 = icmp eq i64 %83, %71
br i1 %84, label %85, label %76
85: ; preds = %76, %70, %63
%86 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([13 x i8], [13 x i8]* @.str.5, i64 0, i64 0), double %68)
%87 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([14 x i8], [14 x i8]* @.str.6, i64 0, i64 0), double %64)
%88 = bitcast [4 x [4 x double]]* %6 to i8*
br i1 %9, label %89, label %114
89: ; preds = %85
%90 = fneg double %68
%91 = zext i32 %8 to i64
br label %92
92: ; preds = %111, %89
%93 = phi i64 [ 0, %89 ], [ %112, %111 ]
%94 = getelementptr inbounds [4 x double], [4 x double]* %5, i64 0, i64 %93
%95 = load double, double* %94
%96 = fmul double %95, %90
%97 = fmul double %95, %90
%98 = tail call double @llvm.fmuladd.f64(double %97, double %95, double 1.000000e+00)
br label %99
99: ; preds = %92, %106
%100 = phi i64 [ 0, %92 ], [ %109, %106 ]
%101 = icmp eq i64 %93, %100
br i1 %101, label %106, label %102
102: ; preds = %99
%103 = getelementptr inbounds [4 x double], [4 x double]* %5, i64 0, i64 %100
%104 = load double, double* %103
%105 = fmul double %96, %104
br label %106
106: ; preds = %99, %102
%107 = phi double [ %105, %102 ], [ %98, %99 ]
%108 = getelementptr inbounds [4 x [4 x double]], [4 x [4 x double]]* %6, i64 0, i64 %93, i64 %100
store double %107, double* %108
%109 = add nuw nsw i64 %100, 1
%110 = icmp eq i64 %109, %91
br i1 %110, label %111, label %99
111: ; preds = %106
%112 = add nuw nsw i64 %93, 1
%113 = icmp eq i64 %112, %91
br i1 %113, label %114, label %92
114: ; preds = %111, %85
%115 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([8 x i8], [8 x i8]* @str, i64 0, i64 0))
%116 = getelementptr inbounds [4 x [4 x double]], [4 x [4 x double]]* %6, i64 0, i64 0, i64 0
%117 = call i32 (double*, i32, i32, ...) bitcast (i32 (...)* @show to i32 (double*, i32, i32, ...)*)(double* %116, i32 4, i32 4)
%118 = icmp sgt i32 %2, 0
%119 = sext i32 %2 to i64
%120 = zext i32 %2 to i64
%121 = getelementptr inbounds double, double* %0, i64 %119
%122 = shl nsw i64 %119, 1
%123 = mul nsw i64 %119, 3
%124 = getelementptr inbounds double, double* %0, i64 %122
%125 = getelementptr inbounds double, double* %0, i64 %123
br label %126
126: ; preds = %114, %170
%127 = phi i64 [ 0, %114 ], [ %178, %170 ]
%128 = phi double [ undef, %114 ], [ %172, %170 ]
br i1 %118, label %129, label %170
129: ; preds = %126
%130 = getelementptr inbounds double, double* %0, i64 %127
%131 = getelementptr inbounds double, double* %121, i64 %127
%132 = getelementptr inbounds double, double* %124, i64 %127
%133 = getelementptr inbounds double, double* %125, i64 %127
br label %134
134: ; preds = %129, %134
%135 = phi i64 [ 0, %129 ], [ %168, %134 ]
%136 = getelementptr inbounds [4 x [4 x double]], [4 x [4 x double]]* %6, i64 0, i64 %135, i64 0
%137 = load double, double* %136
%138 = load double, double* %130
%139 = call double @llvm.fmuladd.f64(double %137, double %138, double 0.000000e+00)
%140 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([22 x i8], [22 x i8]* @.str.8, i64 0, i64 0), double %137)
%141 = load double, double* %130
%142 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([24 x i8], [24 x i8]* @.str.9, i64 0, i64 0), double %141)
%143 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([18 x i8], [18 x i8]* @.str.10, i64 0, i64 0), double %139)
%144 = getelementptr inbounds [4 x [4 x double]], [4 x [4 x double]]* %6, i64 0, i64 %135, i64 1
%145 = load double, double* %144
%146 = load double, double* %131
%147 = call double @llvm.fmuladd.f64(double %145, double %146, double %139)
%148 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([22 x i8], [22 x i8]* @.str.8, i64 0, i64 0), double %145)
%149 = load double, double* %131
%150 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([24 x i8], [24 x i8]* @.str.9, i64 0, i64 0), double %149)
%151 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([18 x i8], [18 x i8]* @.str.10, i64 0, i64 0), double %147)
%152 = getelementptr inbounds [4 x [4 x double]], [4 x [4 x double]]* %6, i64 0, i64 %135, i64 2
%153 = load double, double* %152
%154 = load double, double* %132
%155 = call double @llvm.fmuladd.f64(double %153, double %154, double %147)
%156 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([22 x i8], [22 x i8]* @.str.8, i64 0, i64 0), double %153)
%157 = load double, double* %132
%158 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([24 x i8], [24 x i8]* @.str.9, i64 0, i64 0), double %157)
%159 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([18 x i8], [18 x i8]* @.str.10, i64 0, i64 0), double %155)
%160 = getelementptr inbounds [4 x [4 x double]], [4 x [4 x double]]* %6, i64 0, i64 %135, i64 3
%161 = load double, double* %160
%162 = load double, double* %133
%163 = call double @llvm.fmuladd.f64(double %161, double %162, double %155)
%164 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([22 x i8], [22 x i8]* @.str.8, i64 0, i64 0), double %161)
%165 = load double, double* %133
%166 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([24 x i8], [24 x i8]* @.str.9, i64 0, i64 0), double %165)
%167 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([18 x i8], [18 x i8]* @.str.10, i64 0, i64 0), double %163)
%168 = add nuw nsw i64 %135, 1
%169 = icmp eq i64 %168, %120
br i1 %169, label %170, label %134
170: ; preds = %134, %126
%171 = phi i32 [ 0, %126 ], [ %2, %134 ]
%172 = phi double [ %128, %126 ], [ %163, %134 ]
%173 = mul nsw i32 %171, %2
%174 = trunc i64 %127 to i32
%175 = add nsw i32 %173, %174
%176 = sext i32 %175 to i64
%177 = getelementptr inbounds double, double* %0, i64 %176
store double %172, double* %177
%178 = add nuw nsw i64 %127, 1
%179 = icmp eq i64 %178, 5
br i1 %179, label %180, label %126
180: ; preds = %170
%181 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([12 x i8], [12 x i8]* @str.13, i64 0, i64 0))
%182 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([15 x i8], [15 x i8]* @str.14, i64 0, i64 0))
%183 = call i32 (double*, i32, i32, ...) bitcast (i32 (...)* @show to i32 (double*, i32, i32, ...)*)(double* %0, i32 4, i32 3)
ret double* undef
}